Forwarded mail....

Arpita Mandal arpita at PRL.ERNET.IN
Mon Jun 20 14:43:21 EDT 2005





---------- Forwarded message ----------
Date: Mon, 20 Jun 2005 00:51:03 -0700 (Pacific Daylight Time)
From: arpita at prl.ernet.in
To: GRADSUSR at LIST.CINECA.IT



Dear sir
 i am using the program anal2.f to read the different variables of the
file
anal.1998040100.288 .from GAME  reanalysis site.  sp, slp i could read
correctly and plot in grads. i
am having problem with the u, v,z variables. i can read the data for u,v
wind in binary form and wrote ctl file to plot them. in grads the contours
are coming from 300 to 1000 interval 100. but when i am using the
origininal ctl file anal_125.ctl given in GAME site and plotting u,v, the
values are coming from -20 to 20, -30 to 50, which is correct. the
anal_125.ctl file is provided in the GAME web site.
also i could not set the level.
when i am giving set lev 500
d u
i am getting contours from 300 to 1000,
the same is for all the levels.
kindly help me out
with regards
arpita

-------------- next part --------------
C     Last change:  AM    5 May 2005    2:32 pm
C----------------------------------------------------------------------------------
C    This is a sample fortran program to read 1.25 degree version of GAME-reanalysis
C      Ver.1.5 global objective analysis.
C
C    Variables:
C      slp : Sea Level Pressure
C      sp  : Surface Pressure
C      su  : Surface U-comp
C      sv  : Surface V-comp
C      st  : Surface Temperature
C      ssh : Surface Specific Humidity
C      z   : heights (17 layers)
C      u   : U-comp (17 layers)
C      v   : V-comp (17 layers)
C      t   : Temperature (17 layers)
C      sh  : Specific Humidity (1000-300 hPa)
C
C      17 layers are as follows.   Unit is hPa.
C        1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10
C
C   You should notice that these data were produced on the big_endian machine.
C
C   by K.Takahashi, MRI/JMA on June 2002        
C---------------------------------------------------------------------------------
C     (idim,jdim)=(1,1) and (288,145) correspond to (0E,90N) and (358.75E,90S) respectively.
C
      parameter (idim=288,jdim=145,lyr=17,lytd=8)
      real *4 a(lyr)
      data a / 1000,925,850,700,600,500,400,300,250,200,150,100,70,50,
     #     	 30,20,10/
      real*4    slp(idim,jdim),sp(idim,jdim),su(idim,jdim),
     #          sv(idim,jdim),st(idim,jdim),ssh(idim,jdim),
     #          z(idim,jdim,lyr),u(idim,jdim,lyr),v(idim,jdim,lyr),
     #          t(idim,jdim,lyr),sh(idim,jdim,lytd)
C
      open(10,file='@anal.1998040100.288',recl=idim*jdim*4,
     # form='unformatted', access='direct')
      open (11,FILE='sp.dat',recl=idim*jdim*4,access='direct')
      
      read(10,rec=1) slp
      read(10,rec=2) sp
      read(10,rec=3) su
      read(10,rec=4) sv
      read(10,rec=5) st
      read(10,rec=6) ssh
C
      do l=1,lyr
        read(10,rec=7) ((z(i,j,l),i=1,idim),j=1,jdim)
      enddo
      write(6,*) 'read z'
C
      do  l=1,lyr
        read(10,rec=8)((u(i,j,l),i=1,idim),j=1,jdim)
      enddo
      write(6,*) 'read u'
C
      do  l=1,lyr
        read(10,rec=9) ((v(i,j,l),i=1,idim),j=1,jdim)
      enddo
      write(6,*) 'read v'
C
      do  l=1,lyr
        read(10,rec=10) ((t(i,j,l),i=1,idim),j=1,jdim)
      enddo
      write(6,*) 'read t'
C
      do  l=1,lytd
        read(10,rec=11) ((sh(i,j,l),i=1,idim),j=1,jdim)
      enddo
      write(6,*) 'read sh'
      
      WRITE(11,rec=1) sp
      
      
      end



-------------- next part --------------
dset     ^v.dat
title     sample data set
undef    -9.99E33
options  byteswapped
options  yrev
xdef     288 linear   0.0  1.25
ydef     145 linear   -90  1.25
tdef     856   linear 0z01apr1998 6hr
zdef     17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10
VARS       1
v    17  0  u-comp
endvars
-------------- next part --------------
A non-text attachment was scrubbed...
Name: u.dat
Type: application/octet-stream
Size: 167040 bytes
Desc: dat file
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20050620/d1d8d035/attachment.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: anal_125.ctl
Type: application/octet-stream
Size: 618 bytes
Desc: ctl file
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20050620/d1d8d035/attachment-0001.obj 


More information about the gradsusr mailing list