关于数据转换

haohaoai

    我拿到的ncep的数据，数据读取后的信息如下图


我想把这个nc数据转为grd文件，在gs里编辑如下
'reinit'
'set  gxout  fwrite'
'sdfopen  e:\shuju\vwnd.mon.ltm.nc'
'set  fwrite  e:\shuju\vwnd.grd'
num=1
while(num<=60)
'set t ''num'
'set x 1 145'
'set y 1 73'
'set lev 1000 100'
'd  vwnd'
num=num+1
endwhile
'disable fwrite'
;
然后在grads里运行不了，请帮忙看看是哪儿的问题，谢谢

Aires

t是不是写错了啊？'set t 'num

做个霸气的木头

set t ''num'这句你好好查查，细节问题

残月

'set t ''num'
'set x 1 145'
'set y 1 73'
'set lev 1000 100'   

这个估计还得抽出一维写循环。。
'set t ''num'--------》》‘set t 'num'’

haohaoai

Aires 发表于 2012-10-28 09:51

                               
登录/注册后可看大图

t是不是写错了啊？'set t 'num

我改了，可还是不行，依然显示

兰溪之水

'reinit'
'set  gxout  fwrite'
'sdfopen  e:\shuju\vwnd.mon.ltm.nc'
'set  fwrite  e:\shuju\vwnd.grd'
num=1
while(num<=12)
'set t ''num'
'set x 1 145'
'set y 1 73'
'set lev 1000 100'
'd  vwnd'
num=num+1
endwhile
'disable fwrite'
;
若不行，再添加一个层次循环~

haohaoai

兰溪之水 发表于 2012-10-28 10:13

                               
登录/注册后可看大图

'reinit'
'set  gxout  fwrite'
'sdfopen  e:\shuju\vwnd.mon.ltm.nc'

还是不行，我是这样想的，这是一个四维数据，从1948年的1月开始计算到2017年的12月，每年12个月平均值，然后想把这个数据转换，你所说的再添加一个层次循环，具体要怎么添加，因为新手，对这个了解不多，查看了一些帖子才勉强写出这个转换程序

haohaoai

Aires 发表于 2012-10-28 09:51

                               
登录/注册后可看大图

t是不是写错了啊？'set t 'num

我发现问题了，我把'set lev 1000.100‘换位;set lev1000'就能转了，如果我这样改的话，读取的数据是不是只限定在1000hpa的位置了，其余的数据是不是就抹掉了？是不是需要在镶嵌一个循环语句把其余的读出来

Ivy

把set level也写成一个循环试试
'reinit'
'set  gxout  fwrite'
'sdfopen  e:\shuju\vwnd.mon.ltm.nc'
'set  fwrite  e:\shuju\vwnd.grd'
num=1
while(num<=12)
'set t ''num'
z1=1
while(z1<=*)
'set z 'z1''
'set x ...'
'set y ...'
'd ..'
.....
'set x 1 145'
'set y 1 73'
'set lev 1000 100'
'd  vwnd'
num=num+1
endwhile
'disable fwrite'
;

haohaoai

Ivy 发表于 2012-10-28 10:51

                               
登录/注册后可看大图

把set level也写成一个循环试试
'reinit'
'set  gxout  fwrite'

万分感谢，在你这个基础上我改变成功了