#
#
# ******************
# ECCO cost function
# ******************
#
 &ECCO_COST_NML
#
 data_errfile    = 'data.err',
 tbarfile        = 'tbar',
 sbarfile        = 'sbar',
 psbarfile       = 'psbar',
#
 temp0errfile     = 'input_sigma/sigma_T_atlas2_eccollc.bin',
 salt0errfile     = 'input_sigma/sigma_S_atlas2_eccollc.bin',
#
 whflux0              = 20.0,
 wsflux0              = 3.0E-8,
 wtau0                = 2.0E-2,
#
 watemp0              = 1.0,
 waqh0                = 0.5E-3,
 wprecip0             = 1.5E-8,
 wswflux0             = 15.0,
 wswdown0             = 15.0,
 wwind0               = 1.0,
#
 wkapgm0              = 5.0E3,
 wkapredi0            = 5.0E3,
 wdiffkr0             = 1.0E-4,
 wedtau0              = 1.0E0,
 wbottomdrag0         = 1.0E-5,
#
 wmean_hflux          = 60.0,
 wmean_sflux          = 3.2E-8,
 wmean_tau            = 0.2,
#
 wmean_atemp          = 8.0,
 wmean_aqh            = 4.0E-3,
 wmean_precip         = 9.0E-8,
 wmean_swflux         = 90.0,
 wmean_swdown         = 90.0,
 wmean_wind           = 8.0,
#
 mult_temp0  = 1.,
 mult_salt0  = 1.,
#
 tdatfile        = 'T_OWPv1_M_eccollc_90x50.bin',
 sdatfile        = 'S_OWPv1_M_eccollc_90x50.bin',
 temperrfile     = 'input_sigma/sigma_T_atlas2_eccollc.bin',
 salterrfile     = 'input_sigma/sigma_S_atlas2_eccollc.bin',
 mult_temp   = 1.,
 mult_salt   = 1.,
#
 sststartdate1   = 19920101,
 sststartdate2   = 00000,
 sstdatfile      = 'input_sst_reynolds_v4/reynolds_oiv2_r1',
# tmistartdate1   = 19980101,
 tmistartdate1   = 19920101,
 tmistartdate2   = 00000,
 tmidatfile      = 'input_sst_remss_v4/tmi_amsre_oisst_r1',
 ssterrfile      = 'input_sigma/sigma_half.bin',
 mult_sst    = 1.,
 mult_tmi    = 1.,
#
# mdtdatfile   = 'input_mdt/mdt_pak09.bin',
 mdtdatfile   = 'input_mdt/mdt_max09.bin',
 mdtstartdate1 = 19920101,
 mdtstartdate2 = 00000,
 mdtenddate1 = 20021231,
# mdtenddate1 = 19920129,
 mdtenddate2 = 00000,
 geoid_errfile   = 'input_err/g_TPJ_GRACE_r4.err',
 using_topex     = .TRUE.,
 topexstartdate1 = 19920101,
 topexstartdate2 = 00000,
 topexperiod     = 86400.0,
 topexfile       = 'input_rads_noice_mad/RADS_TJ_v4_noice',
 tp_errfile     = 'input_err/g_PODAAC_TPJ_sd_s_filled_11yr',
 using_ers       = .TRUE.,
 ersstartdate1   = 19920101,
 ersstartdate2   = 00000,
 ersperiod       = 86400.0,
 ersfile         = 'input_rads_noice_mad/RADS_ERS_ENV_v4_noice',
 ers_errfile     = 'input_err/g_PODAAC_TPJ_sd_s_filled_11yr',
 using_gfo       = .TRUE.,
 gfostartdate1   = 19920101,
 gfostartdate2   = 00000,
 gfoperiod       = 86400.0,
 gfofile         = 'input_rads_noice_mad/RADS_GFO_v4_noice',
 gfo_errfile     = 'input_err/g_PODAAC_TPJ_sd_s_filled_11yr',
#
 mult_tp = 1.,
 mult_ers = 1.,
 mult_gfo = 1.,
 mult_sshv4cost = 1.,
#
 cost_iprec  = 32,
 cost_yftype = 'RL',
#
 /
#
 &ECCO_GENCOST_NML
#
 gencost_errfile(1) = 'input_sigma/sigma_half.bin',
 gencost_name(1) = 'sstv4-amsre-lsc',
 gencost_scalefile(1) = 'sshv4_scale_3points.bin',
 mult_gencost(1) = 1.,
 gencost_errfile(2) = 'input_sigma/sigma_half.bin',
 gencost_name(2) = 'sstv4-amsre',
 gencost_datafile(2) = 'input_sst_amsre_daily/tmi_amsre_l2p_r1',
 mult_gencost(2) = 1.,
#
 gencost_errfile(11) = 'input_sigma/sigma_MDT_glob_eccollc.bin',
 gencost_name(11) = 'sshv4-mdt',
 gencost_scalefile(11) = 'sshv4_scale_1p5points.bin',
 mult_gencost(11) = 1.,
 gencost_errfile(12) = 'input_err/slaerr_gridscale_r1.err',
 gencost_name(12) = 'sshv4-tp',
 mult_gencost(12) = 1.,
 gencost_errfile(13) = 'input_err/slaerr_gridscale_r1.err',
 gencost_name(13) = 'sshv4-ers',
 mult_gencost(13) = 1.,
 gencost_errfile(14) = 'input_err/slaerr_gridscale_r1.err',
 gencost_name(14) = 'sshv4-gfo',
 mult_gencost(14) = 1.,
 gencost_errfile(15) = 'input_err/slaerr_largescale_r1.err',
 gencost_name(15) = 'sshv4-lsc',
 gencost_scalefile(15) = 'sshv4_scale_3points.bin',
 mult_gencost(15) = 1.,
#
 /
#