Changeset 1617 for trunk/LMDZ.MARS/libf/phymars/microphys.h

Timestamp:

Oct 10, 2016, 11:26:09 AM (8 years ago)

Author:

jaudouard

Message:

Added modifications for CO2 clouds scheme in physiq_mod.F and added several routines and variables for CO2 microphysics. October 2016 Joachim Audouard

File:

• trunk/LMDZ.MARS/libf/phymars/microphys.h (modified) (2 diffs)

trunk/LMDZ.MARS/libf/phymars/microphys.h

@@ -19 +19 @@
 !     Molecular weight of CO2 (kg.mol-1)
       DOUBLE PRECISION, PARAMETER :: mco2 = 44.d-3
+!     Molecular weight of N2 (kg.mol-1)
+      DOUBLE PRECISION, PARAMETER :: mn2 = 28.01d-3
 !     Effective CO2 gas molecular radius (m)
       DOUBLE PRECISION, PARAMETER :: molco2 = 2.2d-10
@@ -48 +50 @@
 
 
+
+
+!CO2 part
+!      number of bins for nucleation
+      INTEGER, PARAMETER :: nbinco2_cld=40
+!     Surface tension of ice/vapor (J.m-2)
+      DOUBLE PRECISION, PARAMETER :: sigco2 = 0.08
+!     Activation energy for desorption of
+!       water on a dust-like substrate
+!       (J/molecule)
+      DOUBLE PRECISION, PARAMETER :: desorpco2 = 3.25e-20
+!     Jump frequency of a co2 molecule (s-1)
+      DOUBLE PRECISION, PARAMETER :: nusco2 =  2.9e+12
+!     Estimated activation energy for
+!       surface diffusion of co2 molecules
+!       (J/molecule)
+      DOUBLE PRECISION, PARAMETER :: surfdifco2 = desorpco2 / 10.
+!     Weight of a co2 molecule (kg)
+      DOUBLE PRECISION, PARAMETER :: m0co2 = mco2 / nav
+!     Contact parameter ( m=cos(theta) )
+!       (initialized in improvedCO2clouds.F)
+       REAL, parameter :: mtetaco2 = 0.952
+!     Volume of a co2 molecule (m3)
+       DOUBLE PRECISION :: vo1co2
+!     Radius used by the microphysical scheme (m)
+      DOUBLE PRECISION :: rad_cldco2(nbinco2_cld)
+       REAL, parameter :: threshJA = 1.0
+!     COMMON/microphys/vo1co2,rad_cldco2
+
 ! NB: to keep commons aligned: 
 !     split them in groups (reals, integers and characters)
 
-      COMMON/microphys/rad_cld,vo1
-      
-      COMMON/microphys_2/mteta
+      COMMON/microphys/rad_cld,vo1,rad_cldco2,vo1co2
+                  COMMON/microphys_2/mteta
       
 !     EXAMPLE: