program table
*_______________________________________________________________________

* created:  2/17/93

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

* purpose: create a table summarizing results for all cases, all groups

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

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

*_______________________________________________________________________

      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)

      character*4 group

      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))

* ingest spectrally summed results (k=0 only)

      k=0

      read (iounit,'(24x,f16.2)') bandv1(k),bandv2(k)
      read (iounit,'(24x,e16.4)') fluxup(nlevel,k),fluxdn(nlevel,k),
     +                           tropup(       k),tropdn(       k),
     +                           fluxup(     1,k)

      iounit=6

      group=header(1)(1:4)
      if(group.eq.'NASA') group='GLA '
      if(group.eq.'  AT') group='AER '
      if(group.eq.'JUNE') group='KIAE'
      write(iounit,902) group
     +                 ,iphase
     +                 ,icase
     +                 ,atmlbl
     +                 ,pthh2o
     +                 ,int(flgctm)
     +                 ,int(ppmco2)
     +                 ,ptho3
     +                 ,ppmch4
     +                 ,ppmn2o
     +                 ,bandv1(0)
     +                 ,bandv2(0)
     +                 ,fluxup(nlevel,k)
     +                 ,fluxdn(nlevel,k)
     +                 ,tropup(       k)
     +                 ,tropdn(       k)
     +                 ,fluxup(     1,k)
 902  format(1x,a4,2i5,2x,a3,1p,e10.2,2i5,1p,e10.2,0p,
     +       2f5.2,2f10.0,5f10.2)
      stop
      end