| 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>
|
|---|