1 |
# |
2 |
# |
3 |
# ****************** |
4 |
# ECCO cost function |
5 |
# ****************** |
6 |
&ECCO_COST_NML |
7 |
data_errfile = 'data.err', |
8 |
tbarfile = 'tbar', |
9 |
sbarfile = 'sbar', |
10 |
psbarfile = 'psbar', |
11 |
ubarfile = 'ubar', |
12 |
vbarfile = 'vbar', |
13 |
hfluxbarfile = 'hfluxbar', |
14 |
sfluxbarfile = 'sfluxbar', |
15 |
# |
16 |
using_topex = .TRUE., |
17 |
topexstartdate1 = 19920101, |
18 |
topexstartdate2 = 30000, |
19 |
topexperiod = 86400.0, |
20 |
# topexfile = 'mit9yrmn_TPJ_daily_anom', |
21 |
# topexmeanfile = 'mit9yrmn_TP_mean_s', |
22 |
topexfile = 'PODAAC_TPJ_daily_anom_r2', |
23 |
# topexmeanfile = 'PODAAC_TPJ_mean_s_r1_11yr', |
24 |
# topexmeanfile = 'GGM02_DOT_col6', |
25 |
topexmeanfile = 'SM004_GRACE3.bin', |
26 |
#HU! check geoid_errfile!!! |
27 |
# |
28 |
using_ers = .TRUE., |
29 |
ersstartdate1 = 19920101, |
30 |
ersstartdate2 = 30000, |
31 |
ersperiod = 86400.0, |
32 |
ersfile = 'ERS_ENVISAT_daily_r4', |
33 |
ssh_errfile = 'ssh_s.err', |
34 |
# |
35 |
using_gfo = .TRUE., |
36 |
gfostartdate1 = 19920101, |
37 |
gfostartdate2 = 30000, |
38 |
gfoperiod = 86400.0, |
39 |
gfofile = 'GFO_r1', |
40 |
# |
41 |
scatxdatfile = 'QSCAT_ustr_r1', |
42 |
scatydatfile = 'QSCAT_vstr_r1', |
43 |
scatstartdate1 = 19920101, |
44 |
scatstartdate2 = 30000, |
45 |
scatperiod = 2629800, |
46 |
# |
47 |
sststartdate1 = 19920101, |
48 |
sststartdate2 = 30000, |
49 |
sstdatfile = 'SST_monthly', |
50 |
# |
51 |
tmistartdate1 = 19980101, |
52 |
tmistartdate2 = 30000, |
53 |
tmidatfile = 'TMI_monthly', |
54 |
# |
55 |
sssstartdate1 = 19920101, |
56 |
sssstartdate2 = 30000, |
57 |
sssdatfile = 'SSS_tsg_1x1_monthly', |
58 |
# |
59 |
## tdatfile = 'lev01_monthly_ptmp.bin_corK', |
60 |
## sdatfile = 'lev01_monthly_salt.bin_corK', |
61 |
## tdatfile = 'lev_monthly_temp_fl.bin', |
62 |
## sdatfile = 'lev_monthly_salt_fl.bin', |
63 |
## tdatfile = 'LEVITUS_1x1_ptmp_NEW_corK', |
64 |
## sdatfile = 'LEVITUS_1x1_salt_NEW_corK', |
65 |
## tdatfile = 'LEV_wghc_ptmp_corK', |
66 |
## sdatfile = 'LEV_wghc_salt_corK', |
67 |
## tdatfile = 'LEV_wghc_ptmp_r2', |
68 |
## sdatfile = 'LEV_wghc_salt_r2', |
69 |
## tdatfile = 'LEV_wghc_ptmp_r3', |
70 |
## sdatfile = 'LEV_wghc_salt_r3', |
71 |
tdatfile = 'LEV_wghc_ptmp_r4', |
72 |
sdatfile = 'LEV_wghc_salt_r4', |
73 |
# |
74 |
ctdtfile = 'ctdT_r1a', |
75 |
ctdsfile = 'ctdS_r1a', |
76 |
xbtfile = 'XBT_r2', |
77 |
argotfile = 'ARGO_T_r6', |
78 |
argosfile = 'ARGO_S_r6', |
79 |
# |
80 |
udriftfile = 'udrifter_1x1', |
81 |
vdriftfile = 'vdrifter_1x1', |
82 |
# |
83 |
udrifterrfile = 'udrifterror.err', |
84 |
vdrifterrfile = 'vdrifterror.err', |
85 |
# |
86 |
## temperrfile = 'lev01_clim_tmp_sd.bin', |
87 |
## salterrfile = 'lev01_clim_salt_sd.bin', |
88 |
temperrfile = 'insitu_estim_theta_r2.err', |
89 |
salterrfile = 'insitu_estim_salt_r2.err', |
90 |
# |
91 |
temp0errfile = 'errorTtot.bin', |
92 |
salt0errfile = 'errorStot.bin', |
93 |
# |
94 |
## hflux_errfile = 'NCEP_R6_turbflux_sd_1992-2003', |
95 |
## sflux_errfile = 'NCEP_R6_emp_sd_1992-2003', |
96 |
## tauu_errfile = 'NCEP_R6_ustr_anom_sd_1992-2003', |
97 |
## tauv_errfile = 'NCEP_R6_vstr_anom_sd_1992-2003', |
98 |
# hflux_errfile = 'totflux_s.err', |
99 |
# sflux_errfile = 'pme_s.err', |
100 |
# tauu_errfile = 'ustr_ncep_r-qscat_yr2000_s.err', |
101 |
# tauv_errfile = 'vstr_ncep_r-qscat_yr2000_s.err', |
102 |
# |
103 |
atemp_errfile = 'NCEP_R2day_tmp2m_sd_1992-2003', |
104 |
aqh_errfile = 'NCEP_R2day_spfh2m_sd_1992-2003', |
105 |
uwind_errfile = 'NCEP_R2day_u10m_sd_1992-2003', |
106 |
vwind_errfile = 'NCEP_R2day_v10m_sd_1992-2003', |
107 |
precip_errfile = 'NCEP_R2day_rain_sd_1992-2003', |
108 |
swflux_errfile = 'NCEP_R2day_sw_sd_1992-2003', |
109 |
swdown_errfile = 'NCEP_R2day_sw_sd_1992-2003', |
110 |
# |
111 |
whflux0 = 20.0, |
112 |
wsflux0 = 3.0E-8, |
113 |
wtau0 = 2.0E-2, |
114 |
watemp0 = 1.0, |
115 |
waqh0 = 0.5E-3, |
116 |
wprecip0 = 1.5E-8, |
117 |
wswflux0 = 15.0, |
118 |
wswdown0 = 15.0, |
119 |
wwind0 = 1.0, |
120 |
# |
121 |
wmean_hflux = 60.0, |
122 |
wmean_sflux = 3.2E-8, |
123 |
wmean_tau = 0.2, |
124 |
# |
125 |
wmean_atemp = 8.0, |
126 |
wmean_aqh = 4.0E-3, |
127 |
wmean_precip = 9.0E-8, |
128 |
wmean_swflux = 90.0, |
129 |
wmean_swdown = 90.0, |
130 |
wmean_wind = 8.0, |
131 |
# |
132 |
scatx_errfile = 'ustr_ncep_r-qscat_yr2000_s.err', |
133 |
scaty_errfile = 'vstr_ncep_r-qscat_yr2000_s.err', |
134 |
# |
135 |
# geoid_errfile = 'geoid_s.err', |
136 |
# geoid_errfile = '4.5cm.err', |
137 |
# geoid_errfile = 'TPJ_GRACE.err', |
138 |
geoid_errfile = 'TPJ_GRACE_r4.err', |
139 |
### geoid_covariancefile = 'egm96_1x1.pinv', |
140 |
# |
141 |
### mult_hfluxmm = 100000000000., |
142 |
# from iter53: |
143 |
mult_hfluxmm = 1000000000., |
144 |
# |
145 |
### mult_sfluxmm = 10000000000000., |
146 |
# from iter44 |
147 |
mult_sfluxmm = 100000000000 |
148 |
# |
149 |
mult_hflux = 1., |
150 |
mult_sflux = 1., |
151 |
mult_tauu = 1., |
152 |
mult_tauv = 1., |
153 |
# |
154 |
mult_atemp = 0.01, |
155 |
mult_aqh = 0.01, |
156 |
mult_uwind = 0.01, |
157 |
mult_vwind = 0.01, |
158 |
mult_precip = 0.01, |
159 |
mult_swflux = 0.01, |
160 |
mult_swdown = 0.01, |
161 |
# |
162 |
mult_hmean = 1000., |
163 |
mult_h = 1., |
164 |
mult_tp = 1., |
165 |
mult_ers = 1., |
166 |
mult_gfo = 0., |
167 |
mult_temp = 1., |
168 |
mult_salt = 1., |
169 |
mult_temp0 = 0.1, |
170 |
mult_salt0 = 0.1, |
171 |
mult_sst = 1., |
172 |
mult_sss = 0.5, |
173 |
mult_tmi = 0.5, |
174 |
mult_ctdt = 20., |
175 |
mult_ctds = 20., |
176 |
mult_xbt = 20., |
177 |
mult_argot = 20., |
178 |
mult_argos = 20., |
179 |
mult_drift = 1., |
180 |
# from iter53 weaken drift constraints: |
181 |
mult_tdrift = 0.1, |
182 |
mult_sdrift = 0.1, |
183 |
mult_wdrift = 0.1, |
184 |
mult_scatx = 1., |
185 |
mult_scaty = 1., |
186 |
# |
187 |
cost_iprec = 32, |
188 |
cost_yftype = 'RL', |
189 |
# |
190 |
& |