如何将几个不同层次数目（或者不同时次数目）的nc文件的物理量场绘制在一个图里面？

1434272731

如何将几个不同层次数目（或者不同时次数目）的nc文件的物理量场绘制在一个图里面？

需要知道的是，如果几个nc文件的层次数目（或者时次）不同，那么想直接在Grads里面绘制叠加的环流场是不可行的，绘制第二个物理场的时候会报这样的错：

                              

DataRequest Warning:  Request is completelyoutside file limits. Cannot contour grid - all undefined values

在绘制第一个物理量之前，需要设置该物理量所在文件的维度（经纬度、时次、层次），使用'setdfile ....'

例子：现有NCEP资料三种：风场uwnd.nc、vwnd.nc；海平面气压场：ncep2.nc

首先查看两种.netcdf资料的具体信息：

风场信息：这里有1层，时次807个

  

   

海平面气压场信息：这里有17层，时次816个

  

可以看到，风场和海平面气压场文件的时间维总数和层次维总数是不同的，欲把风场和高度场绘制在一张图里面，则应该这样写gs文件：

'reinit'
'sdfopen  ………………uwnd.mon.mean.nc'
'sdfopen  ………………vwnd.mon.mean.nc'
'sdfopen  ………………ncep1.nc'

'set grads off'
'set grid off'
'set mpdset cnworld'
'set dfile 1'
'set dfile 2'
'set lat 20 70'
'set lon 60 140'
'set z 1'
'set t 1'
'set gxout barb'
'd uwnd.1;vwnd.2'
'set dfile 3'
'set lat 20 70'
'set lon 60 140'
'set z 1'
'set t 1'
'set gxout contour'
'set ccolor 1'
'd pres.3/100.'
；

结果贴一个图：

（说明：看不明白的话，可以下载附件中完整地看到图片，能看明白的就不要费金钱啦！）