請教GrADS將nc檔轉dat檔的問題

奕紘

爬文看過論壇上一些相關文章後
自己嘗試寫一個看看
使用資料是1988和1989年的U風場資料(一天四筆)
我想只輸出1988/12、1889/1和1989/2三個月的數值
以下是我的gs檔程式碼 :


'sdfopen uwnd.1988.nc'
'sdfopen uwnd.1989.nc'
'set fwrite 88D89JF.dat'
'set gxout fwrite'
'set x 1 144'
'set y 1 73'
#December
t=1341
while(t<=1464)
z=1
while(z<=17)
'set t 't''
'set z 'z''
'd uwnd.1'
z=z+1
endwhile
t=t+1
endwhile

#January,February
t=1
while(t<=236)
z=1
while(z<=17)
'set t 't''
'set z 'z''
'd uwnd.2'
z=z+1
endwhile
t=t+1
endwhile

但是在執行時，出現錯誤碼
Data Request Warning:  Request is completely outside file limits
意思好像是說超出範圍
但是我的設定都是在資料範圍內呀...
把輸出的資料畫出來看，也只有1988/12有值
麻煩各位前輩幫我看一下了><"

lqouc

看上去感觉没什么问题啊，结尾有写disable fwrite么？

奕紘

lqouc 发表于 2013-2-9 15:15

                               
登录/注册后可看大图

看上去感觉没什么问题啊，结尾有写disable fwrite么？

沒有耶~不過我剛剛加上去後再試跑一次
還是出現了一樣的問題

lqouc

'set fwrite 88D89JF.dat'
'set gxout fwrite'
这两句的顺序换换试试？我也用过你这样的脚本，没出问题啊，真是帮不到你了。
话说你的12月不是写出来了么，实在不行就grads分开写好了，再用fortran写到一起。

奕紘

lqouc 发表于 2013-2-9 15:56

                               
登录/注册后可看大图

'set fwrite 88D89JF.dat'
'set gxout fwrite'
这两句的顺序换换试试？我也用过你这样的脚本，没出问题啊 ...

順序調換後還是一樣呀(sad...)

因為接下來需要同時輸出很多年的資料到一個檔案裡
如果真不行，只後分別寫出再用fortran寫到一起好了..

Mr.治愈系

表示完全看不出来有什么问题- -。。学识不够啊- -