Changeset 4482 for LMDZ6/branches/LMDZ_ECRad/libf/dyn3d/check_isotopes.F90
- Timestamp:
- Mar 29, 2023, 3:14:27 PM (18 months ago)
- Location:
- LMDZ6/branches/LMDZ_ECRad
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ_ECRad
- Property svn:mergeinfo changed
-
LMDZ6/branches/LMDZ_ECRad/libf/dyn3d/check_isotopes.F90
r4143 r4482 2 2 USE strings_mod, ONLY: maxlen, msg, strIdx, strStack, int2str, real2str 3 3 USE infotrac, ONLY: nqtot, niso, nphas, isotope, isoCheck, iqIsoPha, isoSelect, & 4 ntiso, iH2O, nzone, tracers, isoName, itZonIso, tnat4 ntiso, iH2O, nzone, tracers, isoName, itZonIso, getKey 5 5 IMPLICIT NONE 6 6 include "dimensions.h" … … 10 10 CHARACTER(LEN=maxlen) :: modname, msg1, nm(2) 11 11 INTEGER :: ixt, ipha, k, i, iq, iiso, izon, ieau, iqeau, iqpar 12 INTEGER, ALLOCATABLE :: ix(:) 12 INTEGER, ALLOCATABLE :: ix(:) 13 REAL, ALLOCATABLE, SAVE :: tnat(:) 13 14 REAL :: xtractot, xiiso, deltaD, q1, q2 14 15 REAL, PARAMETER :: borne = 1e19, & … … 28 29 IF(niso == 0) RETURN !--- No isotopes => finished 29 30 IF(first) THEN 30 iso_eau = strIdx(isoName,'H2[16]O') 31 iso_HDO = strIdx(isoName,'H[2]HO') 32 iso_O18 = strIdx(isoName,'H2[18]O') 33 iso_O17 = strIdx(isoName,'H2[17]O') 34 iso_HTO = strIdx(isoName,'H[3]HO') 31 iso_eau = strIdx(isoName,'H216O') 32 iso_HDO = strIdx(isoName,'HDO') 33 iso_O18 = strIdx(isoName,'H218O') 34 iso_O17 = strIdx(isoName,'H217O') 35 iso_HTO = strIdx(isoName,'HTO') 36 IF(getKey('tnat', tnat)) CALL abort_gcm(modname, 'missing isotopic parameter', 1) 35 37 first = .FALSE. 36 38 END IF
Note: See TracChangeset
for help on using the changeset viewer.