Changeset 2073 for trunk/LMDZ.GENERIC

Timestamp:

Jan 18, 2019, 10:53:39 PM (6 years ago)

Author:

aslmd

Message:

in largescale enforce type conversion as is done in the other instance of Lcpdqsat_water. corrects an error with picky compilers (e.g. gfortran)

File:

• trunk/LMDZ.GENERIC/libf/phystd/largescale.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/largescale.F90

@@ -151 +151 @@
 !           iterative process to stabilize the scheme when large water amounts JL12
             zcond(i) = 0.0d0
-            Do nn=1,nitermax  
-               call Lcpdqsat_water(zt(i),local_p,psat_tmp,zqs(i),zdqs,dlnpsat_tmp)
+            Do nn=1,nitermax 
+               ! use zqs_temp and not zqs(i) to force type conversion
+               ! -- might not be a good solution, actually
+               ! but this is compliant with "complete cloud cover" case above 
+               call Lcpdqsat_water(zt(i),local_p,psat_tmp,zqs_temp,zdqs,dlnpsat_tmp)
                zcond_iter = MAX(0.0d0,alpha*(zx_q(i)-zqs(i))/(1.d0+zdqs))          
                   !zcond always postive! cannot evaporate clouds!