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

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

File:

• trunk/LMDZ.MARS/libf/phymars/co2cloud.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)