1 |
jmc |
1.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 |
|
|
& |