Changeset 633 for trunk/LMDZ.MARS/libf/phymars/growthrate.F

Timestamp:

Apr 25, 2012, 5:38:49 PM (13 years ago)

Author:

tnavarro

Message:

last scheme in commit r626 led to a wrong physical behaviour. This version uses a new subtimestep for microphysics that should be faster than the last one.

File:

• trunk/LMDZ.MARS/libf/phymars/growthrate.F (modified) (4 diffs)

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

@@ -1 @@
-      subroutine growthrate(temp,pmid,psat,seq,rcrystal,coeff1,coeff2)
+      subroutine growthrate(temp,pmid,psat,rcrystal,res)
 
       IMPLICIT NONE
@@ -31 +31 @@
       REAL temp     ! temperature in the middle of the layer (K)
       REAL pmid     ! pressure in the middle of the layer (K)
-      REAL*8 psat   ! water vapor saturation pressure (Pa) 
-      REAL*8 seq    ! Equilibrium saturation ratio
+      REAL psat   ! water vapor saturation pressure (Pa) 
       REAL rcrystal ! crystal radius before condensation (m)
 
 c     Output
-      REAL coeff1,coeff2     ! coefficients for the analytical relation between time and radius
+      REAL res      ! growth resistance (res=Rk+Rd)
 
 
@@ -42 +41 @@
 c   ------
 
-      REAL*8 k,Lv                 
-      REAL*8 knudsen           ! Knudsen number (gas mean free path/particle radius)
-      REAL*8 afactor,Dv,lambda       ! Intermediate computations for growth rate
-      REAL*8 Rk,Rd
+      REAL k,Lv                 
+      REAL knudsen           ! Knudsen number (gas mean free path/particle radius)
+      REAL afactor,Dv,lambda       ! Intermediate computations for growth rate
+      REAL Rk,Rd
       
       
@@ -86 +85 @@
 
 c     - Compute Rk
-      Rk = Lv**2 * rho_ice * mh2o / (k*rgp*temp*temp)
+      Rk = Lv*Lv* rho_ice * mh2o / (k*rgp*temp*temp)
 c     - Compute Rd
       Rd = rgp * temp *rho_ice / (Dv*psat*mh2o)
       
       
-      coeff1 = real(Rk + Rd*(1. + lambda * knudsen))
-      coeff2 = real(Rk + Rd*(1. - lambda * knudsen))
+      res = Rk + Rd*(1. + lambda * knudsen)
+      
+      !coeff1 = real(Rk + Rd*(1. + lambda * knudsen))
+      !coeff2 = real(Rk + Rd*(1. - lambda * knudsen))
       
 c Below are growth rate used for other schemes