Changeset 3913 for trunk/LMDZ.MARS/libf/phymars/updaterad.F90

Timestamp:

Sep 5, 2025, 3:21:38 PM (2 months ago)

Author:

jbclement

Message:

Mars PCM:
Partial reversion of r3847: removal of the incompatibility test requires 'microphys'?. Therefore, in "updatereffrad_mod.F", without microphysics, we use a typical value for 'tau' (0.2) to assess the number of cloud condensation nuclei in the same way than done before r744. Comments and unused variables are cleaned.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/updaterad.F90 (modified) (6 diffs)

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

@@ -39 +39 @@
 real, parameter :: rccnmax  = 500.e-6
 
-
-
 real, parameter :: ndust_threshold  = 1.
 real, parameter :: qdust_threshold  = 1.e-20
@@ -50 +48 @@
 real, parameter :: rdustmax  = 500.e-6  
 
-
 real, parameter :: rdust0 = 0.8e-6
 
-
-      
 
 contains
@@ -150 +145 @@
 ! Update ice radius from a typical profile if microphys == false
 subroutine updaterice_typ(qice,tau,pzlay,rice)
+
 use tracer_mod, only: rho_ice
-USE comcstfi_h
+USE comcstfi_h, only: pi
+
 implicit none
 
@@ -158 +155 @@
 real, intent(in)  :: pzlay ! altitude at the middle of the layers
 real, intent(out) :: rice
-real rccn,nccn ! radius  and number of ice crystals
+real :: rccn, nccn ! radius  and number of ice crystals
 
 !   Typical CCN profile following Montmessin et al. 2004
-!   (N0=2e6 m-3 has been converted into N0=1.3e8 kg-1, otherwise the equation for rice is not homogeneous...)
-     nccn  = 1.3e+8*max(tau,0.001)/0.1*exp(-pzlay/10000.)
+!   (N0=2e6 m-3 has been converted into N0=1.3e8 kg-1 using a typical atmospheric density,
+!   otherwise the equation for rice is not homogeneous...)
+     nccn = 1.3e+8*max(tau,0.001)/0.1*exp(-pzlay/10000.)
 !   The previously used profile was not correct:
 !   Nccn=( epaisseur/masse ) * 2.e+6/0.1*max(tau,0.001)*exp(-pzlay/10000.)
-     
+
 if (nccn .le. 1) then
 
@@ -174 +172 @@
 ! Typical dust radius profile:
   rccn = max(rdust0*exp(-pzlay/18000.),1.e-9)
-  rice  = qice * 0.75 / pi / rho_ice / nccn + rccn*rccn*rccn
-       
+  rice = qice * 0.75 / pi / rho_ice / nccn + rccn*rccn*rccn
+
   if (rice .le. r3icemin) then
     rice = ricemin
@@ -185 +183 @@
 
 endif
- 
+
 end subroutine updaterice_typ
 !============================================================================