Changeset 5199 for LMDZ6/trunk/libf
- Timestamp:
- Sep 18, 2024, 11:22:37 PM (2 months ago)
- Location:
- LMDZ6/trunk/libf
- Files:
-
- 14 edited
-
dyn3d/dynetat0.F90 (modified) (5 diffs)
-
dyn3d/iniacademic.F90 (modified) (2 diffs)
-
dyn3d/qminimum.F (modified) (1 diff)
-
dyn3d_common/infotrac.F90 (modified) (18 diffs)
-
dyn3dmem/dynetat0_loc.F90 (modified) (4 diffs)
-
dyn3dmem/iniacademic_loc.F90 (modified) (2 diffs)
-
dyn3dmem/qminimum_loc.F (modified) (1 diff)
-
phylmd/infotrac_phy.F90 (modified) (18 diffs)
-
phylmd/phyetat0_mod.F90 (modified) (1 diff)
-
phylmd/physiq_mod.F90 (modified) (1 diff)
-
phylmdiso/isotopes_routines_mod.F90 (modified) (2 diffs)
-
phylmdiso/isotrac_mod.F90 (modified) (1 diff)
-
phylmdiso/phyetat0_mod.F90 (modified) (1 diff)
-
phylmdiso/physiq_mod.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/dyn3d/dynetat0.F90
r5190 r5199 6 6 ! Purpose: Initial state reading. 7 7 !------------------------------------------------------------------------------- 8 USE infotrac, ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName 8 USE infotrac, ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName, & 9 new2oldH2O, newHNO3, oldHNO3, getKey 9 10 USE strings_mod, ONLY: maxlen, msg, strStack, real2str, int2str 10 11 USE netcdf, ONLY: NF90_OPEN, NF90_NOWRITE, NF90_INQ_VARID, & 11 12 NF90_CLOSE, NF90_GET_VAR, NF90_NoErr 12 USE readTracFiles_mod, ONLY: new2oldH2O, newHNO3, oldHNO3, getKey13 13 USE control_mod, ONLY: planet_type 14 14 USE assert_eq_m, ONLY: assert_eq … … 116 116 var="temps" 117 117 IF(NF90_INQ_VARID(fID,var,vID)/=NF90_NoErr) THEN 118 CALL msg(' missing field <temps> ; trying with <Time>', modname)118 CALL msg('Missing field <temps> ; trying with <Time>', modname) 119 119 var="Time" 120 120 CALL err(NF90_INQ_VARID(fID,var,vID),"inq",var) … … 148 148 !-------------------------------------------------------------------------------------------------------------------------- 149 149 ELSE IF(NF90_INQ_VARID(fID, oldVar, vID) == NF90_NoErr) THEN !=== TRY WITH ALTERNATE NAME 150 CALL msg(' Tracer <'//TRIM(var)//'> is missing=> initialized to <'//TRIM(oldVar)//'>', modname)150 CALL msg('Missing tracer <'//TRIM(var)//'> => initialized to <'//TRIM(oldVar)//'>', modname) 151 151 CALL err(NF90_GET_VAR(fID,vID,q(:,:,:,iq)),"get",oldVar) 152 152 !-------------------------------------------------------------------------------------------------------------------------- … … 164 164 CALL abort_gcm(TRIM(modname), 'missing isotopic parameters', 1) 165 165 endif 166 CALL msg(' Tracer <'//TRIM(var)//'> is missing=> initialized with a simplified Rayleigh distillation law.', modname)166 CALL msg('Missing tracer <'//TRIM(var)//'> => initialized with a simplified Rayleigh distillation law.', modname) 167 167 q(:,:,:,iq) = q(:,:,:,iqParent)*tnat*(q(:,:,:,iqParent)/30.e-3)**(alpha_ideal-1.) 168 168 ELSE 169 CALL msg(' Tracer <'//TRIM(var)//'> is missing=> initialized to its parent isotope concentration.', modname)169 CALL msg('Missing tracer <'//TRIM(var)//'> => initialized to its parent isotope concentration.', modname) 170 170 ! Camille 9 mars 2023: attention!! seuls les tags qui correspondent à 171 171 ! izone=izone_init (définie dans isotrac_mod) sont initialisés comme … … 181 181 !-------------------------------------------------------------------------------------------------------------------------- 182 182 ELSE !=== MISSING: SET TO 0 183 CALL msg(' Tracer <'//TRIM(var)//'> is missing=> initialized to zero', modname)183 CALL msg('Missing tracer <'//TRIM(var)//'> => initialized to zero', modname) 184 184 q(:,:,:,iq)=0. 185 185 !-------------------------------------------------------------------------------------------------------------------------- -
LMDZ6/trunk/libf/dyn3d/iniacademic.F90
r5190 r5199 5 5 6 6 USE filtreg_mod, ONLY: inifilr 7 USE infotrac, ONLY: nqtot, niso, iqIsoPha, tracers, getKey, isoName 7 USE infotrac, ONLY: nqtot, niso, iqIsoPha, tracers, getKey, isoName, addPhase 8 8 USE control_mod, ONLY: day_step,planet_type 9 9 use exner_hyb_m, only: exner_hyb … … 21 21 USE temps_mod, ONLY: annee_ref, day_ini, day_ref 22 22 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 23 USE readTracFiles_mod, ONLY: addPhase24 23 use netcdf, only : NF90_NOWRITE,NF90_OPEN,NF90_NOERR,NF90_INQ_VARID 25 24 use netcdf, only : NF90_CLOSE, NF90_GET_VAR -
LMDZ6/trunk/libf/dyn3d/qminimum.F
r5190 r5199 4 4 SUBROUTINE qminimum( q,nqtot,deltap ) 5 5 6 USE infotrac, ONLY: niso, ntiso, iqIsoPha, tracers6 USE infotrac, ONLY: niso, ntiso, iqIsoPha, tracers, addPhase 7 7 USE strings_mod, ONLY: strIdx 8 USE readTracFiles_mod, ONLY: addPhase9 8 IMPLICIT none 10 9 c -
LMDZ6/trunk/libf/dyn3d_common/infotrac.F90
r5190 r5199 7 7 delPhase, niso, getKey, isot_type, processIsotopes, isotope, maxTableWidth, iqIsoPha, nphas, ixIso, isoPhas, & 8 8 addPhase, iH2O, addKey, isoSelect, testTracersFiles, isoKeys, indexUpdate, iqWIsoPha, nbIso, ntiso, isoName, isoCheck 9 USE readTracFiles_mod, ONLY: new2oldH2O, newHNO3, oldHNO3 9 10 IMPLICIT NONE 10 11 … … 16 17 PUBLIC :: nqtot, nbtr, nqo, nqCO2, nqtottr !--- Main dimensions 17 18 PUBLIC :: conv_flg, pbl_flg !--- Convection & boundary layer activation keys 19 PUBLIC :: new2oldH2O, newHNO3, oldHNO3 !--- For backwards compatibility in dynetat0 20 PUBLIC :: addPhase, delPhase !--- Add/remove the phase from the name of a tracer 18 21 19 22 !=== FOR ISOTOPES: General … … 21 24 PUBLIC :: isoSelect, ixIso !--- Isotopes family selection tool + selected family index 22 25 !=== FOR ISOTOPES: Specific to water 23 PUBLIC :: iH2O !--- H2O isotopes class index26 PUBLIC :: iH2O !--- Value of "ixIso" for "H2O" isotopes class 24 27 PUBLIC :: min_qParent, min_qMass, min_ratio !--- Min. values for various isotopic quantities 25 28 !=== FOR ISOTOPES: Depending on the selected isotopes family 26 PUBLIC :: isotope , isoKeys !--- Selected isotopes database + associated keys (cf.getKey)27 PUBLIC :: iso Name, isoZone, isoPhas !--- Isotopes andtagging zones names, phases28 PUBLIC :: niso, nzone, nphas, ntiso !--- " " numbers + isotopes & tagging tracers number29 PUBLIC :: itZonIso !--- i dx "it" (in "isoName(1:niso)") = function(tagging idx, isotope idx)30 PUBLIC :: iqIsoPha !--- i dx "iq" (in "qx") = function(isotope idx, phase idx) + aliases29 PUBLIC :: isotope !--- Selected isotopes database (argument of getKey) 30 PUBLIC :: isoKeys, isoName, isoZone, isoPhas !--- Isotopes keys & names, tagging zones names, phases 31 PUBLIC :: niso, ntiso, nzone, nphas !--- Number of " " 32 PUBLIC :: itZonIso !--- index "it" in "isoName(1:niso)" = f(tagging idx, isotope idx) 33 PUBLIC :: iqIsoPha !--- index "iq" in "qx" = f(isotope idx, phase idx) 31 34 PUBLIC :: isoCheck !--- Run isotopes checking routines 32 35 !=== FOR BOTH TRACERS AND ISOTOPES … … 36 39 ! |--------------------+-----------------------+-----------------+---------------+----------------------------| 37 40 ! | water in different | water tagging | water isotopes | other tracers | additional tracers moments | 38 ! | phases: H2O_[gls ]| isotopes | | | for higher order schemes |41 ! | phases: H2O_[glsrb]| isotopes | | | for higher order schemes | 39 42 ! |--------------------+-----------------------+-----------------+---------------+----------------------------| 40 43 ! | | | | | | … … 61 64 ! |-------------+------------------------------------------------------+-------------+------------------------+ 62 65 ! | name | Name (short) | tname | | 66 ! | keys | key/val pairs accessible with "getKey" routine | / | | 63 67 ! | gen0Name | Name of the 1st generation ancestor | / | | 64 68 ! | parent | Name of the parent | / | | 65 69 ! | longName | Long name (with adv. scheme suffix) for outputs | ttext | | 66 70 ! | type | Type (so far: tracer or tag) | / | tracer,tag | 67 ! | phase | Phases list ("g"as / "l"iquid / "s"olid) | / | [g][l][s] | 71 ! | phase | Phases list ("g"as / "l"iquid / "s"olid | | [g|l|s|r|b] | 72 ! | | "r"(cloud) / "b"lowing) | / | | 68 73 ! | component | Name(s) of the merged/cumulated section(s) | / | coma-separated names | 69 74 ! | iGeneration | Generation (>=1) | / | | … … 72 77 ! | nqDescen | Number of the descendants (all generations) | nqdesc | 1:nqtot | 73 78 ! | nqChildren | Number of childs (1st generation only) | nqfils | 1:nqtot | 74 ! | keys | key/val pairs accessible with "getKey" routine | / | |75 79 ! | iadv | Advection scheme number | iadv | 1,2,10-20(exc.15,19),30| 76 ! | isAdvected | advected tracers flag (.TRUE. if iadv >= 0) | / | nqtrue .TRUE. values |77 ! | isInPhysics | tracers not extracted from the main table in physics | / | nqtottr .TRUE. values |80 ! | isAdvected | Advected tracers flag (.TRUE. if iadv >= 0) | / | nqtrue .TRUE. values | 81 ! | isInPhysics | Tracers not extracted from the main table in physics | / | nqtottr .TRUE. values | 78 82 ! | iso_iGroup | Isotopes group index in isotopes(:) | / | 1:nbIso | 79 83 ! | iso_iName | Isotope name index in isotopes(iso_iGroup)%trac(:) | iso_indnum | 1:niso | … … 91 95 ! | trac | ntiso | Isotopes + tagging tracers list + number | / | ntraciso | | 92 96 ! | zone | nzone | Geographic tagging zones list + number | / | ntraceurs_zone | | 93 ! | phase | nphas | Phases list + number | | [g ][l][s], 1:3|97 ! | phase | nphas | Phases list + number | | [g|l|s|r|b] 1:5 | 94 98 ! | iqIsoPha | Index in "qx" = f(name(1:ntiso)),phas) | iqiso | 1:nqtot | 95 99 ! | itZonIso | Index in "trac(1:ntiso)"= f(zone, name(1:niso)) | index_trac | 1:ntiso | … … 99 103 100 104 !=== DIMENSIONS OF THE TRACERS TABLES AND OTHER SCALAR VARIABLES 101 INTEGER, SAVE :: nqtot, &!--- Tracers nb in dynamics (incl. higher moments + H2O)102 nbtr, &!--- Tracers nb in physics (excl. higher moments + H2O)103 nqo, &!--- Number of water phases105 INTEGER, SAVE :: nqtot, & !--- Tracers nb in dynamics (incl. higher moments + H2O) 106 nbtr, & !--- Tracers nb in physics (excl. higher moments + H2O) 107 nqo, & !--- Number of water phases 104 108 nqtottr, & !--- Number of tracers passed to phytrac (TO BE DELETED ?) 105 109 nqCO2 !--- Number of tracers of CO2 (ThL) 106 CHARACTER(LEN=maxlen), SAVE :: type_trac !--- Keyword for tracers type 110 CHARACTER(LEN=maxlen), SAVE :: type_trac !--- Keyword for tracers type(s) 107 111 108 112 !=== VARIABLES FOR INCA 109 INTEGER, SAVE, ALLOCATABLE :: conv_flg(:), & !--- Convection activation ; needed for INCA (nbtr)110 pbl_flg(:) !--- Boundary layer activation ; needed for INCA(nbtr)113 INTEGER, DIMENSION(:), SAVE, ALLOCATABLE :: & 114 conv_flg, pbl_flg !--- Convection / boundary layer activation (nbtr) 111 115 112 116 CONTAINS … … 115 119 USE control_mod, ONLY: planet_type 116 120 #ifdef REPROBUS 117 USE CHEM_REP, 121 USE CHEM_REP, ONLY: Init_chem_rep_trac 118 122 #endif 119 123 IMPLICIT NONE … … 160 164 TYPE(trac_type), ALLOCATABLE, TARGET :: ttr(:) 161 165 TYPE(trac_type), POINTER :: t1, t(:) 162 CHARACTER(LEN=maxlen), ALLOCATABLE :: types_trac(:) !--- Keyword for tracers type(s), parsed version 163 166 CHARACTER(LEN=maxlen), ALLOCATABLE :: types_trac(:) !--- Keywords for tracers type(s), parsed version 164 167 CHARACTER(LEN=*), PARAMETER :: modname="init_infotrac" 165 168 !------------------------------------------------------------------------------------------------------------------------------ … … 171 174 descrq(10:20) = ['VL1','VLP','FH1','FH2','VLH',' ','PPM','PPS','PPP',' ','SLP'] 172 175 descrq(30) = 'PRA' 173 174 CALL msg('type_trac = "'//TRIM(type_trac)//'"', modname)175 176 176 177 lerr=strParse(type_trac, '|', types_trac, n=nt) 177 178 IF (nt .GT. 1) THEN 178 179 IF (nt .GT. 2) CALL abort_gcm(modname, 'you need to modify type_trac, this version is not supported by lmdz', 1) 179 if(nt .EQ. 2) type_trac=types_trac(2)180 IF (nt .EQ. 2) type_trac=types_trac(2) 180 181 ENDIF 181 182 183 CALL msg('type_trac = "'//TRIM(type_trac)//'"', modname) 182 184 183 185 … … 213 215 214 216 !============================================================================================================================== 215 ! 1) Get the numbers of: true (first order only) tracers "nqtrue", water tracers "nqo" (vapor/liquid/solid)216 !============================================================================================================================== 217 texp = type_trac !=== EXPANDED VERSION OF "type_trac", WITH "|" SEPARATOR217 ! 0) DETERMINE THE TYPE OF THE INPUT TRACERS DESCRIPTION FILE AND READ IT 218 !============================================================================================================================== 219 texp = type_trac !=== EXPANDED (WITH "|" SEPARATOR) "type_trac" 218 220 IF(texp == 'inco') texp = 'co2i|inca' 219 221 IF(texp /= 'lmdz') texp = 'lmdz|'//TRIM(texp) 220 221 !=== DETERMINE THE TYPE OF THE INPUT TRACERS DESCRIPTION FILE222 222 IF(testTracersFiles(modname, texp, fType, .TRUE.)) CALL abort_gcm(modname, 'problem with tracers file(s)',1) 223 223 ttp = type_trac; IF(fType /= 1) ttp = texp 224 225 224 IF(readTracersFiles(ttp, lRepr=type_trac=='repr')) CALL abort_gcm(modname, 'problem with tracers file(s)',1) 226 !--------------------------------------------------------------------------------------------------------------------------- 227 IF(fType == 0) CALL abort_gcm(modname, 'Missing tracers file: "traceur.def", "tracer.def" or "tracer_<keyword>.def file.',1) 228 !--------------------------------------------------------------------------------------------------------------------------- 229 IF(fType == 1 .AND. ANY(['inca','inco']==type_trac)) THEN !=== FOUND OLD STYLE INCA "traceur.def" 225 226 !============================================================================================================================== 227 ! 1) Get various numbers: "nqtrue" (first order only tracers), "nqo" (water phases), 'nbtr' (tracers passed to physics), etc. 228 !============================================================================================================================== 229 !--------------------------------------------------------------------------------------------------------------------------- 230 IF(fType == 0) CALL abort_gcm(modname, 'Missing "traceur.def", "tracer.def" or "tracer_<keyword>.def tracers file.',1) 231 !--------------------------------------------------------------------------------------------------------------------------- 232 IF(fType == 1 .AND. ANY(['inca', 'inco'] == type_trac)) THEN !=== FOUND OLD STYLE INCA "traceur.def" 230 233 !--------------------------------------------------------------------------------------------------------------------------- 231 234 #ifdef INCA … … 264 267 ELSE !=== OTHER CASES (OLD OR NEW FORMAT, NO INCA MODULE) 265 268 !--------------------------------------------------------------------------------------------------------------------------- 266 nqo = COUNT(delPhase(tracers(:)%name) == 'H2O' & 267 .AND. tracers(:)%component == 'lmdz') !--- Number of water phases 268 nqtrue = SIZE(tracers) !--- Total number of "true" tracers 269 nbtr = nqtrue-COUNT(delPhase(tracers(:)%gen0Name) == 'H2O' & 270 .AND. tracers(:)%component == 'lmdz') !--- Number of tracers passed to phytrac 269 nqtrue = SIZE(tracers) !--- "true" tracers 270 nqo = COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%name) == 'H2O') !--- Water phases 271 nbtr = nqtrue-COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%gen0Name) == 'H2O') !--- Passed to phytrac 272 nqCO2 = COUNT( [type_trac == 'inco', type_trac == 'co2i'] ) 271 273 #ifdef INCA 272 nqINCA =COUNT(tracers(:)%component == 'inca')273 #endif 274 lerr = getKey('hadv', hadv, ky=tracers(:)%keys)275 lerr = getKey('vadv', vadv, ky=tracers(:)%keys)274 nqINCA = COUNT(tracers(:)%component == 'inca') 275 #endif 276 IF(getKey('hadv', hadv, ky=tracers(:)%keys)) CALL abort_gcm(modname, 'missing key "hadv"', 1) 277 IF(getKey('vadv', vadv, ky=tracers(:)%keys)) CALL abort_gcm(modname, 'missing key "vadv"', 1) 276 278 !--------------------------------------------------------------------------------------------------------------------------- 277 279 END IF … … 279 281 280 282 #ifdef REPROBUS 281 !--- Transfert the number of tracers to Reprobus282 283 CALL Init_chem_rep_trac(nbtr, nqo, tracers(:)%name) 283 284 #endif 284 #endif 285 285 286 !============================================================================================================================== 286 287 ! 2) Calculate nqtot, number of tracers needed (greater if advection schemes 20 or 30 have been chosen). … … 332 333 IF(iad == -1) CALL abort_gcm(modname, msg1, 1) 333 334 334 !--- SET FIELDS %longName, %iadv, %isAdvected, %isInPhysics335 !--- SET FIELDS longName, iadv, isAdvected, isInPhysics 335 336 t1%longName = t1%name; IF(iad > 0) t1%longName=TRIM(t1%name)//descrq(iad) 336 337 t1%iadv = iad 337 338 t1%isAdvected = iad >= 0 338 t1%isInPhysics= delPhase(t1%gen0Name) /= 'H2O' & 339 .OR. t1%component /= 'lmdz' !=== OTHER EXCEPTIONS TO BE ADDED: CO2i, SURSATURATED WATER CLOUD... 339 t1%isInPhysics= delPhase(t1%gen0Name) /= 'H2O' .OR. t1%component /= 'lmdz' !=== MORE EXCEPTIONS ? CO2i, SURSAT CLOUD H2O 340 340 ttr(iq) = t1 341 341 … … 347 347 ttr(jq+1:jq+nm) = t1 348 348 ttr(jq+1:jq+nm)%name = [ (TRIM(t1%name) //'-'//TRIM(suff(im)), im=1, nm) ] 349 ttr(jq+1:jq+nm)%gen0Name = [ (TRIM(t1%name) //'-'//TRIM(suff(im)), im=1, nm) ] 349 350 ttr(jq+1:jq+nm)%parent = [ (TRIM(t1%parent) //'-'//TRIM(suff(im)), im=1, nm) ] 350 351 ttr(jq+1:jq+nm)%longName = [ (TRIM(t1%longName)//'-'//TRIM(suff(im)), im=1, nm) ] … … 356 357 CALL MOVE_ALLOC(FROM=ttr, TO=tracers) 357 358 358 !--- SET FIELDS %iqParent, %nqChildren, %iGeneration, %iqDescen, %nqDescen359 !--- SET FIELDS iqParent, iqDescen, nqDescen, nqChildren 359 360 IF(indexUpdate(tracers)) CALL abort_gcm(modname, 'problem with tracers indices update', 1) 360 361 361 362 !=== TEST ADVECTION SCHEME 362 DO iq =1,nqtot ; t1 => tracers(iq); iad = t1%iadv363 DO iq = 1, nqtot ; t1 => tracers(iq); iad = t1%iadv 363 364 364 365 !--- ONLY TESTED VALUES FOR TRACERS FOR NOW: iadv = 14, 10 (and 0 for non-transported tracers) … … 405 406 #endif 406 407 t => tracers 407 CALL msg('Information stored in infotrac :', modname) 408 409 IF(dispTable('isssssssssiiiiiiiii', ['iq ', 'name', 'lNam', 'g0Nm', 'prnt', 'type', 'phas', 'comp', & 410 'isPh', 'isAd', 'iadv', 'iGen', 'iqPr', 'nqDe', 'nqCh', 'iGrp', 'iNam', 'iZon', 'iPha'], & 411 cat(t%name, t%longName, t%gen0Name, t%parent, t%type, t%phase, t%component, bool2str(t%isInPhysics), & 412 bool2str(t%isAdvected)), & 408 CALL msg('Information stored in '//TRIM(modname)//': ', modname) 409 IF(dispTable('isssssssiiiiiiiii', ['iq ', 'name', 'lNam', 'g0Nm', 'prnt', 'type', 'phas', 'comp', & 410 'iAdv', 'iGen', 'iqPr', 'nqDe', 'nqCh', 'iGrp', 'iNam', 'iZon', 'iPha'], & 411 cat(t%name, t%longName, t%gen0Name, t%parent, t%type, t%phase, t%component), & 413 412 cat([(iq, iq=1, nqtot)], t%iadv, t%iGeneration, t%iqParent, t%nqDescen, t%nqChildren, t%iso_iGroup, & 414 413 t%iso_iName, t%iso_iZone, t%iso_iPhase), nColMax=maxTableWidth, nHead=2, sub=modname)) & -
LMDZ6/trunk/libf/dyn3dmem/dynetat0_loc.F90
r5190 r5199 7 7 !------------------------------------------------------------------------------- 8 8 USE parallel_lmdz 9 USE infotrac, ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName 9 USE infotrac, ONLY: nqtot, tracers, niso, iqIsoPha, iH2O, isoName, & 10 new2oldH2O, newHNO3, oldHNO3, getKey 10 11 USE strings_mod, ONLY: maxlen, msg, strStack, real2str, int2str, strIdx 11 12 USE netcdf, ONLY: NF90_OPEN, NF90_NOWRITE, NF90_INQUIRE_DIMENSION, NF90_INQ_VARID, & 12 13 NF90_CLOSE, NF90_GET_VAR, NF90_INQUIRE_VARIABLE, NF90_NoErr 13 USE readTracFiles_mod, ONLY: new2oldH2O, newHNO3, oldHNO3, getKey14 14 USE control_mod, ONLY: planet_type 15 15 USE assert_eq_m, ONLY: assert_eq … … 173 173 !-------------------------------------------------------------------------------------------------------------------------- 174 174 ELSE IF(NF90_INQ_VARID(fID, oldVar, vID) == NF90_NoErr) THEN !=== TRY WITH ALTERNATE NAME 175 CALL msg(' Tracer <'//TRIM(var)//'> is missing=> initialized to <'//TRIM(oldVar)//'>', modname)175 CALL msg('Missing tracer <'//TRIM(var)//'> => initialized to <'//TRIM(oldVar)//'>', modname) 176 176 CALL get_var2(oldVar, q_glo); q(ijb_u:ije_u,:,iq)=q_glo(ijb_u:ije_u,:) 177 177 !-------------------------------------------------------------------------------------------------------------------------- … … 189 189 CALL abort_gcm(TRIM(modname), 'missing isotopic parameters', 1) 190 190 endif 191 CALL msg(' Tracer <'//TRIM(var)//'> is missing=> initialized with a simplified Rayleigh distillation law.', modname)191 CALL msg('Missing tracer <'//TRIM(var)//'> => initialized with a simplified Rayleigh distillation law.', modname) 192 192 q(ijb_u:ije_u,:,iq) = q(ijb_u:ije_u,:,iqParent)*tnat*(q(ijb_u:ije_u,:,iqParent)/30.e-3)**(alpha_ideal-1.) 193 193 ! Camille 9 mars 2023: point de vigilence: initialisation incohérente 194 194 ! avec celle de xt_ancien dans la physiq. 195 195 ELSE 196 CALL msg(' Tracer <'//TRIM(var)//'> is missing=> initialized to its parent isotope concentration.', modname)196 CALL msg('Missing tracer <'//TRIM(var)//'> => initialized to its parent isotope concentration.', modname) 197 197 ! Camille 9 mars 2023: attention!! seuls les tags qui correspondent à 198 198 ! izone=izone_init (définie dans isotrac_mod) sont initialisés comme … … 208 208 !-------------------------------------------------------------------------------------------------------------------------- 209 209 ELSE !=== MISSING: SET TO 0 210 CALL msg(' Tracer <'//TRIM(var)//'> is missing=> initialized to zero', modname)210 CALL msg('missing tracer <'//TRIM(var)//'> => initialized to zero', modname) 211 211 q(ijb_u:ije_u,:,iq)=0. 212 212 !-------------------------------------------------------------------------------------------------------------------------- -
LMDZ6/trunk/libf/dyn3dmem/iniacademic_loc.F90
r5190 r5199 5 5 6 6 USE filtreg_mod, ONLY: inifilr 7 USE infotrac, ONLY: nqtot, niso, iqIsoPha, tracers, getKey, isoName 7 USE infotrac, ONLY: nqtot, niso, iqIsoPha, tracers, getKey, isoName, addPhase 8 8 USE control_mod, ONLY: day_step,planet_type 9 9 use exner_hyb_m, only: exner_hyb … … 22 22 USE temps_mod, ONLY: annee_ref, day_ini, day_ref 23 23 USE ener_mod, ONLY: etot0,ptot0,ztot0,stot0,ang0 24 USE readTracFiles_mod, ONLY: addPhase25 24 use netcdf, only : NF90_NOWRITE,NF90_OPEN,NF90_NOERR,NF90_INQ_VARID 26 25 use netcdf, only : NF90_CLOSE, NF90_GET_VAR -
LMDZ6/trunk/libf/dyn3dmem/qminimum_loc.F
r5190 r5199 4 4 SUBROUTINE qminimum_loc( q,nqtot,deltap ) 5 5 USE parallel_lmdz 6 USE infotrac, ONLY: niso, ntiso, iqIsoPha, tracers, 6 USE infotrac, ONLY: niso, ntiso, iqIsoPha, tracers, addPhase, 7 7 & isoCheck, min_qParent 8 8 USE strings_mod, ONLY: strIdx 9 USE readTracFiles_mod, ONLY: addPhase10 9 IMPLICIT none 11 10 c -
LMDZ6/trunk/libf/phylmd/infotrac_phy.F90
r5190 r5199 7 7 delPhase, niso, getKey, isot_type, processIsotopes, isotope, maxTableWidth, iqIsoPha, nphas, ixIso, isoPhas, & 8 8 addPhase, iH2O, addKey, isoSelect, testTracersFiles, isoKeys, indexUpdate, iqWIsoPha, nbIso, ntiso, isoName, isoCheck 9 USE readTracFiles_mod, ONLY: new2oldH2O 9 10 IMPLICIT NONE 10 11 … … 16 17 PUBLIC :: nqtot, nbtr, nqo, nqCO2, nqtottr !--- Main dimensions 17 18 PUBLIC :: conv_flg, pbl_flg !--- Convection & boundary layer activation keys 18 #ifdef CPP_StratAer 19 PUBLIC :: new2oldH2O !--- For backwards compatibility in phyetat0 20 PUBLIC :: addPhase, delPhase !--- Add/remove the phase from the name of a tracer 21 #if defined CPP_StratAer || defined REPROBUS 19 22 PUBLIC :: nbtr_bin, nbtr_sulgas !--- Number of aerosols bins and sulfur gases for StratAer model 20 23 PUBLIC :: id_OCS_strat, id_SO2_strat, id_H2SO4_strat, id_BIN01_strat, id_TEST_strat 21 24 #endif 22 25 23 !=== FOR WATER24 PUBLIC :: ivap, iliq, isol25 26 !=== FOR ISOTOPES: General 26 27 PUBLIC :: isot_type, nbIso !--- Derived type, full isotopes families database + nb of families 27 28 PUBLIC :: isoSelect, ixIso !--- Isotopes family selection tool + selected family index 28 29 !=== FOR ISOTOPES: Specific to water 29 PUBLIC :: iH2O !--- H2O isotopes class index 30 PUBLIC :: iH2O !--- Value of "ixIso" for "H2O" isotopes class 31 PUBLIC :: ivap, iliq, isol 30 32 !=== FOR ISOTOPES: Depending on the selected isotopes family 31 PUBLIC :: isotope, isoKeys !--- Selected isotopes database + associated keys (cf. getKey) 32 PUBLIC :: isoName, isoZone, isoPhas !--- Isotopes and tagging zones names, phases 33 PUBLIC :: niso, nzone, nphas, ntiso !--- " " numbers + isotopes & tagging tracers number 34 PUBLIC :: itZonIso !--- idx "it" (in "isoName(1:niso)") = function(tagging idx, isotope idx) 35 PUBLIC :: iqIsoPha !--- idx "iq" (in "qx") = function(isotope idx, phase idx) + aliases 36 PUBLIC :: iqWIsoPha !--- Same as iqIsoPha but with normal water phases 37 33 PUBLIC :: isotope !--- Selected isotopes database (argument of getKey) 34 PUBLIC :: isoKeys, isoName, isoZone, isoPhas !--- Isotopes keys & names, tagging zones names, phases 35 PUBLIC :: niso, ntiso, nzone, nphas !--- Number of " " 36 PUBLIC :: itZonIso !--- Index "it" in "isoName(1:niso)" = f(tagging idx, isotope idx) 37 PUBLIC :: iqIsoPha !--- Index "iq" in "qx" = f(isotope idx, phase idx) 38 PUBLIC :: iqWIsoPha !--- Same as iqIsoPha but with normal water phases 38 39 PUBLIC :: isoCheck !--- Run isotopes checking routines 39 40 !=== FOR BOTH TRACERS AND ISOTOPES … … 43 44 ! |--------------------+-----------------------+-----------------+---------------+----------------------------| 44 45 ! | water in different | water tagging | water isotopes | other tracers | additional tracers moments | 45 ! | phases: H2O_[gls ]| isotopes | | | for higher order schemes |46 ! | phases: H2O_[glsrb]| isotopes | | | for higher order schemes | 46 47 ! |--------------------+-----------------------+-----------------+---------------+----------------------------| 47 48 ! | | | | | | … … 68 69 ! |-------------+------------------------------------------------------+-------------+------------------------+ 69 70 ! | name | Name (short) | tname | | 71 ! | keys | key/val pairs accessible with "getKey" routine | / | | 70 72 ! | gen0Name | Name of the 1st generation ancestor | / | | 71 73 ! | parent | Name of the parent | / | | 72 74 ! | longName | Long name (with adv. scheme suffix) for outputs | ttext | | 73 75 ! | type | Type (so far: tracer or tag) | / | tracer,tag | 74 ! | phase | Phases list ("g"as / "l"iquid / "s"olid) | / | [g][l][s] | 76 ! | phase | Phases list ("g"as / "l"iquid / "s"olid | | [g|l|s|r|b] | 77 ! | | "r"(cloud) / "b"lowing) | / | | 75 78 ! | component | Name(s) of the merged/cumulated section(s) | / | coma-separated names | 76 79 ! | iGeneration | Generation (>=1) | / | | … … 79 82 ! | nqDescen | Number of the descendants (all generations) | nqdesc | 1:nqtot | 80 83 ! | nqChildren | Number of childs (1st generation only) | nqfils | 1:nqtot | 81 ! | keys | key/val pairs accessible with "getKey" routine | / | | 82 ! | isAdvected | advected tracers flag (.TRUE. if iadv >= 0) | / | nqtrue .TRUE. values | 83 ! | isInPhysics | tracers not extracted from the main table in physics | / | nqtottr .TRUE. values | 84 ! | isAdvected | Advected tracers flag (.TRUE. if iadv >= 0) | / | nqtrue .TRUE. values | 85 ! | isInPhysics | Tracers not extracted from the main table in physics | / | nqtottr .TRUE. values | 84 86 ! | iso_iGroup | Isotopes group index in isotopes(:) | / | 1:nbIso | 85 87 ! | iso_iName | Isotope name index in isotopes(iso_iGroup)%trac(:) | iso_indnum | 1:niso | … … 97 99 ! | trac | ntiso | Isotopes + tagging tracers list + number | / | ntraciso | | 98 100 ! | zone | nzone | Geographic tagging zones list + number | / | ntraceurs_zone | | 99 ! | phase | nphas | Phases list + number | | [g ][l][s], 1:3|101 ! | phase | nphas | Phases list + number | | [g|l|s|r|b] 1:5 | 100 102 ! | iqIsoPha | Index in "qx" = f(name(1:ntiso)),phas) | iqiso | 1:nqtot | 101 ! | iqWIsoPha | Index in "qx" = f(name(1:ntiso)),phas) | iqiso| 1:nqtot |103 ! | iqWIsoPha | Index in "qx" = f(name(1:ntiso+nqo)),phas) | / | 1:nqtot | 102 104 ! | itZonIso | Index in "trac(1:ntiso)"= f(zone, name(1:niso)) | index_trac | 1:ntiso | 103 105 ! +-----------------+--------------------------------------------------+--------------------+-----------------+ 104 106 107 !=== INDICES FOR WATER 108 INTEGER, SAVE :: ivap, iliq, isol 109 !$OMP THREADPRIVATE(ivap, iliq, isol) 110 105 111 !=== DIMENSIONS OF THE TRACERS TABLES AND OTHER SCALAR VARIABLES 106 INTEGER, SAVE :: nqtot, &!--- Tracers nb in dynamics (incl. higher moments + H2O)107 nbtr, &!--- Tracers nb in physics (excl. higher moments + H2O)108 nqo, &!--- Number of water phases112 INTEGER, SAVE :: nqtot, & !--- Tracers nb in dynamics (incl. higher moments + H2O) 113 nbtr, & !--- Tracers nb in physics (excl. higher moments + H2O) 114 nqo, & !--- Number of water phases 109 115 nqtottr, & !--- Number of tracers passed to phytrac (TO BE DELETED ?) 110 116 nqCO2 !--- Number of tracers of CO2 (ThL) … … 112 118 !$OMP THREADPRIVATE(nqtot, nbtr, nqo, nqtottr, nqCO2, type_trac) 113 119 114 !=== INDICES OF WATER115 INTEGER, SAVE :: ivap,iliq,isol ! Indices for vap, liq and ice116 !$OMP THREADPRIVATE(ivap,iliq,isol)117 118 120 !=== VARIABLES FOR INCA 119 INTEGER, SAVE, ALLOCATABLE :: conv_flg(:), & !--- Convection activation ; needed for INCA (nbtr)120 pbl_flg(:) !--- Boundary layer activation ; needed for INCA(nbtr)121 INTEGER, DIMENSION(:), SAVE, ALLOCATABLE :: & 122 conv_flg, pbl_flg !--- Convection / boundary layer activation (nbtr) 121 123 !$OMP THREADPRIVATE(conv_flg, pbl_flg) 122 124 123 #if def CPP_StratAer125 #if defined CPP_StratAer || defined REPROBUS 124 126 !=== SPECIFIC TO STRATOSPHERIC AEROSOLS (CK/OB) 125 127 INTEGER, SAVE :: nbtr_bin, nbtr_sulgas !--- number of aerosols bins and sulfur gases for StratAer model … … 173 175 CHARACTER(LEN=2) :: suff(9) !--- Suffixes for schemes of order 3 or 4 (Prather) 174 176 CHARACTER(LEN=3) :: descrq(30) !--- Advection scheme description tags 175 CHARACTER(LEN=maxlen) :: msg1, texp, ttp !--- String for messages and expanded tracers type177 CHARACTER(LEN=maxlen) :: msg1, texp, ttp !--- Strings for messages and expanded tracers type 176 178 INTEGER :: fType !--- Tracers description file type ; 0: none 177 179 !--- 1/2/3: "traceur.def"/"tracer.def"/"tracer_*.def" 178 180 INTEGER :: nqtrue !--- Tracers nb from tracer.def (no higher order moments) 179 181 INTEGER :: iad !--- Advection scheme number 180 INTEGER :: iq, jq, nt, im, nm , k!--- Indexes and temporary variables182 INTEGER :: iq, jq, nt, im, nm !--- Indexes and temporary variables 181 183 LOGICAL :: lerr, lInit 182 184 TYPE(trac_type), ALLOCATABLE, TARGET :: ttr(:) 183 185 TYPE(trac_type), POINTER :: t1, t(:) 184 CHARACTER(LEN=maxlen), ALLOCATABLE :: types_trac(:) !--- Keyword for tracers type(s), parsed version 185 186 CHARACTER(LEN=maxlen), ALLOCATABLE :: types_trac(:) !--- Keywords for tracers type(s), parsed version 186 187 CHARACTER(LEN=*), PARAMETER :: modname="init_infotrac_phy" 187 188 !------------------------------------------------------------------------------------------------------------------------------ … … 202 203 ENDIF 203 204 204 205 205 CALL msg('type_trac = "'//TRIM(type_trac)//'"', modname) 206 206 lInit = .NOT.ALLOCATED(tracers) … … 243 243 244 244 !============================================================================================================================== 245 ! 1) Get the numbers of: true (first order only) tracers "nqtrue", water tracers "nqo" (vapor/liquid/solid)246 !============================================================================================================================== 247 texp = type_trac !=== EXPANDED VERSION OF "type_trac", WITH "|" SEPARATOR245 ! 0) DETERMINE THE TYPE OF THE INPUT TRACERS DESCRIPTION FILE AND READ IT 246 !============================================================================================================================== 247 texp = type_trac !=== EXPANDED (WITH "|" SEPARATOR) "type_trac" 248 248 IF(texp == 'inco') texp = 'co2i|inca' 249 249 IF(texp /= 'lmdz') texp = 'lmdz|'//TRIM(texp) 250 251 !=== DETERMINE THE TYPE OF THE INPUT TRACERS DESCRIPTION FILE252 250 IF(testTracersFiles(modname, texp, fType, lInit)) CALL abort_physic(modname, 'problem with tracers file(s)',1) 253 254 251 ttp = type_trac; IF(fType /= 1) ttp = texp 255 252 … … 262 259 !############################################################################################################################## 263 260 264 !--------------------------------------------------------------------------------------------------------------------------- 265 IF(fType == 0) CALL abort_physic(modname, 'Missing tracers file: "traceur.def", "tracer.def" or "tracer_<keyword>.def file.',1) 266 !--------------------------------------------------------------------------------------------------------------------------- 267 IF(fType == 1 .AND. ANY(['inca','inco']==type_trac) .AND. lInit) THEN !=== FOUND OLD STYLE INCA "traceur.def" 261 !============================================================================================================================== 262 ! 1) Get various numbers: "nqtrue" (first order only tracers), "nqo" (water phases), 'nbtr' (tracers passed to physics), etc. 263 !============================================================================================================================== 264 !--------------------------------------------------------------------------------------------------------------------------- 265 IF(fType == 0) CALL abort_physic(modname, 'Missing "traceur.def", "tracer.def" or "tracer_<keyword>.def tracers file.',1) 266 !--------------------------------------------------------------------------------------------------------------------------- 267 IF(fType == 1 .AND. ANY(['inca', 'inco'] == type_trac) .AND. lInit) THEN !=== FOUND OLD STYLE INCA "traceur.def" 268 268 !--------------------------------------------------------------------------------------------------------------------------- 269 269 #ifdef INCA … … 281 281 ttr(1+nqo:nqCO2+nqo )%component = 'co2i' 282 282 ttr(1+nqo+nqCO2:nqtrue)%component = 'inca' 283 ttr(1+nqo :nqtrue)%name = [('CO2 ', k=1, nqCO2), solsym_inca]283 ttr(1+nqo :nqtrue)%name = [('CO2 ', iq=1, nqCO2), solsym_inca] 284 284 ttr(1+nqo+nqCO2:nqtrue)%parent = tran0 285 285 ttr(1+nqo+nqCO2:nqtrue)%phase = 'g' … … 302 302 ELSE !=== OTHER CASES (OLD OR NEW FORMAT, NO INCA MODULE) 303 303 !--------------------------------------------------------------------------------------------------------------------------- 304 nqo = COUNT(delPhase(tracers(:)%name) == 'H2O' & 305 .AND. tracers(:)%component == 'lmdz') !--- Number of water phases 306 nqtrue = SIZE(tracers) !--- Total number of "true" tracers 307 nbtr = nqtrue-COUNT(delPhase(tracers(:)%gen0Name) == 'H2O' & 308 .AND. tracers(:)%component == 'lmdz') !--- Number of tracers passed to phytrac 304 nqtrue = SIZE(tracers) !--- "true" tracers 305 nqo = COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%name) == 'H2O') !--- Water phases 306 nbtr = nqtrue-COUNT(tracers(:)%component == 'lmdz' .AND. delPhase(tracers(:)%gen0Name) == 'H2O') !--- Passed to phytrac 307 nqCO2 = COUNT( [type_trac == 'inco', type_trac == 'co2i'] ) 309 308 #ifdef INCA 310 nqINCA = COUNT(tracers(:)%component == 'inca') 311 #endif 312 lerr = getKey('hadv', hadv, ky=tracers(:)%keys) 313 lerr = getKey('vadv', vadv, ky=tracers(:)%keys) 314 !--------------------------------------------------------------------------------------------------------------------------- 315 END IF 316 !--------------------------------------------------------------------------------------------------------------------------- 317 318 !--- Transfert the number of tracers to Reprobus 309 nqINCA = COUNT(tracers(:)%component == 'inca') 310 #endif 311 IF(getKey('hadv', hadv, ky=tracers(:)%keys)) CALL abort_physic(modname, 'missing key "hadv"', 1) 312 IF(getKey('vadv', vadv, ky=tracers(:)%keys)) CALL abort_physic(modname, 'missing key "vadv"', 1) 313 !--------------------------------------------------------------------------------------------------------------------------- 314 END IF 315 !--------------------------------------------------------------------------------------------------------------------------- 316 319 317 #ifdef REPROBUS 320 318 CALL Init_chem_rep_trac(nbtr, nqo, tracers(:)%name) … … 359 357 IF(iad == -1) CALL abort_physic(modname, msg1, 1) 360 358 361 !--- SET FIELDS %longName, %isAdvected, %isInPhysics359 !--- SET FIELDS longName, isAdvected, isInPhysics 362 360 t1%longName = t1%name; IF(iad > 0) t1%longName=TRIM(t1%name)//descrq(iad) 363 361 t1%isAdvected = iad >= 0 364 t1%isInPhysics= delPhase(t1%gen0Name) /= 'H2O' & 365 .OR. t1%component /= 'lmdz' !=== OTHER EXCEPTIONS TO BE ADDED: CO2i, SURSATURATED WATER CLOUD... 362 t1%isInPhysics= delPhase(t1%gen0Name) /= 'H2O' .OR. t1%component /= 'lmdz' !=== MORE EXCEPTIONS ? CO2i, SURSAT CLOUD H2O 366 363 ttr(iq) = t1 367 364 … … 372 369 IF(nm == 0) CYCLE !--- No higher moments 373 370 ttr(jq+1:jq+nm) = t1 374 ttr(jq+1:jq+nm)%name = [(TRIM(t1%name) //'-'//TRIM(suff(im)), im=1, nm) ] 375 ttr(jq+1:jq+nm)%parent = [(TRIM(t1%parent) //'-'//TRIM(suff(im)), im=1, nm) ] 376 ttr(jq+1:jq+nm)%longName = [(TRIM(t1%longName)//'-'//TRIM(suff(im)), im=1, nm) ] 377 ttr(jq+1:jq+nm)%isAdvected = [(.FALSE., im=1, nm) ] 371 ttr(jq+1:jq+nm)%name = [ (TRIM(t1%name) //'-'//TRIM(suff(im)), im=1, nm) ] 372 ttr(jq+1:jq+nm)%gen0Name = [ (TRIM(t1%name) //'-'//TRIM(suff(im)), im=1, nm) ] 373 ttr(jq+1:jq+nm)%parent = [ (TRIM(t1%parent) //'-'//TRIM(suff(im)), im=1, nm) ] 374 ttr(jq+1:jq+nm)%longName = [ (TRIM(t1%longName)//'-'//TRIM(suff(im)), im=1, nm) ] 375 ttr(jq+1:jq+nm)%isAdvected = [ (.FALSE., im=1, nm) ] 376 ttr(jq+1:jq+nm)%isInPhysics = [ (.FALSE., im=1, nm) ] 378 377 jq = jq + nm 379 378 END DO … … 381 380 CALL MOVE_ALLOC(FROM=ttr, TO=tracers) 382 381 383 !--- SET FIELDS %iqParent, %nqChildren, %iGeneration, %iqDescen, %nqDescen384 IF(indexUpdate(tracers)) CALL abort_ gcm(modname, 'problem when processing isotopes parameters', 1)382 !--- SET FIELDS iqParent, iqDescen, nqDescen, nqChildren, iGeneration 383 IF(indexUpdate(tracers)) CALL abort_physic(modname, 'problem with tracers indices update', 1) 385 384 386 385 !############################################################################################################################## … … 404 403 !############################################################################################################################## 405 404 !--- Convection / boundary layer activation for all tracers 406 IF 407 IF (.NOT.ALLOCATED(pbl_flg)) ALLOCATE( pbl_flg(nbtr)); pbl_flg(1:nbtr) = 1405 IF(.NOT.ALLOCATED(conv_flg)) ALLOCATE(conv_flg(nbtr)); conv_flg(1:nbtr) = 1 406 IF(.NOT.ALLOCATED( pbl_flg)) ALLOCATE( pbl_flg(nbtr)); pbl_flg(1:nbtr) = 1 408 407 409 408 !--- Note: nqtottr can differ from nbtr when nmom/=0 … … 424 423 #endif 425 424 t => tracers 426 CALL msg('Information stored in infotrac_phy :', modname) 427 IF(dispTable('issssssssiiiiiiii', ['iq ', 'name', 'lNam', 'g0Nm', 'prnt', 'type', 'phas', 'comp', & 428 'isPh', 'iGen', 'iqPr', 'nqDe', 'nqCh', 'iGrp', 'iNam', 'iZon', 'iPha'], & 429 cat(t%name, t%longName, t%gen0Name, t%parent, t%type, t%phase, t%component, bool2str(t%isInPhysics)),& 425 CALL msg('Information stored in '//TRIM(modname)//': ', modname) 426 IF(dispTable('isssssssssiiiiiiii', ['iq ', 'name', 'lNam', 'g0Nm', 'prnt', 'type', 'phas', 'comp', & 427 'isPh', 'isAd', 'iGen', 'iqPr', 'nqDe', 'nqCh', 'iGrp', 'iNam', 'iZon', 'iPha'], & 428 cat(t%name, t%longName, t%gen0Name, t%parent, t%type, t%phase, t%component, & 429 bool2str(t%isInPhysics), bool2str(t%isAdvected)), & 430 430 cat([(iq, iq=1, nqtot)], t%iGeneration, t%iqParent, t%nqDescen, t%nqChildren, t%iso_iGroup, & 431 431 t%iso_iName, t%iso_iZone, t%iso_iPhase), nColMax=maxTableWidth, nHead=2, sub=modname)) & -
LMDZ6/trunk/libf/phylmd/phyetat0_mod.F90
r5190 r5199 32 32 USE geometry_mod, ONLY: longitude_deg, latitude_deg 33 33 USE iostart, ONLY: close_startphy, get_field, get_var, open_startphy 34 USE infotrac_phy, ONLY: nqtot, nbtr, type_trac, tracers 35 USE readTracFiles_mod,ONLY: maxlen, new2oldH2O34 USE infotrac_phy, ONLY: nqtot, nbtr, type_trac, tracers, new2oldH2O 35 USE strings_mod, ONLY: maxlen 36 36 USE traclmdz_mod, ONLY: traclmdz_from_restart 37 37 USE carbon_cycle_mod, ONLY: carbon_cycle_init, carbon_cycle_cpl, carbon_cycle_tr, carbon_cycle_rad, co2_send, RCO2_glo -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r5190 r5199 39 39 USE ioipsl_getin_p_mod, ONLY : getin_p 40 40 USE indice_sol_mod 41 USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, type_trac 42 USE readTracFiles_mod, ONLY: addPhase 41 USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, type_trac, addPhase 43 42 USE strings_mod, ONLY: strIdx 44 43 USE iophy -
LMDZ6/trunk/libf/phylmdiso/isotopes_routines_mod.F90
r5190 r5199 16419 16419 USE isotopes_mod, ONLY: isoName,iso_HDO,iso_eau 16420 16420 USE phyetat0_get_mod, ONLY: phyetat0_get, phyetat0_srf 16421 USE readTracFiles_mod,ONLY: new2oldH2O16422 USE strings_mod, ONLY: strIdx, str Tail, maxlen, msg, int2str16421 USE infotrac_phy, ONLY: new2oldH2O 16422 USE strings_mod, ONLY: strIdx, strHead, strTail, maxlen, msg, int2str 16423 16423 #ifdef ISOVERIF 16424 16424 USE isotopes_verif_mod … … 16459 16459 outiso = isoName(ixt) 16460 16460 oldIso = strTail(new2oldH2O(outiso), '_') !--- Remove "H2O_" from "H2O_<iso>[_<tag>]" 16461 i = INDEX(outiso, '_', .TRUE.) 16462 oldIso2 = outiso(1:i-1)//outiso(i+1:LEN_TRIM(outiso)) ! CR 2023: on ajoute cette possibilité aussi, elle correspond au cas le plus récent. 16461 oldIso2= TRIM(strHead(outiso,'_'))//strTail(outiso,'_') ! CR 2023: most recent possibility 16463 16462 ! write(*,*) 'tmp 16541:' 16464 16463 ! write(*,*) 'outiso=',outiso -
LMDZ6/trunk/libf/phylmdiso/isotrac_mod.F90
r5190 r5199 3 3 4 4 MODULE isotrac_mod 5 USE infotrac_phy, ONLY: niso, ntiso, nzone 6 USE readTracFiles_mod, ONLY: delPhase 7 USE isotopes_mod, ONLY: ridicule, get_in 5 USE infotrac_phy, ONLY: niso, ntiso, nzone, delPhase 6 USE isotopes_mod, ONLY: ridicule, get_in 8 7 9 8 IMPLICIT NONE -
LMDZ6/trunk/libf/phylmdiso/phyetat0_mod.F90
r5190 r5199 40 40 USE geometry_mod, ONLY: longitude_deg, latitude_deg 41 41 USE iostart, ONLY: close_startphy, get_field, get_var, open_startphy 42 USE infotrac_phy, ONLY: nqtot, nbtr, type_trac, tracers 43 USE readTracFiles_mod,ONLY: maxlen, new2oldH2O42 USE infotrac_phy, ONLY: nqtot, nbtr, type_trac, tracers, new2oldH2O 43 USE strings_mod, ONLY: maxlen 44 44 USE traclmdz_mod, ONLY: traclmdz_from_restart 45 45 USE carbon_cycle_mod, ONLY: carbon_cycle_init, carbon_cycle_cpl, carbon_cycle_tr, carbon_cycle_rad, co2_send, RCO2_glo -
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r5190 r5199 39 39 USE ioipsl_getin_p_mod, ONLY : getin_p 40 40 USE indice_sol_mod 41 USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, type_trac,ivap,iliq,isol 42 USE readTracFiles_mod, ONLY: addPhase 41 USE infotrac_phy, ONLY: nqtot, nbtr, nqo, tracers, type_trac,addPhase, ivap, iliq, isol 43 42 USE strings_mod, ONLY: strIdx 44 43 USE iophy
Note: See TracChangeset
for help on using the changeset viewer.
