GrADS出图

霞丫呀

求助各位大佬，出1000hap天气图时用到的是hgt.1981.nc和air.1981.nc,然后出完后发现显示的数据不应该是1000hpade天气图的数值，想麻烦各位大佬帮忙看看是数据用的不对，还是代码不对。
代码如下：
'c'
'reinit'
'sdfopen E:\hgt.1984.nc'
'sdfopen E:\air.1984.nc'
'set mpdset cnworld'
'set map 15 1 6'
'set grads off'

'set lon 60 200'
'set lat 20 90'
day=275
while(day<306)
'set t 'day''

'set cthick 4'
'set cint 2'
'set ccolor 2'
'd air.2(lev=500)-273.16'

'set cint 2.5'
'set ccolor 4'

if(day-274>=10)
'draw string 4.55 8.1 500hPa 1984 10 'day-274''
'd hgt(lev=500)/10'
'printim E:\198410'day-274'500.png white x1600 y1600'
else
'draw string 4.55 8.1 500hPa 1984 10 'day-274''
'd hgt(lev=500)/10'
'printim E:\1984100'day-274'500.png white x1600 y1600'
endif

'c'
'set cthick 4'
'set cint 2'
'set ccolor 2'
'd air.2(lev=700)-273.16'

'set cint 2.5'
'set ccolor 4'

if(day-274>=10)
'draw string 4.55 8.1 700hPa 1984 10 'day-274''
'd hgt(lev=700)/10'
'printim E:\198410'day-274'700.png white x1600 y1600'
else
'draw string 4.55 8.1 700hPa 1984 10 'day-274''
'd hgt(lev=700)/10'
'printim E:\1984100'day-274'700.png white x1600 y1600'
endif

a1s2d3

没太看懂问题，但根据程序出的图应该是500hPa和700hPa的温度场和高度场，如果题主是想出1000hPa的温度场和高度场的话应该用这两个语句
'd air.2(lev=1000)-273.16'
'd hgt(lev=1000)/10'

霞丫呀

a1s2d3 发表于 2023-1-7 21:34
没太看懂问题，但根据程序出的图应该是500hPa和700hPa的温度场和高度场，如果题主是想出1000hPa的温度场和 ...

谢谢啦