1 | # $Id: $
|
---|
2 |
|
---|
3 | #****c* Configuration Files/cosmic_pp.cf *
|
---|
4 | #
|
---|
5 | # NAME
|
---|
6 | # default_pp.cf - COSMIC data configuration file for pre-processor
|
---|
7 | # implementations in ROPP
|
---|
8 | #
|
---|
9 | # SYNOPSIS
|
---|
10 | # <pp_program> ... -c cosmic_pp.cf ...
|
---|
11 | #
|
---|
12 | # DESCRIPTION
|
---|
13 | # This file reflects the configuration for the PP
|
---|
14 | # implementations within ROPP suitable for use with COSMIC data.
|
---|
15 | #
|
---|
16 | # NOTES
|
---|
17 | #
|
---|
18 | # AUTHOR
|
---|
19 | # Met Office, Exeter, UK.
|
---|
20 | # Any comments on this software should be given via the ROM SAF
|
---|
21 | # Helpdesk at http://www.romsaf.org
|
---|
22 | #
|
---|
23 | # COPYRIGHT
|
---|
24 | # (c) EUMETSAT. All rights reserved.
|
---|
25 | # For further details please refer to the file COPYRIGHT
|
---|
26 | # which you should have received as part of this distribution.
|
---|
27 | #
|
---|
28 | #****
|
---|
29 |
|
---|
30 | #-------------------------------------------------------------------------------
|
---|
31 | # 0. Output options
|
---|
32 | #-------------------------------------------------------------------------------
|
---|
33 | output_tdry = .true. ! Flag to output dry temperature
|
---|
34 |
|
---|
35 | output_diag = .false. ! Flag to output additional diagnostics
|
---|
36 |
|
---|
37 | #-------------------------------------------------------------------------------
|
---|
38 | # 1. Excess phase to bending angle processing
|
---|
39 | #-------------------------------------------------------------------------------
|
---|
40 |
|
---|
41 | # 1.1 Occultation processing method
|
---|
42 | # ---------------------------------
|
---|
43 |
|
---|
44 | # GO - use GEOMETRIC OPTICS processing to derive bending angle as a function of
|
---|
45 | # impact parameter from excess phase as a function of time.
|
---|
46 | # WO - use WAVE OPTICS (CT2 algorithm) processing to derive bending angle as a
|
---|
47 | # function of impact parameter from excess phase as a function of time.
|
---|
48 |
|
---|
49 | occ_method = WO
|
---|
50 |
|
---|
51 | # 1.2 Filtering method
|
---|
52 | # --------------------
|
---|
53 |
|
---|
54 | # optest - use OPTIMAL ESTIMATION: solution of integral equation
|
---|
55 | # slpoly - use SLIDING POLYNOMIAL
|
---|
56 |
|
---|
57 | filter_method = slpoly
|
---|
58 |
|
---|
59 | # 1.3 Smoothing bending angle profile
|
---|
60 | # -----------------------------------
|
---|
61 |
|
---|
62 |
|
---|
63 | fw_go_smooth = 3000.0 # Filter width for smoothed GO bending angles (m)
|
---|
64 |
|
---|
65 | fw_go_full = 3000.0 # Filter width for full resolution GO bending angles (m)
|
---|
66 |
|
---|
67 | fw_wo = 2000.0 # Filter width for wave optics bending angle above 7 km(m)
|
---|
68 |
|
---|
69 | fw_low = -1000.0 # Filter width for wave optics bending angle below 7 km (m)
|
---|
70 |
|
---|
71 | # 1.4 Maximum height for wave optics processing
|
---|
72 | # ---------------------------------------------
|
---|
73 |
|
---|
74 | hmax_wo = 25000.0 # Maximum height for wave optics processing (m)
|
---|
75 |
|
---|
76 | # 1.5 Data cut-off limits
|
---|
77 | # -----------------------
|
---|
78 |
|
---|
79 | Acut = 0.0 # Fractional cut-off limit for amplitude
|
---|
80 |
|
---|
81 | Pcut = -2000.0 # Cut-off limit for impact height
|
---|
82 |
|
---|
83 | Bcut = 0.1 # Cut-off limit for bending angle
|
---|
84 |
|
---|
85 | Hcut = -250000.0 # Cut-off limit for straight-line tangent altitude
|
---|
86 |
|
---|
87 | # 1.6 CT2 options
|
---|
88 | # ---------------
|
---|
89 |
|
---|
90 | CFF = 3 # Complex field filter flag (CFF = 'Pa')
|
---|
91 |
|
---|
92 | dsh = 200.0 # Shadow border width (m)
|
---|
93 |
|
---|
94 | # 1.7 Degraded L2 data flag
|
---|
95 | # -------------------------
|
---|
96 |
|
---|
97 | opt_DL2 = .true.
|
---|
98 |
|
---|
99 | # 1.8 Compute and output spectra flag
|
---|
100 | # -----------------------------------
|
---|
101 |
|
---|
102 | opt_spectra = .false.
|
---|
103 |
|
---|
104 | # 1.9 Paths to EGM96 geoid model coefficients and corrections file
|
---|
105 | # ----------------------------------------------------------------
|
---|
106 |
|
---|
107 | egm96 = ../data/egm96.dat # EGM96 coefficients file
|
---|
108 |
|
---|
109 | corr_egm96 = ../data/corrcoef.dat # Correction coefficients file
|
---|
110 |
|
---|
111 | #-------------------------------------------------------------------------------
|
---|
112 | # 1. Ionospheric correction processing
|
---|
113 | #-------------------------------------------------------------------------------
|
---|
114 |
|
---|
115 | # 1.1 Ionospheric correction method
|
---|
116 | # ---------------------------------
|
---|
117 |
|
---|
118 | # GMSIS - use MSIS climatology bending angle (searching global MSIS profiles
|
---|
119 | # for best fit profile to obs) in ionospheric correction,
|
---|
120 | # statistical optimization and bending angle to refractivity inversion.
|
---|
121 | #
|
---|
122 | # MSIS - use MSIS climatology bending angle in ionospheric correction,
|
---|
123 | # statistical optimization and bending angle to refractivity inversion.
|
---|
124 | #
|
---|
125 | # GBARO - use BAROCLIM bending angle (searching global BAROCLIM profiles
|
---|
126 | # for best fit profile to obs) in ionospheric correction,
|
---|
127 | # statistical optimization and bending angle to refractivity inversion.
|
---|
128 | #
|
---|
129 | # BARO - use BAROCLIM bending angle in ionospheric correction,
|
---|
130 | # statistical optimization and bending angle to refractivity inversion.
|
---|
131 | #
|
---|
132 | # BG - use climatology from a specified input file containing
|
---|
133 | # background temperature, pressure and humidity
|
---|
134 | # (e.g. from an NWP analysis). The input filename can be specified
|
---|
135 | # using the '-bfile' command line argument or setting 'bfile' (see 1.5).
|
---|
136 | #
|
---|
137 | # NONE - linear combination of L1 and L2 bending angles in ionospheric
|
---|
138 | # correction, no additional information above observed profile top
|
---|
139 | # in the inverse Abel to compute refractivity.
|
---|
140 |
|
---|
141 | method = GBARO # Ionospheric correction method
|
---|
142 |
|
---|
143 | # 1.2 Abel integral method
|
---|
144 | # ------------------------
|
---|
145 |
|
---|
146 | # LIN - assume linear variation of bending angle and ln(n) between
|
---|
147 | # observation levels. This algorithm is used in ROM SAF NRT processing
|
---|
148 | #
|
---|
149 | # EXP - assume exponential variation of bending angle and ln(n) between
|
---|
150 | # observation levels. This algorithm is used in ropp_fm module.
|
---|
151 |
|
---|
152 | abel = LIN
|
---|
153 |
|
---|
154 | # 1.3 Statistical optimisation method
|
---|
155 | # -----------------------------------
|
---|
156 |
|
---|
157 | # SO - statistical optimisation.
|
---|
158 | # LCSO - linear combination plus statistical optimisation.
|
---|
159 |
|
---|
160 | so_method = so
|
---|
161 |
|
---|
162 | # 1.4 Climatology model coefficients files
|
---|
163 | # --------------------------------
|
---|
164 |
|
---|
165 | msisfile = MSIS_coeff.nc # MSIS model coefficients file
|
---|
166 | mfile = BAROCLIM_coeff.nc # Model coefficients file for stat.opt.
|
---|
167 |
|
---|
168 | # 1.5 Background model temperature, humidity, pressure file
|
---|
169 | # ---------------------------------------------------------
|
---|
170 |
|
---|
171 | bfile = BG_file.nc # Background meteorology profile file (method=BG)
|
---|
172 |
|
---|
173 | #-------------------------------------------------------------------------------
|
---|
174 | # 2. Impact parameter grid
|
---|
175 | #-------------------------------------------------------------------------------
|
---|
176 |
|
---|
177 | # The ionospheric correction interpolates L1 and L2 bending angle profiles onto a
|
---|
178 | # standard grid.
|
---|
179 |
|
---|
180 | dpi = 100.0 # Step of standard impact parameter grid (m)
|
---|
181 |
|
---|
182 | #-------------------------------------------------------------------------------
|
---|
183 | # 3. Smoothing bending angle profile
|
---|
184 | #-------------------------------------------------------------------------------
|
---|
185 |
|
---|
186 | # A smoothed bending angle profile is derived compute the fit of observed bending
|
---|
187 | # angles to the model bending angle profile.
|
---|
188 |
|
---|
189 | np_smooth = 3 # Polynomial degree for smoothing regression
|
---|
190 |
|
---|
191 | fw_smooth = 1000.0 # Filter width for smoothing profile
|
---|
192 |
|
---|
193 | #-------------------------------------------------------------------------------
|
---|
194 | # 4. Model bending angle profile fit to observations
|
---|
195 | #-------------------------------------------------------------------------------
|
---|
196 |
|
---|
197 | # To avoid systematic deviations from the observed profile with climatology,
|
---|
198 | # the model profile is scaled to the observed profile by a fitting method.
|
---|
199 |
|
---|
200 | sf_method = regular # Search and fit method (convoluted or regular)
|
---|
201 |
|
---|
202 | nparm_fit = 2 # Number of parameters for model fit regression
|
---|
203 |
|
---|
204 | hmin_fit = 20000.0 # Lower limit for model fit regression
|
---|
205 |
|
---|
206 | hmax_fit = 70000.0 # Upper limit for model fit regression
|
---|
207 |
|
---|
208 | omega_fit = 0.3 # A priori standard deviation of regression factor
|
---|
209 |
|
---|
210 | #-------------------------------------------------------------------------------
|
---|
211 | # 5. Ionospheric correction and statistical optimization
|
---|
212 | #-------------------------------------------------------------------------------
|
---|
213 |
|
---|
214 | # The method described by Gorbunov (2002) is implemented to perform ionospheric
|
---|
215 | # correction with statistical optimization.
|
---|
216 |
|
---|
217 | f_width = 2000.0 # Ionospheric correction filter width
|
---|
218 |
|
---|
219 | delta_p = 20.0 # Step of homogeneous impact parameter grid
|
---|
220 |
|
---|
221 | s_smooth = 2000.0 # External ionospheric smoothing scale
|
---|
222 |
|
---|
223 | z_ion = 50000.0 # Lower height limit of ionospheric signal
|
---|
224 |
|
---|
225 | z_str = 35000.0 # Lower height limit of stratospheric signal
|
---|
226 |
|
---|
227 | z_ltr = 12000.0 # Lower height limit of tropospheric signal
|
---|
228 |
|
---|
229 | n_smooth = 11 # Number of points for smoothing (must be odd)
|
---|
230 |
|
---|
231 | model_err = -0.5 # A priori model error std.dev. (dyn.est. if negative)
|
---|
232 |
|
---|
233 | #-------------------------------------------------------------------------------
|
---|
234 | # 6. Bending angle inversion to refractivity
|
---|
235 | #-------------------------------------------------------------------------------
|
---|
236 |
|
---|
237 | # The Abel inversion is computed to retrieve refractivity from corrected
|
---|
238 | # bending angles. The corrected bending angle profile is extended
|
---|
239 | # using MSIS or BAROCLIM data above the observed profile top.
|
---|
240 |
|
---|
241 | ztop_invert = 150000.0 # Height of atmosphere top for inversion
|
---|
242 |
|
---|
243 | dzh_invert = 50.0 # Step of inversion grid above observation top
|
---|
244 |
|
---|
245 | dzr_invert = 20000.0 # Interval for regression in inversion
|
---|
246 |
|
---|