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

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

Location:

trunk/LMDZ.TITAN/libf/muphytitan

Files:

• mm_microphysic.f90 (modified) (1 diff)
• mmp_gcm.f90 (modified) (2 diffs)
• mmp_globals.f90 (modified) (1 diff)
• mmp_optics.f90 (modified) (1 diff)

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

@@ -202 +202 @@
     !! microphysics is disabled (see [[mm_globals(module):mm_w_clouds(variable)]] documentation).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL                 :: aer_prec   !! Aerosols precipitations (both modes) (m).
-    REAL(kind=mm_wp), INTENT(out), OPTIONAL                 :: ccn_prec   !! CCN precipitations (m).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: aer_s_flux !! Spherical aerosol mass flux (\(kg.m^{-2}.s^{-1}\)).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: aer_f_flux !! Fractal aerosol mass flux (\(kg.m^{-2}.s^{-1}\)).
+    REAL(kind=mm_wp), INTENT(out), OPTIONAL                 :: ccn_prec   !! CCN precipitations (m).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: ccn_flux   !! CCN mass flux (\(kg.m^{-2}.s^{-1}\)).
+    REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: ice_prec   !! Ice precipitations (m).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:,:) :: ice_fluxes !! Ice sedimentation fluxes (\(kg.m^{-2}.s^{-1}\)).
     REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:,:) :: gazs_sat   !! Condensible gaz saturation ratios (--).
-    REAL(kind=mm_wp), INTENT(out), OPTIONAL, DIMENSION(:)   :: ice_prec   !! Ice precipitations (m).
 
     IF (PRESENT(aer_prec))   aer_prec   = ABS(mm_aer_prec) 

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 = ''

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

@@ -92 +92 @@
   !> Aerosol electric charge correction control flag.
   LOGICAL, SAVE :: mmp_w_qe   = .true.
-
+  !> Optic look-up table file path.
+  CHARACTER(len=:), ALLOCATABLE, SAVE :: mmp_optic_file
 
   CONTAINS 

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

@@ -62 +62 @@
 
   PRIVATE
+  PUBLIC :: mmp_optic_file ! from mmp_globals :)
   PUBLIC :: mmp_initialize_optics
   PUBLIC :: mmp_sph_optics_vis,mmp_sph_optics_ir