program label
*_______________________________________________________________________

* created:  2/17/93

* author : Bill Ridgway (ridgway@climate.gsfc.nasa.gov, 301-286-9138)

* purpose: create a file of data set labels (all groups)

* input  : standard input, any fm93 file in ../data

* output : standard output, piped to ../doc/labels

*_______________________________________________________________________

      integer iounit,mlayer,mlevel,mbands
      parameter (mlayer=150,mlevel=151,mbands=600)

* variables passed to or created by subroutines:
      character*80 header(10)
      character*3  atmlbl
      integer      iphase,icase,nlayer,nlevel,nbands
      real         psurf,tsurf,ptrop,flgctm
      real         ppmco2,ppmch4,ppmn2o
      real         pthair,pthh2o,pthco2,ptho3,pthch4,pthn2o
      real         bandv1(0:mbands),bandv2(0:mbands)
      real         plevel(mlevel),tlevel(mlevel)
      real         player(mlayer),tlayer(mlayer)
      real         alayer(mlayer),wlayer(mlayer),olayer(mlayer)
      real         fluxup(mlevel,0:mbands)
      real         fluxdn(mlevel,0:mbands)
      real         fluxnt(mlevel,0:mbands)
      real         heatrt(mlayer,0:mbands)
      real         tropup(       0:mbands)
      real         tropdn(       0:mbands)
      real         tropnt(       0:mbands)

      iounit=5

* blank out header initially

      do k=1,10
         write(header(k),'(80a1)') (' ',j=1,80)
      end do

* read 10 lines of header info from fm93 file

      read (iounit,'(a80)') header

* read key parameters which identify this case

      read (iounit,910) iphase,icase,atmlbl,nlayer,nlevel,nbands,
     +                  psurf,tsurf,ptrop,flgctm,
     +                  ppmco2,ppmch4,ppmn2o,
     +                  pthair,pthh2o,pthco2,ptho3,pthch4,pthn2o
 910  format(24x,i16,
     +     /,24x,i16,
     +     /,24x,13x,a3,
     +   3(/,24x,i16),
     +   3(/,24x,f16.2),
     +   1(/,24x,f16.0),
     +   3(/,24x,f16.2),
     +   6(/,24x,f16.4))


      iounit=6

      if(header(1)(1:4).eq.'NASA')
     +   write(iounit,'(1x,a4,2i5,2x,a3,5x,a,2(/,25x,a))')
     +                     'GLA '
     +                    ,iphase
     +                    ,icase
     +                    ,atmlbl
     +                    ,header(5),header(7),header(10)
      if(header(1)(1:4).eq.'  AT')
     +   write(iounit,'(1x,a4,2i5,2x,a3,5x,a,2(/,25x,a))')
     +                     'AER '
     +                    ,iphase
     +                    ,icase
     +                    ,atmlbl
     +                    ,header(5),header(7),header(10)
      if(header(1)(1:4).eq.'GFDL')
     +   write(iounit,'(1x,a4,2i5,2x,a3,5x,a,2(/,25x,a))')
     +                     'GFDL'
     +                    ,iphase
     +                    ,icase
     +                    ,atmlbl
     +                    ,header(3),header(5),header(10)
      if(header(1)(1:4).eq.'JUNE')
     +   write(iounit,'(1x,a4,2i5,2x,a3,5x,a,(/,25x,a))')
     +                     'KIAE'
     +                    ,iphase
     +                    ,icase
     +                    ,atmlbl
     +                    ,header(3),header(10)

 902  format(1x,a4,2i5,2x,a3,5x,a,2(/,25x,a))
      stop
      end