c
c TO READ VIRS DAILY MEAN SST
c
parameter (idim = 2880, jdim = 609)
parameter (imiss = 254, land = 255)
parameter (rmiss = -999.0, rland = -777.0)
real*4 xsst(idim,jdim)
character*1 csst(idim,jdim)
integer*2 isst
character c1(2)*1, c2*2
equivalence (c1,c2)
101 format (a2)
c OPEN DATA FILE
open(10,file='virs_1day.19990101', access='direct',
2 form='unformatted', status='old', recl=idim*jdim, err=99)
c READ COUNT DATA
read(10,rec=1) csst
c CONVERT COUNT TO SST
do 100 i = 1, idim
do 110 j = 1, jdim
c1(2) = csst(i,j)
read(c2,101) isst
if (isst.eq.imiss) then
c MISSING
xsst(i,j) = rmiss
else if (isst.eq.land) then
c LAND
xsst(i,j) = rland
else
xsst(i,j) = real(isst)/10.0 + 10.0
endif
110 continue
100 continue
c CLOSE
99 close(10)
c STOP
stop
end
5. GrADSによる画像表示
このデータは、GrADSの "non float data" のフォーマットで書かれています。
GrADS(The Grid Analysis and Display System)はデスクトップ上でインタラクティブに利用できるツールです。GrADSはほとんどのUNIXワークステーション、マッキントッシュ、DOSまたはLinuxベースのパソコン上で使用可能です。GrADSのソフトウェアはanonymous FTPを通じて、自由に取得することができます。
DSET virs_1day.19990101
TITLE VIRS SST
OPTIONS yrev
UNDEF 254
XDEF 2880 LINEAR 0. 0.125
YDEF 609 LINEAR -38. 0.125
ZDEF 1 LEVELS 1000
TDEF 1 LINEAR 1jan1999 1dy
VARS 1
t1 0 -1,40,1 sst=t1/10+10
ENDVARS
GrADSに関する詳しい情報は、GrADSのホームページのドキュメントを参照下さい。
6. 参考文献
Kachi, M., Imaoka, K., Murakami, H., Nakajima, T. Y., and Shibata, A.
(1999): Preliminary results of TRMM: Part II SST retrieved from TMI
10 GHz and its expected uses. Submitted to Marine Technology
Society Journal.