GrADSコマンドのメモ

コマンドプロンプトで grads と打ってリターンして下さい。 GrADSが起動します。その後、no または yes と打つと、縦長か横長の画面が 現れます。コマンドを打ち込みながら絵を描いていきます。 これらのコマンドを一つのファイルに書いて、"exec hoge.exe" などと 実行することにより、一連の操作をまとめて実行することができます。

以下の実行ファイルはすべて縦長(portrait mode)として実行したものです。 起動時に grads -p とするか、起動した後no としてください。

gradsnc がない場合は、~/.alias に

alias gradsnc='/usr/local/grads-1.9b4/bin/gradsnc'

1. SLPと、風、降水を描いてみる

/home/meer/nishii/gakubuenshu/2008summer/JRA-slp.exe 参照

2. 6月平均の日本付近のSLPと、風、降水を描いてみる

JRA-slp-pr-wind.jun.exe 参照

3. 2007年8月の風分布と高度場分布

/home/meer/nishii/gakubuenshu/2008summer/JRA-u-z.exe 参照

4. 東西平均東西風と温度

/home/meer/nishii/gakubuenshu/2008summer/JRA-u-T-znlmean.exe 参照

5. downstream development

/home/meer/nishii/gakubuenshu/2008summer/JRA-downstreamdevelopment.exe 参照

6. 画像ファイルに出力する

/home/meer/nishii/gakubuenshu/2008summer/figure.exe 参照

/home/meer/nishii/gakubuenshu/2008summer/JRA-u250-jan.exe 参照

/home/meer/nishii/gakubuenshu/2008summer/JRA-v250-var.exe 参照

/home/meer/nishii/gakubuenshu/2008summer/JRA-v250-cor.exe 参照

/home/meer/nishii/gakubuenshu/2008summer/miroc-ua250-tas.exe を各自のディレクトリにコピーした後、"grads -p" で起動し、"exec miroc-ua250-tas.exe" と打ち込んでみて下さい。絵が勝手に描かれて、miroc-ua250-tas.eps とmiroc-ua250-tas.png という画像ファイルができる筈です。

GrADSの本家のページは以下です。
http://www.iges.org/grads

• netCDFのファイルを開く

  > sdfopen /h1a/JRA25/plev2.5/Monthly/nc/z.month-mean.nc
  > sdfopen /h1a/JRA25/plev2.5/Monthly/nc/u.month-mean.nc
  > sdfopen /h1a/JRA25/plev2.5/Monthly/nc/v.month-mean.nc
  

• gradsのマーク?を消す

  > set grads off
  

• 絵を書く領域を指定する。

  > set vpage 0 6 0 4
  

  gradsnc -p で起動した場合は "0 8.5 0 11"、 gradsnc -l で起動した場合は "0 11 0 8.5"がデフォルト。

• 現在開いているファイルの情報を見る

  > q file
  > q file 2
  > q files
  

• 初期化(立ち上げた直後は要らない)

  > reinit
  

• 描く経度を指定する(負の値も可能)。２つ指定する際には小さい値の次に大きな値を書く。

  > set lon 90
  > set lon 0 180
  

• 描く緯度を指定する(負の値は南緯)

  > set lat -45
  > set lat 10 80
  

• 気圧面指定

  > set lev 1000 10
  > set lev 250
  

• 気圧面指定

  > set time 1jan1999
  > set t 200
  

• 平均操作 経度方向

  > ave(z,lon=0,lon=180) 
  > ave(z,x=1,x=100)
  

  緯度方向

  > ave(z,lat=0,lat=80)
  > ave(z,y=10,y=30)
  

  時間平均

  > ave(z,time=1jul,lat=15jul)
  > ave(z,t=1,t=12)
  > ave(z,t-2,t+2)
  

  各年1月の21年間分の平均を平均して zavejan という変数に代入

  >define zavejan = ave(z,time=1jan1979,time=1jan1999,1yr)
  

• 変数に値を代入(今セットされている領域で値が入る)

  > define meanz = ave(z,lon=0,lon=360)
  

• 等値線の太さ(1-10)

  > set cthick 5
  

• 等値線のラベルの色(1-20)、太さ(1-7),大きさ(インチ)

  > set clopts 2 4 0.3
  

• 等値線の間隔

  > set cint 50
  

• 描く等値線の値を決める

  > set clevs 200 400 600 800
  

• 変数の値を表示 (time=,,) で一時的に時間の指定ができる

  > d z
  > d z(time=1jan2000)-zavejan
  

• 矢印でプロットするときに矢の先端の大きさ(インチ)

  > set arrowhead 0.1
  

• 矢印でプロットするときに矢の長さ(0.5インチでどれくらいに相当させるか)

  > set arrscl 0.5 40
  

• 矢印(ベクトル)の表示

  二つ以上のctlファイルを開いている場合は、２つ目、３つ目、、の変数は、.2, .3,,, とつける。

  > d u.2;v.3

• グリッド間隔をあける(vには必要ない)

  > d skip(u.2,4,4);v
  

• タイトルをつける

  > draw title Z wind at 250hPa
  

• 画像ファイル(tmp.gm)に落とす

  > enable print tmp.gm
  > print
  > disable print
  

• shellのコマンドを実行する(ここでは例としてls)

  > !ls
  

• スクリプト(コマンドを順番に列記したテキストファイル)を読み込ませて実行する

  > exec tmp.exe
  
  
  

戻る