垂直螺旋度 gs文件 请问这个到底哪里有问题啊，怎么就是出不来图呢？

流水2012

'reinit'
'clear'
'sdfopen e:\uwnd.2013.nc'
'sdfopen e:\vwnd.2013.nc'
'sdfopen e:\omega.2013.nc'
'set mpdset mres'
'set lon 80 110'
'set lat 30 50'
'set lev 700 200'
'set time 08z08may'
'define wd=hcurl(uwnd.1,vwnd.2)'
'define lx=-(wd*omega.3*10.)/12.64'

'set xlevs 100 101 102 103 104 105'
'set lat 40'
'enable print e:\lxd-4.gmf'
'd lx*1e+6'
'draw title lxd-08-08'
'print'
'disable print'
问题补充：就是我想画高低层的垂直螺旋度，然后，这个运行就是空白图，而且总是说数据要求超出文件限制！！

传说中的谁

你这就算的是单层的物理量，高度设置成多层的了。
PS：提问的时候附上出错信息吧

river

错误信息呢？是要画剖面吗？楼主把自己的问题描述详细一点

nunu18

楼主是画一层的还是某个剖面上的垂直螺旋度？你这个程序自己改一下，若是某一层，z就设成一层，若是剖面，你就该lon~~~

Eegle

这个是我自己画的螺旋度，参考了论坛里的方法，但问题就是为什么输出的图有问题呢？lxd.gmf无法查看，只能在openGrads里面查看？？
'reinit'
'sdfopen E:\zhuang\uwnd.2014.nc'
'sdfopen E:\zhuang\vwnd.2014.nc'
'sdfopen E:\zhuang\omega.2014.nc'
'set grid off'
'set poli on'
'set mpdraw on'
'set mproj scaled'
'set mpdset cnworld'
'set xlab on'
'set ylab on'
'set xlopts 1 4 0.16'
'set ylopts 1 4 0.16'
'set lat 10 35'
'set lon 110 130'
'set lev 700'
'set t 658'
'define vor=hcurl(uwnd.1,vwnd.2)'
'define lx=-(vor*omega.3*10.)/12.64'
'd lx*1e+6'
'set gxout shaded'
'enable print E:\zhuang\picture\lxd.gmf'
'draw title lxd-06-14'
'print'
'disable print'

老A

你前面的那个GS有问题，看看有没有'set time 08z08may'这个时次的资料

subtropical

Eegle 发表于 2014-9-16 10:57
这个是我自己画的螺旋度，参考了论坛里的方法，但问题就是为什么输出的图有问题呢？lxd.gmf无法查看，只能 ...

您好，您的垂直螺旋度感觉很小啊

chuanhao

-(vor*omega.3*10.)/12.64  这里面的负号和除以12.64分别是啥意思啊？