Changeset 1894 for trunk/LMDZ.TITAN/libf/phytitan/tracer_h.F90

Timestamp:

Jan 15, 2018, 12:08:11 PM (7 years ago)

Author:

jvatant

Message:

Making chemistry more flexible - step 3.5
+ Update phyredem and phyetat0 ( with a chem_settings.F90 init routine )
+ Finish the handling of upper chem fields by their names everywhere
( with hardcoded cnames and mmol moved from tracer_h to comchem_h )
--JVO

File:

• trunk/LMDZ.TITAN/libf/phytitan/tracer_h.F90 (modified) (2 diffs)

trunk/LMDZ.TITAN/libf/phytitan/tracer_h.F90

@@ -68 +68 @@
     USE callkeys_mod
     USE comcstfi_mod, only: mugaz
+    USE comchem_h, only: cnames, cmmol
     IMPLICIT NONE
 
@@ -76 +77 @@
     LOGICAL                                      :: verb,found
     CHARACTER(len=20)                            :: str
-    !! Hard-coded chemical species for Titan chemistry
-    CHARACTER(len=10), DIMENSION(44), PARAMETER  :: cnames = &
-      (/"H         ", "H2        ", "CH        ", "CH2s      ", "CH2       ", "CH3       ", &
-        "CH4       ", "C2        ", "C2H       ", "C2H2      ", "C2H3      ", "C2H4      ", &
-        "C2H5      ", "C2H6      ", "C3H3      ", "C3H5      ", "C3H6      ", "C3H7      ", & 
-        "C4H       ", "C4H3      ", "C4H4      ", "C4H2s     ", "CH2CCH2   ", "CH3CCH    ", &
-        "C3H8      ", "C4H2      ", "C4H6      ", "C4H10     ", "AC6H6     ", "C3H2      ", &
-        "C4H5      ", "AC6H5     ", "N2        ", "N4S       ", "CN        ", "HCN       ", &
-        "H2CN      ", "CHCN      ", "CH2CN     ", "CH3CN     ", "C3N       ", "HC3N      ", &
-        "NCCN      ", "C4N2      "/)
-    !! Hard-coded chemical species molar mass (g.mol-1), shares the same indexing than cnames.
-    REAL, DIMENSION(44), PARAMETER               :: cmmol = (/ &
-        1.01   , 2.0158, 13.02, 14.03, 14.03, 15.03, 16.04  , 24.02, 25.03, 26.04  , 27.05  , &
-        28.05  , 29.06 , 30.07, 39.06, 41.07, 42.08, 43.09  , 49.05, 51.07, 52.08  , 50.06  , &
-        40.07  , 40.07 , 44.11, 50.06, 54.09, 58.13, 78.1136, 38.05, 53.07, 77.1136, 28.0134, &
-        14.01  , 26.02 , 27.04, 28.05, 39.05, 40.04, 41.05  , 50.04, 51.05, 52.04  , 76.1   /)
 
     INTEGER :: i,j,n