Changeset 2945 for LMDZ5/trunk/libf/phylmd/fisrtilp.F90

Timestamp:

Jul 12, 2017, 4:20:24 PM (7 years ago)

Author:

jbmadeleine

Message:

• Added a new output called rneblsvol which is the cloud fraction by volume

computed in the thermals (see cloudth_vert in cloudth_mod.F90)

• Added an option called iflag_rain_incloud_vol that computes the conversion

of cloud water to rain using the cloud fraction by volume instead of the cloud
fraction by area, which is larger and otherwise erroneously reduces the in-cloud
water content; iflag_rain_incloud_vol can only be used for iflag_cloudth_vert>=3

File:

• LMDZ5/trunk/libf/phylmd/fisrtilp.F90 (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/fisrtilp.F90

@@ -18 +18 @@
   USE ioipsl_getin_p_mod, ONLY : getin_p
   USE phys_local_var_mod, ONLY: ql_seri,qs_seri
+  USE phys_local_var_mod, ONLY: rneblsvol
   ! flag to include modifications to ensure energy conservation (if flag >0)
   USE add_phys_tend_mod, only : fl_cor_ebil 
@@ -132 +133 @@
   INTEGER ncoreczq  
   REAL ctot(klon,klev)
+  REAL ctot_vol(klon,klev)
   REAL zqs(klon), zdqs(klon), zdelta, zcor, zcvm5  
   REAL zdqsdT_raw(klon)
@@ -217 +219 @@
 !  ice_thermo = iflag_ice_thermo .GE. 1
 
-itap=itap+1
-znebprecip(:)=0.
-
-   ice_thermo = (iflag_ice_thermo .EQ. 1).OR.(iflag_ice_thermo .GE. 3)
+  itap=itap+1
+  znebprecip(:)=0.
+
+  ice_thermo = (iflag_ice_thermo .EQ. 1).OR.(iflag_ice_thermo .GE. 3)
   zdelq=0.0
+  ctot_vol(1:klon,1:klev)=0.0
+  rneblsvol(1:klon,1:klev)=0.0
 
   if (prt_level>9)write(lunout,*)'NUAGES4 A. JAM'
@@ -727 +731 @@
                call cloudth_v3(klon,klev,k,ztv, &
                    zq,zqta,fraca, &
-                   qcloud,ctot,zpspsk,paprs,ztla,zthl, &
+                   qcloud,ctot,ctot_vol,zpspsk,paprs,ztla,zthl, &
                    ratqs,zqs,t)
               endif
               do i=1,klon
                  rneb(i,k)=ctot(i,k)
+                 rneblsvol(i,k)=ctot_vol(i,k)
                  zqn(i)=qcloud(i)
               enddo
@@ -1148 +1153 @@
                          *fallvs(zrhol(i)) * zfice(i)
                  endif
-                 zchau    = zct   *dtime/REAL(ninter) * zoliq(i) &
-                      *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl   )**2)) *(1.-zfice(i))
+
+                 ! si l'heterogeneite verticale est active, on utilise 
+                 ! la fraction volumique "vraie" plutot que la fraction
+                 ! surfacique modifiee, qui est plus grande et reduit
+                 ! sinon l'eau in-cloud de facon artificielle
+                 if ((iflag_cloudth_vert>=3).AND.(iflag_rain_incloud_vol==1)) then
+                    zchau    = zct   *dtime/REAL(ninter) * zoliq(i) &
+                         *(1.0-EXP(-(zoliq(i)/ctot_vol(i,k)/zcl   )**2)) *(1.-zfice(i))
+                 else
+                    zchau    = zct   *dtime/REAL(ninter) * zoliq(i) &
+                         *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl   )**2)) *(1.-zfice(i))
+                 endif
 !AJ<
                  IF (.NOT. ice_thermo) THEN