diff --git a/src/common/turb/mode_turb_ver.F90 b/src/common/turb/mode_turb_ver.F90
index 5b89450718390f304da1b73d36460ef26fe4974f..2ba6afd5310f821681d3cdbd22a876e13c09096b 100644
--- a/src/common/turb/mode_turb_ver.F90
+++ b/src/common/turb/mode_turb_ver.F90
@@ -224,6 +224,7 @@ USE MODE_EMOIST, ONLY: EMOIST
USE MODE_ETHETA, ONLY: ETHETA
USE MODI_GRADIENT_M
USE MODI_GRADIENT_W
+USE MODE_GRADIENT_M_PHY, ONLY : GZ_M_W_PHY
USE MODI_TURB
USE MODE_TURB_VER_THERMO_FLUX, ONLY: TURB_VER_THERMO_FLUX
USE MODE_TURB_VER_THERMO_CORR, ONLY: TURB_VER_THERMO_CORR
@@ -387,8 +388,8 @@ REAL, DIMENSION(D%NIT,D%NJT,D%NKT,KSV) :: &
REAL, DIMENSION(D%NIT,D%NJT,D%NKT) :: ZLM
!
LOGICAL :: GUSERV ! flag to use water vapor
-INTEGER :: IKB,IKE ! index value for the Beginning
- ! and the End of the physical domain for the mass points
+INTEGER :: IKB,IKE,IIB,IIE,IJB,IJE ! index value for the Beginning
+ ! and the End of the physical domain for the mass points
INTEGER :: JSV,JI,JJ,JK ! loop counter
REAL :: ZTIME1
REAL :: ZTIME2
@@ -402,8 +403,12 @@ TYPE(TFIELDDATA) :: TZFIELD
!
IF (LHOOK) CALL DR_HOOK('TURB_VER',0,ZHOOK_HANDLE)
!
-IKB=D%NKB
-IKE=D%NKE
+IKB=D%NKTB
+IKE=D%NKTE
+IIE=D%NIEC
+IIB=D%NIBC
+IJE=D%NJEC
+IJB=D%NJBC
!
!
! 3D Redelsperger numbers
@@ -424,36 +429,38 @@ CALL PRANDTL(D,CST,CSTURB,KRR,KSV,KRRI,OTURB_FLX, &
! Buoyancy coefficient
!
IF (OOCEAN) THEN
- !$mnh_expand_array(JI=1:D%NIT,JJ=1:D%NJT,JK=1:D%NKT)
- ZBETA(:,:,:) = CST%XG*CST%XALPHAOC
- !$mnh_end_expand_array(JI=1:D%NIT,JJ=1:D%NJT,JK=1:D%NKT)
+ !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
+ ZBETA(IIB:IIE,IJB:IJE,1:D%NKT) = CST%XG*CST%XALPHAOC
+ !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
ELSE
- !$mnh_expand_array(JI=1:D%NIT,JJ=1:D%NJT,JK=1:D%NKT)
- ZBETA(:,:,:) = CST%XG/PTHVREF(:,:,:)
- !$mnh_end_expand_array(JI=1:D%NIT,JJ=1:D%NJT,JK=1:D%NKT)
+ !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
+ ZBETA(IIB:IIE,IJB:IJE,1:D%NKT) = CST%XG/PTHVREF(IIB:IIE,IJB:IJE,1:D%NKT)
+ !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
END IF
!
! Square root of Tke
!
-!$mnh_expand_array(JI=1:D%NIT,JJ=1:D%NJT,JK=1:D%NKT)
-ZSQRT_TKE(:,:,:) = SQRT(PTKEM(:,:,:))
-!$mnh_end_expand_array(JI=1:D%NIT,JJ=1:D%NJT,JK=1:D%NKT)
+!$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
+ZSQRT_TKE(IIB:IIE,IJB:IJE,1:D%NKT) = SQRT(PTKEM(IIB:IIE,IJB:IJE,1:D%NKT))
+!$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
!
! gradients of mean quantities at previous time-step
!
-ZDTH_DZ = GZ_M_W(D%NKA, D%NKU, D%NKL,PTHLM(:,:,:),PDZZ)
-ZDR_DZ = 0.
-IF (KRR>0) ZDR_DZ = GZ_M_W(D%NKA, D%NKU, D%NKL,PRM(:,:,:,1),PDZZ)
+CALL GZ_M_W_PHY(D,PTHLM,PDZZ,ZDTH_DZ)
+ZDR_DZ(:,:,:) = 0.
+IF (KRR>0) CALL GZ_M_W_PHY(D,PRM(:,:,:,1),PDZZ,ZDR_DZ)
!
!
! Denominator factor in 3rd order terms
!
IF (.NOT. OHARAT) THEN
- !$mnh_expand_array(JI=1:D%NIT,JJ=1:D%NJT,JK=1:D%NKT)
- ZD(:,:,:) = (1.+ZREDTH1(:,:,:)+ZREDR1(:,:,:)) * (1.+0.5*(ZREDTH1(:,:,:)+ZREDR1(:,:,:)))
- !$mnh_end_expand_array(JI=1:D%NIT,JJ=1:D%NJT,JK=1:D%NKT)
+ !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
+ ZD(IIB:IIE,IJB:IJE,1:D%NKT) = (1.+ZREDTH1(IIB:IIE,IJB:IJE,1:D%NKT)+ZREDR1(IIB:IIE,IJB:IJE,1:D%NKT)) * (1.+0.5*(ZREDTH1(IIB:IIE,IJB:IJE,1:D%NKT)+ZREDR1(IIB:IIE,IJB:IJE,1:D%NKT)))
+ !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
ELSE
- ZD(:,:,:) = 1.
+ !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
+ ZD(IIB:IIE,IJB:IJE,1:D%NKT) = 1.
+ !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
ENDIF
!
! Phi3 and Psi3 Prandtl numbers
@@ -571,7 +578,7 @@ CALL TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,KSV,O2D,OFLAT, &
!
IF (OHARAT) ZLM=PLENGTHH
!
-IF (SIZE(PSVM,4)>0) &
+IF (KSV>0) &
CALL TURB_VER_SV_FLUX(D,CST,CSTURB,ONOMIXLG, &
KSV,KSV_LGBEG,KSV_LGEND, &
OTURB_FLX,HTURBDIM,OHARAT,OBLOWSNOW,OLES_CALL,&
@@ -585,7 +592,7 @@ CALL TURB_VER_SV_FLUX(D,CST,CSTURB,ONOMIXLG, &
PRSVS,PWSV )
!
!
-IF (SIZE(PSVM,4)>0 .AND. OLES_CALL) &
+IF (KSV>0 .AND. OLES_CALL) &
CALL TURB_VER_SV_CORR(D,CST,CSTURB,KRR,KRRL,KRRI,OOCEAN, &
PDZZ,KSV,KSV_LGBEG,KSV_LGEND,ONOMIXLG, &
OBLOWSNOW,OLES_CALL,OCOMPUTE_SRC, &
diff --git a/src/common/turb/mode_turb_ver.f90 b/src/common/turb/mode_turb_ver.f90
deleted file mode 100644
index 2ba6afd5310f821681d3cdbd22a876e13c09096b..0000000000000000000000000000000000000000
--- a/src/common/turb/mode_turb_ver.f90
+++ /dev/null
@@ -1,672 +0,0 @@
-!MNH_LIC Copyright 1994-2022 CNRS, Meteo-France and Universite Paul Sabatier
-!MNH_LIC This is part of the Meso-NH software governed by the CeCILL-C licence
-!MNH_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt
-!MNH_LIC for details. version 1.
-MODULE MODE_TURB_VER
-IMPLICIT NONE
-CONTAINS
-SUBROUTINE TURB_VER(D,CST,CSTURB,TURBN,KRR,KRRL,KRRI, &
- OTURB_FLX,OOCEAN,ODEEPOC,OHARAT,OCOMPUTE_SRC, &
- KSV,KSV_LGBEG,KSV_LGEND, &
- HTURBDIM,HTOM,PIMPL,PEXPL, &
- HPROGRAM, O2D, ONOMIXLG, OFLAT, &
- OLES_CALL,OCOUPLES,OBLOWSNOW, ORMC01, &
- PTSTEP, TPFILE, &
- PDXX,PDYY,PDZZ,PDZX,PDZY,PDIRCOSZW,PZZ, &
- PCOSSLOPE,PSINSLOPE, &
- PRHODJ,PTHVREF, &
- PSFTHM,PSFRM,PSFSVM,PSFTHP,PSFRP,PSFSVP, &
- PCDUEFF,PTAU11M,PTAU12M,PTAU33M, &
- PUM,PVM,PWM,PUSLOPEM,PVSLOPEM,PTHLM,PRM,PSVM, &
- PTKEM,PLM,PLENGTHM,PLENGTHH,PLEPS,MFMOIST, &
- PLOCPEXNM,PATHETA,PAMOIST,PSRCM,PFRAC_ICE, &
- PFWTH,PFWR,PFTH2,PFR2,PFTHR,PBL_DEPTH, &
- PSBL_DEPTH,PLMO, &
- PRUS,PRVS,PRWS,PRTHLS,PRRS,PRSVS, &
- PDP,PTP,PSIGS,PWTH,PWRC,PWSV )
-! ###############################################################
-!
-!
-!!**** *TURB_VER* -compute the source terms due to the vertical turbulent
-!! fluxes.
-!!
-!! PURPOSE
-!! -------
-! The purpose of this routine is to compute the vertical turbulent
-! fluxes of the evolutive variables and give back the source
-! terms to the main program. In the case of large horizontal meshes,
-! the divergence of these vertical turbulent fluxes represent the whole
-! effect of the turbulence but when the three-dimensionnal version of
-! the turbulence scheme is activated (CTURBDIM="3DIM"), these divergences
-! are completed in the next routine TURB_HOR.
-! An arbitrary degree of implicitness has been implemented for the
-! temporal treatment of these diffusion terms.
-! The vertical boundary conditions are as follows:
-! * at the bottom, the surface fluxes are prescribed at the same
-! as the other turbulent fluxes
-! * at the top, the turbulent fluxes are set to 0.
-! It should be noted that the condensation has been implicitely included
-! in this turbulence scheme by using conservative variables and computing
-! the subgrid variance of a statistical variable s indicating the presence
-! or not of condensation in a given mesh.
-!
-!!** METHOD
-!! ------
-!! 1D type calculations are made;
-!! The vertical turbulent fluxes are computed in an off-centered
-!! implicit scheme (a Crank-Nicholson type with coefficients different
-!! than 0.5), which allows to vary the degree of implicitness of the
-!! formulation.
-!! The different prognostic variables are treated one by one.
-!! The contributions of each turbulent fluxes are cumulated into the
-!! tendency PRvarS, and into the dynamic and thermal production of
-!! TKE if necessary.
-!!
-!! In section 2 and 3, the thermodynamical fields are considered.
-!! Only the turbulent fluxes of the conservative variables
-!! (Thetal and Rnp stored in PRx(:,:,:,1)) are computed.
-!! Note that the turbulent fluxes at the vertical
-!! boundaries are given either by the soil scheme for the surface one
-!! ( at the same instant as the others fluxes) and equal to 0 at the
-!! top of the model. The thermal production is computed by vertically
-!! averaging the turbulent flux and multiply this flux at the mass point by
-!! a function ETHETA or EMOIST, which preform the transformation from the
-!! conservative variables to the virtual potential temperature.
-!!
-!! In section 4, the variance of the statistical variable
-!! s indicating presence or not of condensation, is determined in function
-!! of the turbulent moments of the conservative variables and its
-!! squarred root is stored in PSIGS. This information will be completed in
-!! the horizontal turbulence if the turbulence dimensionality is not
-!! equal to "1DIM".
-!!
-!! In section 5, the x component of the stress tensor is computed.
-!! The surface flux <u'w'> is computed from the value of the surface
-!! fluxes computed in axes linked to the orography ( i", j" , k"):
-!! i" is parallel to the surface and in the direction of the maximum
-!! slope
-!! j" is also parallel to the surface and in the normal direction of
-!! the maximum slope
-!! k" is the normal to the surface
-!! In order to prevent numerical instability, the implicit scheme has
-!! been extended to the surface flux regarding to its dependence in
-!! function of U. The dependence in function of the other components
-!! introduced by the different rotations is only explicit.
-!! The turbulent fluxes are used to compute the dynamic production of
-!! TKE. For the last TKE level ( located at PDZZ(:,:,IKB)/2 from the
-!! ground), an harmonic extrapolation from the dynamic production at
-!! PDZZ(:,:,IKB) is used to avoid an evaluation of the gradient of U
-!! in the surface layer.
-!!
-!! In section 6, the same steps are repeated but for the y direction
-!! and in section 7, a diagnostic computation of the W variance is
-!! performed.
-!!
-!! In section 8, the turbulent fluxes for the scalar variables are
-!! computed by the same way as the conservative thermodynamical variables
-!!
-!!
-!! EXTERNAL
-!! --------
-!! GX_U_M, GY_V_M, GZ_W_M : cartesian gradient operators
-!! GX_U_UW,GY_V_VW (X,Y,Z) represent the direction of the gradient
-!! _(M,U,...)_ represent the localization of the
-!! field to be derivated
-!! _(M,UW,...) represent the localization of the
-!! field derivated
-!!
-!! SUBROUTINE TRIDIAG : to compute the split implicit evolution
-!! of a variable located at a mass point
-!!
-!! SUBROUTINE TRIDIAG_WIND: to compute the split implicit evolution
-!! of a variable located at a wind point
-!!
-!! FUNCTIONs ETHETA and EMOIST :
-!! allows to compute:
-!! - the coefficients for the turbulent correlation between
-!! any variable and the virtual potential temperature, of its
-!! correlations with the conservative potential temperature and
-!! the humidity conservative variable:
-!! ------- ------- -------
-!! A' Thv' = ETHETA A' Thl' + EMOIST A' Rnp'
-!!
-!!
-!! IMPLICIT ARGUMENTS
-!! ------------------
-!! Module MODD_CST : contains physical constants
-!!
-!! CST%XG : gravity constant
-!!
-!! Module MODD_CTURB: contains the set of constants for
-!! the turbulence scheme
-!!
-!! CSTURB%XCMFS,XCMFB : cts for the momentum flux
-!! CSTURB%XCSHF : ct for the sensible heat flux
-!! CSTURB%XCHF : ct for the moisture flux
-!! CSTURB%XCTV,CSTURB%XCHV : cts for the T and moisture variances
-!!
-!! Module MODD_PARAMETERS
-!!
-!! JPVEXT_TURB : number of vertical external points
-!!
-!!
-!! REFERENCE
-!! ---------
-!! Book 1 of documentation (Chapter: Turbulence)
-!!
-!! AUTHOR
-!! ------
-!! Joan Cuxart * INM and Meteo-France *
-!!
-!! MODIFICATIONS
-!! -------------
-!! Original August 19, 1994
-!! Modifications: February 14, 1995 (J.Cuxart and J.Stein)
-!! Doctorization and Optimization
-!! Modifications: March 21, 1995 (J.M. Carriere)
-!! Introduction of cloud water
-!! Modifications: June 14, 1995 (J.Cuxart and J. Stein)
-!! Phi3 and Psi3 at w-point + bug in the all
-!! or nothing condens.
-!! Modifications: Sept 15, 1995 (J.Cuxart and J. Stein)
-!! Change the DP computation at the ground
-!! Modifications: October 10, 1995 (J.Cuxart and J. Stein)
-!! Psi for scal var and LES tools
-!! Modifications: November 10, 1995 (J. Stein)
-!! change the surface relations
-!! Modifications: February 20, 1995 (J. Stein) optimization
-!! Modifications: May 21, 1996 (J. Stein)
-!! bug in the vertical flux of the V wind
-!! component for explicit computation
-!! Modifications: May 21, 1996 (N. wood)
-!! modify the computation of the vertical
-!! part or the surface tangential flux
-!! Modifications: May 21, 1996 (P. Jabouille)
-!! same modification in the Y direction
-!!
-!! Modifications: Sept 17, 1996 (J. Stein) change the moist case by using
-!! Pi instead of Piref + use Atheta and Amoist
-!!
-!! Modifications: Nov 24, 1997 (V. Masson) removes the DO loops
-!! Modifications: Mar 31, 1998 (V. Masson) splits the routine TURB_VER
-!! Nov 06, 2002 (V. Masson) LES budgets
-!! Feb 20, 2003 (JP Pinty) Add PFRAC_ICE
-!! July 2005 (S. Tomas, V. Masson)
-!! Add 3rd order moments and
-!! implicitation of PHI3, PSI3
-!! Oct.2009 (C.Lac) Introduction of different PTSTEP according to the
-!! advection schemes
-!! Feb. 2012 (Y. Seity) add possibility to run with
-!! reversed vertical levels
-!! 10/2012 (J.Escobar) Bypass PGI bug , redefine some allocatable array inplace of automatic
-!! 08/2014 (J.Escobar) Bypass PGI memory leak bug , replace IF statement with IF THEN ENDIF
-!! Modifications: July, 2015 (Wim de Rooy) switch for HARATU (Racmo turbulence scheme)
-!! Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
-!! JL Redelsperger 03/2021 : add Ocean LES case
-!!--------------------------------------------------------------------------
-!
-!* 0. DECLARATIONS
-! ------------
-!
-USE PARKIND1, ONLY : JPRB
-USE YOMHOOK , ONLY : LHOOK, DR_HOOK
-!
-USE MODD_CST, ONLY: CST_t
-USE MODD_CTURB, ONLY: CSTURB_t
-USE MODD_DIMPHYEX, ONLY: DIMPHYEX_t
-USE MODD_FIELD, ONLY: TFIELDDATA, TYPEREAL
-USE MODD_IO, ONLY: TFILEDATA
-USE MODD_PARAMETERS, ONLY: JPVEXT_TURB
-USE MODD_LES
-USE MODD_TURB_n, ONLY: TURB_t
-!
-USE MODE_EMOIST, ONLY: EMOIST
-USE MODE_ETHETA, ONLY: ETHETA
-USE MODI_GRADIENT_M
-USE MODI_GRADIENT_W
-USE MODE_GRADIENT_M_PHY, ONLY : GZ_M_W_PHY
-USE MODI_TURB
-USE MODE_TURB_VER_THERMO_FLUX, ONLY: TURB_VER_THERMO_FLUX
-USE MODE_TURB_VER_THERMO_CORR, ONLY: TURB_VER_THERMO_CORR
-USE MODE_TURB_VER_DYN_FLUX, ONLY: TURB_VER_DYN_FLUX
-USE MODE_TURB_VER_SV_FLUX, ONLY: TURB_VER_SV_FLUX
-USE MODE_TURB_VER_SV_CORR, ONLY: TURB_VER_SV_CORR
-USE MODI_LES_MEAN_SUBGRID
-USE MODE_SBL_DEPTH, ONLY: SBL_DEPTH
-USE MODI_SECOND_MNH
-!
-USE MODE_IO_FIELD_WRITE, only: IO_Field_write
-USE MODE_PRANDTL
-!
-IMPLICIT NONE
-!
-!* 0.1 declarations of arguments
-!
-!
-!
-TYPE(DIMPHYEX_t), INTENT(IN) :: D
-TYPE(CST_t), INTENT(IN) :: CST
-TYPE(CSTURB_t), INTENT(IN) :: CSTURB
-TYPE(TURB_t), INTENT(IN) :: TURBN
-INTEGER, INTENT(IN) :: KRR ! number of moist var.
-INTEGER, INTENT(IN) :: KRRL ! number of liquid water var.
-INTEGER, INTENT(IN) :: KRRI ! number of ice water var.
-INTEGER, INTENT(IN) :: KSV, KSV_LGBEG, KSV_LGEND ! number of scalar variables
-LOGICAL, INTENT(IN) :: OTURB_FLX ! switch to write the
- ! turbulent fluxes in the syncronous FM-file
-LOGICAL, INTENT(IN) :: OOCEAN ! switch for Ocean model version
-LOGICAL, INTENT(IN) :: ODEEPOC ! activates sfc forcing for ideal ocean deep conv
-LOGICAL, INTENT(IN) :: OHARAT !
-LOGICAL, INTENT(IN) :: OCOMPUTE_SRC ! flag to define dimensions of SIGS and SRCT variables
-LOGICAL, INTENT(IN) :: OFLAT ! Logical for zero ororography
-LOGICAL, INTENT(IN) :: OLES_CALL ! compute the LES diagnostics at current time-step
-LOGICAL, INTENT(IN) :: OCOUPLES ! switch to activate atmos-ocean LES version
-LOGICAL, INTENT(IN) :: OBLOWSNOW ! switch to activate pronostic blowing snow
-LOGICAL, INTENT(IN) :: ORMC01 ! switch for RMC01 lengths in SBL
-CHARACTER(LEN=6), INTENT(IN) :: HPROGRAM ! HPROGRAM is the program currently running
-LOGICAL, INTENT(IN) :: ONOMIXLG ! to use turbulence for lagrangian variables
-LOGICAL, INTENT(IN) :: O2D ! Logical for 2D model version
-CHARACTER(LEN=4), INTENT(IN) :: HTURBDIM ! dimensionality of the
- ! turbulence scheme
-CHARACTER(LEN=4), INTENT(IN) :: HTOM ! type of Third Order Moment
-REAL, INTENT(IN) :: PIMPL, PEXPL ! Coef. for temporal disc.
-REAL, INTENT(IN) :: PTSTEP ! timestep
-TYPE(TFILEDATA), INTENT(IN) :: TPFILE ! Output file
-!
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PDXX, PDYY, PDZZ, PDZX, PDZY
- ! Metric coefficients
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PDIRCOSZW ! Director Cosinus of the
- ! normal to the ground surface
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PZZ ! altitudes at flux points
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PCOSSLOPE ! cosinus of the angle
- ! between i and the slope vector
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PSINSLOPE ! sinus of the angle
- ! between i and the slope vector
-!
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PRHODJ ! dry density * grid volum
-! MFMOIST used in case of OHARATU
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: MFMOIST ! moist mass flux dual scheme
-
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PTHVREF ! ref. state Virtual
- ! Potential Temperature
-!
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PSFTHM,PSFRM ! surface fluxes at time
-REAL, DIMENSION(D%NIT,D%NJT,KSV), INTENT(IN) :: PSFSVM ! t - deltat
-!
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PSFTHP,PSFRP ! surface fluxes at time
-REAL, DIMENSION(D%NIT,D%NJT,KSV), INTENT(IN) :: PSFSVP ! t + deltat
-!
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PCDUEFF ! Cd * || u || at time t
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PTAU11M ! <uu> in the axes linked
- ! to the maximum slope direction and the surface normal and the binormal
- ! at time t - dt
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PTAU12M ! <uv> in the same axes
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PTAU33M ! <ww> in the same axes
-!
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PUM,PVM,PWM,PTHLM
- ! Wind and potential temperature at t-Delta t
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT,KRR), INTENT(IN) :: PRM ! Mixing ratios
- ! at t-Delta t
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT,KSV), INTENT(IN) :: PSVM ! scalar var. at t-Delta t
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PUSLOPEM ! wind component along the
- ! maximum slope direction
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PVSLOPEM ! wind component along the
- ! direction normal to the maximum slope one
-!
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PTKEM ! TKE at time t
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PLM ! Turb. mixing length
-! PLENGTHM PLENGTHH used in case of OHARATU
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PLENGTHM ! Turb. mixing length momentum
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PLENGTHH ! Turb. mixing length heat/moisture
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PLEPS ! dissipative length
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PLOCPEXNM ! Lv(T)/Cp/Exnref at time t-1
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PATHETA ! coefficients between
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PAMOIST ! s and Thetal and Rnp
-REAL, DIMENSION(MERGE(D%NIT,0,OCOMPUTE_SRC),&
- MERGE(D%NJT,0,OCOMPUTE_SRC),&
- MERGE(D%NKT,0,OCOMPUTE_SRC)), INTENT(IN) :: PSRCM ! normalized
- ! 2nd-order flux s'r'c/2Sigma_s2 at t-1 multiplied by Lambda_3
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PFRAC_ICE ! ri fraction of rc+ri
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PFWTH ! d(w'2th' )/dz
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PFWR ! d(w'2r' )/dz
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PFTH2 ! d(w'th'2 )/dz
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PFR2 ! d(w'r'2 )/dz
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(IN) :: PFTHR ! d(w'th'r')/dz
-REAL, DIMENSION(MERGE(D%NIT,0,HTOM=='TMO6'),&
- MERGE(D%NJT,0,HTOM=='TMO6')),INTENT(INOUT) :: PBL_DEPTH ! BL height for TOMS
-REAL, DIMENSION(MERGE(D%NIT,0,ORMC01),&
- MERGE(D%NJT,0,ORMC01)),INTENT(INOUT) :: PSBL_DEPTH ! SBL depth for RMC01
-REAL, DIMENSION(D%NIT,D%NJT), INTENT(IN) :: PLMO ! Monin-Obukhov length
-!
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(INOUT) :: PRUS, PRVS, PRWS, PRTHLS
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT,KRR), INTENT(INOUT) :: PRRS
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT,KSV), INTENT(INOUT) :: PRSVS
- ! cumulated sources for the prognostic variables
-!
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(OUT) :: PDP,PTP ! Dynamic and thermal
- ! TKE production terms
-REAL, DIMENSION(MERGE(D%NIT,0,OCOMPUTE_SRC),&
- MERGE(D%NJT,0,OCOMPUTE_SRC),&
- MERGE(D%NKT,0,OCOMPUTE_SRC)), INTENT(OUT) :: PSIGS
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(OUT) :: PWTH ! heat flux
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT), INTENT(OUT) :: PWRC ! cloud water flux
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT,KSV),INTENT(OUT) :: PWSV ! scalar flux
-!
-!
-!* 0.2 declaration of local variables
-!
-!JUAN BUG PGI
-!!$REAL, DIMENSION(SIZE(PTHLM,1),SIZE(PTHLM,2),SIZE(PTHLM,3)) :: &
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT) :: &
- ZBETA, & ! buoyancy coefficient
- ZSQRT_TKE,& ! sqrt(e)
- ZDTH_DZ, & ! d(th)/dz
- ZDR_DZ, & ! d(rt)/dz
- ZRED2TH3, & ! 3D Redeslperger number R*2_th
- ZRED2R3, & ! 3D Redeslperger number R*2_r
- ZRED2THR3,& ! 3D Redeslperger number R*2_thr
- ZBLL_O_E, & ! beta * Lk * Leps / tke
- ZETHETA, & ! Coefficient for theta in theta_v computation
- ZEMOIST, & ! Coefficient for r in theta_v computation
- ZREDTH1, & ! 1D Redelsperger number for Th
- ZREDR1, & ! 1D Redelsperger number for r
- ZPHI3, & ! phi3 Prandtl number
- ZPSI3, & ! psi3 Prandtl number for vapor
- ZD, & ! denominator in phi3 terms
- ZWTHV, & ! buoyancy flux
- ZWU, & ! (u'w')
- ZWV, & ! (v'w')
- ZTHLP, & ! guess of potential temperature due to vert. turbulent flux
- ZRP ! guess of total water due to vert. turbulent flux
-
-!!$REAL, DIMENSION(SIZE(PTHLM,1),SIZE(PTHLM,2),SIZE(PTHLM,3),KSV) :: &
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT,KSV) :: &
- ZPSI_SV, & ! Prandtl number for scalars
- ZREDS1, & ! 1D Redelsperger number R_sv
- ZRED2THS, & ! 3D Redelsperger number R*2_thsv
- ZRED2RS ! 3D Redelsperger number R*2_rsv
-REAL, DIMENSION(D%NIT,D%NJT,D%NKT) :: ZLM
-!
-LOGICAL :: GUSERV ! flag to use water vapor
-INTEGER :: IKB,IKE,IIB,IIE,IJB,IJE ! index value for the Beginning
- ! and the End of the physical domain for the mass points
-INTEGER :: JSV,JI,JJ,JK ! loop counter
-REAL :: ZTIME1
-REAL :: ZTIME2
-REAL(KIND=JPRB) :: ZHOOK_HANDLE
-TYPE(TFIELDDATA) :: TZFIELD
-!----------------------------------------------------------------------------
-!----------------------------------------------------------------------------
-!
-!* 1. PRELIMINARIES
-! -------------
-!
-IF (LHOOK) CALL DR_HOOK('TURB_VER',0,ZHOOK_HANDLE)
-!
-IKB=D%NKTB
-IKE=D%NKTE
-IIE=D%NIEC
-IIB=D%NIBC
-IJE=D%NJEC
-IJB=D%NJBC
-!
-!
-! 3D Redelsperger numbers
-!
-!
-CALL PRANDTL(D,CST,CSTURB,KRR,KSV,KRRI,OTURB_FLX, &
- HTURBDIM,OOCEAN,OHARAT,O2D,OCOMPUTE_SRC,&
- TPFILE, &
- PDXX,PDYY,PDZZ,PDZX,PDZY, &
- PTHVREF,PLOCPEXNM,PATHETA,PAMOIST, &
- PLM,PLEPS,PTKEM,PTHLM,PRM,PSVM,PSRCM, &
- ZREDTH1, ZREDR1, &
- ZRED2TH3, ZRED2R3, ZRED2THR3, &
- ZREDS1,ZRED2THS, ZRED2RS, &
- ZBLL_O_E, &
- ZETHETA, ZEMOIST )
-!
-! Buoyancy coefficient
-!
-IF (OOCEAN) THEN
- !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
- ZBETA(IIB:IIE,IJB:IJE,1:D%NKT) = CST%XG*CST%XALPHAOC
- !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
-ELSE
- !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
- ZBETA(IIB:IIE,IJB:IJE,1:D%NKT) = CST%XG/PTHVREF(IIB:IIE,IJB:IJE,1:D%NKT)
- !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
-END IF
-!
-! Square root of Tke
-!
-!$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
-ZSQRT_TKE(IIB:IIE,IJB:IJE,1:D%NKT) = SQRT(PTKEM(IIB:IIE,IJB:IJE,1:D%NKT))
-!$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
-!
-! gradients of mean quantities at previous time-step
-!
-CALL GZ_M_W_PHY(D,PTHLM,PDZZ,ZDTH_DZ)
-ZDR_DZ(:,:,:) = 0.
-IF (KRR>0) CALL GZ_M_W_PHY(D,PRM(:,:,:,1),PDZZ,ZDR_DZ)
-!
-!
-! Denominator factor in 3rd order terms
-!
-IF (.NOT. OHARAT) THEN
- !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
- ZD(IIB:IIE,IJB:IJE,1:D%NKT) = (1.+ZREDTH1(IIB:IIE,IJB:IJE,1:D%NKT)+ZREDR1(IIB:IIE,IJB:IJE,1:D%NKT)) * (1.+0.5*(ZREDTH1(IIB:IIE,IJB:IJE,1:D%NKT)+ZREDR1(IIB:IIE,IJB:IJE,1:D%NKT)))
- !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
-ELSE
- !$mnh_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
- ZD(IIB:IIE,IJB:IJE,1:D%NKT) = 1.
- !$mnh_end_expand_array(JI=IIB:IIE,JJ=IJB:IJE,JK=1:D%NKT)
-ENDIF
-!
-! Phi3 and Psi3 Prandtl numbers
-!
-GUSERV = KRR/=0
-!
-ZPHI3 = PHI3(D,CSTURB,ZREDTH1,ZREDR1,ZRED2TH3,ZRED2R3,ZRED2THR3,HTURBDIM,GUSERV)
-IF(KRR/=0) &
-ZPSI3 = PSI3(D,CSTURB,ZREDR1,ZREDTH1,ZRED2R3,ZRED2TH3,ZRED2THR3,HTURBDIM,GUSERV)
-!
-! Prandtl numbers for scalars
-!
-ZPSI_SV = PSI_SV(D,CSTURB,KSV,ZREDTH1,ZREDR1,ZREDS1,ZRED2THS,ZRED2RS,ZPHI3,ZPSI3)
-!
-! LES diagnostics
-!
-IF (OLES_CALL) THEN
- CALL SECOND_MNH(ZTIME1)
- CALL LES_MEAN_SUBGRID(ZPHI3,X_LES_SUBGRID_PHI3)
- IF(KRR/=0) THEN
- CALL LES_MEAN_SUBGRID(ZPSI3,X_LES_SUBGRID_PSI3)
- END IF
- CALL SECOND_MNH(ZTIME2)
- XTIME_LES = XTIME_LES + ZTIME2 - ZTIME1
-END IF
-!----------------------------------------------------------------------------
-!
-!
-!* 2. SOURCES OF CONSERVATIVE POTENTIAL TEMPERATURE AND
-! PARTIAL THERMAL PRODUCTION
-! ---------------------------------------------------------------
-!
-!* 3. SOURCES OF CONSERVATIVE AND CLOUD MIXING RATIO AND
-! COMPLETE THERMAL PRODUCTION
-! ------------------------------------------------------
-!
-!* 4. TURBULENT CORRELATIONS : <w Rc>, <THl THl>, <THl Rnp>, <Rnp Rnp>
-! ----------------------------------------------------------------
-!
-
-IF (OHARAT) THEN
- ZLM=PLENGTHH
-ELSE
- ZLM=PLM
-ENDIF
-!
- CALL TURB_VER_THERMO_FLUX(D,CST,CSTURB,TURBN, &
- KRR,KRRL,KRRI,KSV, &
- OTURB_FLX,HTURBDIM,HTOM,OOCEAN,ODEEPOC,OHARAT,&
- OCOUPLES,OLES_CALL,OCOMPUTE_SRC, &
- PIMPL,PEXPL,PTSTEP,HPROGRAM,TPFILE, &
- PDXX,PDYY,PDZZ,PDZX,PDZY,PDIRCOSZW,PZZ, &
- PRHODJ,PTHVREF, &
- PSFTHM,PSFRM,PSFTHP,PSFRP, &
- PWM,PTHLM,PRM,PSVM, &
- PTKEM,ZLM,PLEPS, &
- PLOCPEXNM,PATHETA,PAMOIST,PSRCM,PFRAC_ICE, &
- ZBETA, ZSQRT_TKE, ZDTH_DZ, ZDR_DZ, ZRED2TH3, &
- ZRED2R3, ZRED2THR3, ZBLL_O_E, ZETHETA, &
- ZEMOIST, ZREDTH1, ZREDR1, ZPHI3, ZPSI3, ZD, &
- PFWTH,PFWR,PFTH2,PFR2,PFTHR, &
- MFMOIST,PBL_DEPTH,ZWTHV, &
- PRTHLS,PRRS,ZTHLP,ZRP,PTP,PWTH,PWRC )
-!
- CALL TURB_VER_THERMO_CORR(D,CST,CSTURB, &
- KRR,KRRL,KRRI,KSV, &
- OTURB_FLX,HTURBDIM,HTOM, OHARAT,OCOMPUTE_SRC, &
- OCOUPLES,OLES_CALL, &
- PIMPL,PEXPL,TPFILE, &
- PDXX,PDYY,PDZZ,PDZX,PDZY,PDIRCOSZW, &
- PRHODJ,PTHVREF, &
- PSFTHM,PSFRM,PSFTHP,PSFRP, &
- PWM,PTHLM,PRM,PSVM, &
- PTKEM,ZLM,PLEPS, &
- PLOCPEXNM,PATHETA,PAMOIST, &
- ZBETA, ZSQRT_TKE, ZDTH_DZ, ZDR_DZ, ZRED2TH3, &
- ZRED2R3, ZRED2THR3, ZBLL_O_E, ZETHETA, &
- ZEMOIST, ZREDTH1, ZREDR1, ZPHI3, ZPSI3, ZD, &
- PFWTH,PFWR,PFTH2,PFR2,PFTHR, &
- ZTHLP,ZRP,MFMOIST,PSIGS )
-!
-!----------------------------------------------------------------------------
-!
-!
-!
-!* 5. SOURCES OF U,W WIND COMPONENTS AND PARTIAL DYNAMIC PRODUCTION
-! -------------------------------------------------------------
-!
-!* 6. SOURCES OF V,W WIND COMPONENTS AND COMPLETE 1D DYNAMIC PRODUCTION
-! -----------------------------------------------------------------
-!
-!* 7. DIAGNOSTIC COMPUTATION OF THE 1D <W W> VARIANCE
-! -----------------------------------------------
-!
-!
-IF (OHARAT) ZLM=PLENGTHM
-!
-CALL TURB_VER_DYN_FLUX(D,CST,CSTURB,TURBN,KSV,O2D,OFLAT, &
- OTURB_FLX,KRR,OOCEAN,OHARAT,OCOUPLES,OLES_CALL,&
- HTURBDIM,PIMPL,PEXPL,PTSTEP,TPFILE, &
- PDXX,PDYY,PDZZ,PDZX,PDZY,PDIRCOSZW,PZZ, &
- PCOSSLOPE,PSINSLOPE, &
- PRHODJ, &
- PCDUEFF,PTAU11M,PTAU12M,PTAU33M, &
- PTHLM,PRM,PSVM,PUM,PVM,PWM,PUSLOPEM,PVSLOPEM, &
- PTKEM,ZLM,MFMOIST,ZWU,ZWV, &
- PRUS,PRVS,PRWS, &
- PDP,PTP )
-!
-!----------------------------------------------------------------------------
-!
-!
-!* 8. SOURCES OF PASSIVE SCALAR VARIABLES
-! -----------------------------------
-!
-IF (OHARAT) ZLM=PLENGTHH
-!
-IF (KSV>0) &
-CALL TURB_VER_SV_FLUX(D,CST,CSTURB,ONOMIXLG, &
- KSV,KSV_LGBEG,KSV_LGEND, &
- OTURB_FLX,HTURBDIM,OHARAT,OBLOWSNOW,OLES_CALL,&
- PIMPL,PEXPL,PTSTEP, &
- TPFILE, &
- PDZZ,PDIRCOSZW, &
- PRHODJ,PWM, &
- PSFSVM,PSFSVP, &
- PSVM, &
- PTKEM,ZLM,MFMOIST,ZPSI_SV, &
- PRSVS,PWSV )
-!
-!
-IF (KSV>0 .AND. OLES_CALL) &
-CALL TURB_VER_SV_CORR(D,CST,CSTURB,KRR,KRRL,KRRI,OOCEAN, &
- PDZZ,KSV,KSV_LGBEG,KSV_LGEND,ONOMIXLG, &
- OBLOWSNOW,OLES_CALL,OCOMPUTE_SRC, &
- PTHLM,PRM,PTHVREF, &
- PLOCPEXNM,PATHETA,PAMOIST,PSRCM,ZPHI3,ZPSI3, &
- PWM,PSVM, &
- PTKEM,ZLM,PLEPS,ZPSI_SV )
-!
-!
-!----------------------------------------------------------------------------
-!
-!* 9. DIAGNOSTIC OF Surface Boundary Layer Depth
-! ------------------------------------------
-!
-IF (ORMC01) CALL SBL_DEPTH(IKB,IKE,PZZ,ZWU,ZWV,ZWTHV,PLMO,PSBL_DEPTH)
-!
-!----------------------------------------------------------------------------
-!
-!
-!* 10. PRINTS
-! ------
-!
-!
-IF ( OTURB_FLX .AND. TPFILE%LOPENED .AND. .NOT. OHARAT) THEN
-!
-! stores the Turbulent Prandtl number
-!
- TZFIELD%CMNHNAME = 'PHI3'
- TZFIELD%CSTDNAME = ''
- TZFIELD%CLONGNAME = 'PHI3'
- TZFIELD%CUNITS = '1'
- TZFIELD%CDIR = 'XY'
- TZFIELD%CCOMMENT = 'Turbulent Prandtl number'
- TZFIELD%NGRID = 4
- TZFIELD%NTYPE = TYPEREAL
- TZFIELD%NDIMS = 3
- TZFIELD%LTIMEDEP = .TRUE.
- CALL IO_Field_write(TPFILE,TZFIELD,ZPHI3)
-!
-! stores the Turbulent Schmidt number
-!
- TZFIELD%CMNHNAME = 'PSI3'
- TZFIELD%CSTDNAME = ''
- TZFIELD%CLONGNAME = 'PSI3'
- TZFIELD%CUNITS = '1'
- TZFIELD%CDIR = 'XY'
- TZFIELD%CCOMMENT = 'Turbulent Schmidt number'
- TZFIELD%NGRID = 4
- TZFIELD%NTYPE = TYPEREAL
- TZFIELD%NDIMS = 3
- TZFIELD%LTIMEDEP = .TRUE.
- CALL IO_Field_write(TPFILE,TZFIELD,ZPSI3)
-!
-!
-! stores the Turbulent Schmidt number for the scalar variables
-!
- TZFIELD%CSTDNAME = ''
- TZFIELD%CUNITS = '1'
- TZFIELD%CDIR = 'XY'
- TZFIELD%NGRID = 4
- TZFIELD%NTYPE = TYPEREAL
- TZFIELD%NDIMS = 3
- TZFIELD%LTIMEDEP = .TRUE.
- DO JSV=1,KSV
- WRITE(TZFIELD%CMNHNAME, '("PSI_SV_",I3.3)') JSV
- TZFIELD%CLONGNAME = TRIM(TZFIELD%CMNHNAME)
- TZFIELD%CCOMMENT = 'X_Y_Z_'//TRIM(TZFIELD%CMNHNAME)
- CALL IO_Field_write(TPFILE,TZFIELD,ZPSI_SV(:,:,:,JSV))
- END DO
-!
-END IF
-!
-!
-!----------------------------------------------------------------------------
-IF (LHOOK) CALL DR_HOOK('TURB_VER',1,ZHOOK_HANDLE)
-END SUBROUTINE TURB_VER
-END MODULE MODE_TURB_VER