Changeset 1926 for trunk/LMDZ.TITAN/libf/muphytitan/mmp_gcm.f90

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/libf/muphytitan/mmp_gcm.f90 (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/muphytitan/mmp_gcm.f90

@@ -93 +93 @@
     INTEGER                                           :: i
     TYPE(cfgparser)                                   :: cparser
-    CHARACTER(len=st_slen)                            :: spcpath,pssfile,mqfile
+    CHARACTER(len=st_slen)                            :: spcpath,pssfile,mqfile,opt_file
     CHARACTER(len=st_slen), DIMENSION(:), ALLOCATABLE :: species
     REAL(kind=mm_wp), DIMENSION(:), ALLOCATABLE       :: tmp
@@ -138 +138 @@
     err = mm_check_opt(cfg_get_value(cparser,"haze_coag_interactions",coag_choice),coag_choice,7,mm_log)
 
+    ! optic look-up table file path.
+    mmp_optic_file = ''
+    opt_file = ''
+    err = mm_check_opt(cfg_get_value(cparser,"optics/optic_file",opt_file),opt_file,'',mm_log)
+    IF (err /= 0) THEN
+      WRITE(*,'(a)') "Warning: I was unable to retrieve the path of the optic look-up table file:"
+      WRITE(*,'(a)') "  The GCM may abort if it uses YAMMS optical properties calculation module !"
+    ELSE
+      mmp_optic_file = TRIM(opt_file)
+    ENDIF
+
     ! Retrieve clouds species configuration file
     spcpath = ''