Global-merged IR - rffbl気象学的データ解析メモ

Global-merged IR

Globally-merged Full Resolution IR Brightness Temperature Dataもしくは、NCEP/CPC 4km Global (60N - 60S) IR Datasetの取り扱いについて説明する。
これはGOES-8/10、METEOSAT-7/5、GMSで観測された赤外輝度温度のデータを一つのデータセットにまとめたものである。
データの入手: Goddard Earth Sciences Data and Information Services Center (GES DISC)
　　　　　　　Science Portals > Precipitation > Data holdings > Data Product > Ancillary Global-merged IR

このデータセットは、GrADSで読み込むことを想定してctlファイルが用意されている。
Fortranで読み込む場合には1バイト符号なし整数の読み込みを参照してください。

このデータを使うときの注意点

1.衛星の観測になんらかのエラー(例えば、直線状の異常に小さな値)が存する場合がある。
2.必要とする領域の観測がおこなわれていない場合、その領域が一定値で占められているか、遠くの衛星を使って(解像度低い)場合がある。
いずれにしても、 使うデータはすべて一度図化して、問題がないことを確認する 必要がある。

Global-IRで、まとめて図化するためのGrADSスクリプトの例

*draw Tbb figure for each observation time (hourly)
****Parameter settings   Modify this!!!!*****
year='2010'
month='07'
west='120'
east='150'
south='30'
north='45'
***************************

'open sample.ctl'
'set display color white'
'c'
'set xlopts 1 0 0.15'
'set ylopts 1 0 0.15'
'set gxout grfill'
'set mpdset hires'
'set lon 'west' 'east 
'set lat 'south' 'north

***Conversion of month**********
if(month='01');monthEN='JAN';endif
if(month='02');monthEN='FEB';endif
if(month='03');monthEN='MAR';endif
if(month='04');monthEN='APR';endif
if(month='05');monthEN='MAY';endif
if(month='06');monthEN='JUN';endif
if(month='07');monthEN='JUL';endif
if(month='08');monthEN='AUG';endif
if(month='09');monthEN='SEP';endif
if(month='10');monthEN='OCT';endif
if(month='11');monthEN='NOV';endif
if(month='12');monthEN='DEC';endif
********************************

i=1
while (i<=31)

  j=0
  while (j<=23)
    hour=j

    'set grads off'
    if(i<=31)
      day=i
    'set time 'hour'z'day'monthEN'year
    endif

    'tbb = smth9(ch4+75)'
    'color 220 270 5 -kind red->yellow->lime->aqua->blue->white -gxout grfill'
    'd tbb'
    'cbarn'
    'q dims'
    dim_result = result
    dim_time = sublin(dim_result,5)
    time = subwrd(dim_time,6)
    hour = substr(time,1,2)
    day = substr(time,4,2)

    'draw title 'year'/'month'/'day'_'hour
    'printim tbb_snap' year'_'month'_'day'_'hour'.png'
    say year'_'month'_'day'_'hour

    'c'
    j=j+1
  endwhile
  i=i+1
endwhile
'quit'

高い雲の影響などで衛星から見て影になっている部分では値がundefであったり小さくなっているのでsmth9を使って、undefの部分を埋めています。
このスクリプトでは、1ヶ月分作成しますが31*24枚作成するのでかなり時間がかかります。