nc转grd

大果子

我想把nc数据输出成四维grd，可是grads总是说looping dimension does not vary，谁能帮我看看哪里错了？谢谢！


'reinit'
y=2001
while(y<=2012)
'sdfopen f:\nc\vwnd.'y'.nc'
'set fwrite e:\30-60\hl\vmul'y'.grd'
t1=1
while(t1<=365)
'set gxout fwrite'
'set lon 100 130'
'set lat 0 60'
'set lev 1000 100'
'set t 't1''
'd vwnd'
t1=t1+1
endwhile
'close 1'
'disable fwrite'
y=y+1
endwhile
;

随缘

'set lon 100 130'
'set lat 0 60'
'set lev 1000 100'

你的这些维也要变才可以写的，grads不能直接写一个三维的数据块

river

你只循环了时间，高度也是需要用循环的。只有经纬度是可以直接像你那样设置的

sun_shine_Xia

i=1
while(i<789)

'set t 'i''

'set lev 1000'
*'d u'
'd v'
'set lev 925'
'd v'
'set lev 850'
'd v'
'set lev 700'
...
...
i=i+1
endwhile

pip2013

这样直接写是不行的，要用两个while嵌套循环，就是时间就像现在这样用while循环控制，高度层也要用while控制，嵌套在时间维里面，纬度和经度可以直接写

1092995269

两个while嵌套

大果子

sun_shine_Xia 发表于 2013-12-5 12:48

                               
登录/注册后可看大图

i=1
while(i

好的，谢谢啦！

大果子

pip2013 发表于 2013-12-5 14:49

                               
登录/注册后可看大图

这样直接写是不行的，要用两个while嵌套循环，就是时间就像现在这样用while循环控制，高度层也要用while控制 ...

恩，非常感谢！

flyapple1992

{:5_213:}{:5_213:}

flyapple1992

大果子 发表于 2013-12-13 16:16
恩，非常感谢！

你好，看了你的帖子，也想提取四维数据
gs如下
'sdfopen F:/ERAdata/allwinter.nc'

'set gxout fwrite'
'set fwrite F:/t.grd'
'set x 1 144'
'set y 1 37'
q=3
while(q<=110)

'set t'q
w=1
while(w<=37)
'set z' w
'd u'
w=w+1
endwhile
'disable fwrite'
q=q+1
endwhile

；
怎么会出错呢？
百思不得其解啊
望多多指教