Ignore:
Timestamp:
Jan 10, 2018, 12:25:01 PM (7 years ago)
Author:
jvatant
Message:

Making chemistry handling more flexible - step 2.5
+ For more convenience I introduce specific modules
for chemistry stuff specific to start2archive and newstart
and not to pollute main module comchem_h.
--JVO

File:
1 edited

Legend:

Unmodified
Added
Removed

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

    r1889 r1890  
    11MODULE comchem_h
    22
    3 ! Stores data relative to chemistry and routines for chemistry management in startfiles
     3! Stores data relative to upper chemistry in the GCM
     4! For the newstart and satrt2archive, there are specific comchem_***_h modules.
    45
    56! Author : Jan Vatant d'Ollone (2017)
    … …  
    1112!$OMP_THREADPRIVATE(nlaykim_up,nlay_kim_tot)
    1213
    13 ! NB : For the startfile we use nlaykim_up grid (upper atm) and for outputs we use nlaykim_tot grid (all layers)
     14  ! NB : For the startfile we use nlaykim_up grid (upper atm) and for outputs we use nlaykim_tot grid (all layers)
    1415
    1516   REAL,SAVE,ALLOCATABLE,DIMENSION(:) :: preskim  ! Pressure (Pa) of upper chemistry (mid)-layers
    … …  
    2021
    2122
    22   ! Variable and allocatables for regriding chemistry newstart
    23   ! ----------------------------------------------------------
    24    
    25   INTEGER :: nlaykimold ! Number of upper atm. layers for chemistry in the start_archive file
    26   REAL, ALLOCATABLE, DIMENSION(:) :: preskimold ! Pressure grid of upper chemistry in the start_archive file
    27 
    28   ! Allocatable arrays for start2archive
    29   ! ------------------------------------
    30  
    31   ! 1. Grille physique
    32   REAL,ALLOCATABLE :: H(:,:) 
    33   REAL,ALLOCATABLE :: H2(:,:) 
    34   REAL,ALLOCATABLE :: CH(:,:) 
    35   REAL,ALLOCATABLE :: CH2s(:,:) 
    36   REAL,ALLOCATABLE :: CH2(:,:) 
    37   REAL,ALLOCATABLE :: CH3(:,:) 
    38   REAL,ALLOCATABLE :: CH4(:,:) 
    39   REAL,ALLOCATABLE :: C2(:,:) 
    40   REAL,ALLOCATABLE :: C2H(:,:) 
    41   REAL,ALLOCATABLE :: C2H2(:,:) 
    42   REAL,ALLOCATABLE :: C2H3(:,:) 
    43   REAL,ALLOCATABLE :: C2H4(:,:) 
    44   REAL,ALLOCATABLE :: C2H5(:,:) 
    45   REAL,ALLOCATABLE :: C2H6(:,:) 
    46   REAL,ALLOCATABLE :: C3H3(:,:) 
    47   REAL,ALLOCATABLE :: C3H5(:,:) 
    48   REAL,ALLOCATABLE :: C3H6(:,:) 
    49   REAL,ALLOCATABLE :: C3H7(:,:) 
    50   REAL,ALLOCATABLE :: C4H(:,:) 
    51   REAL,ALLOCATABLE :: C4H3(:,:) 
    52   REAL,ALLOCATABLE :: C4H4(:,:) 
    53   REAL,ALLOCATABLE :: C4H2s(:,:) 
    54   REAL,ALLOCATABLE :: CH2CCH2(:,:) 
    55   REAL,ALLOCATABLE :: CH3CCH(:,:) 
    56   REAL,ALLOCATABLE :: C3H8(:,:) 
    57   REAL,ALLOCATABLE :: C4H2(:,:) 
    58   REAL,ALLOCATABLE :: C4H6(:,:) 
    59   REAL,ALLOCATABLE :: C4H10(:,:) 
    60   REAL,ALLOCATABLE :: AC6H6(:,:) 
    61   REAL,ALLOCATABLE :: C3H2(:,:) 
    62   REAL,ALLOCATABLE :: C4H5(:,:) 
    63   REAL,ALLOCATABLE :: AC6H5(:,:) 
    64   REAL,ALLOCATABLE :: N2(:,:) 
    65   REAL,ALLOCATABLE :: N4S(:,:) 
    66   REAL,ALLOCATABLE :: CN(:,:) 
    67   REAL,ALLOCATABLE :: HCN(:,:) 
    68   REAL,ALLOCATABLE :: H2CN(:,:) 
    69   REAL,ALLOCATABLE :: CHCN(:,:) 
    70   REAL,ALLOCATABLE :: CH2CN(:,:) 
    71   REAL,ALLOCATABLE :: CH3CN(:,:) 
    72   REAL,ALLOCATABLE :: C3N(:,:) 
    73   REAL,ALLOCATABLE :: HC3N(:,:) 
    74   REAL,ALLOCATABLE :: NCCN(:,:) 
    75   REAL,ALLOCATABLE :: C4N2(:,:) 
    76  
    77   ! 2. Grille scalaire 
    78   REAL,ALLOCATABLE :: H_S(:,:) 
    79   REAL,ALLOCATABLE :: H2_S(:,:) 
    80   REAL,ALLOCATABLE :: CH_S(:,:) 
    81   REAL,ALLOCATABLE :: CH2s_S(:,:) 
    82   REAL,ALLOCATABLE :: CH2_S(:,:) 
    83   REAL,ALLOCATABLE :: CH3_S(:,:) 
    84   REAL,ALLOCATABLE :: CH4_S(:,:) 
    85   REAL,ALLOCATABLE :: C2_S(:,:) 
    86   REAL,ALLOCATABLE :: C2H_S(:,:) 
    87   REAL,ALLOCATABLE :: C2H2_S(:,:) 
    88   REAL,ALLOCATABLE :: C2H3_S(:,:) 
    89   REAL,ALLOCATABLE :: C2H4_S(:,:) 
    90   REAL,ALLOCATABLE :: C2H5_S(:,:) 
    91   REAL,ALLOCATABLE :: C2H6_S(:,:) 
    92   REAL,ALLOCATABLE :: C3H3_S(:,:) 
    93   REAL,ALLOCATABLE :: C3H5_S(:,:) 
    94   REAL,ALLOCATABLE :: C3H6_S(:,:) 
    95   REAL,ALLOCATABLE :: C3H7_S(:,:) 
    96   REAL,ALLOCATABLE :: C4H_S(:,:) 
    97   REAL,ALLOCATABLE :: C4H3_S(:,:) 
    98   REAL,ALLOCATABLE :: C4H4_S(:,:) 
    99   REAL,ALLOCATABLE :: C4H2s_S(:,:) 
    100   REAL,ALLOCATABLE :: CH2CCH2_S(:,:) 
    101   REAL,ALLOCATABLE :: CH3CCH_S(:,:) 
    102   REAL,ALLOCATABLE :: C3H8_S(:,:) 
    103   REAL,ALLOCATABLE :: C4H2_S(:,:) 
    104   REAL,ALLOCATABLE :: C4H6_S(:,:) 
    105   REAL,ALLOCATABLE :: C4H10_S(:,:) 
    106   REAL,ALLOCATABLE :: AC6H6_S(:,:) 
    107   REAL,ALLOCATABLE :: C3H2_S(:,:) 
    108   REAL,ALLOCATABLE :: C4H5_S(:,:) 
    109   REAL,ALLOCATABLE :: AC6H5_S(:,:) 
    110   REAL,ALLOCATABLE :: N2_S(:,:) 
    111   REAL,ALLOCATABLE :: N4S_S(:,:) 
    112   REAL,ALLOCATABLE :: CN_S(:,:) 
    113   REAL,ALLOCATABLE :: HCN_S(:,:) 
    114   REAL,ALLOCATABLE :: H2CN_S(:,:) 
    115   REAL,ALLOCATABLE :: CHCN_S(:,:) 
    116   REAL,ALLOCATABLE :: CH2CN_S(:,:) 
    117   REAL,ALLOCATABLE :: CH3CN_S(:,:) 
    118   REAL,ALLOCATABLE :: C3N_S(:,:) 
    119   REAL,ALLOCATABLE :: HC3N_S(:,:) 
    120   REAL,ALLOCATABLE :: NCCN_S(:,:) 
    121   REAL,ALLOCATABLE :: C4N2_S(:,:)
    122  
    123  
    124   CONTAINS
    125 
    126 
    127   SUBROUTINE allokim_start2archive(ngridmx,ip1jmp1)
    128  
    129     IMPLICIT NONE
    130 
    131     INTEGER,INTENT(IN) :: ngridmx
    132     INTEGER,INTENT(IN) :: ip1jmp1
    133    
    134     allocate(H(ngridmx,nlaykim_up))
    135     allocate(H_S(ip1jmp1,nlaykim_up))
    136     allocate(H2(ngridmx,nlaykim_up))
    137     allocate(H2_S(ip1jmp1,nlaykim_up))
    138     allocate(CH(ngridmx,nlaykim_up))
    139     allocate(CH_S(ip1jmp1,nlaykim_up))
    140     allocate(CH2s(ngridmx,nlaykim_up))
    141     allocate(CH2s_S(ip1jmp1,nlaykim_up))
    142     allocate(CH2(ngridmx,nlaykim_up))
    143     allocate(CH2_S(ip1jmp1,nlaykim_up))
    144     allocate(CH3(ngridmx,nlaykim_up))
    145     allocate(CH3_S(ip1jmp1,nlaykim_up))
    146     allocate(CH4(ngridmx,nlaykim_up))
    147     allocate(CH4_S(ip1jmp1,nlaykim_up))
    148     allocate(C2(ngridmx,nlaykim_up))
    149     allocate(C2_S(ip1jmp1,nlaykim_up))
    150     allocate(C2H(ngridmx,nlaykim_up))
    151     allocate(C2H_S(ip1jmp1,nlaykim_up))
    152     allocate(C2H2(ngridmx,nlaykim_up))
    153     allocate(C2H2_S(ip1jmp1,nlaykim_up))
    154     allocate(C2H3(ngridmx,nlaykim_up))
    155     allocate(C2H3_S(ip1jmp1,nlaykim_up))
    156     allocate(C2H4(ngridmx,nlaykim_up))
    157     allocate(C2H4_S(ip1jmp1,nlaykim_up))
    158     allocate(C2H5(ngridmx,nlaykim_up))
    159     allocate(C2H5_S(ip1jmp1,nlaykim_up))
    160     allocate(C2H6(ngridmx,nlaykim_up))
    161     allocate(C2H6_S(ip1jmp1,nlaykim_up))
    162     allocate(C3H3(ngridmx,nlaykim_up))
    163     allocate(C3H3_S(ip1jmp1,nlaykim_up))
    164     allocate(C3H5(ngridmx,nlaykim_up))
    165     allocate(C3H5_S(ip1jmp1,nlaykim_up))
    166     allocate(C3H6(ngridmx,nlaykim_up))
    167     allocate(C3H6_S(ip1jmp1,nlaykim_up))
    168     allocate(C3H7(ngridmx,nlaykim_up))
    169     allocate(C3H7_S(ip1jmp1,nlaykim_up))
    170     allocate(C4H(ngridmx,nlaykim_up))
    171     allocate(C4H_S(ip1jmp1,nlaykim_up))
    172     allocate(C4H3(ngridmx,nlaykim_up))
    173     allocate(C4H3_S(ip1jmp1,nlaykim_up))
    174     allocate(C4H4(ngridmx,nlaykim_up))
    175     allocate(C4H4_S(ip1jmp1,nlaykim_up))
    176     allocate(C4H2s(ngridmx,nlaykim_up))
    177     allocate(C4H2s_S(ip1jmp1,nlaykim_up))
    178     allocate(CH2CCH2(ngridmx,nlaykim_up))
    179     allocate(CH2CCH2_S(ip1jmp1,nlaykim_up))
    180     allocate(CH3CCH(ngridmx,nlaykim_up))
    181     allocate(CH3CCH_S(ip1jmp1,nlaykim_up))
    182     allocate(C3H8(ngridmx,nlaykim_up))
    183     allocate(C3H8_S(ip1jmp1,nlaykim_up))
    184     allocate(C4H2(ngridmx,nlaykim_up))
    185     allocate(C4H2_S(ip1jmp1,nlaykim_up))
    186     allocate(C4H6(ngridmx,nlaykim_up))
    187     allocate(C4H6_S(ip1jmp1,nlaykim_up))
    188     allocate(C4H10(ngridmx,nlaykim_up))
    189     allocate(C4H10_S(ip1jmp1,nlaykim_up))
    190     allocate(AC6H6(ngridmx,nlaykim_up))
    191     allocate(AC6H6_S(ip1jmp1,nlaykim_up))
    192     allocate(C3H2(ngridmx,nlaykim_up))
    193     allocate(C3H2_S(ip1jmp1,nlaykim_up))
    194     allocate(C4H5(ngridmx,nlaykim_up))
    195     allocate(C4H5_S(ip1jmp1,nlaykim_up))
    196     allocate(AC6H5(ngridmx,nlaykim_up))
    197     allocate(AC6H5_S(ip1jmp1,nlaykim_up))
    198     allocate(N2(ngridmx,nlaykim_up))
    199     allocate(N2_S(ip1jmp1,nlaykim_up))
    200     allocate(N4S(ngridmx,nlaykim_up))
    201     allocate(N4S_S(ip1jmp1,nlaykim_up))
    202     allocate(CN(ngridmx,nlaykim_up))
    203     allocate(CN_S(ip1jmp1,nlaykim_up))
    204     allocate(HCN(ngridmx,nlaykim_up))
    205     allocate(HCN_S(ip1jmp1,nlaykim_up))
    206     allocate(H2CN(ngridmx,nlaykim_up))
    207     allocate(H2CN_S(ip1jmp1,nlaykim_up))
    208     allocate(CHCN(ngridmx,nlaykim_up))
    209     allocate(CHCN_S(ip1jmp1,nlaykim_up))
    210     allocate(CH2CN(ngridmx,nlaykim_up))
    211     allocate(CH2CN_S(ip1jmp1,nlaykim_up))
    212     allocate(CH3CN(ngridmx,nlaykim_up))
    213     allocate(CH3CN_S(ip1jmp1,nlaykim_up))
    214     allocate(C3N(ngridmx,nlaykim_up))
    215     allocate(C3N_S(ip1jmp1,nlaykim_up))
    216     allocate(HC3N(ngridmx,nlaykim_up))
    217     allocate(HC3N_S(ip1jmp1,nlaykim_up))
    218     allocate(NCCN(ngridmx,nlaykim_up))
    219     allocate(NCCN_S(ip1jmp1,nlaykim_up))
    220     allocate(C4N2(ngridmx,nlaykim_up))
    221     allocate(C4N2_S(ip1jmp1,nlaykim_up))
    222  
    223   END SUBROUTINE allokim_start2archive
    224 
    22523END MODULE comchem_h
Note: See TracChangeset for help on using the changeset viewer.