&H2O

  params deltaO18_oce=0.0        fac_coeff_eq17_liq=0.529                    \
         fac_enrichoce18=0.0005  fac_coeff_eq17_ice=0.529                    \
         talph1_O18=1137.        talps1_O18=11.839       tkcin0_O18=0.006    \
         talph2_O18=-0.4156      talps2_O18=-0.028244    tkcin1_O18=0.000285 \
         talph3_O18=-2.0667E-3   tdifrel_O18=1./0.9723   tkcin2_O18=0.00082  \
         alpha_liq_sol_O18=1.00291

  H216O  alpha=1.000            toce=tnat                            \
         tnat=1.0000            tcorr=1.0                            \
         talph1=0.0             talps1=0.0               tkcin0=0.0  \
         talph2=0.0             talps2=0.0               tkcin1=0.0  \
         talph3=0.0             tdifrel=1.0              tkcin2=0.0  \
         alpha_liq_sol=1.0      Rdefault=tnat*1.0        Rmethox=1.0

  H217O  alpha=1.003            toce=tnat*(1.0+deltaO18_oce/1000.0)^pente_MWL          \
         tnat=40.000e-6         tcorr=1.0+fac_enrichoce18*pente_MWL                    \
         talph1=talph1_O18      talps1=talps1_O18        tkcin0=tkcin0_O18*fac_kcin    \
         talph2=talph2_O18      talps2=talps2_O18        tkcin1=tkcin1_O18*fac_kcin    \
         talph3=talph3_O18      tdifrel=1./0.98555       tkcin2=tkcin2_O18*fac_kcin    \
         alpha_liq_sol=alpha_liq_sol_O18^fac_coeff_eq17_liq                            \
         Rdefault=tnat*(1.0-3.15/1000.)                  Rmethox=tnat*(1.0+230./1000.) \              
         fac_kcin=(tdifrel-1.0)/(tdifrel_O18-1.0)        pente_MWL=0.528

  H218O  alpha=1.006            toce=tnat*(1.0+deltaO18_oce/1000.0)         \
         tnat=2005.2e-6         tcorr=1.0+fac_enrichoce18                   \
         talph1=talph1_O18      talps1=talps1_O18        tkcin0=tkcin0_O18  \
         talph2=talph2_O18      talps2=talps2_O18        tkcin1=tkcin1_O18  \
         talph3=talph3_O18      tdifrel=tdifrel_O18      tkcin2=tkcin2_O18  \
         alpha_liq_sol=alpha_liq_sol_O18                                    \
         Rdefault=tnat*(1.0-6.00/1000.)                  Rmethox=tnat*(1.0+130./1000.)  

  HDO    alpha=1.010            toce=tnat*(1.0+deltaO18_oce/1000.0*pente_MWL)          \
         tnat=155.76e-6         tcorr=1.0+fac_enrichoce18*pente_MWL                    \
         talph1=24844.          talps1=16288.            tkcin0=tkcin0_O18*fac_kcin    \
         talph2=-76.248         talps2=-0.0934           tkcin1=tkcin1_O18*fac_kcin    \
         talph3=52.612E-3       tdifrel=1./0.9755        tkcin2=tkcin2_O18*fac_kcin    \
         alpha_liq_sol=1.0212                                                          \
         Rdefault=tnat*(1.0+(10.0-6.0*pente_MWL)/1000.)  Rmethox=tnat*(1.0-25.0/1000.) \
         fac_kcin=(tdifrel-1.0)/(tdifrel_O18-1.0)        pente_MWL=8.0

  HTO    alpha=1.000            toce=4.0E-19                               \
         tnat=0.0000            tcorr=1.0                                  \
         talph1=46480.          talps1=46480.            tkcin0=0.01056    \
         talph2=-103.87         talps2=-103.87           tkcin1=0.0005016  \
         talph3=0.0             tdifrel=1./0.968         tkcin2=0.0014432  \
         alpha_liq_sol=1.0      Rdefault=0.0             Rmethox=0.0