求助 画剖面图

sharezkg

我自己用fortran写了一个三维的二进制数据
do k=1,15
write(200+wz,rec=k)((month(i,j,k),i=1,192),j=1,105)
enddo   
ctl是自己的写的  不知道有问题没有
DSET   NOX_3d.asc
TITLE  NOX data
UNDEF  -9999
XDEF   192  linear 101.7335 0.0902089
YDEF   105  linear 16.53754 0.08415971
ZDEF   15 LEVELS
1000.00
994.00
987.00
974.00
957.00
940.00
920.00
900.00
877.00
850.00
824.00
795.00
765.00
728.00
690.00
TDEF   1 LINEAR 00z4Jan2005 1hr
VARS   1
NOX 15 9999
ENDVARS
然后想把这个的剖面图画出来
'reinit'
  'open D:\work\biyelunwen\simulationresult\3D\01base\NOX.ctl'
'enable print NOX.gmf'
'set grads off'
'set grid off'
'set x 1'
'set y 1'
'set z 1 15'
lon1=109.758347
lon2=113.75056
lat1=19.494931
lat2=23.742993
lon=lon1
'collect 1 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(NOX,'lon','lat')'
lon=lon+0.1
endwhile

'set lon 'lon1' 'lon2'
'set z 1 15'
'set perea off'
'set grads off'
'set gxout shaded'
'd coll2gr(1)
'cbarn 1 0 6. 0.45'
'print'
'disable print' 层数问题解决了  犯了个低级错误，但是剖面图还没有，我再研究研究
C:\Users\zkg\Desktop\1234.png

sharezkg

图片开始没传上去，画剖面图的时候  就提示变量没有Z维

终古未央

太棒了

18393810639

我和你的问题一样呢。楼主解决了吗？