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 |
# |
14 |
# sshv4cost_errfile(1) = 'SO6input/SOSE_DOT_errorguess.err', |
15 |
# sshv4cost_errfile(2) = 'SO6input/SOSE_SSHanom_error.err', |
16 |
# sshv4cost_errfile(3) = 'SO6input/SOSE_SSHanom_error.err', |
17 |
# sshv4cost_errfile(4) = 'SO6input/SOSE_SSHanom_error.err', |
18 |
# sshv4cost_errfile(5) = 'SO6input/SOSE_SSHanom_error.err', |
19 |
# |
20 |
using_topex = .TRUE., |
21 |
topexstartdate1 = 20050101, |
22 |
topexstartdate2 = 000000, |
23 |
topexperiod = 86400.0, |
24 |
topexfile = 'SO6input/RADS_TPJSN_SO_r3', |
25 |
# topexfile = 'SO6input/PODAAC_TPJ_daily_1_6_anomDT', |
26 |
# topexmeanfile = 'SO6input/GGM02_DOT_col6_SO6.bin', |
27 |
# topexmeanfile = 'SO6input/SM004_GRACE3_SO6.bin', |
28 |
topexmeanfile = 'DOT_egm08.bin', |
29 |
# SO6input/DOT_DNSCMSS08-EGM08_SOSE.bin', |
30 |
# |
31 |
using_ers = .TRUE., |
32 |
ersstartdate1 = 20050101, |
33 |
ersstartdate2 = 00000, |
34 |
ersperiod = 86400.0, |
35 |
ersfile = 'SO6input/RADS_ENVERS_SO_r3', |
36 |
# ersfile = 'SO6input/ERS_ENVISAT_daily_SO6_DT', |
37 |
# ersfile = 'SO6input/ERS_ENVISAT_SO_1_6', |
38 |
|
39 |
ssh_errfile = 'SO6input/SOSE_SSHanom_error.err', |
40 |
# ssh_errfile = 'SO6input/ssh_s_SO6.err', |
41 |
### geoid_errfile = 'SO6input/TPJ_GRACE_r4_SO6.err', |
42 |
### geoid_covariancefile = 'egm96_1x1.pinv', |
43 |
geoid_errfile = 'DOT_StdDev_capped.err', |
44 |
# SO6input/SOSE_DOT_errorguess.err', |
45 |
# |
46 |
using_gfo = .TRUE., |
47 |
gfostartdate1 = 20050101, |
48 |
gfostartdate2 = 00000, |
49 |
gfoperiod = 86400.0, |
50 |
gfofile = 'SO6input/RADS_GFO_SO_r3', |
51 |
# gfofile = 'SO6input/GFO_SO_1_6_DT', |
52 |
# |
53 |
# scatxdatfile = 'QSCAT_ustr_r1', |
54 |
# scatydatfile = 'QSCAT_vstr_r1', |
55 |
# scatstartdate1 = 19920101, |
56 |
# scatstartdate2 = 60000, |
57 |
# scatperiod = 2629800, |
58 |
# |
59 |
# sststartdate1 = 20050101, |
60 |
# sststartdate2 = 00000, |
61 |
# sstdatfile = 'SO6input/sst4-path-eot.mean_SO6', |
62 |
sstdatfile = 'SO6input/TMI_AMSRE_daily_sst_SO6', |
63 |
# |
64 |
# tmistartdate1 = 19980101, |
65 |
# tmistartdate2 = 60000, |
66 |
# tmidatfile = 'SO6input/tmi_amsre.fusion.mean_SO6', |
67 |
# |
68 |
# sssstartdate1 = 19920101, |
69 |
# sssstartdate2 = 60000, |
70 |
# sssdatfile = 'SSS_tsg_1x1_monthly', |
71 |
# |
72 |
tdatfile = 'SO6input/LEV_wghc_ptmp_r4_pc_SO6.data', |
73 |
sdatfile = 'SO6input/LEV_wghc_salt_r4_pc_SO6.data', |
74 |
temp0errfile = 'SO6input/errorTtot_cap_SO6.err', |
75 |
salt0errfile = 'SO6input/errorStot_cap_SO6.err', |
76 |
# |
77 |
# ctdtfile = 'SO6input/ctdT_r1a_1_6', |
78 |
# ctdsfile = 'SO6input/ctdS_r1a_1_6', |
79 |
# xbtfile = 'SO6input/XBT_r2_1_6', |
80 |
# argotfile = 'SO6input/ARGO_T_r6_1_6', |
81 |
# argosfile = 'SO6input/ARGO_S_r6_1_6', |
82 |
# |
83 |
# udriftfile = 'udrifter_1x1', |
84 |
# vdriftfile = 'vdrifter_1x1', |
85 |
# |
86 |
# udrifterrfile = 'udrifterror.err', |
87 |
# vdrifterrfile = 'vdrifterror.err', |
88 |
# |
89 |
temperrfile = 'SO6input/errorTtot_cap_SO6_4SST.err', |
90 |
salterrfile = 'SO6input/errorStot_cap_SO6.err', |
91 |
# |
92 |
atemp_errfile = 'SO6input/NCEP_R6_tmp2m_degC_SO6_sd.err', |
93 |
aqh_errfile = 'SO6input/NCEP_R6_spfh2m_SO6_sd.err', |
94 |
uwind_errfile = 'SO6input/NCEP_R6_u10m_SO6_sd.err', |
95 |
vwind_errfile = 'SO6input/NCEP_R6_v10m_SO6_sd.err', |
96 |
precip_errfile = 'SO6input/NCEP_R6_rain_SO6_sd.err', |
97 |
swdown_errfile = 'SO6input/NCEP_R6_dsw_SO6_sd.err', |
98 |
# |
99 |
# atemp_errfile = 'SO6input/NCEP_R2day_tmp2m_sd_1992-2003_SO6.data', |
100 |
# aqh_errfile = 'SO6input/NCEP_R2day_spfh2m_sd_1992-2003_SO6.data', |
101 |
# uwind_errfile = 'SO6input/NCEP_R2day_u10m_sd_1992-2003_SO6.data', |
102 |
# vwind_errfile = 'SO6input/NCEP_R2day_v10m_sd_1992-2003_SO6.data', |
103 |
# precip_errfile = 'SO6input/NCEP_R2day_rain_sd_1992-2003_SO6.data', |
104 |
# swdown_errfile = 'SO6input/NCEP_R2day_sw_sd_1992-2003_SO6.data', |
105 |
# swflux_errfile = 'SO6input/NCEP_R2day_sw_sd_1992-2003_SO6.data', |
106 |
# |
107 |
whflux0 = 20., |
108 |
wsflux0 = 3.0E-8, |
109 |
wtau0 = 2.0E-2, |
110 |
watemp0 = 1.0, |
111 |
waqh0 = 0.5E-3, |
112 |
wprecip0 = 1.5E-8, |
113 |
wswflux0 = 15.0, |
114 |
wswdown0 = 15.0, |
115 |
wwind0 = 1.0, |
116 |
# |
117 |
wmean_hflux = 30., |
118 |
wmean_sflux = 1.6E-8, |
119 |
wmean_tau = 0.1, |
120 |
# wmean_atemp = 4.0, |
121 |
# wmean_aqh = 2.0E-3, |
122 |
# wmean_precip = 6.0E-8, |
123 |
# wmean_swflux = 60., |
124 |
wmean_swdown = 60., |
125 |
# wmean_wind = 4.0, |
126 |
wmean_atemp = 0.20, |
127 |
wmean_aqh = 1.0E-4, |
128 |
wmean_precip = 3.0E-7, |
129 |
wmean_swdown = 3., |
130 |
wmean_wind = 0.20, |
131 |
# scatx_errfile = 'ustr_ncep_r-qscat_yr2005_s.err', |
132 |
# scaty_errfile = 'vstr_ncep_r-qscat_yr2005_s.err', |
133 |
# |
134 |
mult_hflux = 1., |
135 |
mult_sflux = 1., |
136 |
mult_tauu = 1., |
137 |
mult_tauv = 1., |
138 |
# |
139 |
mult_atemp = 1.0, |
140 |
mult_aqh = 1.0, |
141 |
mult_precip = 1.0, |
142 |
mult_swflux = 1.0, |
143 |
mult_swdown = 1.0, |
144 |
mult_uwind = 1.0, |
145 |
mult_vwind = 1.0, |
146 |
# |
147 |
mult_hmean = 1.0, |
148 |
mult_h = 1., |
149 |
mult_tp = 1.0, |
150 |
mult_ers = 1.0, |
151 |
mult_gfo = 1.0, |
152 |
# mult_hmean = 3.0, |
153 |
# mult_h = 1., |
154 |
# mult_tp = 0.4, |
155 |
# mult_ers = 0.4, |
156 |
# mult_gfo = 0.4, |
157 |
mult_temp = 0.002, |
158 |
mult_salt = 0.002, |
159 |
mult_temp0 = 0.1, |
160 |
mult_salt0 = 0.1, |
161 |
mult_sst = 0.1, |
162 |
mult_sss = 1., |
163 |
mult_tmi = 1.0, |
164 |
mult_ctdt = 1., |
165 |
mult_ctds = 1., |
166 |
mult_xbt = 1., |
167 |
mult_argot = 1., |
168 |
mult_argos = 1., |
169 |
mult_drift = 1., |
170 |
mult_tdrift = 1., |
171 |
mult_sdrift = 1., |
172 |
#mult wdrift is now mult on w variance |
173 |
mult_wdrift = 1., |
174 |
mult_scatx = 1., |
175 |
mult_scaty = 1., |
176 |
mult_obcsn = 1., |
177 |
mult_smooth_ic = 1., |
178 |
mult_smooth_bc = 1., |
179 |
mult_obcsvol = 1., |
180 |
# mult ageo is now weight on w^2 |
181 |
mult_ageos = 0.00001, |
182 |
# |
183 |
cost_iprec = 32, |
184 |
cost_yftype = 'RL', |
185 |
/ |