program rd_t1a_flds
c
c   sample routine to read in 4x daily t1a fields average over eba
c
c   g77 rd_t1a_flds.f -o rd_t1a_flds.x  
c
      parameter (npds = 160, kk=40 )       

      parameter (xmis = -999.)

c   pr  - pressure (hPa)
c   z   - height 
c   u   - zonal-wind component
c   v   - meridional-wind component
c   t   - temperature (C)
c   th  - potential temperature (K)
c   r   - water vapor mixing ratio (g/kg)
c   rh  - relative humidity (%)
c   div  - horizontal divergence (1/s)*10^6
c   vor  - horizontal vorticity (1/s)*10^6
c   w  -  omega (mb/hr)
c   q1  -  apparent heating (K/day) 
c   q2  -  apparent drying (K/day) 
c   ht  -  horizontal advection of T (K/s) 
c   vt  -  vertical  advection of T (K/s) 
c   hq  -  horizontal advection of q ((g of vapor/g of air)/s) 
c   vq  -  vertical  advection of q ((g of vapor/g of air)/s)) 
c   np  -  number of points in EBA average
c
      real z(npds,kk), t(npds,kk), th(npds,kk), r(npds,kk), rh(npds,kk)
      real u(npds,kk), v(npds,kk), w(npds,kk)
      real vor(npds,kk), div(npds,kk)
      real pr(kk), q1(npds,kk), q2(npds,kk)
      real ht(npds,kk), vt(npds,kk), hq(npds,kk), vq(npds,kk)

      integer iy, im, id, ih, np(kk)                        

c   *************************************************************************

c   open files for input

      open(27, file='basic_v3.t1a', status='old')
      open(28, file='q1q2_v3.t1a', status='old')
      open(29, file='lsf_v3.t1a',   status='old')

      do ip=1,npds
         read(27,223) iy, im, id, ih              
         read(28,223) iy, im, id, ih              
         read(29,223) iy, im, id, ih             
         print *, 'reading in data for', iy, im, id, ih

         do  k=1,kk
            read(27,227) pr(k), z(ip,k), u(ip,k), v(ip,k), w(ip,k), 
     2                          t(ip,k), th(ip,k), r(ip,k), rh(ip,k), 
     3                          div(ip,k), vor(ip,k), np(k)
            read(28,228) pr(k), q1(ip,k), q2(ip,k), np(k)
            read(29,229) pr(k), ht(ip,k), vt(ip,k), hq(ip,k), 
     2                          vq(ip,k), np(k)
         enddo
      enddo


  223 format(4i3)
  227 format(11f9.2,i5)
  228 format(3f8.2,i5)
  229 format(f8.2,1p,4e11.3,i5)

      end