- 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/dyn3dmem/check_isotopes_loc.F90
r4143 r4482 3 3 USE strings_mod, ONLY: maxlen, msg, strIdx, strStack, int2str, real2str 4 4 USE infotrac, ONLY: nqtot, niso, nphas, isotope, isoCheck, iqIsoPha, isoSelect, & 5 ntiso, iH2O, nzone, tracers, isoName, itZonIso, tnat5 ntiso, iH2O, nzone, tracers, isoName, itZonIso, getKey 6 6 IMPLICIT NONE 7 7 include "dimensions.h" … … 11 11 CHARACTER(LEN=maxlen) :: modname, msg1, nm(2) 12 12 INTEGER :: ixt, ipha, k, i, iq, iiso, izon, ieau, iqeau, iqpar 13 INTEGER, ALLOCATABLE :: ix(:) 13 INTEGER, ALLOCATABLE :: ix(:) 14 REAL, ALLOCATABLE, SAVE :: tnat(:) !--- OpenMP shared variable 14 15 REAL :: xtractot, xiiso, deltaD, q1, q2 15 16 REAL, PARAMETER :: borne = 1e19, & … … 31 32 IF(first) THEN 32 33 !$OMP MASTER 33 iso_eau = strIdx(isoName,'H2[16]O') 34 iso_HDO = strIdx(isoName,'H[2]HO') 35 iso_O18 = strIdx(isoName,'H2[18]O') 36 iso_O17 = strIdx(isoName,'H2[17]O') 37 iso_HTO = strIdx(isoName,'H[3]HO') 34 iso_eau = strIdx(isoName,'H216O') 35 iso_HDO = strIdx(isoName,'HDO') 36 iso_O18 = strIdx(isoName,'H218O') 37 iso_O17 = strIdx(isoName,'H217O') 38 iso_HTO = strIdx(isoName,'HTO') 39 IF(getKey('tnat', tnat)) CALL abort_gcm(modname, 'missing isotopic parameter', 1) 38 40 !$OMP END MASTER 39 41 !$OMP BARRIER
Note: See TracChangeset
for help on using the changeset viewer.