Changeset 2685

Timestamp:

Oct 26, 2016, 5:12:55 PM (8 years ago)

Author:

musat

Message:

Go back to the 2016' summer versions of cloudth.F90 for the
6.0.6. tuning.

Replace cloudth.F90 by a module containing two versions
of cloudth.F90 and cloudth_vert.F90 routines.

• cloudth.F90 and cloudth_vert.F90 correspond to summer of 2016.
• cloudth_v3.F90 and cloudth_vert_v3.F90 correspond to autumn 2016.

Note that the older versions are called if iflag_cloudth_vert=2
and the newer versions are called if iflag_cloudth_vert=3.

This can be usefull if we want to make sensitivity tests to the newer
versions without recompiling the code.

Location:

LMDZ5/branches/testing/libf/phylmd

Files:

• cloudth_mod.F90 (added)
• fisrtilp.F90 (modified) (2 diffs)

LMDZ5/branches/testing/libf/phylmd/fisrtilp.F90

@@ -15 +15 @@
   USE icefrac_lsc_mod ! compute ice fraction (JBM 3/14)
   USE print_control_mod, ONLY: prt_level, lunout
+  USE cloudth_mod 
   IMPLICIT none
   !======================================================================
@@ -574 +575 @@
            if (iflag_cld_th>=5) then
 
-              call cloudth(klon,klev,k,ztv, &
+              if (iflag_cloudth_vert.EQ.2) then
+               call cloudth(klon,klev,k,ztv, &
                    zq,zqta,fraca, &
                    qcloud,ctot,zpspsk,paprs,ztla,zthl, &
                    ratqs,zqs,t)
-
+              elseif (iflag_cloudth_vert.EQ.3) then
+               call cloudth_v3(klon,klev,k,ztv, &
+                   zq,zqta,fraca, &
+                   qcloud,ctot,zpspsk,paprs,ztla,zthl, &
+                   ratqs,zqs,t)
+              endif
               do i=1,klon
                  rneb(i,k)=ctot(i,k)