No subject

[Arpita Mandal]

Hi I have a dat file which i am attaching herewith. i have gridded it
using IDL. when i write the ctl file and plot tpw, the contours are coming
in inclined manner. kindly suggest help. the data is from GAME reanalysis
site. tpw(288,145) is the array of total ppt water for the entire globe of
0-360 lon and -90 90 lat. the grid is 1.25*1.25. i am using a fortran prog
to read the data from the file at 2dmon.1998040100.288.
please suggest some help
arpita



-------------- next part --------------
C     Last change:  AM   25 May 2005    3:56 pm
C----------------------------------------------------------------------------
C   This is a sample fortran program to read 1.25 dgree version of GAME-reanalysis
C     ver.1.5 physical monitor.
C
C   Variables are as follows.
C
C       rsub: Upward shortwave flux at bottom[W/m**2]
C       rsdb: Downward shortwave flux at bottom[W/m**2]
C       rlub: Upward longwave flux at bottom[W/m**2]
C       rldb: Downward longwave flux at bottom[W/m**2]
C       rsut: Upward shortwave flux at top[W/m**2]
C       rsdt: Downward shortwave flux at top[W/m**2]
C       rlut: Upward longwave flux at top[W/m**2]
C       cla : Total cloud amount
C       rcsb: Cloudforcing of shortwave at bottom[W/m**2]
C       rclb: Cloudforcing of longwave at bottom[W/m**2]
C       rof : Run off(m/day)
C       fllh: Latent heat flux[W/m**2]
C       flsh: Sensible heat flux[W/m**2]
C       ltrs: Transpiration[W/m**2]
C       lint: Interception[W/m**2]
C       rrc : Convective precipitation[mm/day]
C       rrl : Largescale precipitation[mm/day]
C       tpw : Total precipitable water[kg/m**2]
C       uwv : Water vapor flux[m/s]
C       vwv : Water vapor flux[m/s]
C
C   You should notice that these data were produced on the big_endian machine.
C
C   by K.Takahashi, MRI/JMA on June 27, 2002
C------------------------------------------------------------------------------
C     (idim,jdim)=(1,1) and (288,145) correspond to (0E,90N) and (358.75W,90S) respectively.
C
      parameter (idim=288, jdim=145)
      real*4 rsub(idim,jdim),rsdb(idim,jdim),rlub(idim,jdim),
     #       rldb(idim,jdim),rsut(idim,jdim),rsdt(idim,jdim),
     #       rlut(idim,jdim),cla(idim,jdim),rcsb(idim,jdim),
     #       rclb(idim,jdim),rof(idim,jdim),fllh(idim,jdim),
     #       flsh(idim,jdim),ltrs(idim,jdim),lint(idim,jdim),
     #       rrc(idim,jdim),rrl(idim,jdim),tpw(idim,jdim),
     #       uwv(idim,jdim),vwv(idim,jdim)
      OPEN(10,FILE='./@2dmon.1998040100.288',recl=idim*jdim*4,
     #       access='direct')
      write(6,*)'file read'
      OPEN(11,file = '2tpw0.dat',recl=idim*jdim*4,FORM='formatted')

      read(10,rec=18)tpw

      write(6,*)'tpw read'

      WRITE(11,5)((tpw(i,j),i=1,idim),j=1,jdim)
5     FORMAT(145(F12.6,4X))
      END






-------------- next part --------------
A non-text attachment was scrubbed...
Name: 2tpw0.xls
Type: application/vnd.ms-excel
Size: 775680 bytes
Desc: excel file
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20050609/88bef8e2/attachment.xls 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 2tpw0.dat
Type: application/octet-stream
Size: 667295 bytes
Desc: dat file
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20050609/88bef8e2/attachment.obj 
-------------- next part --------------
openr, 1, 'drh%96.dat'
a=fltarr(75,12)
readf, 1, a
close, 1
b=transpose(a)
openw, 1, 'drh%96.grd'
writeu, 1, b
close, 1
end













-------------- next part --------------
dset      2tpw1.grd
undef    -9.99e3
title     sample data set
options   little_endian yrev
xdef     288 linear   0  1.25
ydef     145 linear   -90 1.25
zdef      1 levels 1000
tdef     856 linear 00z01apr1998 6hr
VARS      1
tpw  0 0 total ppt water
endvars