推 erdosl :恩,抱歉,濾是在IDL濾,濾之前,要先在grads擷取出130天 09/10 20:38
推 erdosl :好,我試試看這樣行不行,謝謝 09/10 20:40
推 erdosl :感謝m大,問題出在ctl要改成17才能畫z=2以上,解決了 09/10 21:20
※ 引述《erdosl (stupid)》之銘言：
: 小弟想用GrADS自.nc檔擷取多層z層的氣象資料(u,v,...etc),因為用IDL濾波,但是不知道
: 我的"山寨版"IDL是否真的濾對(因為我跟老師說,我的電腦濾很快,他不相信,要我檢查,囧)
: 我的想法是,固定某一地理位置(x,y固定),挑選某一層(z固定),時間選一段(t=1~16),產生
: 時間序列,濾前和濾後圖型比較;如果不一樣,不見得代表濾成功了,但是,如果一樣,就根本
: 沒濾到了.
我不太懂你所謂的「濾」是要做啥，似乎你上面寫的跟你 gs 檔的內容不太一樣？
根據你提供的 fwrite-130d.gs 檔案
看起來會是 fwrite 出：
j=223～352（130個j），x=1～144，y=1～73，z=1～17
: 可是,我碰到一個詭異的問題,就是,我寫.gs檔擷取後產生的.dat檔,附上.ctl說明檔後,用
: GrADS繪圖,竟然畫不出來z = 2以上!!出現以下錯誤訊息:
: Data Request Warning: Request beyond file limits
: Entire grid contents are set to missing data
: Dimension ranges are X = 100 144 Y = 50 73 Z = 2 1 T = 1 130 E = 1 1
: Warning issued for variable = air
: Cannot plot data - all undefined values
: 請問這是犯了甚麼愚蠢的錯誤嗎??
: ps :
: 1. 這是set x 100 y 50 z 2 t 1 16; d air所產生的錯誤訊息
: 2. 我的.gs檔(fwrite) :http://0rz.tw/p8HvM
: .ctl檔 :http://0rz.tw/pZQdK
ctl 檔裡面，第九行應該是：
"air 17 99 T"，而不是"air 1 99 T"
先試試看這樣改行不行
如果不行，我會猜是 "Entire grid contents are set to missing data"
也就是先不要 set gxout fwrite
先 set gxout stat
（stat 表示要 GrADS 輸出資料的統計結果）
看看你設定的範圍是否都是 missing data
另外，set gxout stat 的結果，如果你用 gs 檔連續輸出，螢幕不會顯示
要一個時間一個時間先檢查看看
Mark
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.67.14