GrADS講座

MENU

1. GrADSって?
2. データセットを作る
3. コントロールファイルを作る
4. GrADSを使う
5. おまけコーナー

• Fortranダイレクトアクセス形式(デフォルト)
• netCDF形式
• GRIB形式

---------------------------------------------------------
       real z(100,50,10)
             .
             .
       open (1,file='test.dat',form='unformatted',
      &      access='direct',recl=100*50*4)
       irec=1
       do 10 k=1,10
         write(1,rec=irec) ((z(i,j,k),i=1,100),j=1,50)
         irec=irec+1
10     continue
             .
             .
       stop
       end
----------------------------------------------------------

次に実際に私が使っているGANALデータ用のコントロールファイルを例に、コ ントロールファイルの内容を説明します。 ----------------------------------------------------------------

DSET	/home/yoshida/ganal/data/d%y2%m2%d2.%h2g
TITLE	GANAL DATA
UNDEF	-9.99E33
OPTIONS	template
XDEF	192 LINEAR 0.00 1.875
YDEF	49 LINEAR 0.00 1.875
ZDEF	9 LEVELS 1000 850 700 500 400 300 250 200 150
TDEF	1000 LINEAR 00z01mar1994 12hr
VARS	17
zps	0 99 sea surface pressure
ups	0 99 u winds at surface
vps	0 99 v winds at surface
tms	0 99 temperature at surface
qs	0 99 specific humidity at surface
tps	0 99 potential temperature at surface
tes	0 99 equivalent potential temperature at surface
zp	9 99 geopotential height
up	9 99 u winds
vp	9 99 v winds
tm	9 99 temperature
q	9 99 specific humidity
tp	9 99 potantial temperature
te	9 99 equivalent potential temperature
dtdp	9 99 dtdp
vort	9 99 vorticity
pvort	9 99 potential vorticity
ENDVARS

----------------------------------------------------------------

各行は次のような意味です。

DSET データファイル名

データの入っているデータファイルの名前です。ここでは、templateオ プションによって、/home/yoshida/ganal/data/の中の「d年2桁月2桁日2桁.時2桁g」の名前のついたファイルを 扱うようになっています。

TITLE タイトル

このデータのタイトルです。

UNDEF 値

データのない格子点を示す値です。この値が入っている格子点は表示な どの時に無視します。

OPTION template

ファイルの種類や扱い方によって、いくつかのオプションが使えます。 ここでは、ファイル名に任意の日付を用いて複数のファイルを1つのコントロー ルファイルで扱うことができる「template」オプションを使っています。

XDEF

X方向(東西方向)の格子の並び方を示します。表し方には「LINEAR」と 「LEVELS」の2通りの方法があります。どちらの場合も西から東向きに進みま す。

• LINEAR 始点の経度 格子間隔(経度)
• LEVELS 最初の点の経度 2番目の点の経度 3番目の点の経度 ．．．

この例では、「東経0度から1.875度間隔で192格子のデータ」であることを表してい ます。
（注:「LEVELS」は格子数が2以上の時に使えます。それ以外の時は「LINEAR」を使って下さい。）

YDEF

Y方向の格子の並び方を示します。XDEFと同様に「LINEAR」と「LEVELS」 があります。南から北へ向かって進みます。

ZDEF

Z方向の格子の並び方を示します。基本的に気圧で表します。同様に 「LINEAR」、「LEVELS」があります。下から上に向かって進みます。 この例では、「1000hPa、850hPa、700hPa、500hPa、400hPa、300hPa、 250hPa、200hPa、150hPaの9高度のデータ」であることを表します。

TDEF 時刻数 LINEAR 開始時刻 時間間隔

時間に関する表記です。この例では、「1994年3月01日0UTCから12時間間 隔1000時刻のデータ」であることを表しています。

VARS 変数数

入っている変数の数です。ここでは、17個の変数を使っています。

略称 入っている高度数 ユニット 変数の意味

入っている変数についての情報です。

略称	GrADSの中でこの変数を使う時の名前です。アルファベットで始まってアルファベットか数字で12文字以内で決めてください。大文字、小文字の区別はありません。
入っている高度数	その変数が入っている高度の数です。地上データなどの高度を気圧で表せない変数は0にしておきます。
ユニット	GRIBデータなど特殊なデータフォーマットを使 う時に使います。通常は0か99にしておきます。
変数の意味	その変数の情報です。40文字以内でご自由に お書き下さい。

ENDVARS

コントロールファイルの終りを示します。この行から後ろには何も書か ないで下さい。

/home/yoshida/@arika#grads

とコマンドgradsを打ち込みます。 すると、

Grid Analysis and Display System (GrADS) Version 1.7Beta9
Copyright (c) 1988-1997 by Brian Doty 
Center for Ocean-Land-Atmosphere Studies
Institute for Global Environment and Society 
All Rights Reserved

Config: v1.7Beta9 32-bit little-endian readline sdf/xdf 
hdf-sds lats athena im/image-output

Issue 'q config' command for more information.

Landscape mode? (no for portrait):  

と出てきます。

ここでそのままリターンキーを押すと横長の画面、「no」と打ち込んでリター ンキーを押すと縦長の画面が現れます。 そして、コンソール画面は

ga->

となります。これがGrADSのコマンド入力画面です。
まず、GrADSではコントロールファイルを開きます。 そのためのコマンドが「open」です。

ga->open ganal.ctl

ga->set time 0z11sep94
ga->set lon 125 150
ga->set lat 30 50

ga->d zps

---------------------------------------------------------------------
set lon 125 150
set lat 30 50
set time 0z11sep94
d zps
---------------------------------------------------------------------

そして

ga->exec test.ex