高度-经度剖面图问题

坏人！

求大家帮我看看，为啥我'set lev 10 1000‘能画出来，'set lev 1000 10’grads就会崩溃呀{:eb303:}

gs:
'reinit'
'clear'
'sdfopen C:\Users\lalala\Desktop\nwp_homework1\data_read\NetCDF\mLat.198101-201012.clt.nc'
latn=45
pi=3.1415926
r=6371000
omega=2*pi/24/3600
g=9.8
f=2*omega*math_sin(latn/180*pi)
'set lon -180 180'
'set lev 1000 10'
'set t 1'
'define dx=cdiff(lon,x)*'pi'/180*'r'*cos('latn'*'pi'/180)'
'define dh=cdiff(hgt,x)'
'define pressure=-'g'*dh/dx*1e5'
'define coriolis='f'*vwnd*1e5'
'set gxout shaded'
'set timelab off'
'd coriolis'
'cbarn 0.5 1'
'draw title height-longtitude section of coriolis term'
'set strsiz 0.15'
'draw string 3.3 0.2 latitude:'latn'N'
'set strsiz 0.12'
'draw string 8 0.2 e-5 ms'
'set strsiz 0.08'
'draw string 8.8 0.3 -2'
;

报错（grads直接崩溃了。。。）：
Define Error:  Unable to allocate data memory
  Size of request was -2175 grid elements
1855129216 [sig] opengrads 8612 open_stackdumpfile: Dumping stack trace to opengrads.exe.stackdump

river

你display之前设置里都没固定纬度啊，先把纬度固定了看看

坏人！

river 发表于 2016-3-26 10:57
你display之前设置里都没固定纬度啊，先把纬度固定了看看

忘说了。。这个纬度是固定的。。。只有一维的数据。。。我错了。。

补充一下q file
File 1 :
  Descriptor: C:\Users\lalala\Desktop\nwp_homework1\data_read\NetCDF\mLat.198101-201012.clt.nc
  Binary: C:\Users\lalala\Desktop\nwp_homework1\data_read\NetCDF\mLat.198101-201012.clt.nc
  Type = Gridded
  Xsize = 144  Ysize = 1  Zsize = 17  Tsize = 12  Esize = 1
  Number of Variables = 4
     uwnd  17  t,z,y,x  U wind
     vwnd  17  t,z,y,x  V wind
     hgt  17  t,z,y,x  geopotential height
     air  17  t,z,y,x  air temperature

river

q ctlinfo 看看

坏人！

river 发表于 2016-3-26 11:24
q ctlinfo 看看

q ctlinfo
dset C:\Users\lalala\Desktop\nwp_homework1\data_read\NetCDF\mLat.198101-201012.clt.nc
title
undef -32767
dtype netcdf
xdef 144 linear 0 2.5
ydef 1 linear 0 1
zdef 17 levels 10 20 30 50 70 100 150 200
250 300 400 500 600 700 850 925 1000
tdef 12 linear 00Z01JAN1981 1mo
vars 4
uwnd=>uwnd  17  t,z,y,x  U wind
vwnd=>vwnd  17  t,z,y,x  V wind
hgt=>hgt  17  t,z,y,x  geopotential height
air=>air  17  t,z,y,x  air temperature
endvars

坏人！

river 发表于 2016-3-26 11:24
q ctlinfo 看看

关键是我倒过来‘set lev 10 1000'就画得出来。。。

康理想

Grads的确存在这样的问题，建议你把数据提取出来，按照自己想要的顺序存储，重新编写ctl绘图

river

坏人！ 发表于 2016-3-26 11:54
关键是我倒过来‘set lev 10 1000'就画得出来。。。

这个，你自己看啊，ctl上面就已经告诉你这个资料是从高到低排列的啊。那你自己加一个Y轴反向命令不就好了吗

坏人！

river 发表于 2016-3-26 21:50
这个，你自己看啊，ctl上面就已经告诉你这个资料是从高到低排列的啊。那你自己加一个Y轴反向命令不就好了 ...

受教了~！谢谢~~~

坏人！

康理想 发表于 2016-3-26 15:30
Grads的确存在这样的问题，建议你把数据提取出来，按照自己想要的顺序存储，重新编写ctl绘图

按楼下的方法比较简单呢~谢谢~~