Changeset 5995 for LMDZ6/trunk/libf/phylmd/lmdz_ratqs_multi.f90

Timestamp:

Jan 5, 2026, 5:09:54 PM (8 days ago)

Author:

evignon

Message:

nettoyage d'une option obsolete de ratqs_multi

File:

• LMDZ6/trunk/libf/phylmd/lmdz_ratqs_multi.f90 (modified) (5 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_ratqs_multi.f90

@@ -38 +38 @@
 
 !========================================================================
-! L. d'Alen??on, 25/02/2021
+! L. d'Alencon, 25/02/2021
 ! Cette subroutine calcule une valeur de ratqsbas interactive  
 ! Elle est appel??e par la subroutine ratqs lorsque iflag_ratqs = 11. 
@@ -84 +84 @@
       enddo
 !-------------------------------------------------------------------------
-!  Caclul du terme de d??trainement de la variance pour les thermiques
+!  Caclul du terme de detrainement de la variance pour les thermiques
 !-------------------------------------------------------------------------
 
@@ -102 +102 @@
       detr0(:,:) = detr_therm(:,:) 
 
-! calcul du carr?? de l'humidit?? sp??cifique et circulation dans les thermiques
+! calcul du carre de l'humidite specifique et circulation dans les thermiques
       po(:,:) = q_seri(:,:)
       call thermcell_dq(klon,klev,dqimpl,pdtphys,fm0,entr0,zmasse,  &
@@ -115 +115 @@
 
 
-! variance de l'humidit?? sp??cifique totale dans les thermiques     
+! variance de l'humidit?? specifique totale dans les thermiques     
       do k=1,klev
          do i=1,klon      
@@ -364 +364 @@
 !=============================================
 
-SUBROUTINE ratqs_tke(klon,klev,pdtphys,temp,q,qsat,pplay,paprs,omega,tke,tke_dissip,lmix,wprime,ratqs_tke_)
-
-! References:
-!
-! Etienne Vignon: effect of subgrid turbulence on ratqs
-!
-! Field, P.R., Hill, A., Furtado, K., Korolev, A., 2014b. Mixed-phase clouds in a turbulent environment. Part
-! 2: analytic treatment. Q. J. R. Meteorol. Soc. 21, 2651???2663. https://doi.org/10.1002/qj.2175.
-! 
-! Furtado, K., Field, P.R., Boutle, I.A., Morcrette, C.R., Wilkinson, J., 2016. A physically-based, subgrid
-! parametrization for the production and maintenance of mixed-phase clouds in a general circulation
-! model. J. Atmos. Sci. 73, 279???291. https://doi.org/10.1175/JAS-D-15-0021.
-
-USE lmdz_ratqs_ini, ONLY : RG,RV,RD,RCPD,RLSTT,RLVTT,RTT
-
-IMPLICIT NONE
-
-! INPUTS
-
-INTEGER, INTENT(IN) :: klon                             ! number of horizontal grid points
-INTEGER, INTENT(IN) :: klev                             ! number of vertical layers
-REAL, INTENT(IN) :: pdtphys                             ! physics time step [s]
-REAL, DIMENSION(klon,klev), INTENT(IN) :: temp          ! air temperature [K]
-REAL, DIMENSION(klon,klev), INTENT(IN) :: q             ! specific humidity [kg/kg]
-REAL, DIMENSION(klon,klev), INTENT(IN) :: qsat          ! saturation specific humidity [kg/kg]
-REAL, DIMENSION(klon,klev), INTENT(IN) :: pplay         ! air pressure, layer's center [Pa]
-REAL, DIMENSION(klon,klev+1), INTENT(IN) :: paprs       ! air pressure, lower inteface [Pa]
-REAL, DIMENSION(klon,klev), INTENT(IN) :: omega       ! air pressure, lower inteface [Pa]
-REAL, DIMENSION(klon,klev+1), INTENT(IN) :: tke         ! Turbulent Kinetic Energy [m2/s2]
-REAL, DIMENSION(klon,klev+1), INTENT(IN) :: tke_dissip  ! Turbulent Kinetic Energy Dissipation rate [m2/s3]
-REAL, DIMENSION(klon,klev+1), INTENT(IN) :: lmix  ! Turbulent mixing length
-REAL, DIMENSION(klon,klev+1), INTENT(IN) :: wprime      ! Turbulent vertical velocity scale [m/s]
-
-! OUTPUTS
-
-REAL, DIMENSION(klon,klev), INTENT(out) :: ratqs_tke_  ! ratsq profile due to subgrid TKE
-
-! LOCAL
-INTEGER :: i, k
-REAL :: AA, DD, NW, AAprime, VARLOG,rho,Lvs,taue,lhomo,dissmin,maxvarlog
-REAL, DIMENSION(klon,klev) :: sigmaw,w
-REAL, PARAMETER :: C0=10.0
-REAL, PARAMETER :: lmin=0.001
-REAL, PARAMETER :: ratqsmin=1E-6
-REAL, PARAMETER :: ratqsmax=0.5
-
-
-! Calculation of large scale and turbulent vertical velocities
-!---------------------------------------------------------------
-
-DO k=1,klev
-    DO i=1,klon
-        rho=pplay(i,k)/temp(i,k)/RD
-        w(i,k)=-rho*RG*omega(i,k)
-        sigmaw(i,k)=0.5*(wprime(i,k+1)+wprime(i,k)) ! turbulent vertical velocity at the middle of model layers. 
-    END DO
-END DO
-
-! Calculation of ratqs
-!---------------------------------------------------------------
-ratqs_tke_(:,1)=ratqsmin ! set to a very low value to avoid division by 0 in order parts
-                        ! of the code
-DO k=2,klev ! we start from second model level since TKE is not defined at k=1
-    DO i=1,klon
-
-       IF (temp(i,k) .LT. RTT) THEN
-           Lvs=RLSTT
-       ELSE
-           Lvs=RLVTT
-       ENDIF
-       dissmin=0.01*(0.5*(tke(i,k)+tke(i,k+1))/pdtphys)
-       maxvarlog=LOG(1.0+ratqsmax**2)! to prevent ratqs from exceeding an arbitrary threshold value 
-       AA=RG*(Lvs/(RCPD*temp(i,k)*temp(i,k)*RV) - 1./(RD*temp(i,k)))
-       lhomo=MAX(0.5*(lmix(i,k)+lmix(i,k+1)),lmin)
-       taue=(lhomo*lhomo/MAX(0.5*(tke_dissip(i,k)+tke_dissip(i,k+1)),dissmin))**(1./3) ! Fields et al. 2014
-       DD=1.0/taue
-       NW=(sigmaw(i,k)**2)*SQRT(2./(C0*MAX(0.5*(tke_dissip(i,k)+tke_dissip(i,k+1)),dissmin)))
-       AAprime=AA*NW
-       VARLOG=AAprime/2./DD
-       VARLOG=MIN(VARLOG,maxvarlog)
-       ratqs_tke_(i,k)=SQRT(MAX(EXP(VARLOG)-1.0,ratqsmin))
-       END DO
-END DO
-END SUBROUTINE ratqs_tke
-
 END MODULE lmdz_ratqs_multi