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

Timestamp:

Nov 3, 2017, 1:57:01 PM (7 years ago)

Author:

jaudouard

Message:

Final update from J.A. about the CO2 clouds scheme for the LMDZ.MARS GCM

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• callkeys.h (modified) (2 diffs)
• co2cloud.F (modified) (3 diffs)
• conf_phys.F (modified) (1 diff)
• improvedCO2clouds.F (modified) (2 diffs)
• nucleaCO2.F (modified) (1 diff)
• physiq_mod.F (modified) (1 diff)

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

@@ -14 +14 @@
      &   ,activice,water,tifeedback,microphys,supersat,caps,photochem   &
      &   ,calltherm,callrichsl,callslope,tituscap,callyamada4,co2clouds,&
-     &   co2useh2o,meteo_flux,CLFvaryingCO2,spantCO2,CLFvarying
+     &   co2useh2o,meteo_flux,CLFvaryingCO2,spantCO2,CLFvarying,        &
+     &   satindexco2
      
       COMMON/callkeys_i/iradia,iaervar,iddist,ilwd,ilwb,ilwn,ncouche    &
@@ -60 +61 @@
       logical sedimentation
       logical activice,tifeedback,supersat,caps
-      logical co2clouds,co2useh2o,meteo_flux,CLFvaryingCO2
+      logical co2clouds,co2useh2o,meteo_flux,CLFvaryingCO2,satindexco2
       integer spantCO2
       logical CLFvarying

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

@@ -44 +44 @@
 c 07/2017 J.Audouard
 c Several logicals and integer must be set to .true. in callphys.def 
-c uf not, default values are .false.
+c if not, default values are .false.
 c co2clouds=.true. call this routine
 c co2useh2o=.true. allow the use of water ice particles as CCN for CO2
 c meteo_flux=.true. supply meteoritic particles 
 c CLFvaryingCO2=.true. allows a subgrid temperature distribution
-c of amplitude spantCO2(=integer in callphys.def)
+c of amplitude spantCO2(=integer in callphys.def, typically 3)
+c satindexco2=.true. allows the filtering out of the sub-grid T distribution
+c                    if the GW saturates in the column. Based on Spiga et al
+c                    2012
+c                    An index is computed for the column, and the sub-grid T 
+c                    distribution is applied if the index remains < 0.1
+c                    setting to .false. applies the sub-grid T everywhere.
+c                    default value is .true., only applies if 
+c                    CLFvaryingCO2=.true. anyway.
 c imicroco2=50 
 c
@@ -369 +377 @@
      &        pplay)
 
+         if (satindexco2) then !logical in callphys.def
          DO l=12,26
             ! layers 12 --> 26 ~ 12->85 km
@@ -391 +400 @@
          call WRITEDIAGFI(ngrid,"SatIndexmap","SatIndexmap","km",2,
      &        SatIndexmap)
+         else
+            DO ig=1,ngrid
+               SatIndexmap(ig)=0.05 !maxval(SatIndex(ig,12:26))
+            enddo
+         endif
 
 !Modulate the DeltaT by GW propagation index : 

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

@@ -472 +472 @@
          call getin("spantCO2",spantCO2)
          write(*,*)  " spantCO2 = ",spantCO2
+!Do you want to filter the sub-grid T distribution by a Saturation index?
+         write(*,*) "filter sub-grid temperature by Saturation index?"
+         satindexco2=.true.
+         call getin("satindexco2",satindexco2)
+         write(*,*)  " satindexco2 = ",satindexco2
+
 
 ! thermal inertia feedback

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

@@ -29 +29 @@
 c When CO2 ice is deposited onto a water ice particles, the particle is
 c removed from the water tracers.
+c Memory of the origin of the co2 particles is kept and thus the 
+c water cycle shouldn't be modified by this.
 cWARNING: no sedimentation of the water ice origin is performed
 c in the microphysical timestep in co2cloud.F. 
@@ -39 +41 @@
 c Adaptation for CO2 clouds by Joachim Audouard (09/16), based on the work
 c of Constantino Listowski 
+c There is an energy limit to how much co2 can sublimate/condensate. It is 
+c defined by the difference of the GCM temperature with the co2 condensation 
+c temperature. 
 c Warning:
 c If meteoritic particles are activated and turn into co2 ice particles,

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

@@ -18 +18 @@
 *   CO2 nucleation routine dev. by Constantino        *
 *     Listowski and Joachim Audouard (2016-2017),     *
-*     adapted from the water ice nucleation           *
+*     adapted from the water ice nucleation     
+* It computes two different nucleation rates : one 
+* on the dust CCN distribution and the other one on
+* the water ice particles distribution
 *******************************************************
  ! nucrate = output

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

@@ -1211 +1211 @@
 
 c   6a bis. CO2 clouds (CL & JA)
-
 c        ---------------------------------------
 c        CO2 ice cloud condensation in the atmosphere
 c        ----------------------------------------
-
+c flag needed in callphys.def:
+c               co2clouds=.true. is mandatory (default is .false.)
+c               co2useh2o=.true. if you want to allow co2 condensation 
+c                                on water ice particles 
+c               meteo_flux=.true. if you want to add a meteoritic
+c                                 supply of CCN
+c               CLFvaryingCO2=.true. if you want to have a sub-grid
+c                                    temperature distribution 
+c               spantCO2=integer (i.e. 3) amplitude of the sub-grid T disti
+c               nuiceco2_sed=0.2 variance of the size distribution for the 
+c                                sedimentation
+c               nuiceco2_ref=0.2 variance of the size distribution for the 
+c                                nucleation
+c               imicroco2=50     micro-timestep is 1/50 of physical timestep
          
          IF (co2clouds ) THEN