1 |
ou.wang |
1.1 |
# |
2 |
|
|
# |
3 |
|
|
# ****************** |
4 |
|
|
# ECCO cost function |
5 |
|
|
# ****************** |
6 |
|
|
# |
7 |
|
|
&ECCO_COST_NML |
8 |
|
|
# |
9 |
|
|
data_errfile = 'data.err', |
10 |
|
|
temp0errfile = 'Theta_sigma_smoothed_method_02_masked_merged_capped_areascaled.bin', |
11 |
|
|
salt0errfile = 'Salt_sigma_smoothed_method_02_masked_merged_capped_areascaled.bin', |
12 |
|
|
cost_iprec = 32, |
13 |
|
|
cost_yftype = 'RL', |
14 |
|
|
# |
15 |
|
|
/ |
16 |
|
|
# |
17 |
|
|
&ECCO_GENCOST_NML |
18 |
|
|
# |
19 |
|
|
gencost_avgperiod(1) = 'day', |
20 |
|
|
gencost_barfile(1) = 'gbar_area', |
21 |
|
|
gencost_datafile(1) = 'NOAA_NSIDC_DAILY_v3b_MAPPED_TO_LLC90', |
22 |
|
|
gencost_errfile(1) = 'sigma_iceconc_eccollc_coastalmask_areascaled.bin', |
23 |
|
|
gencost_name(1) = 'siv4-conc', |
24 |
|
|
gencost_spmin(1) = -999., |
25 |
|
|
gencost_spmax(1) = 999., |
26 |
|
|
gencost_spzero(1) = -9999., |
27 |
|
|
gencost_outputlevel(1)=1, |
28 |
|
|
mult_gencost(1) = 0., |
29 |
|
|
# |
30 |
|
|
gencost_avgperiod(2) = 'day', |
31 |
|
|
gencost_barfile(2) = 'gbar_deconc', |
32 |
|
|
gencost_name(2) = 'siv4-deconc', |
33 |
|
|
gencost_datafile(2) = 'NOAA_NSIDC_DAILY_v3b_MAPPED_TO_LLC90', |
34 |
|
|
gencost_errfile(2) = 'seaiceproxy_sigma_coastalmask_areascaled.bin', |
35 |
|
|
gencost_spmin(2) = -999., |
36 |
|
|
gencost_spmax(2) = 999., |
37 |
|
|
gencost_spzero(2) = -9999., |
38 |
|
|
gencost_outputlevel(2)=1, |
39 |
|
|
mult_gencost(2) = 0.071, |
40 |
|
|
# |
41 |
|
|
gencost_avgperiod(3) = 'day', |
42 |
|
|
gencost_barfile(3) = 'gbar_exconc', |
43 |
|
|
gencost_name(3) = 'siv4-exconc', |
44 |
|
|
gencost_datafile(3) = 'NOAA_NSIDC_DAILY_v3b_MAPPED_TO_LLC90', |
45 |
|
|
gencost_errfile(3) = 'seaiceproxy_sigma_coastalmask_areascaled.bin', |
46 |
|
|
gencost_spmin(3) = -999., |
47 |
|
|
gencost_spmax(3) = 999., |
48 |
|
|
gencost_spzero(3) = -9999., |
49 |
|
|
gencost_outputlevel(3)=1, |
50 |
|
|
mult_gencost(3) = 0., |
51 |
|
|
# |
52 |
|
|
gencost_avgperiod(4) = 'month', |
53 |
|
|
gencost_barfile(4) = 'm_sstmon', |
54 |
|
|
gencost_datafile(4) = 'reynolds_oiv2_r1_mar2018_simasked', |
55 |
|
|
gencost_errfile(4) = 'sigma_half_areascaled.bin', |
56 |
|
|
gencost_name(4) = 'sst-reynoldsmean', |
57 |
|
|
gencost_spmin(4) = -1.8, |
58 |
|
|
gencost_spmax(4) = 40., |
59 |
|
|
gencost_spzero(4) = 0., |
60 |
|
|
gencost_preproc(1,4)='mean', |
61 |
|
|
gencost_outputlevel(4)=1, |
62 |
|
|
mult_gencost(4) = 1., |
63 |
|
|
# |
64 |
|
|
gencost_avgperiod(5) = 'month', |
65 |
|
|
gencost_barfile(5) = 'm_sstmon', |
66 |
|
|
gencost_datafile(5) = 'reynolds_oiv2_r1_mar2018_simasked', |
67 |
|
|
gencost_errfile(5) = 'sigma_half_areascaled.bin', |
68 |
|
|
gencost_name(5) = 'sst-reynoldsanom', |
69 |
|
|
gencost_spmin(5) = -1.8, |
70 |
|
|
gencost_spmax(5) = 40., |
71 |
|
|
gencost_spzero(5) = 0., |
72 |
|
|
gencost_preproc(1,5)='anom', |
73 |
|
|
gencost_outputlevel(5)=1, |
74 |
|
|
mult_gencost(5) = 1., |
75 |
|
|
|
76 |
|
|
# |
77 |
|
|
# gencost_avgperiod(5) = 'month', |
78 |
|
|
# gencost_barfile(5) = 'sstmon', |
79 |
|
|
# gencost_datafile(5) = 'tmi_amsre_oisst_r1', |
80 |
|
|
# gencost_errfile(5) = 'sigma_half_areascaled.bin', |
81 |
|
|
# gencost_name(5) = 'sst-tmi-amsre', |
82 |
|
|
# gencost_startdate1(5) = 19980101, |
83 |
|
|
# gencost_startdate2(5)= 00000, |
84 |
|
|
# gencost_spmin(5) = -1.8, |
85 |
|
|
# gencost_spmax(5) = 40., |
86 |
|
|
# gencost_spzero(5) = 0., |
87 |
|
|
# gencost_outputlevel(5)=1, |
88 |
|
|
# mult_gencost(5) = 0., |
89 |
|
|
# |
90 |
|
|
gencost_avgperiod(6) = 'month', |
91 |
|
|
gencost_barfile(6) = 'm_bp_nopabarmon', |
92 |
|
|
gencost_datafile(6) = 'GRACE_jpl_rl05m_20180122', |
93 |
|
|
gencost_errfile(6) = 'GRACE_jpl_rl05m_20180122_err_areascaled', |
94 |
|
|
gencost_name(6) = 'bpv4-grace', |
95 |
|
|
gencost_outputlevel(6)=1, |
96 |
|
|
mult_gencost(6) = 1., |
97 |
|
|
# |
98 |
|
|
gencost_avgperiod(7) = 'month', |
99 |
|
|
gencost_barfile(7) = 'm_thetamon', |
100 |
|
|
gencost_datafile(7) = 'T_monthly_woa09_highlatnearsurface_masked', |
101 |
|
|
gencost_errfile(7) = 'Theta_sigma_smoothed_method_02_masked_merged_capped_areascaled.bin', |
102 |
|
|
# gencost_name(7) = 'thetaclim', |
103 |
|
|
gencost_name(7) = 'thetamean', |
104 |
|
|
gencost_spmin(7) = -1.8, |
105 |
|
|
gencost_spmax(7) = 40., |
106 |
|
|
gencost_spzero(7) = 0., |
107 |
|
|
gencost_is3d(7)=.TRUE., |
108 |
|
|
gencost_preproc(1,7)='mean', |
109 |
|
|
gencost_outputlevel(7)=1, |
110 |
|
|
mult_gencost(7) = 1., |
111 |
|
|
# |
112 |
|
|
gencost_avgperiod(8) = 'month', |
113 |
|
|
gencost_barfile(8) = 'm_saltmon', |
114 |
|
|
gencost_datafile(8) = 'S_monthly_woa09_highlatnearsurface_masked', |
115 |
|
|
gencost_errfile(8) = 'Salt_sigma_smoothed_method_02_masked_merged_capped_areascaled.bin', |
116 |
|
|
# gencost_name(8) = 'saltclim', |
117 |
|
|
gencost_name(8) = 'saltmean', |
118 |
|
|
gencost_spmin(8) = 25., |
119 |
|
|
gencost_spmax(8) = 40., |
120 |
|
|
gencost_spzero(8) = 0., |
121 |
|
|
gencost_is3d(8)=.TRUE., |
122 |
|
|
gencost_preproc(1,8)='mean', |
123 |
|
|
gencost_outputlevel(8)=1, |
124 |
|
|
mult_gencost(8) = 1., |
125 |
|
|
# |
126 |
|
|
gencost_datafile(11) = 'mdt2_dtu17.bin', |
127 |
|
|
gencost_startdate1(11) = 19930101, |
128 |
|
|
gencost_startdate2(11) = 00000, |
129 |
|
|
gencost_enddate1(11) = 20171231, |
130 |
|
|
gencost_enddate2(11) = 00000, |
131 |
|
|
gencost_errfile(11) = 'mdt_dtu13_error_m_areascaled.bin', |
132 |
|
|
gencost_name(11) = 'sshv4-mdt', |
133 |
|
|
gencost_posproc(1,11)='smooth', |
134 |
|
|
gencost_posproc_c(1,11)='sshv4_scale_1p5points.bin', |
135 |
|
|
gencost_posproc_i(1,11)=300, |
136 |
|
|
# gencost_scalefile(11) = 'sshv4_scale_1p5points.bin', |
137 |
|
|
gencost_outputlevel(11)=1, |
138 |
|
|
mult_gencost(11) = 1., |
139 |
|
|
# |
140 |
|
|
gencost_barfile(12) = 'm_eta_dynday', |
141 |
|
|
gencost_datafile(12) = 'RADS_TJ_mar2018', |
142 |
|
|
gencost_startdate1(12) = 19920101, |
143 |
|
|
gencost_startdate2(12) = 00000, |
144 |
|
|
gencost_avgperiod(12) = 'day', |
145 |
|
|
gencost_errfile(12) = 'slaerr_gridscale_r5_areascaled.err', |
146 |
|
|
gencost_name(12) = 'sshv4-tp', |
147 |
|
|
gencost_outputlevel(12)=1, |
148 |
|
|
mult_gencost(12) = 0., |
149 |
|
|
# |
150 |
|
|
gencost_datafile(13) = 'RADS_ERS_ENV_SA_mar2018', |
151 |
|
|
gencost_startdate1(13) = 19920101, |
152 |
|
|
gencost_startdate2(13) = 00000, |
153 |
|
|
gencost_avgperiod(13) = 'day', |
154 |
|
|
gencost_errfile(13) = 'slaerr_gridscale_r5_areascaled.err', |
155 |
|
|
gencost_name(13) = 'sshv4-ers', |
156 |
|
|
gencost_outputlevel(13)=1, |
157 |
|
|
mult_gencost(13) = 0., |
158 |
|
|
# |
159 |
|
|
gencost_datafile(14) = 'RADS_GFO_C2_mar2018', |
160 |
|
|
gencost_startdate1(14) = 19920101, |
161 |
|
|
gencost_startdate2(14) = 00000, |
162 |
|
|
gencost_avgperiod(14) = 'day', |
163 |
|
|
gencost_errfile(14) = 'slaerr_gridscale_r5_areascaled.err', |
164 |
|
|
gencost_name(14) = 'sshv4-gfo', |
165 |
|
|
gencost_outputlevel(14)=1, |
166 |
|
|
mult_gencost(14) = 0., |
167 |
|
|
# |
168 |
|
|
gencost_errfile(15) = 'slaerr_largescale_r5_areascaled.err', |
169 |
|
|
gencost_name(15) = 'sshv4-lsc', |
170 |
|
|
#gencost_scalefile(15) = 'sshv4_scale_3points.bin', |
171 |
|
|
gencost_posproc(1,15)='smooth', |
172 |
|
|
gencost_posproc_c(1,15)='sshv4_scale_3points.bin', |
173 |
|
|
gencost_posproc_i(1,15)=300, |
174 |
|
|
gencost_outputlevel(15)=1, |
175 |
|
|
mult_gencost(15) = 0.001, |
176 |
|
|
# |
177 |
|
|
#gencost_name(16) = 'sshv4-gmsl', |
178 |
|
|
#gencost_outputlevel(16)=1, |
179 |
|
|
#mult_gencost(16) = 0., |
180 |
|
|
# |
181 |
|
|
gencost_avgperiod(17) = 'month', |
182 |
|
|
gencost_barfile(17) = 'm_thetamon', |
183 |
|
|
gencost_datafile(17) = 'T_monthly_woa09_highlatnearsurface_masked', |
184 |
|
|
gencost_errfile(17) = 'Theta_sigma_smoothed_method_02_masked_merged_capped_areascaled.bin', |
185 |
|
|
# gencost_name(17) = 'thetaclim', |
186 |
|
|
gencost_name(17) = 'thetaanom', |
187 |
|
|
gencost_spmin(17) = -1.8, |
188 |
|
|
gencost_spmax(17) = 40., |
189 |
|
|
gencost_spzero(17) = 0., |
190 |
|
|
gencost_is3d(17)=.TRUE., |
191 |
|
|
gencost_preproc(1,17)='clim', |
192 |
|
|
gencost_preproc_i(1,17)=12, |
193 |
|
|
gencost_preproc(2,17)='anom', |
194 |
|
|
gencost_outputlevel(17)=1, |
195 |
|
|
mult_gencost(17) = 1., |
196 |
|
|
# |
197 |
|
|
gencost_avgperiod(18) = 'month', |
198 |
|
|
gencost_barfile(18) = 'm_saltmon', |
199 |
|
|
gencost_datafile(18) = 'S_monthly_woa09_highlatnearsurface_masked', |
200 |
|
|
gencost_errfile(18) = 'Salt_sigma_smoothed_method_02_masked_merged_capped_areascaled.bin', |
201 |
|
|
# gencost_name(18) = 'saltclim', |
202 |
|
|
gencost_name(18) = 'saltanom', |
203 |
|
|
gencost_spmin(18) = 25., |
204 |
|
|
gencost_spmax(18) = 40., |
205 |
|
|
gencost_spzero(18) = 0., |
206 |
|
|
gencost_is3d(18)=.TRUE., |
207 |
|
|
gencost_preproc(1,18)='clim', |
208 |
|
|
gencost_preproc_i(1,18)=12, |
209 |
|
|
gencost_preproc(2,18)='anom', |
210 |
|
|
gencost_outputlevel(18)=1, |
211 |
|
|
mult_gencost(18) = 1., |
212 |
|
|
# |
213 |
|
|
gencost_avgperiod(19) = 'month', |
214 |
|
|
gencost_barfile(19) = 'm_sssmon', |
215 |
|
|
gencost_datafile(19) = 'monthlyL3v5_20180129', |
216 |
|
|
gencost_errfile(19) = 'errStdAQ_v3_llc_090_floored_areascaled', |
217 |
|
|
gencost_name(19) = 'sss_AQanom', |
218 |
|
|
gencost_spmin(19) = 25., |
219 |
|
|
gencost_spmax(19) = 40., |
220 |
|
|
gencost_spzero(19) = 0., |
221 |
|
|
gencost_preproc(1,19) = 'anom', |
222 |
|
|
# gencost_is3d(19)=.TRUE., |
223 |
|
|
gencost_outputlevel(19)=1, |
224 |
|
|
mult_gencost(19) = 1., |
225 |
|
|
# |
226 |
|
|
gencost_avgperiod(20) = 'month', |
227 |
|
|
gencost_barfile(20) = 'm_sssmon', |
228 |
|
|
gencost_datafile(20) = 'monthlyL3v5_20180129', |
229 |
|
|
gencost_errfile(20) = 'biasStdAQ_v3_llc_090_floored_areascaled', |
230 |
|
|
gencost_name(20) = 'sss_AQmean', |
231 |
|
|
gencost_spmin(20) = 25., |
232 |
|
|
gencost_spmax(20) = 40., |
233 |
|
|
gencost_spzero(20) = 0., |
234 |
|
|
gencost_preproc(1,20) = 'mean', |
235 |
|
|
# gencost_is3d(20)=.TRUE., |
236 |
|
|
gencost_outputlevel(20)=1, |
237 |
|
|
mult_gencost(20) = 1., |
238 |
|
|
# |
239 |
|
|
gencost_avgperiod(21) = 'month', |
240 |
|
|
gencost_barfile(21) = 'm_bp_nopabarmon', |
241 |
|
|
gencost_datafile(21) = 'GRACE_jpl_rl05m_20180122_SpatialMean.asc', |
242 |
|
|
gencost_name(21) = 'gmbp', |
243 |
|
|
gencost_wei1d(21) = 0.1583, |
244 |
|
|
gencost_spmin(21) = -990., |
245 |
|
|
gencost_spmax(21) = 990., |
246 |
|
|
gencost_spzero(21) = 0., |
247 |
|
|
gencost_is1d(21)=.TRUE., |
248 |
|
|
gencost_outputlevel(21)=1, |
249 |
|
|
mult_gencost(21) = 1., |
250 |
|
|
# |
251 |
|
|
gencost_avgperiod(22) = 'month', |
252 |
|
|
gencost_barfile(22) = 'm_eta_dynmon', |
253 |
|
|
gencost_datafile(22) = '20180228_ensemble_average_gmsl_ref.asc', |
254 |
|
|
gencost_name(22) = 'gmsl', |
255 |
|
|
gencost_wei1d(22) = 0.00382, |
256 |
|
|
gencost_spmin(22) = -990., |
257 |
|
|
gencost_spmax(22) = 990., |
258 |
|
|
gencost_spzero(22) = 0., |
259 |
|
|
gencost_is1d(22)=.TRUE., |
260 |
|
|
gencost_outputlevel(22)=1, |
261 |
|
|
mult_gencost(22) = 1., |
262 |
|
|
# |
263 |
|
|
gencost_barfile(27) = 'm_etaday', |
264 |
|
|
gencost_datafile(27) = 'RADS_TJ_mar2018', |
265 |
|
|
gencost_startdate1(27) = 19920101, |
266 |
|
|
gencost_startdate2(27) = 00000, |
267 |
|
|
gencost_avgperiod(27) = 'day', |
268 |
|
|
gencost_errfile(27) = 'slaerr_gridscale_r5_areascaled.err', |
269 |
|
|
gencost_spmin(27) = -990., |
270 |
|
|
gencost_spmax(27) = 990., |
271 |
|
|
gencost_spzero(27) = 0., |
272 |
|
|
gencost_outputlevel(27)=1, |
273 |
|
|
mult_gencost(27) = 0., |
274 |
|
|
# |
275 |
|
|
gencost_avgperiod(28) = 'month', |
276 |
|
|
gencost_barfile(28) = 'm_etamon', |
277 |
|
|
gencost_datafile(28) = '20180228_ensemble_average_gmsl_ref.asc', |
278 |
|
|
gencost_wei1d(28) = 0.00382, |
279 |
|
|
gencost_spmin(28) = -990., |
280 |
|
|
gencost_spmax(28) = 990., |
281 |
|
|
gencost_spzero(28) = 0., |
282 |
|
|
gencost_is1d(28)=.TRUE., |
283 |
|
|
gencost_outputlevel(28)=1, |
284 |
|
|
mult_gencost(28) = 0., |
285 |
|
|
# |
286 |
|
|
gencost_avgperiod(29) = 'month', |
287 |
|
|
gencost_barfile(29) = 'm_bpmon', |
288 |
|
|
gencost_datafile(29) = 'GRACE_jpl_rl05m_20180122_SpatialMean.asc', |
289 |
|
|
gencost_wei1d(29) = 0.1583, |
290 |
|
|
gencost_spmin(29) = -990., |
291 |
|
|
gencost_spmax(29) = 990., |
292 |
|
|
gencost_spzero(29) = 0., |
293 |
|
|
gencost_is1d(29)=.TRUE., |
294 |
|
|
gencost_outputlevel(29)=1, |
295 |
|
|
mult_gencost(29) = 0., |
296 |
|
|
# |
297 |
|
|
gencost_avgperiod(30) = 'day', |
298 |
|
|
gencost_barfile(30) = 'm_bpday', |
299 |
|
|
gencost_datafile(30) = 'GRACE_fake_daily', |
300 |
|
|
gencost_errfile(30) = 'GRACE_jpl_rl05m_20180122_err_areascaled', |
301 |
|
|
gencost_spmin(30) = -990., |
302 |
|
|
gencost_spmax(30) = 990., |
303 |
|
|
gencost_spzero(30) = 0., |
304 |
|
|
gencost_outputlevel(30)=1, |
305 |
|
|
mult_gencost(30) = 0., |
306 |
|
|
# |
307 |
|
|
gencost_avgperiod(31) = 'day', |
308 |
|
|
gencost_barfile(31) = 'm_bp_nopabarday', |
309 |
|
|
gencost_datafile(31) = 'GRACE_fake_daily' |
310 |
|
|
gencost_errfile(31) = 'GRACE_jpl_rl05m_20180122_err_areascaled', |
311 |
|
|
gencost_spmin(31) = -990., |
312 |
|
|
gencost_spmax(31) = 990., |
313 |
|
|
gencost_spzero(31) = 0., |
314 |
|
|
gencost_outputlevel(31)=1, |
315 |
|
|
mult_gencost(31) = 0., |
316 |
|
|
/ |
317 |
|
|
# |