1 | <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
|
---|
2 | <html>
|
---|
3 | <head>
|
---|
4 |
|
---|
5 |
|
---|
6 |
|
---|
7 | <meta content="text/html; charset=ISO-8859-1" http-equiv="content-type">
|
---|
8 |
|
---|
9 |
|
---|
10 | <title>grib2txt converter</title>
|
---|
11 | </head>
|
---|
12 |
|
---|
13 |
|
---|
14 | <body style="background-color: rgb(255, 255, 204); color: rgb(0, 0, 0);" alink="#ee0000" link="#0000ee" vlink="#551a8b">
|
---|
15 |
|
---|
16 |
|
---|
17 | <h1 style="background-color: rgb(255, 255, 204); color: rgb(51, 51, 255);">GRIB
|
---|
18 | to TXT converter</h1>
|
---|
19 |
|
---|
20 |
|
---|
21 | This page descibes the logic of the grib2txt routine, written by Hans
|
---|
22 | Gleisner (DMI), in<br>
|
---|
23 |
|
---|
24 |
|
---|
25 | <br>
|
---|
26 |
|
---|
27 |
|
---|
28 | <span style="font-weight: bold;"> /data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/hw_grib2ropp/hgl_script/grib2txt.</span><br>
|
---|
29 |
|
---|
30 |
|
---|
31 | <br>
|
---|
32 |
|
---|
33 |
|
---|
34 | <br>
|
---|
35 |
|
---|
36 |
|
---|
37 | <font size="+1"><span style="font-weight: bold; color: rgb(51, 51, 255);">Interface:</span></font><br>
|
---|
38 |
|
---|
39 |
|
---|
40 | <br style="font-family: Courier;">
|
---|
41 |
|
---|
42 |
|
---|
43 | <font style="font-family: Courier;" size="+1">!==============================================================================!<br>
|
---|
44 | !
|
---|
45 | !<br>
|
---|
46 | !
|
---|
47 | NAME
|
---|
48 | !<br>
|
---|
49 | !
|
---|
50 | get_ecmwf_profs
|
---|
51 | !<br>
|
---|
52 | !
|
---|
53 | !<br>
|
---|
54 | !
|
---|
55 | DESCRIPTION
|
---|
56 | !<br>
|
---|
57 | ! Retrieve vertical profiles from ECMWF GRIB
|
---|
58 | files.
|
---|
59 | !<br>
|
---|
60 | !
|
---|
61 | !<br>
|
---|
62 | !
|
---|
63 | USAGE
|
---|
64 | !<br>
|
---|
65 | ! get_ecmwf_profs -gribroot
|
---|
66 | dir
|
---|
67 | &
|
---|
68 | !<br>
|
---|
69 | !
|
---|
70 | -anfc
|
---|
71 | afstrng
|
---|
72 | &
|
---|
73 | !<br>
|
---|
74 | !
|
---|
75 | -date dstrng -time
|
---|
76 | tstrng
|
---|
77 | &
|
---|
78 | !<br>
|
---|
79 | !
|
---|
80 | -lat lat -lon
|
---|
81 | lon
|
---|
82 | &
|
---|
83 | !<br>
|
---|
84 | !
|
---|
85 | -outfile
|
---|
86 | outfile
|
---|
87 | &
|
---|
88 | !<br>
|
---|
89 | !
|
---|
90 | -datdir
|
---|
91 | datadir
|
---|
92 | &
|
---|
93 | !<br>
|
---|
94 | !
|
---|
95 | [-ascii |
|
---|
96 | -ropp]
|
---|
97 | &
|
---|
98 | !<br>
|
---|
99 | !
|
---|
100 | [-sil]
|
---|
101 | &
|
---|
102 | !<br>
|
---|
103 | !
|
---|
104 | [-h]
|
---|
105 | !<br>
|
---|
106 | !
|
---|
107 | !<br>
|
---|
108 | ! get_ecmwf_profs -gribroot
|
---|
109 | dir
|
---|
110 | &
|
---|
111 | !<br>
|
---|
112 | !
|
---|
113 | -anfc
|
---|
114 | afstrng
|
---|
115 | &
|
---|
116 | !<br>
|
---|
117 | !
|
---|
118 | -infile
|
---|
119 | infile
|
---|
120 | &
|
---|
121 | !<br>
|
---|
122 | !
|
---|
123 | -datdir
|
---|
124 | datadir
|
---|
125 | &
|
---|
126 | !<br>
|
---|
127 | !
|
---|
128 | [-ascii |
|
---|
129 | -ropp]
|
---|
130 | &
|
---|
131 | !<br>
|
---|
132 | !
|
---|
133 | [-sil]
|
---|
134 | &
|
---|
135 | !<br>
|
---|
136 | !
|
---|
137 | [-h]
|
---|
138 | !<br>
|
---|
139 | !
|
---|
140 | !<br>
|
---|
141 | !
|
---|
142 | ARGUMENTS
|
---|
143 | !<br>
|
---|
144 |
|
---|
145 |
|
---|
146 | ! -gribroot
|
---|
147 | <dir> Root
|
---|
148 | directory for the GRIB file
|
---|
149 | tree.
|
---|
150 | !<br>
|
---|
151 |
|
---|
152 |
|
---|
153 | !
|
---|
154 | -anfc
|
---|
155 | <anfc> Analysis or
|
---|
156 | forecast: 'AN' or
|
---|
157 | 'FC'
|
---|
158 | !<br>
|
---|
159 | !
|
---|
160 | -date
|
---|
161 | <dstrng>
|
---|
162 | Date string, e.g.
|
---|
163 | '20080915'.
|
---|
164 | !<br>
|
---|
165 |
|
---|
166 |
|
---|
167 | !
|
---|
168 | -time
|
---|
169 | <tstrng>
|
---|
170 | Time string, e.g. '0937' or
|
---|
171 | '122534'.
|
---|
172 | !<br>
|
---|
173 | !
|
---|
174 | -lat
|
---|
175 | <lat> Latitude:
|
---|
176 | -90.0 to
|
---|
177 | +90.0.
|
---|
178 | !<br>
|
---|
179 | !
|
---|
180 | -lon
|
---|
181 | <lon> Longitude:
|
---|
182 | -180.0 to
|
---|
183 | +180.0.
|
---|
184 | !<br>
|
---|
185 | ! -outfile
|
---|
186 | <outfile> Name of output
|
---|
187 | file
|
---|
188 | !<br>
|
---|
189 |
|
---|
190 |
|
---|
191 | ! -datdir
|
---|
192 | <datadir> Root directory
|
---|
193 | to the GSM-96 data
|
---|
194 | directory.
|
---|
195 | !<br>
|
---|
196 |
|
---|
197 |
|
---|
198 | ! -infile
|
---|
199 | <infile> Name of
|
---|
200 | input file (each line contains date,
|
---|
201 | time, !<br>
|
---|
202 | !
|
---|
203 | lat, lon,
|
---|
204 | outfile).
|
---|
205 | !<br>
|
---|
206 | !
|
---|
207 | -ascii
|
---|
208 | Output stored in a ROPP
|
---|
209 | file.
|
---|
210 | !<br>
|
---|
211 | !
|
---|
212 | -ropp
|
---|
213 | Output stored in an ASCII
|
---|
214 | file.
|
---|
215 | !<br>
|
---|
216 | !
|
---|
217 | -sil
|
---|
218 | Execute
|
---|
219 | silently.
|
---|
220 | !<br>
|
---|
221 |
|
---|
222 |
|
---|
223 | !
|
---|
224 | -h
|
---|
225 | Causes a summary help to be
|
---|
226 | output.
|
---|
227 | !<br>
|
---|
228 | !
|
---|
229 | !<br>
|
---|
230 | !
|
---|
231 | DEPENDENCIES
|
---|
232 | !<br>
|
---|
233 | ! module
|
---|
234 | rom_const
|
---|
235 | !<br>
|
---|
236 | ! module
|
---|
237 | rom_utils
|
---|
238 | !<br>
|
---|
239 | ! module
|
---|
240 | rom_num
|
---|
241 | !<br>
|
---|
242 | ! module
|
---|
243 | rom_io
|
---|
244 | !<br>
|
---|
245 | ! module
|
---|
246 | grib_api
|
---|
247 | !<br>
|
---|
248 | ! module
|
---|
249 | ncdf
|
---|
250 | !<br>
|
---|
251 | !
|
---|
252 | !<br>
|
---|
253 | !
|
---|
254 | NOTES
|
---|
255 | !<br>
|
---|
256 | ! Already existing output files will be
|
---|
257 | overwritten.
|
---|
258 | !<br>
|
---|
259 | !
|
---|
260 | !<br>
|
---|
261 | !
|
---|
262 | AUTHOR
|
---|
263 | !<br>
|
---|
264 | ! Hans Gleisner
|
---|
265 | (hgl@dmi.dk)
|
---|
266 | !<br>
|
---|
267 | !
|
---|
268 | !<br>
|
---|
269 | !
|
---|
270 | VERSION
|
---|
271 | !<br>
|
---|
272 | !
|
---|
273 | 10-Jul-2012
|
---|
274 | !<br>
|
---|
275 | !
|
---|
276 | !<br>
|
---|
277 | ! CODE
|
---|
278 | DESCRIPTION:
|
---|
279 | !<br>
|
---|
280 | ! Language: Fortran
|
---|
281 | 90
|
---|
282 | !<br>
|
---|
283 | !
|
---|
284 | !<br>
|
---|
285 |
|
---|
286 |
|
---|
287 | !------------------------------------------------------------------------------!</font><br style="font-family: Courier;">
|
---|
288 |
|
---|
289 |
|
---|
290 | <br>
|
---|
291 |
|
---|
292 |
|
---|
293 | <font size="+1"><span style="font-weight: bold; color: rgb(51, 51, 255);">Steps:</span></font><br>
|
---|
294 |
|
---|
295 |
|
---|
296 | <br>
|
---|
297 |
|
---|
298 |
|
---|
299 | 1) Read required occ dates and times.<br>
|
---|
300 |
|
---|
301 |
|
---|
302 | <br>
|
---|
303 |
|
---|
304 |
|
---|
305 | 2) Choose corresponding GRIB file. <br>
|
---|
306 |
|
---|
307 |
|
---|
308 | If anfc = 'AN', this is the file with validity time within
|
---|
309 | ±3h of the
|
---|
310 | occ time, in the subdirectory of gribroot called
|
---|
311 | yyyymmddhh/000. <br style="font-family: Courier;">
|
---|
312 |
|
---|
313 |
|
---|
314 | If anfc = 'FC', this is the file with validity time within
|
---|
315 | ±3h of the
|
---|
316 | occ time, in the subdirectory of gribroot called yyyymmddhh/HHH, which
|
---|
317 | contains the short range forecast from the most
|
---|
318 | recent analysis, which is valid at either 00Z or 12Z.
|
---|
319 | <br>
|
---|
320 |
|
---|
321 |
|
---|
322 | <br>
|
---|
323 |
|
---|
324 |
|
---|
325 | Thus, for example, specifying an occultation time of 0212Z on 20110508
|
---|
326 | would result in the extraction of a background from the 00Z analysis
|
---|
327 | (ie 2011050800/000) if anfc='AN', but a background from the T+12
|
---|
328 | forecast starting from the previous day's 12Z analysis (ie
|
---|
329 | 2011050712/012) if
|
---|
330 | anfc='FC'.<br>
|
---|
331 |
|
---|
332 |
|
---|
333 | <br>
|
---|
334 |
|
---|
335 |
|
---|
336 | <span style="color: rgb(51, 51, 255);">[Explanation from
|
---|
337 | Hallgeir of DMI's ECMWF gribfile archive structure:</span><br style="color: rgb(51, 51, 255);">
|
---|
338 |
|
---|
339 |
|
---|
340 | <br style="color: rgb(51, 51, 255);">
|
---|
341 |
|
---|
342 |
|
---|
343 | <span style="color: rgb(51, 51, 255);">You may know this
|
---|
344 | already, as it is probably standard, but here comes a quick explanation
|
---|
345 | of the files in:
|
---|
346 | </span><br style="color: rgb(51, 51, 255);">
|
---|
347 |
|
---|
348 |
|
---|
349 | <span style="color: rgb(51, 51, 255);">https://svn.grassaf.org/ropp/ropp_src/branches/dev/Share/hw_grib2ropp/grib_files
|
---|
350 | </span><br style="color: rgb(51, 51, 255);">
|
---|
351 |
|
---|
352 |
|
---|
353 | <span style="color: rgb(51, 51, 255);"> </span><br style="color: rgb(51, 51, 255);">
|
---|
354 |
|
---|
355 |
|
---|
356 | <span style="color: rgb(51, 51, 255);">The 000 files are
|
---|
357 | analysis
|
---|
358 | files. The other files (006, 012, 018) are forecasts based on the 000
|
---|
359 | analysis in that directory. The numbers mean 006 hour-, 012 hour- and
|
---|
360 | 018 hour forecasts. :)
|
---|
361 | </span><br style="color: rgb(51, 51, 255);">
|
---|
362 |
|
---|
363 |
|
---|
364 | <span style="color: rgb(51, 51, 255);"> </span><br style="color: rgb(51, 51, 255);">
|
---|
365 |
|
---|
366 |
|
---|
367 | <span style="color: rgb(51, 51, 255);">We get forecasts at
|
---|
368 | 00 and 12
|
---|
369 | hours every day, and analysis at 00, 06, 12, 18 hours every day. Our
|
---|
370 | file structure for a certain day, here January 2nd 2012, therefore
|
---|
371 | looks like this:
|
---|
372 | </span><br style="color: rgb(51, 51, 255);">
|
---|
373 |
|
---|
374 |
|
---|
375 | <span style="color: rgb(51, 51, 255);"> </span><br style="color: rgb(51, 51, 255);">
|
---|
376 |
|
---|
377 |
|
---|
378 | <span style="color: rgb(51, 51, 255);">00 hour:
|
---|
379 | </span><br style="color: rgb(51, 51, 255);">
|
---|
380 |
|
---|
381 |
|
---|
382 | <span style="color: rgb(51, 51, 255);">2012010200/000
|
---|
383 | </span><br style="color: rgb(51, 51, 255);">
|
---|
384 |
|
---|
385 |
|
---|
386 | <span style="color: rgb(51, 51, 255);">2012010200/006
|
---|
387 | </span><br style="color: rgb(51, 51, 255);">
|
---|
388 |
|
---|
389 |
|
---|
390 | <span style="color: rgb(51, 51, 255);">2012010200/012
|
---|
391 | </span><br style="color: rgb(51, 51, 255);">
|
---|
392 |
|
---|
393 |
|
---|
394 | <span style="color: rgb(51, 51, 255);">2012010200/018
|
---|
395 | </span><br style="color: rgb(51, 51, 255);">
|
---|
396 |
|
---|
397 |
|
---|
398 | <span style="color: rgb(51, 51, 255);"> </span><br style="color: rgb(51, 51, 255);">
|
---|
399 |
|
---|
400 |
|
---|
401 | <span style="color: rgb(51, 51, 255);">06 hour:
|
---|
402 | </span><br style="color: rgb(51, 51, 255);">
|
---|
403 |
|
---|
404 |
|
---|
405 | <span style="color: rgb(51, 51, 255);">2012010206/000
|
---|
406 | </span><br style="color: rgb(51, 51, 255);">
|
---|
407 |
|
---|
408 |
|
---|
409 | <span style="color: rgb(51, 51, 255);"> </span><br style="color: rgb(51, 51, 255);">
|
---|
410 |
|
---|
411 |
|
---|
412 | <span style="color: rgb(51, 51, 255);">12 hour:
|
---|
413 | </span><br style="color: rgb(51, 51, 255);">
|
---|
414 |
|
---|
415 |
|
---|
416 | <span style="color: rgb(51, 51, 255);">2012010212/000
|
---|
417 | </span><br style="color: rgb(51, 51, 255);">
|
---|
418 |
|
---|
419 |
|
---|
420 | <span style="color: rgb(51, 51, 255);">2012010212/006
|
---|
421 | </span><br style="color: rgb(51, 51, 255);">
|
---|
422 |
|
---|
423 |
|
---|
424 | <span style="color: rgb(51, 51, 255);">2012010212/012
|
---|
425 | </span><br style="color: rgb(51, 51, 255);">
|
---|
426 |
|
---|
427 |
|
---|
428 | <span style="color: rgb(51, 51, 255);">2012010212/018
|
---|
429 | </span><br style="color: rgb(51, 51, 255);">
|
---|
430 |
|
---|
431 |
|
---|
432 | <span style="color: rgb(51, 51, 255);"> </span><br style="color: rgb(51, 51, 255);">
|
---|
433 |
|
---|
434 |
|
---|
435 | <span style="color: rgb(51, 51, 255);">18 hour:
|
---|
436 | </span><br style="color: rgb(51, 51, 255);">
|
---|
437 |
|
---|
438 |
|
---|
439 | <span style="color: rgb(51, 51, 255);">2012010218/000
|
---|
440 | </span><br style="color: rgb(51, 51, 255);">
|
---|
441 |
|
---|
442 |
|
---|
443 | <span style="color: rgb(51, 51, 255);"> </span><br style="color: rgb(51, 51, 255);">
|
---|
444 |
|
---|
445 |
|
---|
446 | <span style="color: rgb(51, 51, 255);">In other words, the
|
---|
447 | files you have gotten in svn is only based on one single analysis...]</span><br>
|
---|
448 |
|
---|
449 |
|
---|
450 | <br>
|
---|
451 |
|
---|
452 |
|
---|
453 | All we have are:<br>
|
---|
454 |
|
---|
455 |
|
---|
456 | <br>
|
---|
457 |
|
---|
458 |
|
---|
459 | idculv@eld037:> ls -ltR
|
---|
460 | /net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files<br>
|
---|
461 |
|
---|
462 |
|
---|
463 | /net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files:<br>
|
---|
464 |
|
---|
465 |
|
---|
466 | total 6<br>
|
---|
467 |
|
---|
468 |
|
---|
469 | drwxr-xr-x 4 idculv satsense 50 Sep 17 12:15 grib2<br>
|
---|
470 |
|
---|
471 |
|
---|
472 | drwxr-xr-x 4 idculv satsense 50 Sep 17 12:09 grib1<br>
|
---|
473 |
|
---|
474 |
|
---|
475 | drwxr-xr-x 4 idculv satsense 50 Sep 17 12:02 grib2_double<br>
|
---|
476 |
|
---|
477 |
|
---|
478 | <br>
|
---|
479 |
|
---|
480 |
|
---|
481 | /net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib2:<br>
|
---|
482 |
|
---|
483 |
|
---|
484 | total 2<br>
|
---|
485 |
|
---|
486 |
|
---|
487 | drwxr-xr-x 3 idculv satsense 106 Sep 17 12:22 2011050800<br>
|
---|
488 |
|
---|
489 |
|
---|
490 | <br>
|
---|
491 |
|
---|
492 |
|
---|
493 | /net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib2/2011050800:<br>
|
---|
494 |
|
---|
495 |
|
---|
496 | total 89390<br>
|
---|
497 |
|
---|
498 |
|
---|
499 | -rw-r--r-- 1 idculv satsense 18056976 Sep 17 12:22 018<br>
|
---|
500 |
|
---|
501 |
|
---|
502 | -rw-r--r-- 1 idculv satsense 18056976 Sep 17 12:22 006<br>
|
---|
503 |
|
---|
504 |
|
---|
505 | -rw-r--r-- 1 idculv satsense 18056976 Sep 17 12:22 012<br>
|
---|
506 |
|
---|
507 |
|
---|
508 | -rw-r--r-- 1 idculv satsense 18188228 Sep 17 12:22 000<br>
|
---|
509 |
|
---|
510 |
|
---|
511 | <br>
|
---|
512 |
|
---|
513 |
|
---|
514 | /net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib1:<br>
|
---|
515 |
|
---|
516 |
|
---|
517 | total 2<br>
|
---|
518 |
|
---|
519 |
|
---|
520 | drwxr-xr-x 3 idculv satsense 106 Sep 17 12:15 2011050800<br>
|
---|
521 |
|
---|
522 |
|
---|
523 | <br>
|
---|
524 |
|
---|
525 |
|
---|
526 | /net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib1/2011050800:<br>
|
---|
527 |
|
---|
528 |
|
---|
529 | total 89294<br>
|
---|
530 |
|
---|
531 |
|
---|
532 | -rw-r--r-- 1 idculv satsense 18036480 Sep 17 12:15 018<br>
|
---|
533 |
|
---|
534 |
|
---|
535 | -rw-r--r-- 1 idculv satsense 18036480 Sep 17 12:15 006<br>
|
---|
536 |
|
---|
537 |
|
---|
538 | -rw-r--r-- 1 idculv satsense 18036480 Sep 17 12:15 012<br>
|
---|
539 |
|
---|
540 |
|
---|
541 | -rw-r--r-- 1 idculv satsense 18167620 Sep 17 12:15 000<br>
|
---|
542 |
|
---|
543 |
|
---|
544 | <br>
|
---|
545 |
|
---|
546 |
|
---|
547 | /net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib2_double:<br>
|
---|
548 |
|
---|
549 |
|
---|
550 | total 2<br>
|
---|
551 |
|
---|
552 |
|
---|
553 | drwxr-xr-x 3 idculv satsense 106 Sep 17 12:09 2011070200<br>
|
---|
554 |
|
---|
555 |
|
---|
556 | <br>
|
---|
557 |
|
---|
558 |
|
---|
559 | /net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib2_double/2011070200:<br>
|
---|
560 |
|
---|
561 |
|
---|
562 | total 102086<br>
|
---|
563 |
|
---|
564 |
|
---|
565 | -rw-r--r-- 1 idculv satsense 18056976 Sep 17 12:09 018<br>
|
---|
566 |
|
---|
567 |
|
---|
568 | -rw-r--r-- 1 idculv satsense 28417536 Sep 17 12:09 006<br>
|
---|
569 |
|
---|
570 |
|
---|
571 | -rw-r--r-- 1 idculv satsense 18056976 Sep 17 12:09 012<br>
|
---|
572 |
|
---|
573 |
|
---|
574 | -rw-r--r-- 1 idculv satsense 18188228 Sep 17 12:09 000<br>
|
---|
575 |
|
---|
576 |
|
---|
577 | <br>
|
---|
578 |
|
---|
579 |
|
---|
580 | 3) Call read_ECMWF_profs. This redoes (?) the calculation of
|
---|
581 | the gribfilename, and checks for its existence. It then opens
|
---|
582 | gribfile (if available), which contains {lnsp,Φ0, T(i), q(i)},
|
---|
583 | ie 2N_lev+2 fields. Currently there is a check that it only
|
---|
584 | contains this many fields, so we'd need to switch that off for ROPP,
|
---|
585 | which should presumably be able to handle data from a standard ECMWF
|
---|
586 | GRIB file. <br>
|
---|
587 |
|
---|
588 |
|
---|
589 | <br>
|
---|
590 |
|
---|
591 |
|
---|
592 | 4) This calls read_ECMWF_3Dgrid, which redoes (?) the check that the
|
---|
593 | variable 'pv' is of size 2N_lev+2, and if so, it extracts the level
|
---|
594 | coeffs Ak and Bk from the 1st and 2nd halves of this array (in reverse
|
---|
595 | order). It then reads {lnsp, Φ0, T(i),
|
---|
596 | q(i)}, converts to {P0=exp(lnsp), Z0=Φ0/g_wmo, T, Q(g/kg)=q(kg/kg)*1000}, and currently fails if it comes across any
|
---|
597 | other fields - again, we may want to switch this off in the ROPP
|
---|
598 | implementation. It also seems to extract the lat and the lon every
|
---|
599 | time? Could be taken out of loop, since they're the same
|
---|
600 | here. <br>
|
---|
601 |
|
---|
602 |
|
---|
603 | <br>
|
---|
604 |
|
---|
605 |
|
---|
606 | 5) 3D pressure field on full levels, P, calculated according to<br>
|
---|
607 |
|
---|
608 |
|
---|
609 | P(ilat,ilon,ilev) = 0.5d0*(ak(ilev) +
|
---|
610 | bk(ilev)*P0(ilat,ilon)) + 0.5d0*(ak(ilev+1) + bk(ilev+1)*P0(ilat,ilon))<br>
|
---|
611 |
|
---|
612 |
|
---|
613 | <br>
|
---|
614 |
|
---|
615 |
|
---|
616 | 6) 3D GPH field on full levels, Z, calculated by integrating R_dry
|
---|
617 | * Tvirt / g_wmo from Z0. If it didn't find Z0, then it
|
---|
618 | returns Z0=Z=mdi. There's some complicated business involving
|
---|
619 | call read_ECMWF_orography, which as far as I can tell is just looking
|
---|
620 | for Z0 again. <br>
|
---|
621 |
|
---|
622 | <br>
|
---|
623 |
|
---|
624 | 7) It also recalcalutes Z_grb if anfc='FC', upon return to subroutine read_ECMWF_profs. <br>
|
---|
625 |
|
---|
626 | <br>
|
---|
627 |
|
---|
628 | 8) It locates the target point's 4 nearest neighbours (in lat & lon) and bilinearly interpolates {Z0, P0, Z(i), P(i), T(i),
|
---|
629 | Q(i)} to it.
|
---|
630 | <br>
|
---|
631 |
|
---|
632 | <br>
|
---|
633 |
|
---|
634 | 9) It computes orthometric (mean-sea level) heights, H, for the profiles. <br>
|
---|
635 |
|
---|
636 | <br>
|
---|
637 | a) u(iocc) = GeoidHeight(datadir, lat(iocc), lon(iocc)) =
|
---|
638 | the height of the EGM96 geoid above the WGS84 ellipsoid (ie, the
|
---|
639 | undulation). The location of the dir containing the data is under user
|
---|
640 | control. We could put this somewhere central in ROPP.
|
---|
641 | Copyright?<br>
|
---|
642 | <br>
|
---|
643 | b) ! call curvature(lat(iocc), lon(iocc), az, Rcoc, Rc) is the
|
---|
644 | ROPP routine, with az=0. Why? - it can make a huge difference.
|
---|
645 | Presumably that's why it's <span style="color: rgb(255, 0, 0);">not currently called.</span> If it were, it would return Rcoc and Rc. <br>
|
---|
646 | <br>
|
---|
647 | c) call geopot2msl_ropp(lat(iocc), u(iocc), Z(iocc,:),
|
---|
648 | H(iocc,:)). Currently this is basically the ROPP routine H =
|
---|
649 | geopotential2geometric(lat, Z), ie the geoid height isn't used. So we
|
---|
650 | convert Z wrt geoid to H wrt geoid (the MSL one), using Somigliana's
|
---|
651 | R_eff and g.<br>
|
---|
652 | <br>
|
---|
653 | 10) It then writes a profile in ascii format:<br>
|
---|
654 | <br>
|
---|
655 | write(9,'(A)') '# Profile pressure (hPa), temperature (K), and specific humidity (g/kg).'<br>
|
---|
656 | do ilev=1,Nlev<br>
|
---|
657 | write(9,'(F8.3,1X,F7.2,1X,F8.5)') Pecm(i,ilev)/100.0d0, Tecm(i,ilev), Qecm(i,ilev)<br>
|
---|
658 | enddo<br>
|
---|
659 | <br>
|
---|
660 | (So H isn't used at present).<br>
|
---|
661 | <br>
|
---|
662 | <br>
|
---|
663 | <span style="font-weight: bold; color: rgb(51, 51, 255);">Example:</span><br>
|
---|
664 | <br>
|
---|
665 | cd /data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/hgl_script/grib2txt<br>
|
---|
666 | <br>
|
---|
667 | ./compile_getprofs.bash<br>
|
---|
668 | <br>
|
---|
669 | ./get_ecmwf_profs.x \<br>
|
---|
670 | -gribroot '/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib2' \<br>
|
---|
671 | -datdir '/net/data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/hgl_script/ancillary' \<br>
|
---|
672 | -anfc 'AN' \<br>
|
---|
673 | -date '20110508' \<br>
|
---|
674 | -time '0212' \<br>
|
---|
675 | -lat 50.0 \<br>
|
---|
676 | -lon -3.0 \<br>
|
---|
677 | -ascii \<br>
|
---|
678 | -outfile 'example1.out'<br>
|
---|
679 | <br>
|
---|
680 | produces:<br>
|
---|
681 | <br>
|
---|
682 | <br>
|
---|
683 | #
|
---|
684 | /data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib2/2011050800/000
|
---|
685 | at LAT: 50.000 LON: -3.000<br>
|
---|
686 | # Surface geopotential (m), pressure (hPa), LS mask, and Nlev.<br>
|
---|
687 | -1.43 1007.129 -1 91<br>
|
---|
688 | # Profile pressure (hPa), temperature (K), and specific humidity (g/kg).<br>
|
---|
689 | 1005.936 286.18 8.83387<br>
|
---|
690 | 1003.050 286.13 8.92123<br>
|
---|
691 | 999.052 286.00 8.99149<br>
|
---|
692 | 993.804 285.83 9.07474<br>
|
---|
693 | 987.317 285.79 9.11437<br>
|
---|
694 | 979.704 285.56 9.12646<br>
|
---|
695 | 970.818 285.31 9.04204<br>
|
---|
696 | 960.400 285.02 8.94932<br>
|
---|
697 | 948.421 284.64 8.79558<br>
|
---|
698 | 934.885 284.54 8.47089<br>
|
---|
699 | 919.795 283.73 8.19792<br>
|
---|
700 | 903.156 282.82 7.83847<br>
|
---|
701 | 885.007 282.13 6.77560<br>
|
---|
702 | 865.423 280.86 6.09711<br>
|
---|
703 | 844.477 279.39 5.66150<br>
|
---|
704 | 822.233 277.96 6.10021<br>
|
---|
705 | 798.780 275.75 5.28384<br>
|
---|
706 | 774.238 274.14 5.06247<br>
|
---|
707 | 748.725 273.02 4.77243<br>
|
---|
708 | 722.344 271.41 4.45200<br>
|
---|
709 | 695.212 269.52 4.02497<br>
|
---|
710 | 667.472 267.62 3.54029<br>
|
---|
711 | 639.265 265.44 3.09754<br>
|
---|
712 | 610.708 263.43 2.81931<br>
|
---|
713 | 581.927 261.62 2.55990<br>
|
---|
714 | 553.069 259.44 2.14411<br>
|
---|
715 | 524.355 256.91 1.67860<br>
|
---|
716 | 496.066 254.23 1.12169<br>
|
---|
717 | 468.478 250.35 0.85643<br>
|
---|
718 | 441.848 247.22 0.45206<br>
|
---|
719 | 416.415 244.90 0.18884<br>
|
---|
720 | 392.282 242.20 0.07058<br>
|
---|
721 | 369.411 238.82 0.04959<br>
|
---|
722 | 347.751 235.26 0.05150<br>
|
---|
723 | 327.245 231.91 0.04244<br>
|
---|
724 | 307.838 228.77 0.03338<br>
|
---|
725 | 289.477 225.85 0.02766<br>
|
---|
726 | 272.111 223.53 0.02314<br>
|
---|
727 | 255.693 221.78 0.02003<br>
|
---|
728 | 240.176 220.33 0.01669<br>
|
---|
729 | 225.517 218.46 0.01454<br>
|
---|
730 | 211.671 216.09 0.01329<br>
|
---|
731 | 198.600 214.65 0.01097<br>
|
---|
732 | 186.264 216.03 0.00703<br>
|
---|
733 | 174.627 217.57 0.00450<br>
|
---|
734 | 163.652 220.72 0.00328<br>
|
---|
735 | 153.307 222.11 0.00311<br>
|
---|
736 | 143.559 220.27 0.00294<br>
|
---|
737 | 134.377 218.90 0.00288<br>
|
---|
738 | 125.732 218.71 0.00289<br>
|
---|
739 | 117.586 218.55 0.00290<br>
|
---|
740 | 109.887 217.88 0.00291<br>
|
---|
741 | 102.579 216.87 0.00288<br>
|
---|
742 | 95.616 216.75 0.00282<br>
|
---|
743 | 88.958 217.58 0.00281<br>
|
---|
744 | 82.572 218.56 0.00282<br>
|
---|
745 | 76.429 218.52 0.00284<br>
|
---|
746 | 70.506 217.73 0.00285<br>
|
---|
747 | 64.798 217.80 0.00286<br>
|
---|
748 | 59.315 217.86 0.00288<br>
|
---|
749 | 54.070 218.37 0.00291<br>
|
---|
750 | 49.073 219.28 0.00293<br>
|
---|
751 | 44.331 220.09 0.00294<br>
|
---|
752 | 39.848 219.71 0.00294<br>
|
---|
753 | 35.628 218.64 0.00295<br>
|
---|
754 | 31.674 218.01 0.00296<br>
|
---|
755 | 27.986 218.12 0.00299<br>
|
---|
756 | 24.565 218.94 0.00302<br>
|
---|
757 | 21.410 220.96 0.00308<br>
|
---|
758 | 18.518 223.01 0.00319<br>
|
---|
759 | 15.884 223.78 0.00326<br>
|
---|
760 | 13.504 224.73 0.00326<br>
|
---|
761 | 11.369 226.88 0.00326<br>
|
---|
762 | 9.471 231.47 0.00326<br>
|
---|
763 | 7.799 236.67 0.00330<br>
|
---|
764 | 6.342 238.83 0.00337<br>
|
---|
765 | 5.086 244.75 0.00346<br>
|
---|
766 | 4.018 249.20 0.00354<br>
|
---|
767 | 3.121 256.42 0.00362<br>
|
---|
768 | 2.380 263.96 0.00373<br>
|
---|
769 | 1.778 268.75 0.00383<br>
|
---|
770 | 1.299 271.71 0.00391<br>
|
---|
771 | 0.924 272.85 0.00396<br>
|
---|
772 | 0.640 270.36 0.00398<br>
|
---|
773 | 0.428 265.95 0.00400<br>
|
---|
774 | 0.277 260.09 0.00402<br>
|
---|
775 | 0.172 251.76 0.00401<br>
|
---|
776 | 0.101 237.05 0.00406<br>
|
---|
777 | 0.057 217.36 0.00339<br>
|
---|
778 | 0.030 201.11 0.00241<br>
|
---|
779 | 0.010 187.84 0.00156<br>
|
---|
780 | <br>
|
---|
781 | This agrees pretty much with:<br>
|
---|
782 | <br>
|
---|
783 | /data/nwp1/idculv/GRIB/grib_api-1.9.9_dir/bin<br>
|
---|
784 | idculv@eld037:> grib_ls -l50,-3,1 $file |sed
|
---|
785 | -es/"2
|
---|
786 | ecmf
|
---|
787 | 20110508
|
---|
788 | an
|
---|
789 | regular_ll
|
---|
790 | 0
|
---|
791 | hybrid"/""/ |sort -nr<br>
|
---|
792 | <br>
|
---|
793 | 184 of 184 total grib messages in 1 files<br>
|
---|
794 | 184 of 184 grib messages in /data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib2/2011050800/000<br>
|
---|
795 |
|
---|
796 | 91
|
---|
797 | t
|
---|
798 | grid_simple 286.182 <br>
|
---|
799 |
|
---|
800 | 91
|
---|
801 | q
|
---|
802 | grid_simple 0.00883387 <br>
|
---|
803 |
|
---|
804 | 90
|
---|
805 | t
|
---|
806 | grid_simple 286.133 <br>
|
---|
807 |
|
---|
808 | 90
|
---|
809 | q
|
---|
810 | grid_simple 0.00892123 <br>
|
---|
811 |
|
---|
812 | 89
|
---|
813 | t
|
---|
814 | grid_simple 286.004 <br>
|
---|
815 |
|
---|
816 | 89
|
---|
817 | q
|
---|
818 | grid_simple 0.00899149 <br>
|
---|
819 |
|
---|
820 | 88
|
---|
821 | t
|
---|
822 | grid_simple 285.828 <br>
|
---|
823 |
|
---|
824 | 88
|
---|
825 | q
|
---|
826 | grid_simple 0.00907474 <br>
|
---|
827 |
|
---|
828 | 87
|
---|
829 | t
|
---|
830 | grid_simple 285.791 <br>
|
---|
831 |
|
---|
832 | 87
|
---|
833 | q
|
---|
834 | grid_simple 0.00911437 <br>
|
---|
835 |
|
---|
836 | 86
|
---|
837 | t
|
---|
838 | grid_simple 285.562 <br>
|
---|
839 |
|
---|
840 | 86
|
---|
841 | q
|
---|
842 | grid_simple 0.00912646 <br>
|
---|
843 |
|
---|
844 | 85
|
---|
845 | t
|
---|
846 | grid_simple 285.31 <br>
|
---|
847 |
|
---|
848 | 85
|
---|
849 | q
|
---|
850 | grid_simple 0.00904204 <br>
|
---|
851 |
|
---|
852 | 84
|
---|
853 | t
|
---|
854 | grid_simple 285.02 <br>
|
---|
855 |
|
---|
856 | 84
|
---|
857 | q
|
---|
858 | grid_simple 0.00894932 <br>
|
---|
859 |
|
---|
860 | 83
|
---|
861 | t
|
---|
862 | grid_simple 284.644 <br>
|
---|
863 |
|
---|
864 | 83
|
---|
865 | q
|
---|
866 | grid_simple 0.00879558 <br>
|
---|
867 |
|
---|
868 | 82
|
---|
869 | t
|
---|
870 | grid_simple 284.54 <br>
|
---|
871 |
|
---|
872 | 82
|
---|
873 | q
|
---|
874 | grid_simple 0.00847089 <br>
|
---|
875 |
|
---|
876 | 81
|
---|
877 | t
|
---|
878 | grid_simple 283.734 <br>
|
---|
879 |
|
---|
880 | 81
|
---|
881 | q
|
---|
882 | grid_simple 0.00819792 <br>
|
---|
883 |
|
---|
884 | 80
|
---|
885 | t
|
---|
886 | grid_simple 282.816 <br>
|
---|
887 |
|
---|
888 | 80
|
---|
889 | q
|
---|
890 | grid_simple 0.00783847 <br>
|
---|
891 |
|
---|
892 | 79
|
---|
893 | t
|
---|
894 | grid_simple 282.131 <br>
|
---|
895 |
|
---|
896 | 79
|
---|
897 | q
|
---|
898 | grid_simple 0.0067756 <br>
|
---|
899 |
|
---|
900 | 78
|
---|
901 | t
|
---|
902 | grid_simple 280.865 <br>
|
---|
903 |
|
---|
904 | 78
|
---|
905 | q
|
---|
906 | grid_simple 0.00609711 <br>
|
---|
907 |
|
---|
908 | 77
|
---|
909 | t
|
---|
910 | grid_simple 279.389 <br>
|
---|
911 |
|
---|
912 | 77
|
---|
913 | q
|
---|
914 | grid_simple 0.0056615 <br>
|
---|
915 |
|
---|
916 | 76
|
---|
917 | t
|
---|
918 | grid_simple 277.959 <br>
|
---|
919 |
|
---|
920 | 76
|
---|
921 | q
|
---|
922 | grid_simple 0.00610021 <br>
|
---|
923 |
|
---|
924 | 75
|
---|
925 | t
|
---|
926 | grid_simple 275.75 <br>
|
---|
927 |
|
---|
928 | 75
|
---|
929 | q
|
---|
930 | grid_simple 0.00528384 <br>
|
---|
931 |
|
---|
932 | 74
|
---|
933 | t
|
---|
934 | grid_simple 274.142 <br>
|
---|
935 |
|
---|
936 | 74
|
---|
937 | q
|
---|
938 | grid_simple 0.00506247 <br>
|
---|
939 |
|
---|
940 | 73
|
---|
941 | t
|
---|
942 | grid_simple 273.019 <br>
|
---|
943 |
|
---|
944 | 73
|
---|
945 | q
|
---|
946 | grid_simple 0.00477243 <br>
|
---|
947 |
|
---|
948 | 72
|
---|
949 | t
|
---|
950 | grid_simple 271.411 <br>
|
---|
951 |
|
---|
952 | 72
|
---|
953 | q
|
---|
954 | grid_simple 0.004452 <br>
|
---|
955 |
|
---|
956 | 71
|
---|
957 | t
|
---|
958 | grid_simple 269.523 <br>
|
---|
959 |
|
---|
960 | 71
|
---|
961 | q
|
---|
962 | grid_simple 0.00402497 <br>
|
---|
963 |
|
---|
964 | 70
|
---|
965 | t
|
---|
966 | grid_simple 267.619 <br>
|
---|
967 |
|
---|
968 | 70
|
---|
969 | q
|
---|
970 | grid_simple 0.00354029 <br>
|
---|
971 |
|
---|
972 | 69
|
---|
973 | t
|
---|
974 | grid_simple 265.44 <br>
|
---|
975 |
|
---|
976 | 69
|
---|
977 | q
|
---|
978 | grid_simple 0.00309754 <br>
|
---|
979 |
|
---|
980 | 68
|
---|
981 | t
|
---|
982 | grid_simple 263.427 <br>
|
---|
983 |
|
---|
984 | 68
|
---|
985 | q
|
---|
986 | grid_simple 0.00281931 <br>
|
---|
987 |
|
---|
988 | 67
|
---|
989 | t
|
---|
990 | grid_simple 261.623 <br>
|
---|
991 |
|
---|
992 | 67
|
---|
993 | q
|
---|
994 | grid_simple 0.0025599 <br>
|
---|
995 |
|
---|
996 | 66
|
---|
997 | t
|
---|
998 | grid_simple 259.445 <br>
|
---|
999 |
|
---|
1000 | 66
|
---|
1001 | q
|
---|
1002 | grid_simple 0.00214411 <br>
|
---|
1003 |
|
---|
1004 | 65
|
---|
1005 | t
|
---|
1006 | grid_simple 256.911 <br>
|
---|
1007 |
|
---|
1008 | 65
|
---|
1009 | q
|
---|
1010 | grid_simple 0.0016786 <br>
|
---|
1011 |
|
---|
1012 | 64
|
---|
1013 | t
|
---|
1014 | grid_simple 254.227 <br>
|
---|
1015 |
|
---|
1016 | 64
|
---|
1017 | q
|
---|
1018 | grid_simple 0.00112169 <br>
|
---|
1019 |
|
---|
1020 | 63
|
---|
1021 | t
|
---|
1022 | grid_simple 250.347 <br>
|
---|
1023 |
|
---|
1024 | 63
|
---|
1025 | q
|
---|
1026 | grid_simple 0.00085643 <br>
|
---|
1027 |
|
---|
1028 | 62
|
---|
1029 | t
|
---|
1030 | grid_simple 247.216 <br>
|
---|
1031 |
|
---|
1032 | 62
|
---|
1033 | q
|
---|
1034 | grid_simple 0.000452057 <br>
|
---|
1035 |
|
---|
1036 | 61
|
---|
1037 | t
|
---|
1038 | grid_simple 244.898 <br>
|
---|
1039 |
|
---|
1040 | 61
|
---|
1041 | q
|
---|
1042 | grid_simple 0.000188837 <br>
|
---|
1043 |
|
---|
1044 | 60
|
---|
1045 | t
|
---|
1046 | grid_simple 242.201 <br>
|
---|
1047 |
|
---|
1048 | 60
|
---|
1049 | q
|
---|
1050 | grid_simple 7.05798e-05 <br>
|
---|
1051 |
|
---|
1052 | 59
|
---|
1053 | t
|
---|
1054 | grid_simple 238.817 <br>
|
---|
1055 |
|
---|
1056 | 59
|
---|
1057 | q
|
---|
1058 | grid_simple 4.95934e-05 <br>
|
---|
1059 |
|
---|
1060 | 58
|
---|
1061 | t
|
---|
1062 | grid_simple 235.257 <br>
|
---|
1063 |
|
---|
1064 | 58
|
---|
1065 | q
|
---|
1066 | grid_simple 5.14984e-05 <br>
|
---|
1067 |
|
---|
1068 | 57
|
---|
1069 | t
|
---|
1070 | grid_simple 231.915 <br>
|
---|
1071 |
|
---|
1072 | 57
|
---|
1073 | q
|
---|
1074 | grid_simple 4.24385e-05 <br>
|
---|
1075 |
|
---|
1076 | 56
|
---|
1077 | t
|
---|
1078 | grid_simple 228.768 <br>
|
---|
1079 |
|
---|
1080 | 56
|
---|
1081 | q
|
---|
1082 | grid_simple 3.33786e-05 <br>
|
---|
1083 |
|
---|
1084 | 55
|
---|
1085 | t
|
---|
1086 | grid_simple 225.854 <br>
|
---|
1087 |
|
---|
1088 | 55
|
---|
1089 | q
|
---|
1090 | grid_simple 2.766e-05 <br>
|
---|
1091 |
|
---|
1092 | 54
|
---|
1093 | t
|
---|
1094 | grid_simple 223.533 <br>
|
---|
1095 |
|
---|
1096 | 54
|
---|
1097 | q
|
---|
1098 | grid_simple 2.31432e-05 <br>
|
---|
1099 |
|
---|
1100 | 53
|
---|
1101 | t
|
---|
1102 | grid_simple 221.777 <br>
|
---|
1103 |
|
---|
1104 | 53
|
---|
1105 | q
|
---|
1106 | grid_simple 2.00323e-05 <br>
|
---|
1107 |
|
---|
1108 | 52
|
---|
1109 | t
|
---|
1110 | grid_simple 220.329 <br>
|
---|
1111 |
|
---|
1112 | 52
|
---|
1113 | q
|
---|
1114 | grid_simple 1.66893e-05 <br>
|
---|
1115 |
|
---|
1116 | 51
|
---|
1117 | t
|
---|
1118 | grid_simple 218.458 <br>
|
---|
1119 |
|
---|
1120 | 51
|
---|
1121 | q
|
---|
1122 | grid_simple 1.45435e-05 <br>
|
---|
1123 |
|
---|
1124 | 50
|
---|
1125 | t
|
---|
1126 | grid_simple 216.088 <br>
|
---|
1127 |
|
---|
1128 | 50
|
---|
1129 | q
|
---|
1130 | grid_simple 1.32918e-05 <br>
|
---|
1131 |
|
---|
1132 | 49
|
---|
1133 | t
|
---|
1134 | grid_simple 214.647 <br>
|
---|
1135 |
|
---|
1136 | 49
|
---|
1137 | q
|
---|
1138 | grid_simple 1.09673e-05 <br>
|
---|
1139 |
|
---|
1140 | 48
|
---|
1141 | t
|
---|
1142 | grid_simple 216.025 <br>
|
---|
1143 |
|
---|
1144 | 48
|
---|
1145 | q
|
---|
1146 | grid_simple 7.03335e-06 <br>
|
---|
1147 |
|
---|
1148 | 47
|
---|
1149 | t
|
---|
1150 | grid_simple 217.567 <br>
|
---|
1151 |
|
---|
1152 | 47
|
---|
1153 | q
|
---|
1154 | grid_simple 4.50015e-06 <br>
|
---|
1155 |
|
---|
1156 | 46
|
---|
1157 | t
|
---|
1158 | grid_simple 220.724 <br>
|
---|
1159 |
|
---|
1160 | 46
|
---|
1161 | q
|
---|
1162 | grid_simple 3.27826e-06 <br>
|
---|
1163 |
|
---|
1164 | 45
|
---|
1165 | t
|
---|
1166 | grid_simple 222.106 <br>
|
---|
1167 |
|
---|
1168 | 45
|
---|
1169 | q
|
---|
1170 | grid_simple 3.10689e-06 <br>
|
---|
1171 |
|
---|
1172 | 44
|
---|
1173 | t
|
---|
1174 | grid_simple 220.271 <br>
|
---|
1175 |
|
---|
1176 | 44
|
---|
1177 | q
|
---|
1178 | grid_simple 2.93553e-06 <br>
|
---|
1179 |
|
---|
1180 | 43
|
---|
1181 | t
|
---|
1182 | grid_simple 218.904 <br>
|
---|
1183 |
|
---|
1184 | 43
|
---|
1185 | q
|
---|
1186 | grid_simple 2.87592e-06 <br>
|
---|
1187 |
|
---|
1188 | 42
|
---|
1189 | t
|
---|
1190 | grid_simple 218.715 <br>
|
---|
1191 |
|
---|
1192 | 42
|
---|
1193 | q
|
---|
1194 | grid_simple 2.8871e-06 <br>
|
---|
1195 |
|
---|
1196 | 41
|
---|
1197 | t
|
---|
1198 | grid_simple 218.554 <br>
|
---|
1199 |
|
---|
1200 | 41
|
---|
1201 | q
|
---|
1202 | grid_simple 2.90386e-06 <br>
|
---|
1203 |
|
---|
1204 | 40
|
---|
1205 | t
|
---|
1206 | grid_simple 217.875 <br>
|
---|
1207 |
|
---|
1208 | 40
|
---|
1209 | q
|
---|
1210 | grid_simple 2.90945e-06 <br>
|
---|
1211 |
|
---|
1212 | 39
|
---|
1213 | t
|
---|
1214 | grid_simple 216.871 <br>
|
---|
1215 |
|
---|
1216 | 39
|
---|
1217 | q
|
---|
1218 | grid_simple 2.87965e-06 <br>
|
---|
1219 |
|
---|
1220 | 38
|
---|
1221 | t
|
---|
1222 | grid_simple 216.75 <br>
|
---|
1223 |
|
---|
1224 | 38
|
---|
1225 | q
|
---|
1226 | grid_simple 2.81602e-06 <br>
|
---|
1227 |
|
---|
1228 | 37
|
---|
1229 | t
|
---|
1230 | grid_simple 217.58 <br>
|
---|
1231 |
|
---|
1232 | 37
|
---|
1233 | q
|
---|
1234 | grid_simple 2.80837e-06 <br>
|
---|
1235 |
|
---|
1236 | 36
|
---|
1237 | t
|
---|
1238 | grid_simple 218.555 <br>
|
---|
1239 |
|
---|
1240 | 36
|
---|
1241 | q
|
---|
1242 | grid_simple 2.82151e-06 <br>
|
---|
1243 |
|
---|
1244 | 35
|
---|
1245 | t
|
---|
1246 | grid_simple 218.522 <br>
|
---|
1247 |
|
---|
1248 | 35
|
---|
1249 | q
|
---|
1250 | grid_simple 2.83583e-06 <br>
|
---|
1251 |
|
---|
1252 | 34
|
---|
1253 | t
|
---|
1254 | grid_simple 217.725 <br>
|
---|
1255 |
|
---|
1256 | 34
|
---|
1257 | q
|
---|
1258 | grid_simple 2.84932e-06 <br>
|
---|
1259 |
|
---|
1260 | 33
|
---|
1261 | t
|
---|
1262 | grid_simple 217.795 <br>
|
---|
1263 |
|
---|
1264 | 33
|
---|
1265 | q
|
---|
1266 | grid_simple 2.86236e-06 <br>
|
---|
1267 |
|
---|
1268 | 32
|
---|
1269 | t
|
---|
1270 | grid_simple 217.858 <br>
|
---|
1271 |
|
---|
1272 | 32
|
---|
1273 | q
|
---|
1274 | grid_simple 2.88346e-06 <br>
|
---|
1275 |
|
---|
1276 | 31
|
---|
1277 | t
|
---|
1278 | grid_simple 218.372 <br>
|
---|
1279 |
|
---|
1280 | 31
|
---|
1281 | q
|
---|
1282 | grid_simple 2.90988e-06 <br>
|
---|
1283 |
|
---|
1284 | 30
|
---|
1285 | t
|
---|
1286 | grid_simple 219.283 <br>
|
---|
1287 |
|
---|
1288 | 30
|
---|
1289 | q
|
---|
1290 | grid_simple 2.92827e-06 <br>
|
---|
1291 |
|
---|
1292 | 29
|
---|
1293 | t
|
---|
1294 | grid_simple 220.089 <br>
|
---|
1295 |
|
---|
1296 | 29
|
---|
1297 | q
|
---|
1298 | grid_simple 2.93849e-06 <br>
|
---|
1299 |
|
---|
1300 | 28
|
---|
1301 | t
|
---|
1302 | grid_simple 219.708 <br>
|
---|
1303 |
|
---|
1304 | 28
|
---|
1305 | q
|
---|
1306 | grid_simple 2.94274e-06 <br>
|
---|
1307 |
|
---|
1308 | 27
|
---|
1309 | t
|
---|
1310 | grid_simple 218.645 <br>
|
---|
1311 |
|
---|
1312 | 27
|
---|
1313 | q
|
---|
1314 | grid_simple 2.94806e-06 <br>
|
---|
1315 |
|
---|
1316 | 26
|
---|
1317 | t
|
---|
1318 | grid_simple 218.014 <br>
|
---|
1319 |
|
---|
1320 | 26
|
---|
1321 | q
|
---|
1322 | grid_simple 2.96126e-06 <br>
|
---|
1323 |
|
---|
1324 | 25
|
---|
1325 | t
|
---|
1326 | grid_simple 218.125 <br>
|
---|
1327 |
|
---|
1328 | 25
|
---|
1329 | q
|
---|
1330 | grid_simple 2.98719e-06 <br>
|
---|
1331 |
|
---|
1332 | 24
|
---|
1333 | t
|
---|
1334 | grid_simple 218.94 <br>
|
---|
1335 |
|
---|
1336 | 24
|
---|
1337 | q
|
---|
1338 | grid_simple 3.02117e-06 <br>
|
---|
1339 |
|
---|
1340 | 23
|
---|
1341 | t
|
---|
1342 | grid_simple 220.96 <br>
|
---|
1343 |
|
---|
1344 | 23
|
---|
1345 | q
|
---|
1346 | grid_simple 3.08349e-06 <br>
|
---|
1347 |
|
---|
1348 | 22
|
---|
1349 | t
|
---|
1350 | grid_simple 223.01 <br>
|
---|
1351 |
|
---|
1352 | 22
|
---|
1353 | q
|
---|
1354 | grid_simple 3.18947e-06 <br>
|
---|
1355 |
|
---|
1356 | 21
|
---|
1357 | t
|
---|
1358 | grid_simple 223.778 <br>
|
---|
1359 |
|
---|
1360 | 21
|
---|
1361 | q
|
---|
1362 | grid_simple 3.26119e-06 <br>
|
---|
1363 |
|
---|
1364 | 20
|
---|
1365 | t
|
---|
1366 | grid_simple 224.728 <br>
|
---|
1367 |
|
---|
1368 | 20
|
---|
1369 | q
|
---|
1370 | grid_simple 3.25947e-06 <br>
|
---|
1371 |
|
---|
1372 | 19
|
---|
1373 | t
|
---|
1374 | grid_simple 226.885 <br>
|
---|
1375 |
|
---|
1376 | 19
|
---|
1377 | q
|
---|
1378 | grid_simple 3.26296e-06 <br>
|
---|
1379 |
|
---|
1380 | 18
|
---|
1381 | t
|
---|
1382 | grid_simple 231.473 <br>
|
---|
1383 |
|
---|
1384 | 18
|
---|
1385 | q
|
---|
1386 | grid_simple 3.25986e-06 <br>
|
---|
1387 |
|
---|
1388 | 17
|
---|
1389 | t
|
---|
1390 | grid_simple 236.668 <br>
|
---|
1391 |
|
---|
1392 | 17
|
---|
1393 | q
|
---|
1394 | grid_simple 3.3009e-06 <br>
|
---|
1395 |
|
---|
1396 | 16
|
---|
1397 | t
|
---|
1398 | grid_simple 238.828 <br>
|
---|
1399 |
|
---|
1400 | 16
|
---|
1401 | q
|
---|
1402 | grid_simple 3.37098e-06 <br>
|
---|
1403 |
|
---|
1404 | 15
|
---|
1405 | t
|
---|
1406 | grid_simple 244.755 <br>
|
---|
1407 |
|
---|
1408 | 15
|
---|
1409 | q
|
---|
1410 | grid_simple 3.45707e-06 <br>
|
---|
1411 |
|
---|
1412 | 14
|
---|
1413 | t
|
---|
1414 | grid_simple 249.2 <br>
|
---|
1415 |
|
---|
1416 | 14
|
---|
1417 | q
|
---|
1418 | grid_simple 3.53864e-06 <br>
|
---|
1419 |
|
---|
1420 | 13
|
---|
1421 | t
|
---|
1422 | grid_simple 256.419 <br>
|
---|
1423 |
|
---|
1424 | 13
|
---|
1425 | q
|
---|
1426 | grid_simple 3.6194e-06 <br>
|
---|
1427 |
|
---|
1428 | 12
|
---|
1429 | t
|
---|
1430 | grid_simple 263.959 <br>
|
---|
1431 |
|
---|
1432 | 12
|
---|
1433 | q
|
---|
1434 | grid_simple 3.73099e-06 <br>
|
---|
1435 |
|
---|
1436 | 11
|
---|
1437 | t
|
---|
1438 | grid_simple 268.754 <br>
|
---|
1439 |
|
---|
1440 | 11
|
---|
1441 | q
|
---|
1442 | grid_simple 3.82798e-06 <br>
|
---|
1443 |
|
---|
1444 | 10
|
---|
1445 | t
|
---|
1446 | grid_simple 271.708 <br>
|
---|
1447 |
|
---|
1448 | 10
|
---|
1449 | q
|
---|
1450 | grid_simple 3.9072e-06 <br>
|
---|
1451 |
|
---|
1452 | 9
|
---|
1453 | t
|
---|
1454 | grid_simple 272.846 <br>
|
---|
1455 |
|
---|
1456 | 9
|
---|
1457 | q
|
---|
1458 | grid_simple 3.95614e-06 <br>
|
---|
1459 |
|
---|
1460 | 8
|
---|
1461 | t
|
---|
1462 | grid_simple 270.363 <br>
|
---|
1463 |
|
---|
1464 | 8
|
---|
1465 | q
|
---|
1466 | grid_simple 3.97863e-06 <br>
|
---|
1467 |
|
---|
1468 | 7
|
---|
1469 | t
|
---|
1470 | grid_simple 265.951 <br>
|
---|
1471 |
|
---|
1472 | 7
|
---|
1473 | q
|
---|
1474 | grid_simple 3.99851e-06 <br>
|
---|
1475 |
|
---|
1476 | 6
|
---|
1477 | t
|
---|
1478 | grid_simple 260.086 <br>
|
---|
1479 |
|
---|
1480 | 6
|
---|
1481 | q
|
---|
1482 | grid_simple 4.01725e-06 <br>
|
---|
1483 |
|
---|
1484 | 5
|
---|
1485 | t
|
---|
1486 | grid_simple 251.757 <br>
|
---|
1487 |
|
---|
1488 | 5
|
---|
1489 | q
|
---|
1490 | grid_simple 4.01412e-06 <br>
|
---|
1491 |
|
---|
1492 | 4
|
---|
1493 | t
|
---|
1494 | grid_simple 237.045 <br>
|
---|
1495 |
|
---|
1496 | 4
|
---|
1497 | q
|
---|
1498 | grid_simple 4.05898e-06 <br>
|
---|
1499 |
|
---|
1500 | 3
|
---|
1501 | t
|
---|
1502 | grid_simple 217.361 <br>
|
---|
1503 |
|
---|
1504 | 3
|
---|
1505 | q
|
---|
1506 | grid_simple 3.39441e-06 <br>
|
---|
1507 |
|
---|
1508 | 2
|
---|
1509 | t
|
---|
1510 | grid_simple 201.113 <br>
|
---|
1511 |
|
---|
1512 | 2
|
---|
1513 | q
|
---|
1514 | grid_simple 2.40852e-06 <br>
|
---|
1515 |
|
---|
1516 | 1
|
---|
1517 | z
|
---|
1518 | grid_simple -14.0557 = <span style="color: rgb(51, 51, 255);">9.80665 * -1.4332825 </span><br>
|
---|
1519 |
|
---|
1520 | 1
|
---|
1521 | t
|
---|
1522 | grid_simple 187.843 <br>
|
---|
1523 |
|
---|
1524 | 1
|
---|
1525 | q
|
---|
1526 | grid_simple 1.56344e-06 <br>
|
---|
1527 |
|
---|
1528 | 1
|
---|
1529 | lnsp grid_simple
|
---|
1530 | 11.52 <span style="color: rgb(51, 51, 255);">= ln(100710.)</span><br>
|
---|
1531 | edition
|
---|
1532 | centre
|
---|
1533 | date
|
---|
1534 | dataType gridType
|
---|
1535 | stepRange typeOfLevel
|
---|
1536 | level
|
---|
1537 | shortName packingType value <br>
|
---|
1538 | Other grid Points<br>
|
---|
1539 | Input Point: latitude=50.00 longitude=-3.00<br>
|
---|
1540 | Grid Point chosen #2 index=14577 latitude=50.00 longitude=-3.00 distance=0.00 (Km)<br>
|
---|
1541 | /data/nwp1/idculv/ROPP/ropp_src/branches/dev/Share/ic_grib2ropp/grib_files/grib2/2011050800/000<br>
|
---|
1542 | - 4 - index=14217 latitude=51.00 longitude=-3.00 distance=111.20 (Km)<br>
|
---|
1543 | - 3 - index=14218 latitude=51.00 longitude=-2.00 distance=131.79 (Km)<br>
|
---|
1544 | - 2 - index=14577 latitude=50.00 longitude=-3.00 distance=0.00 (Km)<br>
|
---|
1545 | - 1 - index=14578 latitude=50.00 longitude=-2.00 distance=71.48 (Km)<br>
|
---|
1546 | <br>
|
---|
1547 | <br>
|
---|
1548 | <br>
|
---|
1549 | <span style="font-weight: bold; color: rgb(51, 51, 255);">Aspects relevant to ROPP implementation:</span><br>
|
---|
1550 | <br>
|
---|
1551 | 1) Need to be able to read from general ECMWF grib2 (and grib1?) dump file.<br>
|
---|
1552 | 2) Need to include GRIB API as an external library, like netCDF or
|
---|
1553 | BUFR. Copyright? Currently I use grib_api_1.9.9, as 1.9.16 has
|
---|
1554 | problems with DMI's example file. Under discussion with ECMWF.<br>
|
---|
1555 | 3) Need to clean up the reading/checking (don't generate lat and lon
|
---|
1556 | for every call: bracket this with an (IF Nmess ==1) control statement.
|
---|
1557 | What's the need for call_read_ECMWF_orography? What's the
|
---|
1558 | need for the Geoid Height and Curvature business? Forget about them
|
---|
1559 | intially.<br>
|
---|
1560 | 4) Try to write with the flexibility to read two grib files and
|
---|
1561 | interpolate (linearly) in time between them (eventually). Ensure
|
---|
1562 | we can handle GRIB1 and GRIB2 (and GRIB2_double?). Keep ability
|
---|
1563 | to produce multiple bgr files per call. I suggest we don't bother
|
---|
1564 | with Hallgeir's suggestion to allow forward modelling as part of the
|
---|
1565 | script. <br>
|
---|
1566 | 5) Natural place to put it: ropp_io/tools/grib2text.<br>
|
---|
1567 | <br>
|
---|
1568 | <br>
|
---|
1569 | <br>
|
---|
1570 |
|
---|
1571 |
|
---|
1572 | <br>
|
---|
1573 |
|
---|
1574 |
|
---|
1575 | <hr style="width: 100%; height: 2px;"><span style="font-style: italic;">Last updated </span><span style="font-weight: bold; font-style: italic;">24th
|
---|
1576 | September 2012</span><span style="font-style: italic;">
|
---|
1577 | by </span><a style="font-style: italic;" href="mailto:ian.culverwell@metoffice.gov.uk">ian.culverwell@metoffice.gov.uk</a><br>
|
---|
1578 |
|
---|
1579 |
|
---|
1580 | <br>
|
---|
1581 |
|
---|
1582 |
|
---|
1583 | </body>
|
---|
1584 | </html>
|
---|