[2320] | 1 | |
---|
| 2 | ! $Id: $ |
---|
| 3 | |
---|
| 4 | MODULE infotrac_phy |
---|
| 5 | |
---|
[4170] | 6 | USE strings_mod, ONLY: msg, maxlen, strStack, strHead, strParse, strIdx, int2str |
---|
[4120] | 7 | USE readTracFiles_mod, ONLY: trac_type, isot_type, keys_type, delPhase, getKey, tnom_iso => newH2OIso |
---|
[2320] | 8 | |
---|
[4120] | 9 | IMPLICIT NONE |
---|
[4046] | 10 | |
---|
[4120] | 11 | PRIVATE |
---|
[2320] | 12 | |
---|
[4120] | 13 | !=== FOR TRACERS: |
---|
| 14 | PUBLIC :: init_infotrac_phy !--- Initialization of the tracers |
---|
[4170] | 15 | PUBLIC :: tracers, type_trac, types_trac !--- Full tracers database, tracers type keyword |
---|
[4120] | 16 | PUBLIC :: nqtot, nbtr, nqo, nqCO2, nqtottr !--- Main dimensions |
---|
[4124] | 17 | PUBLIC :: conv_flg, pbl_flg !--- Convection & boundary layer activation keys |
---|
[2320] | 18 | |
---|
[4120] | 19 | !=== FOR ISOTOPES: General |
---|
| 20 | PUBLIC :: isotopes, nbIso !--- Derived type, full isotopes families database + nb of families |
---|
| 21 | PUBLIC :: isoSelect, ixIso !--- Isotopes family selection tool + selected family index |
---|
| 22 | !=== FOR ISOTOPES: Specific to water |
---|
| 23 | PUBLIC :: iH2O !--- H2O isotopes index |
---|
| 24 | !=== FOR ISOTOPES: Depending on the selected isotopes family |
---|
| 25 | PUBLIC :: isotope, isoKeys !--- Selected isotopes database + associated keys (cf. getKey) |
---|
| 26 | PUBLIC :: isoName, isoZone, isoPhas !--- Isotopes and tagging zones names, phases |
---|
| 27 | PUBLIC :: niso, nzone, nphas, ntiso !--- " " numbers + isotopes & tagging tracers number |
---|
| 28 | PUBLIC :: itZonIso !--- iq = function(tagging zone idx, isotope idx) |
---|
[4143] | 29 | PUBLIC :: iqIsoPha !--- idx of tagging tracer in iName = function(isotope idx, phase idx) |
---|
[4120] | 30 | PUBLIC :: isoCheck !--- Run isotopes checking routines |
---|
| 31 | !=== FOR BOTH TRACERS AND ISOTOPES |
---|
| 32 | PUBLIC :: getKey !--- Get a key from "tracers" or "isotope" |
---|
[2320] | 33 | |
---|
[4120] | 34 | INTERFACE isoSelect; MODULE PROCEDURE isoSelectByIndex, isoSelectByName; END INTERFACE isoSelect |
---|
[3865] | 35 | |
---|
[4120] | 36 | !=== CONVENTIONS FOR TRACERS NUMBERS: |
---|
| 37 | ! |--------------------+-----------------------+-----------------+---------------+----------------------------| |
---|
| 38 | ! | water in different | water tagging | water isotopes | other tracers | additional tracers moments | |
---|
| 39 | ! | phases: H2O_[gls] | isotopes | | | for higher order schemes | |
---|
| 40 | ! |--------------------+-----------------------+-----------------+---------------+----------------------------| |
---|
| 41 | ! | | | | | | |
---|
| 42 | ! |<-- nqo -->|<-- nqo*niso* nzone -->|<-- nqo*niso -->|<-- nbtr -->|<-- (nmom) -->| |
---|
| 43 | ! | | | | |
---|
| 44 | ! | |<-- nqo*niso*(nzone+1) = nqo*ntiso -->|<-- nqtottr = nbtr + nmom -->| |
---|
| 45 | ! | = nqtot - nqo*(ntiso+1) | |
---|
| 46 | ! | | |
---|
| 47 | ! |<-- nqtrue = nbtr + nqo*(ntiso+1) -->| | |
---|
| 48 | ! | | |
---|
| 49 | ! |<-- nqtot = nqtrue + nmom -->| |
---|
| 50 | ! | | |
---|
| 51 | ! |-----------------------------------------------------------------------------------------------------------| |
---|
| 52 | ! NOTES FOR THIS TABLE: |
---|
| 53 | ! * Used "niso", "nzone" and "ntiso" are components of "isotopes(ip)" for water (isotopes(ip)%parent == 'H2O'), |
---|
| 54 | ! since water is so far the sole tracers family, except passive CO2, removed from the main tracers table. |
---|
| 55 | ! * For water, "nqo" is equal to the more general field "isotopes(ip)%nphas". |
---|
| 56 | ! * "niso", "nzone", "ntiso", "nphas" are defined for other isotopic tracers families, if any. |
---|
| 57 | ! |
---|
| 58 | !=== DERIVED TYPE EMBEDDING MOST OF THE TRACERS-RELATED QUANTITIES (LENGTH: nqtot) |
---|
| 59 | ! Each entry is accessible using "%" sign. |
---|
| 60 | ! |-------------+------------------------------------------------------+-------------+------------------------+ |
---|
| 61 | ! | entry | Meaning | Former name | Possible values | |
---|
| 62 | ! |-------------+------------------------------------------------------+-------------+------------------------+ |
---|
| 63 | ! | name | Name (short) | tname | | |
---|
| 64 | ! | gen0Name | Name of the 1st generation ancestor | / | | |
---|
| 65 | ! | parent | Name of the parent | / | | |
---|
| 66 | ! | longName | Long name (with adv. scheme suffix) for outputs | ttext | | |
---|
| 67 | ! | type | Type (so far: tracer or tag) | / | tracer,tag | |
---|
| 68 | ! | phase | Phases list ("g"as / "l"iquid / "s"olid) | / | [g][l][s] | |
---|
| 69 | ! | component | Name(s) of the merged/cumulated section(s) | / | coma-separated names | |
---|
| 70 | ! | iadv | Advection scheme number | iadv | 1-20,30 exc. 3-9,15,19 | |
---|
| 71 | ! | iGeneration | Generation (>=1) | / | | |
---|
| 72 | ! | isAdvected | advected tracers flag (.TRUE. if iadv >= 0) | / | nqtrue .TRUE. values | |
---|
| 73 | ! | isInPhysics | tracers not extracted from the main table in physics | / | nqtottr .TRUE. values | |
---|
| 74 | ! | iqParent | Index of the parent tracer | iqpere | 1:nqtot | |
---|
| 75 | ! | iqDescen | Indexes of the childs (all generations) | iqfils | 1:nqtot | |
---|
| 76 | ! | nqDescen | Number of the descendants (all generations) | nqdesc | 1:nqtot | |
---|
| 77 | ! | nqChilds | Number of childs (1st generation only) | nqfils | 1:nqtot | |
---|
| 78 | ! | iso_iGroup | Isotopes group index in isotopes(:) | / | 1:nbIso | |
---|
| 79 | ! | iso_iName | Isotope name index in isotopes(iso_iGroup)%trac(:) | iso_indnum | 1:niso | |
---|
| 80 | ! | iso_iZone | Isotope zone index in isotopes(iso_iGroup)%zone(:) | zone_num | 1:nzone | |
---|
| 81 | ! | iso_iPhas | Isotope phase index in isotopes(iso_iGroup)%phas(:) | phase_num | 1:nphas | |
---|
| 82 | ! | keys | key/val pairs accessible with "getKey" routine | / | | |
---|
| 83 | ! +-------------+------------------------------------------------------+-------------+------------------------+ |
---|
| 84 | ! |
---|
| 85 | !=== DERIVED TYPE EMBEDDING MOST OF THE ISOTOPES-RELATED QUANTITIES (LENGTH: nbIso, NUMBER OF ISOTOPES FAMILIES) |
---|
| 86 | ! Each entry is accessible using "%" sign. |
---|
| 87 | ! |-----------------+--------------------------------------------------+--------------------+-----------------+ |
---|
| 88 | ! | entry | length | Meaning | Former name | Possible values | |
---|
| 89 | ! |-----------------+--------------------------------------------------+--------------------+-----------------+ |
---|
| 90 | ! | parent | Parent tracer (isotopes family name) | | | |
---|
| 91 | ! | keys | niso | Isotopes keys/values pairs list + number | | | |
---|
| 92 | ! | trac | ntiso | Isotopes + tagging tracers list + number | / | ntraciso | | |
---|
| 93 | ! | zone | nzone | Geographic tagging zones list + number | / | ntraceurs_zone | | |
---|
| 94 | ! | phase | nphas | Phases list + number | | [g][l][s], 1:3 | |
---|
[4143] | 95 | ! | iqIsoPha | Index in "qx" = f(name(1:ntiso)),phas) | iqiso | 1:nqtot | |
---|
[4120] | 96 | ! | itZonIso | Index in "trac(1:ntiso)"= f(zone, name(1:niso)) | index_trac | 1:ntiso | |
---|
| 97 | ! +-----------------+--------------------------------------------------+--------------------+-----------------+ |
---|
| 98 | |
---|
| 99 | !=== DIMENSIONS OF THE TRACERS TABLES AND OTHER SCALAR VARIABLES |
---|
| 100 | INTEGER, SAVE :: nqtot, & !--- Tracers nb in dynamics (incl. higher moments + H2O) |
---|
| 101 | nbtr, & !--- Tracers nb in physics (excl. higher moments + H2O) |
---|
| 102 | nqo, & !--- Number of water phases |
---|
| 103 | nbIso, & !--- Number of available isotopes family |
---|
| 104 | nqtottr, & !--- Number of tracers passed to phytrac (TO BE DELETED ?) |
---|
| 105 | nqCO2 !--- Number of tracers of CO2 (ThL) |
---|
| 106 | CHARACTER(LEN=maxlen), SAVE :: type_trac !--- Keyword for tracers type |
---|
[4170] | 107 | CHARACTER(LEN=maxlen), SAVE, ALLOCATABLE :: types_trac(:) !--- Keyword for tracers type |
---|
| 108 | !$OMP THREADPRIVATE(nqtot, nbtr, nqo, nbIso, nqtottr, nqCO2, type_trac, types_trac) |
---|
| 109 | |
---|
[4120] | 110 | !=== DERIVED TYPES EMBEDDING MOST INFORMATIONS ABOUT TRACERS AND ISOTOPES FAMILIES |
---|
| 111 | TYPE(trac_type), TARGET, SAVE, ALLOCATABLE :: tracers(:) !=== TRACERS DESCRIPTORS VECTOR |
---|
| 112 | TYPE(isot_type), TARGET, SAVE, ALLOCATABLE :: isotopes(:) !=== ISOTOPES PARAMETERS VECTOR |
---|
| 113 | !$OMP THREADPRIVATE(tracers, isotopes) |
---|
| 114 | |
---|
| 115 | !=== ALIASES FOR CURRENTLY SELECTED ISOTOPES FAMILY OF VARIABLES EMBEDDED IN THE VECTOR "isotopes" |
---|
| 116 | TYPE(isot_type), SAVE, POINTER :: isotope !--- CURRENTLY SELECTED ISOTOPES FAMILY DESCRIPTOR |
---|
| 117 | INTEGER, SAVE :: ixIso, iH2O !--- Index of the selected isotopes family and H2O family |
---|
[4154] | 118 | LOGICAL, SAVE :: isoCheck !--- Flag to trigger the checking routines |
---|
[4120] | 119 | TYPE(keys_type), SAVE, POINTER :: isoKeys(:) !--- ONE SET OF KEYS FOR EACH ISOTOPE (LISTED IN isoName) |
---|
| 120 | CHARACTER(LEN=maxlen), SAVE, POINTER :: isoName(:), & !--- ISOTOPES NAMES FOR THE CURRENTLY SELECTED FAMILY |
---|
| 121 | isoZone(:), & !--- TAGGING ZONES FOR THE CURRENTLY SELECTED FAMILY |
---|
| 122 | isoPhas !--- USED PHASES FOR THE CURRENTLY SELECTED FAMILY |
---|
[4154] | 123 | INTEGER, SAVE :: niso, nzone, & !--- NUMBER OF ISOTOPES, TAGGING ZONES AND PHASES |
---|
| 124 | nphas, ntiso !--- NUMBER OF PHASES AND ISOTOPES + ISOTOPIC TAGGING TRACERS |
---|
| 125 | INTEGER, SAVE, POINTER ::itZonIso(:,:), & !--- INDEX IN "isoTrac" AS f(tagging zone idx, isotope idx) |
---|
[4143] | 126 | iqIsoPha(:,:) !--- INDEX IN "qx" AS f(isotopic tracer idx, phase idx) |
---|
| 127 | !$OMP THREADPRIVATE(isotope, ixIso,iH2O, isoCheck, isoKeys, isoName,isoZone,isoPhas, niso,nzone,nphas,ntiso, itZonIso,iqIsoPha) |
---|
[4120] | 128 | |
---|
| 129 | !=== VARIABLES FOR ISOTOPES INITIALIZATION AND FOR INCA |
---|
| 130 | INTEGER, SAVE, ALLOCATABLE ::conv_flg(:), & !--- Convection activation ; needed for INCA (nbtr) |
---|
| 131 | pbl_flg(:) !--- Boundary layer activation ; needed for INCA (nbtr) |
---|
[4124] | 132 | !$OMP THREADPRIVATE(conv_flg, pbl_flg) |
---|
[4120] | 133 | |
---|
[3677] | 134 | #ifdef CPP_StratAer |
---|
[4056] | 135 | !=== SPECIFIC TO STRATOSPHERIC AEROSOLS (CK/OB) |
---|
| 136 | INTEGER, SAVE :: nbtr_bin, nbtr_sulgas !--- number of aerosols bins and sulfur gases for StratAer model |
---|
| 137 | !$OMP THREADPRIVATE(nbtr_bin, nbtr_sulgas) |
---|
| 138 | INTEGER, SAVE :: id_OCS_strat, id_SO2_strat, id_H2SO4_strat, id_BIN01_strat, id_TEST_strat |
---|
| 139 | !$OMP THREADPRIVATE(id_OCS_strat, id_SO2_strat, id_H2SO4_strat, id_BIN01_strat, id_TEST_strat) |
---|
[3677] | 140 | #endif |
---|
| 141 | |
---|
[4120] | 142 | CONTAINS |
---|
[2320] | 143 | |
---|
[4124] | 144 | SUBROUTINE init_infotrac_phy(type_trac_, tracers_, isotopes_, nqtottr_, nqCO2_, pbl_flg_, conv_flg_) |
---|
[2320] | 145 | |
---|
[4120] | 146 | USE print_control_mod, ONLY: prt_level, lunout |
---|
[2320] | 147 | |
---|
[4120] | 148 | IMPLICIT NONE |
---|
| 149 | CHARACTER(LEN=*),INTENT(IN) :: type_trac_ |
---|
| 150 | TYPE(trac_type), INTENT(IN) :: tracers_(:) |
---|
| 151 | TYPE(isot_type), INTENT(IN) :: isotopes_(:) |
---|
| 152 | INTEGER, INTENT(IN) :: nqtottr_ |
---|
| 153 | INTEGER, INTENT(IN) :: nqCO2_ |
---|
| 154 | INTEGER, INTENT(IN) :: conv_flg_(:) |
---|
| 155 | INTEGER, INTENT(IN) :: pbl_flg_(:) |
---|
[3924] | 156 | |
---|
[4120] | 157 | INTEGER :: iq, ixt |
---|
| 158 | #ifdef CPP_StratAer |
---|
| 159 | CHARACTER(LEN=maxlen), ALLOCATABLE :: tnames(:) |
---|
| 160 | #endif |
---|
| 161 | CHARACTER(LEN=maxlen) :: modname="init_infotrac_phy" |
---|
[2320] | 162 | |
---|
[4120] | 163 | type_trac = type_trac_ |
---|
[4191] | 164 | IF(strParse(type_trac, '|', types_trac)) CALL abort_physic(modname,'can''t parse "type_trac = '//TRIM(type_trac)//'"',1) |
---|
[4120] | 165 | tracers = tracers_ |
---|
| 166 | isotopes = isotopes_ |
---|
| 167 | nqtottr = nqtottr_ |
---|
| 168 | nqCO2 = nqCO2_ |
---|
| 169 | pbl_flg = pbl_flg_ |
---|
| 170 | conv_flg = conv_flg_ |
---|
| 171 | nqtot = SIZE(tracers_) |
---|
[4138] | 172 | nqo = COUNT(delPhase(tracers%name)=='H2O' .AND. tracers%iGeneration==0 .AND. tracers%component=='lmdz') |
---|
[4120] | 173 | nbtr = SIZE(conv_flg) |
---|
| 174 | nbIso = SIZE(isotopes_) |
---|
[2690] | 175 | |
---|
[4120] | 176 | !=== Determine selected isotopes class related quantities: |
---|
[4143] | 177 | ! ixIso, isotope, niso,isoKeys, ntiso,isoName, nzone,isoZone, nphas,isoPhas, itZonIso, iqIsoPha, isoCheck |
---|
[4120] | 178 | IF(.NOT.isoSelect('H2O')) iH2O = ixIso |
---|
| 179 | IF(prt_level > 1) THEN |
---|
| 180 | CALL msg('nqtot = '//TRIM(int2str(nqtot)), modname) |
---|
| 181 | CALL msg('nbtr = '//TRIM(int2str(nbtr )), modname) |
---|
| 182 | CALL msg('nqo = '//TRIM(int2str(nqo )), modname) |
---|
| 183 | CALL msg('niso = '//TRIM(int2str(niso )), modname) |
---|
| 184 | CALL msg('ntiso = '//TRIM(int2str(ntiso)), modname) |
---|
| 185 | CALL msg('nqtottr = '//TRIM(int2str(nqtottr)), modname) |
---|
| 186 | CALL msg('nqCO2 = '//TRIM(int2str(nqCO2)), modname) |
---|
| 187 | END IF |
---|
[2690] | 188 | |
---|
[3677] | 189 | #ifdef CPP_StratAer |
---|
[4170] | 190 | IF (ANY(types_trac == 'coag')) THEN |
---|
[4056] | 191 | nbtr_bin = COUNT([(tracers(iq)%name(1:3)=='BIN', iq=1, nqtot)]) |
---|
| 192 | nbtr_sulgas = COUNT([(tracers(iq)%name(1:3)=='GAS', iq=1, nqtot)]) |
---|
[4071] | 193 | tnames = PACK(tracers(:)%name, MASK=tracers(:)%isInPhysics) |
---|
[4056] | 194 | id_BIN01_strat = strIdx(tnames, 'BIN01' ) |
---|
[4232] | 195 | id_OCS_strat = strIdx(tnames, 'GASOCS' ) |
---|
[4056] | 196 | id_SO2_strat = strIdx(tnames, 'GASSO2' ) |
---|
| 197 | id_H2SO4_strat = strIdx(tnames, 'GASH2SO4') |
---|
| 198 | id_TEST_strat = strIdx(tnames, 'GASTEST' ) |
---|
[4120] | 199 | CALL msg('nbtr_bin ='//TRIM(int2str(nbtr_bin )), modname) |
---|
| 200 | CALL msg('nbtr_sulgas ='//TRIM(int2str(nbtr_sulgas )), modname) |
---|
| 201 | CALL msg('id_BIN01_strat ='//TRIM(int2str(id_BIN01_strat)), modname) |
---|
| 202 | CALL msg('id_OCS_strat ='//TRIM(int2str(id_OCS_strat )), modname) |
---|
| 203 | CALL msg('id_SO2_strat ='//TRIM(int2str(id_SO2_strat )), modname) |
---|
| 204 | CALL msg('id_H2SO4_strat ='//TRIM(int2str(id_H2SO4_strat)), modname) |
---|
| 205 | CALL msg('id_TEST_strat ='//TRIM(int2str(id_TEST_strat )), modname) |
---|
| 206 | END IF |
---|
[3677] | 207 | #endif |
---|
[4056] | 208 | #ifdef ISOVERIF |
---|
[4120] | 209 | CALL msg('iso_iName = '//strStack(int2str(PACK(tracers(:)%iso_iName, MASK=tracers(:)%iso_iGroup==iH2O))), modname) |
---|
[4056] | 210 | #endif |
---|
| 211 | |
---|
[4120] | 212 | END SUBROUTINE init_infotrac_phy |
---|
[2320] | 213 | |
---|
[4120] | 214 | |
---|
| 215 | !============================================================================================================================== |
---|
| 216 | !=== THE ROUTINE isoSelect IS USED TO SWITCH FROM AN ISOTOPE FAMILY TO ANOTHER: ISOTOPES DEPENDENT PARAMETERS ARE UPDATED |
---|
| 217 | ! Single generic "isoSelect" routine, using the predefined index of the parent (fast version) or its name (first call). |
---|
| 218 | !============================================================================================================================== |
---|
| 219 | LOGICAL FUNCTION isoSelectByName(iName, lVerbose) RESULT(lerr) |
---|
| 220 | IMPLICIT NONE |
---|
| 221 | CHARACTER(LEN=*), INTENT(IN) :: iName |
---|
| 222 | LOGICAL, OPTIONAL, INTENT(IN) :: lVerbose |
---|
| 223 | INTEGER :: iIso |
---|
| 224 | LOGICAL :: lV |
---|
| 225 | lV = .FALSE.; IF(PRESENT(lVerbose)) lV = lVerbose |
---|
| 226 | iIso = strIdx(isotopes(:)%parent, iName) |
---|
| 227 | lerr = iIso == 0 |
---|
[4154] | 228 | IF(lerr) THEN |
---|
| 229 | niso = 0; ntiso = 0; nzone=0; nphas=nqo; isoCheck=.FALSE. |
---|
| 230 | CALL msg('no isotope family named "'//TRIM(iName)//'"', ll=lV) |
---|
| 231 | RETURN |
---|
| 232 | END IF |
---|
[4120] | 233 | lerr = isoSelectByIndex(iIso, lV) |
---|
| 234 | END FUNCTION isoSelectByName |
---|
| 235 | !============================================================================================================================== |
---|
| 236 | LOGICAL FUNCTION isoSelectByIndex(iIso, lVerbose) RESULT(lerr) |
---|
| 237 | IMPLICIT NONE |
---|
| 238 | INTEGER, INTENT(IN) :: iIso |
---|
| 239 | LOGICAL, OPTIONAL, INTENT(IN) :: lVerbose |
---|
| 240 | LOGICAL :: lv |
---|
| 241 | lv = .FALSE.; IF(PRESENT(lVerbose)) lv = lVerbose |
---|
| 242 | lerr = .FALSE. |
---|
| 243 | IF(iIso == ixIso) RETURN !--- Nothing to do if the index is already OK |
---|
| 244 | lerr = iIso<=0 .OR. iIso>nbIso |
---|
| 245 | CALL msg('Inconsistent isotopes family index '//TRIM(int2str(iIso))//': should be > 0 and <= '//TRIM(int2str(nbIso))//'"',& |
---|
| 246 | ll=lerr .AND. lV) |
---|
| 247 | IF(lerr) RETURN |
---|
| 248 | ixIso = iIso !--- Update currently selected family index |
---|
| 249 | isotope => isotopes(ixIso) !--- Select corresponding component |
---|
[4154] | 250 | isoKeys => isotope%keys; niso = isotope%niso |
---|
| 251 | isoName => isotope%trac; ntiso = isotope%ntiso |
---|
| 252 | isoZone => isotope%zone; nzone = isotope%nzone |
---|
| 253 | isoPhas => isotope%phase; nphas = isotope%nphas |
---|
| 254 | itZonIso => isotope%itZonIso; isoCheck = isotope%check |
---|
[4143] | 255 | iqIsoPha => isotope%iqIsoPha |
---|
[4120] | 256 | END FUNCTION isoSelectByIndex |
---|
| 257 | !============================================================================================================================== |
---|
| 258 | |
---|
| 259 | |
---|
[2320] | 260 | END MODULE infotrac_phy |
---|