Changeset 4399
- Timestamp:
- Jan 27, 2023, 12:35:36 AM (2 years ago)
- Location:
- LMDZ6/trunk
- Files:
-
- 6 edited
-
DefLists/isotopes_params.def (modified) (1 diff)
-
DefLists/tracer_RN_PB.def_iso (modified) (1 diff)
-
libf/dyn3d/check_isotopes.F90 (modified) (1 diff)
-
libf/dyn3dmem/check_isotopes_loc.F90 (modified) (1 diff)
-
libf/misc/readTracFiles_mod.f90 (modified) (1 diff)
-
libf/phylmdiso/isotopes_mod.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/DefLists/isotopes_params.def
r4193 r4399 1 1 &H2O 2 2 3 params 4 5 6 7 8 3 params deltaO18_oce=0.0 fac_coeff_eq17_liq=0.529 \ 4 fac_enrichoce18=0.0005 fac_coeff_eq17_ice=0.529 \ 5 talph1_O18=1137. talps1_O18=11.839 tkcin0_O18=0.006 \ 6 talph2_O18=-0.4156 talps2_O18=-0.028244 tkcin1_O18=0.000285 \ 7 talph3_O18=-2.0667E-3 tdifrel_O18=1./0.9723 tkcin2_O18=0.00082 \ 8 alpha_liq_sol_O18=1.00291 9 9 10 H2 [16]O alpha=1.000 toce=tnat \11 12 13 14 15 10 H216O alpha=1.000 toce=tnat \ 11 tnat=1.0000 tcorr=1.0 \ 12 talph1=0.0 talps1=0.0 tkcin0=0.0 \ 13 talph2=0.0 talps2=0.0 tkcin1=0.0 \ 14 talph3=0.0 tdifrel=1.0 tkcin2=0.0 \ 15 alpha_liq_sol=1.0 Rdefault=tnat*1.0 Rmethox=1.0 16 16 17 H2 [17]O alpha=1.003 toce=tnat*(1.0+deltaO18_oce/1000.0)^pente_MWL \18 19 20 21 22 23 24 17 H217O alpha=1.003 toce=tnat*(1.0+deltaO18_oce/1000.0)^pente_MWL \ 18 tnat=40.000e-6 tcorr=1.0+fac_enrichoce18*pente_MWL \ 19 talph1=talph1_O18 talps1=talps1_O18 tkcin0=tkcin0_O18*fac_kcin \ 20 talph2=talph2_O18 talps2=talps2_O18 tkcin1=tkcin1_O18*fac_kcin \ 21 talph3=talph3_O18 tdifrel=1./0.98555 tkcin2=tkcin2_O18*fac_kcin \ 22 alpha_liq_sol=alpha_liq_sol_O18^fac_coeff_eq17_liq \ 23 Rdefault=tnat*(1.0-3.15/1000.) Rmethox=tnat*(1.0+230./1000.) \ 24 fac_kcin=(tdifrel-1.0)/(tdifrel_O18-1.0) pente_MWL=0.528 25 25 26 H2 [18]O alpha=1.006 toce=tnat*(1.0+deltaO18_oce/1000.0) \27 28 29 30 31 32 26 H218O alpha=1.006 toce=tnat*(1.0+deltaO18_oce/1000.0) \ 27 tnat=2005.2e-6 tcorr=1.0+fac_enrichoce18 \ 28 talph1=talph1_O18 talps1=talps1_O18 tkcin0=tkcin0_O18 \ 29 talph2=talph2_O18 talps2=talps2_O18 tkcin1=tkcin1_O18 \ 30 talph3=talph3_O18 tdifrel=tdifrel_O18 tkcin2=tkcin2_O18 \ 31 alpha_liq_sol=alpha_liq_sol_O18 \ 32 Rdefault=tnat*(1.0-6.00/1000.) Rmethox=tnat*(1.0+130./1000.) 33 33 34 H [2]HOalpha=1.010 toce=tnat*(1.0+deltaO18_oce/1000.0*pente_MWL) \35 36 37 38 39 40 41 34 HDO alpha=1.010 toce=tnat*(1.0+deltaO18_oce/1000.0*pente_MWL) \ 35 tnat=155.76e-6 tcorr=1.0+fac_enrichoce18*pente_MWL \ 36 talph1=24844. talps1=16288. tkcin0=tkcin0_O18*fac_kcin \ 37 talph2=-76.248 talps2=-0.0934 tkcin1=tkcin1_O18*fac_kcin \ 38 talph3=52.612E-3 tdifrel=1./0.9755 tkcin2=tkcin2_O18*fac_kcin \ 39 alpha_liq_sol=1.0212 \ 40 Rdefault=tnat*(1.0+(10.0-6.0*pente_MWL)/1000.) Rmethox=tnat*(1.0-25.0/1000.) \ 41 fac_kcin=(tdifrel-1.0)/(tdifrel_O18-1.0) pente_MWL=8.0 42 42 43 H [3]HOalpha=1.000 toce=4.0E-19 \44 45 46 47 48 43 HTO alpha=1.000 toce=4.0E-19 \ 44 tnat=0.0000 tcorr=1.0 \ 45 talph1=46480. talps1=46480. tkcin0=0.01056 \ 46 talph2=-103.87 talps2=-103.87 tkcin1=0.0005016 \ 47 talph3=0.0 tdifrel=1./0.968 tkcin2=0.0014432 \ 48 alpha_liq_sol=1.0 Rdefault=0.0 Rmethox=0.0 49 49 -
LMDZ6/trunk/DefLists/tracer_RN_PB.def_iso
r4158 r4399 1 1 &version=1.0 2 2 &lmdz 3 default 4 H2O 5 H2O 6 H2 [18]O,H[2]HO,H2[16]Ophases=gls parent=H2O3 default type=tracer phases=g hadv=10 vadv=10 parent=air 4 H2O hadv=14 vadv=14 5 H2O phases=ls 6 H218O,HDO,H216O phases=gls parent=H2O 7 7 RN,PB -
LMDZ6/trunk/libf/dyn3d/check_isotopes.F90
r4367 r4399 29 29 IF(niso == 0) RETURN !--- No isotopes => finished 30 30 IF(first) THEN 31 iso_eau = strIdx(isoName,'H2 [16]O')32 iso_HDO = strIdx(isoName,'H [2]HO')33 iso_O18 = strIdx(isoName,'H2 [18]O')34 iso_O17 = strIdx(isoName,'H2 [17]O')35 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 36 IF(getKey('tnat', tnat)) CALL abort_gcm(modname, 'missing isotopic parameter', 1) 37 37 first = .FALSE. -
LMDZ6/trunk/libf/dyn3dmem/check_isotopes_loc.F90
r4367 r4399 32 32 IF(first) THEN 33 33 !$OMP MASTER 34 iso_eau = strIdx(isoName,'H2 [16]O')35 iso_HDO = strIdx(isoName,'H [2]HO')36 iso_O18 = strIdx(isoName,'H2 [18]O')37 iso_O17 = strIdx(isoName,'H2 [17]O')38 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 39 IF(getKey('tnat', tnat)) CALL abort_gcm(modname, 'missing isotopic parameter', 1) 40 40 !$OMP END MASTER -
LMDZ6/trunk/libf/misc/readTracFiles_mod.f90
r4394 r4399 125 125 126 126 !--- CORRESPONDANCE BETWEEN OLD AND NEW WATER NAMES 127 CHARACTER(LEN=maxlen), SAVE :: oldH2OIso(5) = ['eau', 'HDO', 'O18', 'O17', 'HTO']128 CHARACTER(LEN=maxlen), SAVE :: newH2OIso(5) = ['H2 [16]O', 'H[2]HO ', 'H2[18]O', 'H2[17]O', 'H[3]HO']127 CHARACTER(LEN=maxlen), SAVE :: oldH2OIso(5) = ['eau', 'HDO', 'O18', 'O17', 'HTO' ] 128 CHARACTER(LEN=maxlen), SAVE :: newH2OIso(5) = ['H216O', 'HDO ', 'H218O', 'H217O', 'HTO '] 129 129 130 130 !--- CORRESPONDANCE BETWEEN OLD AND NEW HNO3 RELATED SPECIES NAMES -
LMDZ6/trunk/libf/phylmdiso/isotopes_mod.F90
r4325 r4399 174 174 175 175 !--- Type of water isotopes: 176 iso_eau = strIdx(isoName, 'H2 [16]O'); CALL msg('iso_eau='//int2str(iso_eau), modname)177 iso_HDO = strIdx(isoName, 'H [2]HO');CALL msg('iso_HDO='//int2str(iso_HDO), modname)178 iso_O18 = strIdx(isoName, 'H2 [18]O'); CALL msg('iso_O18='//int2str(iso_O18), modname)179 iso_O17 = strIdx(isoName, 'H2 [17]O'); CALL msg('iso_O17='//int2str(iso_O17), modname)180 iso_HTO = strIdx(isoName, 'H [3]HO');CALL msg('iso_HTO='//int2str(iso_HTO), modname)176 iso_eau = strIdx(isoName, 'H216O'); CALL msg('iso_eau='//int2str(iso_eau), modname) 177 iso_HDO = strIdx(isoName, 'HDO'); CALL msg('iso_HDO='//int2str(iso_HDO), modname) 178 iso_O18 = strIdx(isoName, 'H218O'); CALL msg('iso_O18='//int2str(iso_O18), modname) 179 iso_O17 = strIdx(isoName, 'H217O'); CALL msg('iso_O17='//int2str(iso_O17), modname) 180 iso_HTO = strIdx(isoName, 'HTO'); CALL msg('iso_HTO='//int2str(iso_HTO), modname) 181 181 182 182 !--- Initialiaation: reading the isotopic parameters. … … 246 246 ! bugs quand temperature dans ascendances convs est mal calculee 247 247 CALL get_in('cond_temp_env', cond_temp_env, .FALSE.) 248 IF(ANY(isoName == 'H [3]HO')) &248 IF(ANY(isoName == 'HTO')) & 249 249 CALL get_in('ok_prod_nucl_tritium', ok_prod_nucl_tritium, .FALSE., .FALSE.) 250 250
Note: See TracChangeset
for help on using the changeset viewer.
