SUBROUTINE SUMMARY

      INCLUDE 'PARAMETERS.INC'


      CHARACTER*8 AMOL(MB) , APHOTO(MB) , APCE(MB) , AREACT(IB,6)
      CHARACTER SYMBOL(IB,6) , ACOEFF(IB,6)

      COMMON

     &      /ONED17/  NTRANS

     &      /ONED40/  MTRANS(MB)

     &      /ONED43/  AMOL

     &      /BLOCK5/  ACOEFF , SYMBOL , APHOTO , APCE ,  NPHOTO

     &      /ONED28/  DEN(NALT)   !!!!!!!!!mes 3/14/97

     &      /BLOCK37/ AREACT

     &      /BLOCK38/ NREACT

     &      /BLOCK7/  RATE(IB,NALT)

     &      /DEBUG/   RPROD(MB,NALT)

     &      /ONED59/  RLOSS(MB,NALT)

     &      /ONED60/  BNDRYFLX(MB,NALT)

      DIMENSION RPCOLMN(MB) , RLCOLMN(MB) , RCOLUMN(IB)


C THIS SUBROUTINE PRODUCES A SUMMARY FILE WHICH CONSISTS OF COLUMN INTEGRATED
C NET PRODUCTION RATES, NET LOSS RATES AND TOP AND BOTTOM FLUXES FOR EACH
C SPECIES. FOR EACH REACTION, THE COLUMN RATE IS OUTPUT AS WELL.

C*********************************************************************
      print * , 'Subroutine SUMMARY called'

C		First get Total Column Rates and Fluxes
      deltax = 2.0e+5


C CALCULATE COLUMN PRODUCTION/LOSS RATES

      DO M = 1 , NTRANS
        MVARY = MTRANS(M)
        RP = 0.0
        RL = 0.0
        DO K = 1 , NALT
          RP = RP + RPROD(MVARY,K)
          RL = RL + RLOSS(MVARY,K)
        ENDDO
        RLCOLMN(MVARY) = RL*deltax
        RPCOLMN(MVARY) = RP*deltax
      ENDDO

C CALCULATE COLUMN REACTION RATES

      DO I = 1 , NREACT
        R = 0.0
        DO K = 1 , NALT
          R = R + RATE(I,K)
        ENDDO
        RCOLUMN(I) = R*deltax
      ENDDO

C PRODUCE SUMMARY FILE

      WRITE(10,100)
      DO M = 1 , NTRANS
        MVARY = MTRANS(M)
        WRITE(10,110) MVARY , AMOL(MVARY) , RPCOLMN(MVARY) ,
     &                RLCOLMN(MVARY) , BNDRYFLX(M,1) ,
     &                BNDRYFLX(M,NALT)
      ENDDO

      WRITE(10,120)
      DO I = 1 , NREACT
        WRITE(10,130) I , (ACOEFF(I,J),AREACT(I,J),SYMBOL(I,J),J=1,6) ,
     &                RCOLUMN(I)
      ENDDO


100   FORMAT('SUMMARY FILE FOR CHEM1D',//,'INDEX  SPECIES',4X,
     &       'PRODUCTION',3X,'   LOSS',12X,'FLUX(B)',9X,'FLUX(T)',//)
110   FORMAT(I3,4X,A8,2X,4(1PE10.3,6X))
120   FORMAT(//,'COLUMN REACTION RATES',/)
130   FORMAT(I3,2X,6(A1,A8,A1),2X,1PE9.3)

C************************************************************************

C	Now get Column Rates and Fluxes for only the 50-80 km region

      LGP = 26
      IGP = 41


C CALCULATE COLUMN PRODUCTION/LOSS RATES

      DO M = 1 , NTRANS
        MVARY = MTRANS(M)
        RP = 0.0
        RL = 0.0
        DO K = LGP , IGP
          RP = RP + RPROD(MVARY,K)
          RL = RL + RLOSS(MVARY,K)
        ENDDO
        RLCOLMN(MVARY) = RL*deltax
        RPCOLMN(MVARY) = RP*deltax
      ENDDO

C CALCULATE COLUMN REACTION RATES

      DO I = 1 , NREACT
        R = 0.0
        DO K = LGP , IGP
          R = R + RATE(I,K)
        ENDDO
        RCOLUMN(I) = R*deltax
      ENDDO

C PRODUCE SUMMARY FILE

      WRITE(10,140)
140   FORMAT(//,'SUMMARY FOR 50-80 KM REGION',//,'INDEX  SPECIES',4X,
     &       'PRODUCTION',3X,'   LOSS',12X,'FLUX(B)',9X,'FLUX(T)',//)

C      WRITE(10,100)
      DO M = 1 , NTRANS
        MVARY = MTRANS(M)
        WRITE(10,110) MVARY , AMOL(MVARY) , RPCOLMN(MVARY) ,
     &                RLCOLMN(MVARY) , BNDRYFLX(M,LGP) ,
     &                BNDRYFLX(M,IGP)
      ENDDO

      WRITE(10,150)
      DO I = 1 , NREACT
        WRITE(10,130) I , (ACOEFF(I,J),AREACT(I,J),SYMBOL(I,J),J=1,6) ,
     &                RCOLUMN(I)
      ENDDO

150   FORMAT(//,'50-80KM COLUMN REACTION RATES',/)


C***********************************************************************

C	Now get Column Rates and Fluxes for only the 80-130 km region

      LGP = 41
      IGP = 66


C CALCULATE COLUMN PRODUCTION/LOSS RATES

      DO M = 1 , NTRANS
        MVARY = MTRANS(M)
        RP = 0.0
        RL = 0.0
        DO K = LGP , IGP
          RP = RP + RPROD(MVARY,K)
          RL = RL + RLOSS(MVARY,K)
        ENDDO
        RLCOLMN(MVARY) = RL*deltax
        RPCOLMN(MVARY) = RP*deltax
      ENDDO

C CALCULATE COLUMN REACTION RATES

      DO I = 1 , NREACT
        R = 0.0
        DO K = LGP , IGP
          R = R + RATE(I,K)
        ENDDO
        RCOLUMN(I) = R*deltax
      ENDDO

C PRODUCE SUMMARY FILE

      WRITE(10,141)
141   FORMAT(//,'SUMMARY FOR 80-130 KM REGION',//,'INDEX  SPECIES',4X,
     &       'PRODUCTION',3X,'   LOSS',12X,'FLUX(B)',9X,'FLUX(T)',//)

C      WRITE(10,100)
      DO M = 1 , NTRANS
        MVARY = MTRANS(M)
        WRITE(10,110) MVARY , AMOL(MVARY) , RPCOLMN(MVARY) ,
     &                RLCOLMN(MVARY) , BNDRYFLX(M,LGP) ,
     &                BNDRYFLX(M,IGP)
      ENDDO

      WRITE(10,151)
      DO I = 1 , NREACT
        WRITE(10,130) I , (ACOEFF(I,J),AREACT(I,J),SYMBOL(I,J),J=1,6) ,
     &                RCOLUMN(I)
      ENDDO

151   FORMAT(//,'80-130KM COLUMN REACTION RATES',/)

C---------------------------------------------------------------------

      WRITE (10,209)
209   FORMAT(//,2X,'FLUX OF H2O VS ALTITUDE',/,2X,'ALTITUDE   FLUX',/)
      ZALT = 0
      DO I = 1 , NALT
        WRITE (10,201) ZALT, BNDRYFLX(17,I)
        ZALT = ZALT + 2.
      ENDDO
201   FORMAT(1F6.2,1PE10.2)

C***********************************************************************

      RETURN
      END