画纬度-时间剖面度程序问题(更新计算的GS）

影形

各位大大好，初学grads，很多东西还在模仿阶段，自己找不到错误，请多指正：
我的ctl文件：
dset E:\China_dqv.grd
title WCD
undef -9.99E33
xdef 1 linear 110.0 2.5
ydef 18 linear 15.0 2.5
zdef 1 levels 700
tdef 64 linear JAN1960 1mo
vars 1
dqv 0 99 div
endvars
画图的gs：
'reinit'
'open E:\div.ctl'
'enable print E:\div.gmf'
'set xlopts 1 5 0.16'
'set ylopts 1 5 0.16'
'c'
'set grads off'
'set grid off'
'set t 1 648'
'set cthick 6'
'set xyrev on'
'd dqv'
'print'
'disable print'
;
现在的问题是run gs的时候提示：
The "enable print "command is no longer valid,use "gxprint" instead
Non-terminated constant
然而我改成gxprint以后却变成了：
Unable to determine output format
求帮助。。。是gs出问题还是我前期grd数据有问题

计算GRD的gs如下

'reinit'

'sdfopen E:\NCEP\pres.mon.mean.nc'
'sdfopen E:\NCEP\hgt.mon.mean.nc'
'sdfopen E:\NCEP\uwnd.mon.mean.nc'
'sdfopen E:\NCEP\vwnd.mon.mean.nc'
'sdfopen E:\NCEP\shum.mon.mean.nc'

'set mpdset cnworld'

'set gxout fwrite'
      
'set fwrite E:\F\essay\data\data\NCEP\div\China_dqv.grd'

year = 1960
while(year<=2013)
    mon = 1
    while(mon<=12)
        it=(year-1960)*12+mon+144
'set dfile 1'
'set x 29 57'
'set y 43 60'
'set t 'it''
'define pres1=pres'

'set dfile 2'
'set x 29 57'shum
'set y 43 60'
'set z 4'
'set t 'it''
'define hgt1=hgt'

'set dfile 3'
'set x 29 57'
'set y 43 60'
'set z 1 8'
'set t 'it''
'define U=uwnd'

'set dfile 4'
'set x 29 57'
'set y 43 60'
'set z 1 8'
'set t 'it''
'define V=vwnd'

'set dfile 5'
'set x 29 57'
'set y 43 60'
'set z 1 8'
'set t 'it''
'define q=shum'

'define qu=U*q/1000'
'define qv=V*q/1000'
'define uqs=vint(pres1,qu,850)'
'define vqs=vint(pres1,qv,850)'
'define qs=vint(pres1,q,850)/1000'

'define dqv = hdivg(uqs,vqs)'
'set z 1'

'set lon 110'

'd dqv'
mon=mon+1
endwhile
year=year+1
endwhile
'disable fwrite'
这个grd用MATLAB打开看只有17以前有数值，以后都是-9999

爱糖的孩子lhk

把中间的c去掉然后再把set x set y都加上呢？

river

'reinit'
'open E:\div.ctl'
'gxprint E:\div.png white'
'set xlopts 1 5 0.16'
'set ylopts 1 5 0.16'
'set grads off'
'set grid off'
'set t 1 648'
'set cthick 6'
'set xyrev on'
'd dqv'
'gxprint E:\div.png white'
;

影形

river 发表于 2015-8-8 09:54
'reinit'
'open E:\div.ctl'
'gxprint E:\div.png white'

谢谢您~
程序没有再报错了，但是图也仅显示了坐标，没有内容。。。

river

影形 发表于 2015-8-8 12:08
谢谢您~
程序没有再报错了，但是图也仅显示了坐标，没有内容。。。

那就检查你的grd资料吧

影形

river 发表于 2015-8-8 17:42
那就检查你的grd资料吧

'reinit'

'sdfopen E:\NCEP\pres.mon.mean.nc'
'sdfopen E:\NCEP\hgt.mon.mean.nc'
'sdfopen E:\NCEP\uwnd.mon.mean.nc'
'sdfopen E:\NCEP\vwnd.mon.mean.nc'
'sdfopen E:\NCEP\shum.mon.mean.nc'

'set mpdset cnworld'

'set gxout fwrite'
      
'set fwrite E:\F\essay\data\data\NCEP\div\China_dqv.grd'

year = 1960
while(year<=2013)
    mon = 1
    while(mon<=12)
        it=(year-1960)*12+mon+144
'set dfile 1'
'set x 29 57'
'set y 43 60'
'set t 'it''
'define pres1=pres'

'set dfile 2'
'set x 29 57'shum
'set y 43 60'
'set z 4'
'set t 'it''
'define hgt1=hgt'

'set dfile 3'
'set x 29 57'
'set y 43 60'
'set z 1 8'
'set t 'it''
'define U=uwnd'

'set dfile 4'
'set x 29 57'
'set y 43 60'
'set z 1 8'
'set t 'it''
'define V=vwnd'

'set dfile 5'
'set x 29 57'
'set y 43 60'
'set z 1 8'
'set t 'it''
'define q=shum'

'define qu=U*q/1000'
'define qv=V*q/1000'
'define uqs=vint(pres1,qu,850)'
'define vqs=vint(pres1,qv,850)'
'define qs=vint(pres1,q,850)/1000'

'define dqv = hdivg(uqs,vqs)'
'set z 1'

'set lon 110'

'd dqv'
mon=mon+1
endwhile
year=year+1
endwhile
'disable fwrite'
您好~这是我的gs。。。麻烦您再帮我看一眼，非常感谢！

river

影形 发表于 2015-8-9 08:55
'reinit'

'sdfopen E:\NCEP\pres.mon.mean.nc'

你这个没办法改，明明可以一步到位的，没必要这么多步骤。建议你看一下论坛里关于剖面图作法的帖子