Changeset 4049 for trunk/LMDZ.TITAN/libf


Ignore:
Timestamp:
Feb 5, 2026, 5:50:40 PM (2 weeks ago)
Author:
aslmd
Message:

Titan:

  • add computation of molar ratio for CH4 in tracer_h (rat_mmol), for moist dynamics without callchim
  • aesthetic changes in comchem_h

EMo

Location:
trunk/LMDZ.TITAN/libf/phytitan
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.TITAN/libf/phytitan/comchem_h.F90

    r2793 r4049  
    1414!
    1515! NB : A given order is assumed for the 44 chemistry tracers :
    16 !      H, H2, CH, CH2s, CH2, CH3, CH4, C2, C2H, C2H2, C2H3, C2H4, C2H5,
    17 !      C2H6, C3H3, C3H5, C3H6, C3H7, C4H, C4H3, C4H4, C4H2s, CH2CCH2,
    18 !      CH3CCH, C3H8, C4H2, C4H6, C4H10, AC6H6, C3H2, C4H5, AC6H5, N2,
    19 !      N4S, CN, HCN, H2CN, CHCN, CH2CN, CH3CN, C3N, HC3N, NCCN, C4N2
     16!      H,       H2,      CH,      CH2s,      CH2,      CH3,
     17!      CH4,     C2,      C2H,     C2H2,      C2H3,     C2H4,
     18!      C2H5,    C2H6,    C3H3,    C3H5,      C3H6,     C3H7,
     19!      C4H,     C4H3,    C4H4,    CH2CCH2,   CH3CCH,   C3H8,
     20!      C4H2,    C4H6,    C4H10,   AC6H6,     C3H2,     C4H5,
     21!      AC6H5,   N2,      N4S,     CN,        HCN,      H2CN,
     22!      CHCN,    CH2CN,   CH3CN,   C3N,       HC3N,     NCCN,
     23!      C4N2
    2024!
    2125!       
    … …  
    4751   !! Hard-coded chemical species molar mass (g.mol-1), shares the same indexing than cnames.
    4852   REAL, DIMENSION(nkim), PARAMETER               :: cmmol = (/ &
    49        1.01   , 2.0158, 13.02, 14.03, 14.03, 15.03, 16.04  , 24.02, 25.03, 26.04  , 27.05  , &
    50        28.05  , 29.06 , 30.07, 39.06, 41.07, 42.08, 43.09  , 49.05, 51.07, 52.08  , 40.07  , &
    51        40.07 , 44.11, 50.06, 54.09, 58.13, 78.1136, 38.05, 53.07, 77.1136, 28.0134, 14.01  , &
    52        26.02 , 27.04, 28.05, 39.05, 40.04, 41.05  , 50.04, 51.05, 52.04  , 76.1   /)
     53       1.01        , 2.0158      , 13.02       , 14.03       , 14.03       , 15.03       , &
     54       16.04       , 24.02       , 25.03       , 26.04       , 27.05       , 28.05       , &
     55       29.06       , 30.07       , 39.06       , 41.07       , 42.08       , 43.09       , &
     56       49.05       , 51.07       , 52.08       , 40.07       , 40.07       , 44.11       , &
     57       50.06       , 54.09       , 58.13       , 78.1136     , 38.05       , 53.07       , &
     58       77.1136     , 28.0134     , 14.01       , 26.02       , 27.04       , 28.05       , &
     59       39.05       , 40.04       , 41.05       , 50.04       , 51.05       , 52.04       , &
     60       76.1        /)
    5361
    5462   !! Hard-coded molar fraction of surface methane

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

    r3090 r4049  
    105105    mmol(:)  = 0.0
    106106    rat_mmol(:) = 1.0
     107    ! moist dynamics without callchim
     108    DO i=1,nq
     109      IF (TRIM(noms(i)) == 'CH4') THEN
     110        rat_mmol(i) = cmmol(7)/mugaz
     111      ENDIF
     112    ENDDO
    107113
    108114    ! Compute number of microphysics tracers:
Note: See TracChangeset for help on using the changeset viewer.