在一个二维(高度-时间)图中渐变来自不同文件(不同时间点)的打印数据



我从美国gfs天气模型下载了数据ftp://ftp.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2014101812/master,如果阅读时间超过1个月,请编辑日期)。那些是grib文件。

我使用梯度来绘制数据,如果我在地图上绘制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'

我使用的是opengrades-网格分析和显示系统(GrADS)2.0.2.oga.2版本。

有谁能帮助我将不同时间步长(=不同文件)的数据绘制到一个图中?

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

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

基本上,在数据集中,你会放一个类似的文件名

dset gfs.t12z.mastergrb2f%h2.grb

其中%h2表示2位小时。

最新更新