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

Timestamp:

Sep 5, 2019, 4:30:23 PM (5 years ago)

Author:

aslmd

Message:

in improvedCO2clouds_mod, changed dlog to log, dexp to exp, and a dble instance so that the model can work in simple precision e.g. mesoscale

File:

• trunk/LMDZ.MARS/libf/phymars/improvedCO2clouds_mod.F (modified) (8 diffs)

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

@@ -198 +198 @@
 c       dr_cld is the width of each rad_cldco2 bin.
 
-        vrat_cld = dlog(rmax_cld/rmin_cld) / float(nbinco2_cld-1) *3.
-        vrat_cld = dexp(vrat_cld)
+        vrat_cld = log(rmax_cld/rmin_cld) / float(nbinco2_cld-1) *3.
+        vrat_cld = exp(vrat_cld)
         rb_cldco2(1)  = rbmin_cld
         rad_cldco2(1) = rmin_cld
@@ -226 +226 @@
         print*,'-----------------------------------'
         do i=1,nbinco2_cld+1
-            rb_cldco2(i) = dlog(rb_cldco2(i))  !! we save that so that it is not computed
+            rb_cldco2(i) = log(rb_cldco2(i))  !! we save that so that it is not computed
                                          !! at each timestep and gridpoint
         enddo
@@ -402 +402 @@
 
               Mo =4.*pi*rho_dust*No*rdust(ig,l)**(3.)
-     &           *dexp(9.*nuiceco2_ref/2.)/3. ! in Madeleine et al 2011
+     &           *exp(9.*nuiceco2_ref/2.)/3. ! in Madeleine et al 2011
      
               No = zq(ig,l,igcm_dust_number)* tauscaling(ig)+1.e-30
@@ -410 +410 @@
 
               Rn = rdust(ig,l)
-              Rn = -dlog(Rn) 
+              Rn = -log(Rn) 
               Rm = Rn - 3. * sigma_iceco2*sigma_iceco2  
               n_derf = derf( (rb_cldco2(1)+Rn) *dev2)
@@ -448 +448 @@
                 No = zq(ig,l,igcm_ccn_number)* tauscaling(ig) + 1.e-30
                 Rn = rice(ig,l)
-                Rn = -dlog(Rn) 
+                Rn = -log(Rn) 
                 Rm = Rn - 3. * sigma_ice*sigma_ice  
                 n_derf = derf( (rb_cldco2(1)+Rn) *dev3)
@@ -473 +473 @@
               dMh2o = 0.
               do i = 1, nbinco2_cld
-                 Proba    =1.0-dexp(-1.*microtimestep*rate(i))
+                 Proba    =1.0-exp(-1.*microtimestep*rate(i))
                  dN       = dN + n_aer(i) * Proba
                  dM       = dM + m_aer(i) * Proba             
@@ -479 +479 @@
               if (co2useh2o) then
                  do i = 1, nbinco2_cld
-                    Probah2o = 1.0-dexp(-1.*microtimestep*rateh2o(i))
+                    Probah2o = 1.0-exp(-1.*microtimestep*rateh2o(i))
                     dNh2o    = dNh2o + n_aer_h2oice(i) * Probah2o
                     dMh2o    = dMh2o + m_aer_h2oice(i) * Probah2o
@@ -535 +535 @@
      &         * tauscaling(ig)+1.e-30.ge. 1)THEN   ! we trigger crystal growth
               
-              call updaterice_microco2(zq(ig,l,igcm_co2_ice),
-     &            zq(ig,l,igcm_ccnco2_mass),zq(ig,l,igcm_ccnco2_number),
+              call updaterice_microco2(dble(zq(ig,l,igcm_co2_ice)),
+     &            dble(zq(ig,l,igcm_ccnco2_mass)),
+     &            dble(zq(ig,l,igcm_ccnco2_number)),
      &            tauscaling(ig),riceco2(ig,l),rhocloudco2(ig,l))