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

Timestamp:

Jul 9, 2014, 7:26:01 PM (11 years ago)

Author:

tnavarro

Message:

Changed two microphysics thresholds for timesteps used in Mesoscale and LES

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• improvedclouds.F (modified) (4 diffs)
• updaterad.F90 (modified) (2 diffs)

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

@@ -149 +149 @@
 
 c       Volume ratio between two adjacent bins
-   !     vrat_cld = log(rmax_cld/rmin_cld) / float(nbin_cld-1) *3.
-   !     vrat_cld = exp(vrat_cld)
+!        vrat_cld = log(rmax_cld/rmin_cld) / float(nbin_cld-1) *3.
+!        vrat_cld = exp(vrat_cld)
         vrat_cld = dlog(rmax_cld/rmin_cld) / float(nbin_cld-1) *3.
         vrat_cld = dexp(vrat_cld)
@@ -158 +158 @@
         rad_cld(1) = rmin_cld
         vol_cld(1) = 4./3. * dble(pi) * rmin_cld*rmin_cld*rmin_cld
-   !     vol_cld(1) = 4./3. * pi * rmin_cld*rmin_cld*rmin_cld
+!        vol_cld(1) = 4./3. * pi * rmin_cld*rmin_cld*rmin_cld
 
         do i=1,nbin_cld-1
@@ -224 +224 @@
       pdtcloud(1:ngrid,1:nlay)=0
     
-c     Initialize the tendencies
-      pdqcloud(1:ngrid,1:nlay,1:nq)=0
-      pdtcloud(1:ngrid,1:nlay)=0
-      
-c     Initialize the tendencies
-      pdqcloud(1:ngrid,1:nlay,1:nq)=0
-      pdtcloud(1:ngrid,1:nlay)=0
-      
       
       zt(1:ngrid,1:nlay) = 
@@ -403 +395 @@
 c         If all the ice particles sublimate, all the condensation
 c         nuclei are released:
-          if (zq(ig,l,igcm_h2o_ice).le.1.e-8) then
+          if (zq(ig,l,igcm_h2o_ice).le.1.e-28) then
           
 c           Water

trunk/LMDZ.MARS/libf/phymars/updaterad.F90

@@ -23 +23 @@
 real, parameter :: r3icemax = 125.e-12 ! ie ricemax  = 500 microns
 real, parameter :: ricemax  = 500.e-6 
-
-real, parameter :: qice_threshold  = 1.e-15 ! 1.e-10
-
 
 
@@ -79 +76 @@
 !! Indeed, one can obtain ricemin for small but non negligible qice values, and therefore hugely opaque clouds.
  
-if (qice .ge. qice_threshold) then
-
-  rhocloud = (qice*rho_ice + qccn_true*rho_dust) / (qice + qccn_true)
-  rhocloud = min(max(rhocloud,rho_ice),rho_dust)
-
-  rice = (qice + qccn_true) * 0.75 / pi / rhocloud / nccn_true
+
+rhocloud = (qice*rho_ice + qccn_true*rho_dust) / (qice + qccn_true)
+rhocloud = min(max(rhocloud,rho_ice),rho_dust)
+
+rice = (qice + qccn_true) * 0.75 / pi / rhocloud / nccn_true
   
-  if (rice .le. r3icemin) then
-    rice = ricemin
-  else if (rice .ge. r3icemax) then
-    rice = ricemax
-  else
-    rice = rice**(1./3.) ! here rice is always positive
-  endif
-
-else
-
+if (rice .le. r3icemin) then
   rice = ricemin
-  rhocloud = rho_dust
-
-endif
+else if (rice .ge. r3icemax) then
+  rice = ricemax
+else
+  rice = rice**(1./3.) ! here rice is always positive
+endif
+
 
 end subroutine updaterice_micro