c     ===================================================
c
c     example routine to read the ape lat-alt o3 vmr data
c
c     ===================================================
c
c     resolution: t42
c
c     data basis: annual mean amip-2 ozone climatology
c                 by Liang and Wang, symmetrized about
c                 the equator
c
c     date of release: february 26, 2003
c
c-----references:
c     liang, x.-z., and w.-c. wang, 1996: atmospheric ozone climatology
c            for use in gcms.
c     wang, w.-c., x.-z. liang, m.p. dudek, d. pollard, and
c            s.l. thompson, 1995: atmospheric ozone as a climate gas.
c            atmos. res., 37, 247-256.
c
c     -------------------
c     -------------------
      subroutine read_ape
c     -------------------
c     -------------------
c
      parameter (jg=64,kg=59,lg=kg+1)
      dimension glat(jg),galt(kg),gprs(kg),gdat(jg,kg)
     &                  ,gali(lg),gpri(lg)
      character finame*80,label*120
c
      lato3d  = jg  ! number of data latitudinal grid
      layo3d  = kg  ! number of data vertical layers
      lvlo3d  = lg  ! number of data vertical layer interfaces
c     o3vubc  =     upper boundary o3 vmr (ppmv)
c     prsubc  =     upper boundary pressure (mb)
c     altubc  =     upper boundary altitude (km)
c
      finame      ='APE-Ozone.T42-Lat-Alt'
      open (01,file=finame,status='old')

      read (01,'(i3)') lato3d

      read (01,'(i3)') layo3d
      lvlo3d = 1 + layo3d

      read (01,'(3(1pe12.5))') o3vubc,prsubc,altubc

      read (01,'(a)') label   ! grid latitudes (deg) --------->
      read (01,'(10(1pe12.5))') (glat(j),j=1,lato3d)

      read (01,'(a)') label   ! layer pressure (mb) ---------->
      read (01,'(10(1pe12.5))') (gprs(k),k=1,layo3d)

      read (01,'(a)') label   ! layer altitude (km) ---------->
      read (01,'(10(1pe12.5))') (galt(k),k=1,layo3d)

      read (01,'(a)') label   ! interface pressure (mb) ------>
      read (01,'(10(1pe12.5))') (gpri(k),k=1,lvlo3d)

      read (01,'(a)') label   ! interface altitude (km) ------>
      read (01,'(10(1pe12.5))') (gali(k),k=1,lvlo3d)

      read (01,'(a)') label   ! o3 vmr (ppmv) lat-alt
      read (01,'(10(1pe12.5))')
     &          ((gdat(j,k),j=1,lato3d),k=1,layo3d)
      close(01)

      return
      end