# # # ****************** # ECCO cost function # ****************** &ECCO_COST_NML data_errfile = 'data.err', tbarfile = 'tbar', sbarfile = 'sbar', psbarfile = 'psbar', ubarfile = 'ubar', vbarfile = 'vbar', # # sshv4cost_errfile(1) = 'SO6input/SOSE_DOT_errorguess.err', # sshv4cost_errfile(2) = 'SO6input/SOSE_SSHanom_error.err', # sshv4cost_errfile(3) = 'SO6input/SOSE_SSHanom_error.err', # sshv4cost_errfile(4) = 'SO6input/SOSE_SSHanom_error.err', # sshv4cost_errfile(5) = 'SO6input/SOSE_SSHanom_error.err', # using_topex = .TRUE., topexstartdate1 = 20050101, topexstartdate2 = 000000, topexperiod = 86400.0, topexfile = 'SO6input/RADS_TPJSN_SO_r3', # topexfile = 'SO6input/PODAAC_TPJ_daily_1_6_anomDT', # topexmeanfile = 'SO6input/GGM02_DOT_col6_SO6.bin', # topexmeanfile = 'SO6input/SM004_GRACE3_SO6.bin', topexmeanfile = 'DOT_egm08.bin', # SO6input/DOT_DNSCMSS08-EGM08_SOSE.bin', # using_ers = .TRUE., ersstartdate1 = 20050101, ersstartdate2 = 00000, ersperiod = 86400.0, ersfile = 'SO6input/RADS_ENVERS_SO_r3', # ersfile = 'SO6input/ERS_ENVISAT_daily_SO6_DT', # ersfile = 'SO6input/ERS_ENVISAT_SO_1_6', ssh_errfile = 'SO6input/SOSE_SSHanom_error.err', # ssh_errfile = 'SO6input/ssh_s_SO6.err', ### geoid_errfile = 'SO6input/TPJ_GRACE_r4_SO6.err', ### geoid_covariancefile = 'egm96_1x1.pinv', geoid_errfile = 'DOT_StdDev_capped.err', # SO6input/SOSE_DOT_errorguess.err', # using_gfo = .TRUE., gfostartdate1 = 20050101, gfostartdate2 = 00000, gfoperiod = 86400.0, gfofile = 'SO6input/RADS_GFO_SO_r3', # gfofile = 'SO6input/GFO_SO_1_6_DT', # # scatxdatfile = 'QSCAT_ustr_r1', # scatydatfile = 'QSCAT_vstr_r1', # scatstartdate1 = 19920101, # scatstartdate2 = 60000, # scatperiod = 2629800, # # sststartdate1 = 20050101, # sststartdate2 = 00000, # sstdatfile = 'SO6input/sst4-path-eot.mean_SO6', sstdatfile = 'SO6input/TMI_AMSRE_daily_sst_SO6', # # tmistartdate1 = 19980101, # tmistartdate2 = 60000, # tmidatfile = 'SO6input/tmi_amsre.fusion.mean_SO6', # # sssstartdate1 = 19920101, # sssstartdate2 = 60000, # sssdatfile = 'SSS_tsg_1x1_monthly', # tdatfile = 'SO6input/LEV_wghc_ptmp_r4_pc_SO6.data', sdatfile = 'SO6input/LEV_wghc_salt_r4_pc_SO6.data', temp0errfile = 'SO6input/errorTtot_cap_SO6.err', salt0errfile = 'SO6input/errorStot_cap_SO6.err', # # ctdtfile = 'SO6input/ctdT_r1a_1_6', # ctdsfile = 'SO6input/ctdS_r1a_1_6', # xbtfile = 'SO6input/XBT_r2_1_6', # argotfile = 'SO6input/ARGO_T_r6_1_6', # argosfile = 'SO6input/ARGO_S_r6_1_6', # # udriftfile = 'udrifter_1x1', # vdriftfile = 'vdrifter_1x1', # # udrifterrfile = 'udrifterror.err', # vdrifterrfile = 'vdrifterror.err', # temperrfile = 'SO6input/errorTtot_cap_SO6_4SST.err', salterrfile = 'SO6input/errorStot_cap_SO6.err', # atemp_errfile = 'SO6input/NCEP_R6_tmp2m_degC_SO6_sd.err', aqh_errfile = 'SO6input/NCEP_R6_spfh2m_SO6_sd.err', uwind_errfile = 'SO6input/NCEP_R6_u10m_SO6_sd.err', vwind_errfile = 'SO6input/NCEP_R6_v10m_SO6_sd.err', precip_errfile = 'SO6input/NCEP_R6_rain_SO6_sd.err', swdown_errfile = 'SO6input/NCEP_R6_dsw_SO6_sd.err', # # atemp_errfile = 'SO6input/NCEP_R2day_tmp2m_sd_1992-2003_SO6.data', # aqh_errfile = 'SO6input/NCEP_R2day_spfh2m_sd_1992-2003_SO6.data', # uwind_errfile = 'SO6input/NCEP_R2day_u10m_sd_1992-2003_SO6.data', # vwind_errfile = 'SO6input/NCEP_R2day_v10m_sd_1992-2003_SO6.data', # precip_errfile = 'SO6input/NCEP_R2day_rain_sd_1992-2003_SO6.data', # swdown_errfile = 'SO6input/NCEP_R2day_sw_sd_1992-2003_SO6.data', # swflux_errfile = 'SO6input/NCEP_R2day_sw_sd_1992-2003_SO6.data', # whflux0 = 20., 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, # wmean_hflux = 30., wmean_sflux = 1.6E-8, wmean_tau = 0.1, # wmean_atemp = 4.0, # wmean_aqh = 2.0E-3, # wmean_precip = 6.0E-8, # wmean_swflux = 60., wmean_swdown = 60., # wmean_wind = 4.0, wmean_atemp = 0.20, wmean_aqh = 1.0E-4, wmean_precip = 3.0E-7, wmean_swdown = 3., wmean_wind = 0.20, # scatx_errfile = 'ustr_ncep_r-qscat_yr2005_s.err', # scaty_errfile = 'vstr_ncep_r-qscat_yr2005_s.err', # mult_hflux = 1., mult_sflux = 1., mult_tauu = 1., mult_tauv = 1., # mult_atemp = 1.0, mult_aqh = 1.0, mult_precip = 1.0, mult_swflux = 1.0, mult_swdown = 1.0, mult_uwind = 1.0, mult_vwind = 1.0, # mult_hmean = 1.0, mult_h = 1., mult_tp = 1.0, mult_ers = 1.0, mult_gfo = 1.0, # mult_hmean = 3.0, # mult_h = 1., # mult_tp = 0.4, # mult_ers = 0.4, # mult_gfo = 0.4, mult_temp = 0.002, mult_salt = 0.002, mult_temp0 = 0.1, mult_salt0 = 0.1, mult_sst = 0.1, mult_sss = 1., mult_tmi = 1.0, mult_ctdt = 1., mult_ctds = 1., mult_xbt = 1., mult_argot = 1., mult_argos = 1., mult_drift = 1., mult_tdrift = 1., mult_sdrift = 1., #mult wdrift is now mult on w variance mult_wdrift = 1., mult_scatx = 1., mult_scaty = 1., mult_obcsn = 1., mult_smooth_ic = 1., mult_smooth_bc = 1., mult_obcsvol = 1., # mult ageo is now weight on w^2 mult_ageos = 0.00001, # cost_iprec = 32, cost_yftype = 'RL', /