c      
c       This dataset is from NCAR achived 
c       ECMWF TOGA 2.5 degree Global Surface and Upper Air Analyses
c       http://dss.ucar.edu/datasets/ds111.2/

        program output
        parameter (lons=144,lats=73,levs=14,kvar=6)
        real*4 data(lons,lats,levs)

        open(8,file='198801.dat',
     &  form='unformatted',access='sequential')



        k=31     ! or 28,or 29, or 30 
        do 110 idy=1,kdy*2

        do 115 ivar=1,kvar
c       for ivar= 
c                1 Temperature [K]
C                2 Geopotential [m**2 s**-2]
c                3 Vertical velocity [Pa/s]
c                4 Relative humidity [%]
c                5 u-velocity [m/s]
c                6 v-velocity [m/s]
        do 120 L=1,levs
        read(8) ((data(i,j,L),i=1,lons),j=1,lats)
  120   continue

  115   continue
        read(8) ((data(i,j,1),i=1,lons),j=1,lats)
c       Read Mean sea level pressure [Pa]

  110   continue
 
        stop
        end