!
      PROGRAM P72
!
!     DEFINE 3 MATRICES
!     AND AUXILLIARY ARRAYS
!
      IMPLICIT NONE
      CHARACTER (LEN=3) :: SEASON(4,3)
      INTEGER :: NUMBER(4,3),NCOL(3),CMEAN(3)
      INTEGER :: J,KS,M,NTOTAL,NAVE,NGRAND,NGAVE,KOL,KROW
      REAL :: REVENU(4,3),RCOL(3),DMEAN(3),RTOTAL,RAVE,RGRAND,RGAVE
!
!
      PRINT *, 'This is Program >> P72  - Matrix summary'
!
!     Tell program where data for  READ *  is coming from
      OPEN(UNIT=5, FILE='P72.DAT')      ! UNIT=5 is the default input
!
!
!     READ IN DATA
!     12 LINES IN MONTH ORDER
!
L1:    DO  KS=1,4
L2:      DO  M=1,3
            READ 127,SEASON(KS,M),  &
            NUMBER(KS,M),REVENU(KS,M)
         END DO L2
       END DO L1
 127  FORMAT(A3,I7,F5.1)
      PRINT 1
   1  FORMAT(/' SEASONAL SUMMARY')
!
!     DO 4 QUARTER AVERAGE + TOTALS
!
      NGRAND=0
      RGRAND=0.0
L3:   DO  KS=1,4
         NTOTAL=0
         RTOTAL=0.0
L4:         DO  J=1,3
              NTOTAL=NTOTAL+NUMBER(KS,J)
              RTOTAL=RTOTAL+REVENU(KS,J)
            END DO L4
         NGRAND=NGRAND+NTOTAL
         RGRAND=RGRAND+RTOTAL
         NAVE=NTOTAL/3
         RAVE=RTOTAL/3.0
         PRINT 220,(SEASON(KS,J),J=1,3)
         PRINT 230,(NUMBER(KS,J),J=1,3),  &
                    NTOTAL,NAVE
         PRINT 240,(REVENU(KS,J),J=1,3),  &
                    RTOTAL,RAVE
     END DO L3
 220  FORMAT(/ 3('     ',A3),'   ',  &
             'TOTAL','     ','AVE')
 230  FORMAT(5I8)
 240  FORMAT('  ',5(F6.1,'M$'))
!
!     DO SAME CALCULATIONS
!     FOR EACH COLUMN
!
L5:   DO  KOL=1,3
           NCOL(KOL)=0
           RCOL(KOL)=0.0
L6:        DO  KROW=1,4
               NCOL(KOL)=NCOL(KOL)+NUMBER(KROW,KOL)
               RCOL(KOL)=RCOL(KOL)+REVENU(KROW,KOL)
           END DO L6
           CMEAN(KOL)=NCOL(KOL)/4
           DMEAN(KOL)=RCOL(KOL)/4.0
      END DO L5
      NGAVE=NGRAND/12
      RGAVE=RGRAND/12.0
!
!     OUTPUT COLUMN STATISTICS
!
      PRINT 410,NCOL,NGRAND,NGAVE
      PRINT 420,RCOL,RGRAND,RGAVE
      PRINT 430,CMEAN,DMEAN
 410  FORMAT(/ 40('-')/  &
               5I8)
 420  FORMAT('  ',5(F6.1,'M$')/)
 430  FORMAT(/3I8,'    COLUMN AVERAGE'/  &
            '  ',3(F6.1,'M$')// 40('=')//)
      STOP
      END PROGRAM P72