GrADSでは、コントロールファイルの設定によってメソ領域のデータを描画することができます(右上図)。しかしながら、ランベルト座標系のまま直接描画することはできず、等緯度経度座標に変換されて出力されます。 これをランベルト座標系のまま長方形の領域で描画し、地図を重ねることができるようにしました(右下図)。 入力は、テキスト(経度 緯度 値)またはグリッドデータが可能です。 なお、描画できるのは、ランベルト座標(正確には割円錐型のランベルト正角円錐座標)で基準緯度は北緯30度と60度、基準経度は東経140度の場合のみです。いわゆる MSM 領域で、2013年の領域拡張の前後に対応しています。 |
環境変数GASCRPに設定されているディレクトリに GrADS_Scripts.tar.gz を展開する。(右クリックして保存。うまく展開できない場合)
環境変数GASCRPを設定していない場合 >>
テキストで [ 経度 緯度 値 ] をタブ区切りor空白区切りで1行ずつ並べます。カンマ区切り(CSV形式)のデータは扱えませんので要注意!
並びが異なる場合(例えば [ 緯度 経度 値 ] の並びの場合)は、描画時に -xyz 2 1 3
などとオプションをつけることで読み込みのカラムの位置(列)を設定できます。
改行コードは「LF
」にしてください。(「CR+LF
」の場合は行末にスペースを入れると上手く動くかも?)
使用したいファイル名が tcpw.txt だとします。
と入力すると、単にforground colorでプロットされます。 |
!注意 デフォルトでは縦横比が正しくありません。 デフォルトの領域設定では68:55なので、
横の長さが 8.036、縦の長さが 6.5 で、約 68:55 となっています。 |
色々設定して、値に応じて色が変化するように描画したりなどもできます。 詳細は dtextの解説ページ をご参照ください。
|
領域を拡張前の旧領域にするには、 縦横比が異なる(5:4)ので、
以下、 |
全領域ではなく部分領域を指定するときには、 たとえば、デフォルト領域の南西4分の1を指定する場合は、以下のようになります(ついでに点の大きさを大きくしています)。
|
経緯線のスタイルを指定するときには、
さらに、1度毎の経緯線を追加するには、
|
海岸線のスタイルを指定するときには、 1つ目の引数で、解像度を指定します。 2〜4番目の引数で、線の色、スタイル、太さを指定します。
|
ここでは 15km の解像度の場合(273×221)で考えます。
ランベルト座標のグリッドのまま、X×Y(×Z×要素×時刻)で出力させれば問題ありません。
Y座標は北から並んでいるのをデフォルトとします。そうでない場合にも出力を並び替える必要は無く、ctlファイルで options yrev
としておけばよいです。
フツーに元の座標のまま、
xdef 273 linear 1 1
ydef 221 linear 1 1
とすればよいです。273,221 の値は適宜変えてください。他の部分も適宜作成してください。
(参考)通常のGrADSのpdefを使う場合は例えば以下のようになります。
pdef 273 221 lcc 30 140 189 73 30 60 140 15000 15000
xdef 571 linear 102 0.1
ydef 351 linear 16 0.1
使用するファイル名が rain.ctl、変数名が rain だとします。
ファイルを毎回閉じないと、地図を正しく描画できないのでご注意ください。
|
|
部分領域で描画するときには、
さらにテキストデータを重ねる場合は、
|
|