Changeset 6007 for LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

Timestamp:

Jan 7, 2026, 10:29:50 PM (2 weeks ago)

Author:

evignon

Message:

modif routines ratqs et lscp, fin des modifs avant le regroupement
dans un call unique

File:

• LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

@@ -108 +108 @@
     USE lmdz_thermcell_dtke, ONLY : thermcell_dtke
     USE lmdz_blowing_snow_ini, ONLY : blowing_snow_ini , qbst_bs
-    USE lmdz_lscp_ini, ONLY : lscp_ini
+    USE lmdz_lscp_ini, ONLY : lscp_ini, ratqsbas
     USE lmdz_lscp_subgridvarq, ONLY : ratqs_main, ratqs_main_first
     USE lmdz_cloud_optics_prop_ini, ONLY : cloud_optics_prop_ini
@@ -4222 +4222 @@
           enddo
        ENDIF
-
-       ! =================================================================== c
-       ! Calcul des proprietes des nuages convectifs
-       !
-
-       !   calcul des proprietes des nuages convectifs
-       clwcon0(:,:)=fact_cldcon*clwcon0(:,:)
-       IF (iflag_cld_cv == 0) THEN
-          CALL clouds_gno &
-               (klon,klev,q_seri,zqsat,clwcon0,ptconv,ratqsc,rnebcon0)
-       ELSE
-          CALL clouds_bigauss &
-               (klon,klev,q_seri,zqsat,qtc_cv,sigt_cv,ptconv,ratqsc,rnebcon0)
-       ENDIF
-
-
-       ! =================================================================== c
+         
+          ! =========================================================================
+          ! Calculation of deep convective clouds properties
+          ! that is, the convective ratqs (raqtsc) and the convective cloud fraction
+          ! =========================================================================
+
+          clwcon0(:,:)=fact_cldcon*clwcon0(:,:)
+          IF (iflag_cld_cv == 0) THEN
+             CALL clouds_gno &
+                  (klon,klev,q_seri,zqsat,clwcon0,ptconv,ratqsc,rnebcon0)
+          ELSE
+             CALL clouds_bigauss &
+                  (klon,klev,q_seri,zqsat,qtc_cv,sigt_cv,ptconv,ratqsc,rnebcon0)
+          ENDIF
+
+          ! old-style convective ratqs computation as a function of q(z=0)-q / q
+          ! the ratqsc computed by clouds_gno is replaced
+          if (iflag_cld_th.eq.1) then
+            do k=1,klev
+               do i=1,klon
+                  if (ptconv(i,k)) then
+                      ratqsc(i,k)=ratqsbas &
+                      +fact_cldcon*(q_seri(i,1)-q_seri(i,k))/q_seri(i,k)
+                  else
+                      ratqsc(i,k)=0.
+                  endif
+               enddo
+             enddo
+          ! through log-normal distribution inversion
+          else if (iflag_cld_th.eq.4) then
+             ptconvth(:,:)=.false.
+             ratqsc(:,:)=0.
+             call clouds_gno &
+             (klon,klev,q_seri,zqsat,clwcon0th,ptconvth,ratqsc,rnebcon0th)
+          endif
+
+          if (iflag_cld_th.eq.1 .or.iflag_cld_th.eq.2.or.iflag_cld_th.eq.4) then
+          ! We add a small constant value to ratqsc*2 to account for small-scale fluctuations
+            do k=1,klev
+               do i=1,klon
+                  if ((fm_therm(i,k)>1.e-10)) then
+                     ratqsc(i,k)=sqrt(ratqsc(i,k)**2+0.05**2)
+                  endif
+               enddo
+            enddo
+          endif
+
+          ! ======================================================================= 
+
 
        DO i = 1, klon
@@ -5096 +5128 @@
     CALL ratqs_main_first(klon, cell_area)
     CALL ratqs_main(klon,klev,nbsrf,is_ter,is_lic,  &
-         iflag_cld_th,fact_cldcon,pdtphys,          &
+         iflag_cld_th,pdtphys,          &
          pctsrf,s_pblh,zstd, wake_s, wake_deltaq,   &
-         ptconv,clwcon0th, rnebcon0th,     &
          paprs,pplay,t_seri,q_seri, &
          qtc_cv, sigt_cv,detrain_cv,fm_cv,fqd,fqcomp,sigd,zqsat, &
          fm_therm,entr_therm,detr_therm,cell_area, &
-         ratqs,ratqsc,ratqs_inter_,sigma_qtherm)
+         ratqss,ratqs,ratqs_inter_,sigma_qtherm)
 
     !