Changeset 1926 for trunk/LMDZ.TITAN/deftank

Timestamp:

May 28, 2018, 11:47:03 PM (7 years ago)

Author:

jvatant

Message:

1) Microphysics diags / outputs :

---

+ Add supplementary diagnostics outputs for microphysics ( precip, flux, rc ... ) ( new muphy_diag.F90 module )
+ Correct the outputs of microphys tracers to be in X/m-3 to be comparable to "standard values"

+ Also update the deftank callphys.def with latest revs modifs for microphysics

2) Condensation / chemistry updates :

---

+ Moved chemistry AFTER microphysics

• To have mufi condensation before photochem
• Chemistry called last coherent with the fact that it brings back fields to equilibrium

+ If 2D chemistry, make zonally averaged fields go through mufi and chem condensation

to have non saturated profiles in input of photochemistry
( other 'short' processes neglected as 2D -> no diurnal cycle, just seasonal evolution )

+ Also corrected the positivity check ( took Mars GCM syntax ) after chemistry ( could previously lead to negs )

3) Noticed a weird behaviour ( bug? ) :

---

+ If generalize the use of arrays *_indx for tracers, to get rid of ugly "iq+nmicro",

it ends up with weird results / crash in optim mode ( ok in debug ) but didn't find out why ...

--JVO

File:

• trunk/LMDZ.TITAN/deftank/callphys.def (modified) (2 diffs)

trunk/LMDZ.TITAN/deftank/callphys.def

@@ -77 +77 @@
 # call chemistry ? (must have tracer=true)
 callchim = .false.
-# Activate zonal mean ? (must be true if callchim)
+# Activate zonal mean ? 
 moyzon_ch = .false.
 # chemistry is computed every "ichim" phys. timestep
@@ -86 +86 @@
 # call microphysics ? (must have tracer=true)
 callmufi = .false.
-# Activate zonal mean ? (must be true if callmufi by now)
-moyzon_mu = .false.
-## If yes, clouds ?
+## If yes, compute clouds ?
 callclouds = .false.
 ### If yes, number of ices ? (must be compatible with traceur.def AND microphysical model )
 nices = 3
+# Desactivate radiative coupling of microphysics ( i.e. use mean opacity profile )
+uncoupl_optic_haze = .true.
 # ~~~~~ Parameters for microphysics ~~~~~
 # Path to microphys. config file ?
 config_mufi = ../../datagcm/microphysics/config.ne15.cfg
-# Fractal dimension ?
-df_mufi = 2.0
-# Monomer radius (m) ?
-rm_mufi = 6.66e-08
-# Aerosol density (kg.m-3) ?
-rho_aer_mufi = 1.e3
 # Pressure level of aer. production (Pa) ?
 p_prod = 1.0