Changeset 2261 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Mar 16, 2020, 2:22:13 PM (5 years ago)

Author:

emillour

Message:

Mars GCM:
Revert co2cloud.F and improvedco2clouds_mod.F to what they were in r2257 (they were accidentaly reverted to earlier version by r2260).
EM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• co2cloud.F (modified) (1 diff)
• improvedco2clouds_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/co2cloud.F

@@ -840 +840 @@
                 No=Nccnco2*tauscaling(ig)
                 Rn=-dlog(riceco2(ig,l))
-                n_derf = derf( (rb_cldco2(1)+Rn) *dev2)
+                n_derf = erf( (rb_cldco2(1)+Rn) *dev2)
                 Qext1bins2(ig,l)=0.
                 do i = 1, nbinco2_cld 
                  n_aer(i) = -0.5 * No * n_derf !! this ith previously computed
-                 n_derf = derf((rb_cldco2(i+1)+Rn) *dev2)
+                 n_derf = erf((rb_cldco2(i+1)+Rn) *dev2)
                  n_aer(i) = n_aer(i) + 0.5 * No * n_derf
                  Qext1bins2(ig,l)=Qext1bins2(ig,l)+Qext1bins(i)*n_aer(i)

trunk/LMDZ.MARS/libf/phymars/improvedco2clouds_mod.F

@@ -552 +552 @@
 c Ic_rice Mass transfer rate (kg/s) for a rice particle >0 si croissance ! 
              
-              if (isnan(Ic_rice) .or. Ic_rice .eq. 0.) then 
+              if ( (Ic_rice.ne.Ic_rice) ! will be true if it is Nan  
+     &               .or. (Ic_rice.eq.0.) ) then
                  Ic_rice=0.
                  subpdtcloudco2(ig,l)=-sum_subpdt(ig,l)