Changeset 1818 for trunk

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

Files:

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

trunk/LMDZ.MARS/README

@@ -2470 +2470 @@
 - Some cleaning of the CO2 ice clouds routine has been done. Not perfect yet!
  
-== 02/10/2017 == JA
+== 02/11/2017 == JA
 - Update for CO2 clouds scheme.
 - Several bugs have been corrected, further cleaning performed.
 - The main routine of co2 clouds is co2cloud.F. Please read its comments if interested.
 
+== 03/11/2017 == JA
+- Further cleaning of the CO2 ice clouds code
+- We can now have different micro-timesteps bewteen the water ice clouds and the co2 ice clouds
+- callphys.def needs imicroco2 = integer (typically 50) (water microtimestep is imicro=integer)
+- A filter based on Spiga et al 2012 Saturation index has been implemented. It is supposed to represent
+  if a gravity wave propagates to the mesosphere or not (wheter the wave saturates on its way or not). 
+  The sub-grid temperature distribution is cancelled if the saturation index value of the column (between 12 and 80 km) is > 0.1, and is applied (+/- 3K) otherwise. if the keyword satindexco2=.true. in callphys.def (only applies if CLFvaryingCO2=.true. anyway). If set to .false., deactivate this filter and the sub grid T distribution applies everywhere and everytime. See comments in co2clouds.F
+- All that you need to launch a GCM run with co2 clouds is described in co2clouds.F comments. Ehouarn Millour and Anni Määttänen have the files. 
  

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