Grads在一个2D(高度 - 时间)图中绘制来自不同文件(不同时间点)的数据

时间:2014-10-18 19:15:20

标签: plot noaa grib gfs grads

我从美国gfs天气模式下载数据(来自ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2014101812/master,如果您阅读超过1个月后编辑日期)。那些是夹子文件。

我正在使用梯度来绘制数据,如果我在地图上绘制1个时间步长(= 1个文件)的数据,它可以正常工作。 所以我的第一次尝试是打开前两个文件并制作一个阴影图,显示相同风速的相同颜色的区域。我得到的是一个正确轴的空图。

'open gfs.t12z.mastergrb2f00.ctl'

'open gfs.t12z.mastergrb2f03.ctl'

'set dfile 1'

'set time 12Z17OCT2014'
'set gxout shaded'
'set lat 50'
'set lon 0'
'set lev 1000 200'
'set t 1 2'
'd sqrt(ugrdprs*ugrdprs+vgrdprs*vgrdprs)*1.9438'

'set dfile 2'
'set time 15Z17OCT2014'
'set t 1 2'
'd sqrt(ugrdprs*ugrdprs+vgrdprs*vgrdprs)*1.9438'

我正在使用opengrads - 网格分析和显示系统(GrADS)版本2.0.2.oga.2。

任何可以帮助我将不同时间步(=不同文件)的数据绘制成一个图的人?

1 个答案:

答案 0 :(得分:1)

您可以使用一个.ctl文件来绘制多个数据文件。 我用.gdat文件完成了这个,所以我猜你可以用GRIB做到这一点。 您应该使用grads模板创建新的.ctl文件。

http://www.iges.org/grads/gadoc/templates.html

基本上,在dset行中,您将放置一个类似于

的文件名
dset gfs.t12z.mastergrb2f%h2.grb

其中%h2表示2位数小时。