aste_90x150x60/input_ad/data.ecco

#
#
# ******************
# 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.,
#
 /
#
1	#
2	#
3	# ******************
4	# ECCO cost function
5	# ******************
6	#
7	&ECCO_COST_NML
8	#
9	data_errfile = 'data.err',
10	tbarfile = 'tbar',
11	sbarfile = 'sbar',
12	psbarfile = 'psbar',
13	#
14	temp0errfile = 'input_sigma/sigma_T_atlas2_eccollc.bin',
15	salt0errfile = 'input_sigma/sigma_S_atlas2_eccollc.bin',
16	#
17	whflux0 = 20.0,
18	wsflux0 = 3.0E-8,
19	wtau0 = 2.0E-2,
20	#
21	watemp0 = 1.0,
22	waqh0 = 0.5E-3,
23	wprecip0 = 1.5E-8,
24	wswflux0 = 15.0,
25	wswdown0 = 15.0,
26	wwind0 = 1.0,
27	#
28	wkapgm0 = 5.0E3,
29	wkapredi0 = 5.0E3,
30	wdiffkr0 = 1.0E-4,
31	wedtau0 = 1.0E0,
32	wbottomdrag0 = 1.0E-5,
33	#
34	wmean_hflux = 60.0,
35	wmean_sflux = 3.2E-8,
36	wmean_tau = 0.2,
37	#
38	wmean_atemp = 8.0,
39	wmean_aqh = 4.0E-3,
40	wmean_precip = 9.0E-8,
41	wmean_swflux = 90.0,
42	wmean_swdown = 90.0,
43	wmean_wind = 8.0,
44	#
45	mult_temp0 = 1.,
46	mult_salt0 = 1.,
47	#
48	tdatfile = 'T_OWPv1_M_eccollc_90x50.bin',
49	sdatfile = 'S_OWPv1_M_eccollc_90x50.bin',
50	temperrfile = 'input_sigma/sigma_T_atlas2_eccollc.bin',
51	salterrfile = 'input_sigma/sigma_S_atlas2_eccollc.bin',
52	mult_temp = 1.,
53	mult_salt = 1.,
54	#
55	sststartdate1 = 19920101,
56	sststartdate2 = 00000,
57	sstdatfile = 'input_sst_reynolds_v4/reynolds_oiv2_r1',
58	# tmistartdate1 = 19980101,
59	tmistartdate1 = 19920101,
60	tmistartdate2 = 00000,
61	tmidatfile = 'input_sst_remss_v4/tmi_amsre_oisst_r1',
62	ssterrfile = 'input_sigma/sigma_half.bin',
63	mult_sst = 1.,
64	mult_tmi = 1.,
65	#
66	# mdtdatfile = 'input_mdt/mdt_pak09.bin',
67	mdtdatfile = 'input_mdt/mdt_max09.bin',
68	mdtstartdate1 = 19920101,
69	mdtstartdate2 = 00000,
70	mdtenddate1 = 20021231,
71	# mdtenddate1 = 19920129,
72	mdtenddate2 = 00000,
73	geoid_errfile = 'input_err/g_TPJ_GRACE_r4.err',
74	using_topex = .TRUE.,
75	topexstartdate1 = 19920101,
76	topexstartdate2 = 00000,
77	topexperiod = 86400.0,
78	topexfile = 'input_rads_noice_mad/RADS_TJ_v4_noice',
79	tp_errfile = 'input_err/g_PODAAC_TPJ_sd_s_filled_11yr',
80	using_ers = .TRUE.,
81	ersstartdate1 = 19920101,
82	ersstartdate2 = 00000,
83	ersperiod = 86400.0,
84	ersfile = 'input_rads_noice_mad/RADS_ERS_ENV_v4_noice',
85	ers_errfile = 'input_err/g_PODAAC_TPJ_sd_s_filled_11yr',
86	using_gfo = .TRUE.,
87	gfostartdate1 = 19920101,
88	gfostartdate2 = 00000,
89	gfoperiod = 86400.0,
90	gfofile = 'input_rads_noice_mad/RADS_GFO_v4_noice',
91	gfo_errfile = 'input_err/g_PODAAC_TPJ_sd_s_filled_11yr',
92	#
93	mult_tp = 1.,
94	mult_ers = 1.,
95	mult_gfo = 1.,
96	mult_sshv4cost = 1.,
97	#
98	cost_iprec = 32,
99	cost_yftype = 'RL',
100	#
101	/
102	#
103	&ECCO_GENCOST_NML
104	#
105	gencost_errfile(1) = 'input_sigma/sigma_half.bin',
106	gencost_name(1) = 'sstv4-amsre-lsc',
107	gencost_scalefile(1) = 'sshv4_scale_3points.bin',
108	mult_gencost(1) = 1.,
109	gencost_errfile(2) = 'input_sigma/sigma_half.bin',
110	gencost_name(2) = 'sstv4-amsre',
111	gencost_datafile(2) = 'input_sst_amsre_daily/tmi_amsre_l2p_r1',
112	mult_gencost(2) = 1.,
113	#
114	gencost_errfile(11) = 'input_sigma/sigma_MDT_glob_eccollc.bin',
115	gencost_name(11) = 'sshv4-mdt',
116	gencost_scalefile(11) = 'sshv4_scale_1p5points.bin',
117	mult_gencost(11) = 1.,
118	gencost_errfile(12) = 'input_err/slaerr_gridscale_r1.err',
119	gencost_name(12) = 'sshv4-tp',
120	mult_gencost(12) = 1.,
121	gencost_errfile(13) = 'input_err/slaerr_gridscale_r1.err',
122	gencost_name(13) = 'sshv4-ers',
123	mult_gencost(13) = 1.,
124	gencost_errfile(14) = 'input_err/slaerr_gridscale_r1.err',
125	gencost_name(14) = 'sshv4-gfo',
126	mult_gencost(14) = 1.,
127	gencost_errfile(15) = 'input_err/slaerr_largescale_r1.err',
128	gencost_name(15) = 'sshv4-lsc',
129	gencost_scalefile(15) = 'sshv4_scale_3points.bin',
130	mult_gencost(15) = 1.,
131	#
132	/
133	#