Changeset 3570 for LMDZ6/trunk/libf/phylmd/cloudth_mod.F90

Timestamp:

Sep 27, 2019, 11:17:19 PM (5 years ago)

Author:

fhourdin

Message:

Sortie de nouveaux parametres pour controler la distribution des nuages.

File:

• LMDZ6/trunk/libf/phylmd/cloudth_mod.F90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/cloudth_mod.F90

@@ -885 +885 @@
       REAL,SAVE :: sigma1s_factor=1.1
       REAL,SAVE :: sigma1s_power=0.6
+      REAL,SAVE :: sigma2s_factor=0.09
+      REAL,SAVE :: sigma2s_power=0.5
       REAL,SAVE :: cloudth_ratqsmin=-1.
-      !$OMP THREADPRIVATE(sigma1s_factor,sigma1s_power,cloudth_ratqsmin)
+      !$OMP THREADPRIVATE(sigma1s_factor,sigma1s_power,sigma2s_factor,sigma2s_power,cloudth_ratqsmin)
       INTEGER, SAVE :: iflag_cloudth_vert_noratqs=0
       !$OMP THREADPRIVATE(iflag_cloudth_vert_noratqs)
@@ -951 +953 @@
         CALL getin_p('cloudth_sigma1s_power',sigma1s_power)
         WRITE(*,*) 'cloudth_sigma1s_power = ', sigma1s_power
+        ! Factor used in the calculation of sigma2s
+        CALL getin_p('cloudth_sigma2s_factor',sigma2s_factor)
+        WRITE(*,*) 'cloudth_sigma2s_factor = ', sigma2s_factor
+        ! Power used in the calculation of sigma2s
+        CALL getin_p('cloudth_sigma2s_power',sigma2s_power)
+        WRITE(*,*) 'cloudth_sigma2s_power = ', sigma2s_power
         ! Minimum value for the environmental air subgrid water distrib
         CALL getin_p('cloudth_ratqsmin',cloudth_ratqsmin)
@@ -1019 +1027 @@
       ENDIF
       sigma1s = sigma1s_fraca + sigma1s_ratqs
-      sigma2s=(0.09*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**0.5))+0.002*zqta(ind1,ind2)
+      sigma2s=(sigma2s_factor*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2)
 !      tests
 !      sigma1s=(0.92**0.5)*(fraca(ind1,ind2)**0.5)/(1-fraca(ind1,ind2))*((sth-senv)**2)**0.5+ratqs(ind1,ind2)*po(ind1)