...
Code Block | ||||||
---|---|---|---|---|---|---|
| ||||||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() server.retrieve({ "class": "ei", "dataset": "interim", "date": "2015-08-01", "expver": "1", "grid": "0.75/0.75", "levtype": "ml", "levelist": "all", "param": "t/q", "step": "0", "stream": "oper", "time": "0", "type": "an", "format": "netcdf" "target": "tq_ml_20150801_00.nc", }) |
Code Block | ||||
---|---|---|---|---|
| ||||
#!/usr/bin/env python from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() server.retrieve({ |
...
"class": "ei", |
...
"dataset": "interim", |
...
|
...
"date": "2016-12-01", |
...
"expver": "1", |
...
"grid": "0.75/0.75", |
...
|
...
"levelist": "1", |
...
"levtype": "ml", |
...
"param": "lnsp", |
...
|
...
"step": "0", |
...
|
...
"stream": "oper", |
...
"time": "0", |
...
|
...
"type": "an", |
...
"format": "netcdf" |
...
...
"target": "lnsp_ml_20150801_00.nc", |
...
}) |
Code Block | ||||||
---|---|---|---|---|---|---|
| ||||||
#!/usr/bin/env python
|
...
from ecmwfapi import ECMWFDataServer server = ECMWFDataServer() server.retrieve({ |
...
"class": "ei", |
...
|
...
"dataset": "interim", |
...
|
...
"date": "1989-01-01", |
...
"expver": "1", |
...
|
...
"grid": "0.75/0.75", |
...
"levtype": "sfc", |
...
"param": "z", |
...
|
...
"step": "0", |
...
|
...
"stream": "oper", |
...
"time": "12:00:00", |
...
|
...
"type": "an", |
...
|
...
"format": "netcdf" |
...
...
"target": "zlnsp_ml_20150801_00.nc", }) |
Compute geopotential on model levels
- The script is based on the Metview macro mvl_geopotential_on_ml, code from Nils Wedi, IFS documentation:
- CY41R1 Official IFS Documentation
- part III. Dynamics and numerical procedures
- CY41R1 Official IFS Documentation
Download
Python program | compute_geopotential_on_ml.py |
Wrapper script | compute_geopotential_on_ml.ksh |
Examples
This example will compute the geopotential on the 2015-10-08 time 00 operational analysis model levels (137). Below you can see the MARS user documentation request used to retrive both files. You can set a different class/stream/type for the input data. The gribType and resolution can also be changed.
python compute_geopotential_on_ml.py tq_ml_20151008_00.grib zlnsp_ml_20151008_00.grib python compute_geopotential_on_ml.py tq_ml_20151008_00.grib zlnsp_ml_20151008_00.grib -o my_grib.grib |
tq_ml_20151008_00.grib
zlnsp_ml_20151008_00.grib