このwiki内では、グリッドファイル(GrADS 4byteバイナリ、特に断らなければdirect書き込み)を〜.bin、コントロールファイルを〜.ctl、変数名をvarと表記します。
'set gxout shaded'
open 〜 sh d 変数名
unzip gradslib.zip
export GASCRP=/・・・/gradslib #gradslibへの絶対パス
wget ftp://cola.gmu.edu/grads/scripts/lpoly_lowres.asc wget ftp://cola.gmu.edu/grads/scripts/lpoly_mres.asc wget ftp://cola.gmu.edu/grads/scripts/lpoly_hires.asc wget ftp://cola.gmu.edu/grads/scripts/opoly_lowres.asc wget ftp://cola.gmu.edu/grads/scripts/opoly_mres.asc wget ftp://cola.gmu.edu/grads/scripts/opoly_hires.asc wget ftp://cola.gmu.edu/grads/scripts/lpoly_US.asc
auxlinex 1.0 10.0 2.0 -c 1 -th 4 -d 5 -n 6 -s 0.05
auxlinex 1.0 10.0 2.0 -c 1 -th 4 -d 1 -n 6 -s 0.07 -st 2
'reinit' 'set display color white' 'c' 'sdfopen precip.mon.ltm.nc' #CMAP降水量気候値 'set mproj scaled' #x軸、y軸の位置がずれないように設定 'set lon 0 360' 'set lat -90 90' 'set xlint 60' 'set ylint 30' 'set parea 1 10 2 7' 'set xlopts 1 6 0.22' 'set ylopts 1 6 0.22' 'set gxout shaded' 'color 0 10 0.2 -kind white->blue->lime' 'set grads off' 'd ave(precip,t=1,t=12)' 'auxlinex 1.0 10.0 2.0 -c 1 -th 4 -d 5 -n 6 -s 0.05' #南x軸に補助目盛りを表示 'auxliney 1.0 2.0 7.0 -c 1 -th 4 -d 5 -n 6 -s 0.05' #西y軸に補助目盛りを表示 'auxlinex 1.0 10.0 7.0 -th 6 -d 1 -s -0.07' #北x軸に目盛りを表示 'auxlinex 1.0 10.0 7.0 -th 4 -d 5 -s -0.05' #北x軸に補助目盛りを表示 'auxliney 10.0 2.0 7.0 -th 6 -d 1 -s -0.08' #東y軸に目盛りを表示 'auxliney 10.0 2.0 7.0 -th 4 -d 5 -s -0.05' #東y軸に補助目盛りを表示 'xcbar 4.5 6.5 1.0 1.4 -fwidth 0.18 -fheight 0.19 -edge box -line off -fstep 25'
d var basemap L/O
basemap L 色番号1 色番号2
wget ftp://cola.gmu.edu/grads/scripts/lpoly_lowres.asc wget ftp://cola.gmu.edu/grads/scripts/lpoly_mres.asc wget ftp://cola.gmu.edu/grads/scripts/lpoly_hires.asc wget ftp://cola.gmu.edu/grads/scripts/opoly_lowres.asc wget ftp://cola.gmu.edu/grads/scripts/opoly_mres.asc wget ftp://cola.gmu.edu/grads/scripts/opoly_hires.asc wget ftp://cola.gmu.edu/grads/scripts/lpoly_US.asc
basemap L/O 色番号1 色番号2 L/M/H
cbar_line -x 3.0 -y 6.0 -c 1 -m 0 -l 1 -t "line1"
cbar_line -x 3.0 -y 6.0 -c 1 2 3 -m 0 -l 1 -t "line1" "line2" "line3"
cbar_line_set -x 3.0 -y 6.0 -n 6 -l 1 -c 1 -m 1 -z 0.13 -sn 4 -sc 1 -sz 0.17 -t "line1"
cbar_line_set
cbar_line_set -x 3.0 -y 6.0 -r 0.08 -n 6 -l 1 -c 1 -m 1 -z 0.13 -sn 4 -sc 1 -sz 0.17 -t "line1"
cbarn
cbarn 1 0 5.5 0.18
cbarn scale num xmid ymid
set parea 0.8 10.2 1.0 8.0 set lon 0 360 set lat -90 90 set xlint 30 set ylint 30 set xlopts 1 6 0.2 set ylopts 1 6 0.2 set gxout shaded color 〜 set grads off d var cbarn 1.2 0 5.5 0.3
set mproj nps/sps d var circlat 30 90 0.18
circlat 30 90 0.18 1
set mproj nps/sps d var circlon 30
circlon 30 0.18 1
color 最小値 最大値 間隔 d var
color 最小値 最大値 間隔 -kind saddlebrown->white->green d var
color 最小値 最大値 間隔 -kind white->green d var
color ...
set ccols 16 17 0 0 20 21
darkblue->blue->dodgerblue->skyblue->white->white->gold->darkorange->red->darkred乾燥から湿潤
maroon->saddlebrown->darkgoldenrod->khaki->white->white->palegreen->lightgreen->limegreen->darkgreen一定値以上とそれ以下を強調(降水量や水蒸気量の気候値など)
white->antiquewhite->mistyrose->lightpink->mediumvioletred->navy->darkblue->blue->dodgerblue->aqua
colormaps -l 265 310 1 -map ocean
colormaps -flipped -l 265 310 1 -map ocean
colorcb 塗り分け数 タイプ set clevs -3 -2 … d var
colorcb 9 c set clevs -4 -3 -2 -1 1 2 3 4 d var
d var hatch var 3 5
hatch var 3 5 -angle 45 hatch var 3 5 -angle 135
hatch var 3 5 -density 0.1 -int 0.5
hatch var -stlye 1 -color 1 -thick 3
map type d var
markplot 140.0 40.0 -m 3 -s 0.2 -c 1
'set lon 120 152' 'set lat 20 48' 'set mpdset hires' 'set clevs -1e+30' 'd var' 'markplot 127.41 26.13 -c 1 -m 2 -s 0.25' 'markplot 135.29 34.41 -c 2 -m 3 -s 0.25' 'markplot 139.45 35.41 -c 3 -m 4 -s 0.25' 'markplot 141.21 43.03 -c 4 -m 5 -s 0.25'
mul 2 2 1 1 d var(t=1) mul 2 2 1 2 d var(t=2) mul 2 2 2 1 d var(t=3) mul 2 2 2 2 d var(t=4)
set gxout scatter set t 1 30 d var1;var2 regline var1 var2 -c 1 -l 1 -t 5
'open hoge.ctl' 'set parea 1 5 1 5' 'set line 1 1 5' 'set strsiz 0.15 0.16' 'taylor -s 0.8 1.3 1.7 -r 0.7 0.95 0.9 -t A B C -z 0.25 -c 1 2 4 -m 2 6 8'
trackplot 120.0 20.0 140.0 40.0 -c 1 -l 1 -t 6
#markplotと同じ工程# 'd var' 'markplot 127.41 26.13 -c 1 -m 2 -s 0.25' 'markplot 135.29 34.41 -c 2 -m 3 -s 0.25' 'trackplot 127.41 26.13 135.29 34.41 -c 2 -l 1 -t 6' 'markplot 139.45 35.41 -c 3 -m 4 -s 0.25' 'trackplot 135.29 34.41 139.45 35.41 -c 2 -l 1 -t 6' 'markplot 141.21 43.03 -c 4 -m 5 -s 0.25' 'trackplot 139.45 35.41 141.21 43.03 -c 2 -l 1 -t 6'
set line 1 1 5 set cthick 5 vec var1;var2 -SCL 1 4 -P 3 0.5
set t 最小値 最大値 xamin var -pause
set t 最小値 最大値 xamin var -sec 1
set t 最小値 最大値 xamin var -sec 間隔
xcbar 2.0 9.0 0.6 1.0 -fwidth 0.24 -fheight 0.26 -edge triangle -fstep 2 -line on
-line on -lc 0
set line 1 1 6 # 線の太さを6に xcbar 2.0 9.0 0.6 1.0 -line on
set ccolor 1 d var yrmask 1981 1982 1983 1985 1986 1988 1990 1994 1995 1996 1999 2000 2003 2005 2009 2010 2012 a # 時系列データaとして、指定した年は1、それ以外は欠損値のデータをdefine set ccolor 2 d var*a # 指定した年のみ描く
set cthick 3 set ccolor 1 set cmark 0 d const(lev,0) set cthick 5
set vrange -3 3 zeroline
zeroline color mark linetype thickness value
zeroline -help
set t 1 100 fcorr var1 var2 set t 1 d fcorrout
set t 1 100 ltrend var a b c d a # 線形回帰直線の表示(1次元の場合のみ) set t 1 d b # 傾きの表示 d c # 標準偏差の表示
set t 1 34 ltrend hoge a b c criteria=b*57.2*sqrt(34-2)/sqrt(34)/c
monclim var clim 1 240 set t 1 12 d clim
function clim(args) var=subwrd(args,1) varclim=subwrd(args,2) tst=subwrd(args,3) ten=subwrd(args,4) tst2=tst+12 'set t 'tst' 'tst2'' 'define 'varclim'=ave('var',t+0,t='ten',1yr)' 'modify 'varclim' seasonal' 'set t 'tst''
cp (GASCRPで指定しているディレクトリ)/pinterp.gs (〃)/pinterp.gsf
rc=gsfallow("on") # GrADSスクリプト関数を使用 'open var.ctl' 'set x 1 ???' 'set y 1 ???' 'set lev 1000 700' # 800hPaを挟む二つのレベルを含むように範囲を設定 'define a='pinterp(var, lev, 800) # 800hPa面に内挿したvarをaと定義 'set lev 800' # 適当に高さを一つに設定 'd a'
d var # 変数を表示 pstat var # 画面上部とターミナルウィンドウに統計量(最小値、最大値、 平均値、標準偏差)を表示
textout num hoge.txt var
$GADDIR
cp (GASCRPで指定しているディレクトリ)/zinterp.gs (〃)/zinterp.gsf
rc=gsfallow("on") # GrADSスクリプト関数を使用 'open var.ctl' 'set x 1 ???' 'set y 1 ???' 'set lev 0 2000' # 1,000mを挟む二つのレベルを含むように範囲を設定 'define a='zinterp(var, lev, 1000) # 1,000mの高度に内挿したvarをaと定義 'set lev 1000' # 適当に高さを一つに設定 'd a'