#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 2 ]
then
  echo " "
  echo "            *****     AVG     *****      "
  echo " "
  echo "   The program calculates the average field (file2)"
  echo "      over all record in the input file (file1)"
  echo " "
  echo "           Syntax: avg <file1> <file2>"
  echo " "
  exit 0
else
 echo " "
 echo "            *****     AVG     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 cp $1 fort.1
fi
cat > avg.f <<EOF
      PROGRAM PAVG
      IMPLICIT REAL (A-H,O-Z)
      
      PRINT *,'0**** AVG (,OUTPUT,TAPE1,TAPE2)****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL AVG(NDIM)     
      STOP 
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
C
      INTEGER I1,I2,I3
      REAL F

      NSET=0

      REWIND NCHAN
      READ(NCHAN,END=10) ITIME,INAME,ILEV,ISIZE
      GOTO 20
   10 PRINT *,'INPUT FILE ON TAPE',NCHAN,' IS EMPTY'
      STOP 'MISTAKE'
   20 REWIND NCHAN

      NDIM=0  
    1 READ (NCHAN,END=100) I1,I2,I3,I4
      NDIM=MAX0(NDIM,I4)
      READ (NCHAN,END=101) 
      NSET=NSET+1
      GOTO1
      
  101 PRINT *,' MORE HEADER THEN DATA '     
  100 PRINT 1000,NCHAN,NDIM,NSET
      RETURN    
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5)
      END
C*****************************************************************
      SUBROUTINE AVG(IE)
      IMPLICIT REAL (A-H,O-Z)
C
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 DELETE THE
C THE FOLLOWING LINES
C      PARAMETER(IE77=9E12, UNDEF=UVALUE)
C      DIMENSION X(IE),XBAR(IE),INUM(IE)
C      POINTER (NX,X),(NXBAR,XBAR),(NINUM,INUM)
C      ALLOCATE (NX,IE)
C      ALLOCATE (NXBAR,IE)
C      ALLOCATE (NINUM,IE)
C AND ACTIVATE THE FOLLOWING TWO LINES
       PARAMETER(IE77=20000, UNDEF=0.9E+10)
       DIMENSION X(IE77),XBAR(IE77),INUM(IE77)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      DATA ONE/1.0/
      UNDEFE = 0.99*UNDEF
C
C*  CHECK THE INPUT FILE.
C
      ICNT=0
      REWIND 1
      READ(1,END=10) ITIME,INAME,ILEV,ISIZE
      REWIND 1
      WRITE(6,6020) ITIME,INAME,ILEV,ISIZE
      IF(ISIZE.GT.IE77)GOTO 20
      MAX=IE
C
C*  MAIN LOOP -- FIRST INITIALIZE THE ACCUMULATION VECTOR.
C
      DO 30 J=1,ISIZE
         XBAR(J)=0.0
         INUM(J)=0
   30 CONTINUE
  100 CONTINUE
         READ(1,END=500) ITIME,INAME,ILEV,ISIZE
         ICNT=ICNT+1
         IF(ISIZE.GT.MAX) GOTO 21
         READ(1,END=40) (X(J),J=1,ISIZE)
         DO 150 J=1,ISIZE
           IF(X(J).LT.UNDEFE)THEN
            XBAR(J)=XBAR(J)+X(J)
            INUM(J)=INUM(J)+1
           ENDIF
  150    CONTINUE
      GOTO 100
C
C*  SCALE THE ACCUMULATED SUM AND OUTPUT.
C
  500 CONTINUE
      DO 600 J=1,ISIZE
        IF (INUM(J).GT.0.1)THEN
         XBAR(J)=XBAR(J)/INUM(J)
        ELSE
         XBAR(J)=UNDEF
        ENDIF
  600 CONTINUE
      WRITE(6,6021) ITIME,INAME,ILEV,ISIZE
      ITIME=ICNT
      WRITE(2) ITIME,INAME,ILEV,ISIZE
      WRITE(2) (XBAR(J),J=1,ISIZE)
      WRITE(6,6040) ICNT
C
C*  NORMAL TERMINATION.
C
         STOP 'O.K.'
   10 WRITE(6,6010)
      STOP 'MISTAKE'
   20 WRITE(6,6025)
      STOP 'MISTAKE'
   21 WRITE(6,6026) ICNT
      STOP 'MISTAKE'
   40 WRITE(6,6030) ICNT
      STOP 'MISTAKE'
 6010 FORMAT(' *AVG* -- THE INPUT FILE IS EMPTY.')
 6020 FORMAT(' *AVG* -- LABEL OF FIRST RECORD:', 5I10)
 6021 FORMAT(' *AVG* -- LABEL OF LAST RECORD: ', 5I10)
 6025 FORMAT(' *AVG* -- THE INPUT GRID IS TOO LARGE.')
 6026 FORMAT(' *AVG* -- AT RECORD ',I5,' IS THE INPUT ',
     + 'GRID TOO LARGE.')
 6030 FORMAT(' *AVG* -- UNEXPECTED EOF READING RECORD ',I10)
 6040 FORMAT(' *AVG* -- AVERAGED ',I10,' RECORDS.')
      END
EOF
f77 avg.f -o avg.x
avg.x $1 $2
cp fort.2 $2
rm fort.1 fort.2
rm avg.x avg.f
exit