|
Last change
on this file since 4145 was
3985,
checked in by dcugnet, 3 years ago
|
- fix of the delPhase function.
- getvar1 and getvar2 fixed and modified to avoid the usage of files with
several time records and make the calls rather short.
- works again with iadv==0
- no more issues with tracers numbers (nqo, nqtot, etc.)
- fixes in the algebrical reduction routine used for "isotopes_parems.def"
(containing simple expressions with variables that have to be substituted).
- still to be validated numerically
|
|
File size:
2.8 KB
|
| Line | |
|---|
| 1 | &H2O |
|---|
| 2 | params deltaO18_oce=0.0 fac_enrichoce18=0.0005 talph1_O18=1137. talph2_O18=-0.4156 talph3_O18=-2.0667E-3 talps1_O18=11.839 talps2_O18=-0.028244 tkcin0_O18=0.006 tkcin1_O18=0.000285 tkcin2_O18=0.00082 tdifrel_O18=1./0.9723 alpha_liq_sol_O18=1.00291 fac_coeff_eq17_liq=0.529 fac_coeff_eq17_ice=0.529 |
|---|
| 3 | H2[16]O alpha=1.000 tnat=1.0000 toce=tnat tcorr=1.0 talph1=0.0 talph2=0.0 talph3=0.0 talps1=0.0 talps2=0.0 tkcin0=0.0 tkcin1=0.0 tkcin2=0.0 tdifrel=1.0 alpha_liq_sol=1.0 Rdefault=tnat*1.0 Rmethox=1.0 |
|---|
| 4 | H2[17]O alpha=1.003 tnat=40.000e-6 toce=tnat*(1.0+deltaO18_oce/1000.0)^pente_MWL tcorr=1.0+fac_enrichoce18*pente_MWL talph1=talph1_O18 talph2=talph2_O18 talph3=talph3_O18 talps1=talps1_O18 talps2=talps2_O18 tkcin0=tkcin0_O18*fac_kcin tkcin1=tkcin1_O18*fac_kcin tkcin2=tkcin2_O18*fac_kcin tdifrel=1./0.98555 alpha_liq_sol=alpha_liq_sol_O18^fac_coeff_eq17_liq Rdefault=tnat*(1.0-3.15/1000.) Rmethox=tnat*(1.0+230./1000.) pente_MWL=0.528 fac_kcin=(tdifrel-1.0)/(tdifrel_O18-1.0) |
|---|
| 5 | H2[18]O alpha=1.006 tnat=2005.2e-6 toce=tnat*(1.0+deltaO18_oce/1000.0) tcorr=1.0+fac_enrichoce18 talph1=talph1_O18 talph2=talph2_O18 talph3=talph3_O18 talps1=talps1_O18 talps2=talps2_O18 tkcin0=tkcin0_O18 tkcin1=tkcin1_O18 tkcin2=tkcin2_O18 tdifrel=tdifrel_O18 alpha_liq_sol=alpha_liq_sol_O18 Rdefault=tnat*(1.0-6.00/1000.) Rmethox=tnat*(1.0+130./1000.) |
|---|
| 6 | H[2]HO alpha=1.010 tnat=155.78e-6 toce=tnat*(1.0+deltaO18_oce/1000.0*pente_MWL) tcorr=1.0+fac_enrichoce18*pente_MWL talph1=24844. talph2=-76.248 talph3=52.612E-3 talps1=16288. talps2=-0.0934 tkcin0=tkcin0_O18*fac_kcin tkcin1=tkcin1_O18*fac_kcin tkcin2=tkcin2_O18*fac_kcin tdifrel=1./0.9755 alpha_liq_sol=1.0212 Rdefault=tnat*(1.0+(10.0-60.0*pente_MWL)/1000.) Rmethox=tnat*(1.0-25.0/1000.) pente_MWL=8.0 fac_kcin=(tdifrel-1.0)/(tdifrel_O18-1.0) |
|---|
| 7 | H[3]HO alpha=1.000 tnat=0.0000 toce=4.0E-19 tcorr=1.0 talph1=46480. talph2=-103.87 talph3=0.0 talps1=46480. talps2=-103.87 tkcin0=0.01056 tkcin1=0.0005016 tkcin2=0.0014432 tdifrel=1./0.968 alpha_liq_sol=1.0 Rdefault=0.0 Rmethox=0.0 |
|---|
| 8 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
