Section | |||||
---|---|---|---|---|---|
|
Example Metview Macro
The following macro retrieves several time steps of data from MARS, plots them to a PostScript file and uses the convert
command to generate an animated gif. The result is shown on this page.
Code Block | ||
---|---|---|
| ||
# Metview Macro # **************************** LICENSE START *********************************** # # Copyright 2014- ECMWF. This software is distributed under the terms # of the Apache License version 2.0. In applying this license, ECMWF does not # waive the privileges and immunities granted to it by virtue of its status as # an Intergovernmental Organization or submit itself to any jurisdiction. # # ***************************** LICENSE END ************************************ # retrieve some data from MARS t2m_fc = retrieve( type : "fc", levtype : "sfc", param : "2t", time : 00, step : [00,"to",72,"by",6], grid : [1.5,1.5] ) # define our plotting attributes t_shade = mcont( legend : "on", contour_automatic_setting : "ecchart" ) view = geoview( map_area_definition : "corners", area : [30.62,-25.4,70.12,40.36] ) # plot the data outdir = getenv('SCRATCH') outbasename = "t2m_fc" ps = ps_output(output_name : outdir & "/" & outbasename) setoutput(ps) plot(view, t2m_fc, t_shade) # force Macro to wait for the plot to be generated, then convert # to animated gif setoutput(ps) # setoutput()wait forcesfor the plot file to be finishgenerated shell('convert -delay 100 -rotate "90<" ' & outdir & '/' & outbasename & '.ps' & ' -delay 100 -rotate "90<" ' & outdir & '/' & outbasename & '.gif') |
...