GrADS經緯度比例方法
*取消經緯比例自動調整
'set mproj scaled'
設置繪圖區域,使得一致
'set vpage 1 7 1 7'
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 風羽朝向 set hempref auto|nhem|shem
默認的是坐標在北半球,則風羽朝向風的左側,在南半球,則風羽朝向風的右側。一般情況下是沒問題的。但是,如果畫的是有特殊含義的圖(坐標經緯度是相對的經緯度,那么就需要修改這一選項)
NCL執行出現問題 a warning of 'Get System Resource'
沒有設置NCARG_ROOT環境變量
NCL讀取地圖數據
https://www.ncl.ucar.edu/Applications/topo.shtml
利用grib2ctl對ERA5預報場數據生成ctl
grib2ctl -verf 10fg.grib
