1. 程式人生 > >GrADS,NCL一些經驗

GrADS,NCL一些經驗

.gz ... 變化 如果 無法使用 更改 2個 open sea

GrADS畫特定經線方法

set clevs 23.5; d lat set clevs 120;d lon GrADS 查看多個打開文件的ctl q ctlinfo 1 q ctlinfo 2 ... 2個不同的缺測值 fwrite更改默認缺測值 GrADS用fwrite輸出時,默認的缺測值是-9.9900000E+08,如果輸入源的缺測值與之不一致,則會出現兩種缺測值的情況。解決的辦法,可以在fwrite 聲明後面設置默認缺測值設置。比如 ‘reinit‘ ‘open J:/ec79-15/ansnablameansst79_15.ctl‘ ‘set t 1 12‘ ‘save=ave(absnablamsst,t+0,t=444,12)‘ ‘modify save seasonal‘ ‘set gxout fwrite‘ ‘set undef 9.999E+20‘ ‘set x 1 480‘ ‘set y 1 89‘ ‘set fwrite J:/ec79-15/seasonal.grd‘ ‘d save‘ ‘disable fwrite‘ ‘reinit‘ ; (註:modify save seasonal命令意為將save變量修改成季節變量,方便與其它不同大小的變量計算(比如求季節距平),此外modify還有另外一個選項diurnal日變化) NCL使用自定義色標
清風給的調色盤確實好用,生成的rgb格式也可以。但是最近再次遇到了ncl無法正常讀取colormap錯誤的問題,參考了一下官網。最後發現是這樣的。官網給出的文件有“ncolors=8”這一行,用NCL調色盤生成的腳本也有這一句話。但是在我的ncl裏面無法使用。最後刪了這一句話之後就好了。另外,註意前兩行是用來定義底色和填色的。需要手動添加。比如加的是“255 255 255”和“0 0 0”這兩行。最後能正常生成了。(一句話,用GrADS調色盤生成的rgb文件需要刪掉nclolors這一行,並在文件最前面加上“255 255 255”和“0 0 0”這兩行(定義底色)) (註:今天發現rgb文件最後一行上也要加,貌似是這樣的:第一行定義,第二行定義填充線色,最後一行表示背景色) GrADS使用NCL色標
(20171102註:GrADS使用colormap命令使用NCL色標時候,出現了一個問題,發現colormap打開的相關rgb文件在ncolors那一行後面還要有 # r g b 這一行,否則可能無法正常使用。 Cressman插值(oacre函數) oacres函數實際上就是對數據進行多次cressman插值,其插值的影響半徑是逐漸減少的,默認的影響半徑(單位是格點)是10,7,4,2,1。如果某格點影響半徑範圍內的站點數目小於3個,那麽該格點的值設置為缺測值。(參考http://cola.gmu.edu/grads/gadoc/gradfuncoacres.html) oacres函數源代碼 ftp://cola.gmu.edu/grads/2.1/grads-2.1.1.b0-src.tar.gz,解壓後src目錄下gafunc.c文件,第5313-5602行 GrADS 坐標
調字體 set xlopts 4 1 0.1 設置X軸標記數字的顏色、粗細和大小(英寸) set ylopts 4 1 0.1 設置Y軸標記數字的顏色、粗細和大小(英寸) 調坐標間隔 set xlint interval set ylint interval set xlevs/ylevs 再控制標記的分布

GrADS,NCL一些經驗