【精细版】用GrADS转换nc数据

cxyykl

刚好要用到，已存，谢

饭团TT

nc转grd一直不成功，没解决，继续学

lehuoxiaoyu

学习学习{:eb502:}{:eb502:}

xmgcqlmw

楼主你好，跟你的方法一样，想输出12层的资料，但是只要第一层是对的，后面的都是错的，想了很久实在不知道哪里出了问题，我取的第一个时次，程序如下：
'reinit'      
    'sdfopen D:\sy\1980EC\nc\1000-100hpa\1980.6.nc'
    'set fwrite D:\sy\1980EC\1000-100hpa\1980.6grd\1.grd'
    'set gxout fwrite'
     'set t 1'
    i=1
     while(i<=12)
     'set z 'i''
     'set lat 10 60'
     'set lon 80 140'
     'd t'
     'd w'
     'd vo'
     'd d'
     'd r'
     i=i+1
    endwhile
'disable fwrite'
    'reinit'
   ;

对应的ctl：
dset D:\sy\1980EC\1000-100hpa\1980.6grd\1.grd
title ec
undef -9.99E+33
xdef   61  linear    80.0000    1
ydef   51  linear    10.0000    1
zdef   12    levels    1000 925 850 700 600 500 400 300 250 200 150 100
tdef 1 linear 00z01JUN1980    360mn
vars 5
t 12 99 temperature
w 12 99 vertical velocity
vo  12 99  vorticiy
d 12 99 divergence
r 12 99 relative humidity
endvars

xmgcqlmw

求解，为什么display单个变量没有问题，同时display多个变量数据就不对了
'set lev 1000'
     'd t'
'd w'
     'set lev 925'
     'd t'
    'd w'
'set lev 850'
     'd t'
    'd w'
'set lev 700'
     'd t'
    'd w'