#! /bin/sh
#
#     Last change:  Mihai Dima  2001
#
if [ $# -lt 1 ]
then
  echo " "
  echo "    *****     GGSTAT     *****    "
  echo " "
  echo "  The program provides information "
  echo "      about the file1's contents "
  echo " "
  echo "      Syntax: ggstat <file1>"
  echo " "
  exit 0
else
 echo " "
 echo "    *****     GGSTAT     *****    "
 echo " "
 echo "file1="$1
 cp $1 fort.1
fi
cat > ggstat.f <<EOF
      PROGRAM  PGGSTAT

      IMPLICIT REAL (A-H,O-Z)
      
      PRINT *,'****     GGSTAT     ****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)

      CALL GGSTAT (NDIM,NSET)

      STOP 
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
C
C  SUBROUTINE FOR READING THE DIMENSION PARAMETER OF THE FIELDS ON
C
      INTEGER I1,I2,I3

      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
      REWIND NCHAN
      RETURN    
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5)
      END
C*****************************************************************
      SUBROUTINE GGSTAT (IE,NS)
      IMPLICIT REAL (A-H,O-Z)
C
      PARAMETER (UNDEF=9E9,UNDEFE=0.9*UNDEF)
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 DELETE THE
C THE FOLLOWING TWO LINES
C      PARAMETER (IE77=9E12)
C      DIMENSION Y(IE)
C      POINTER (NPY,Y)
C      ALLOCATE (NPY,IE)
C AND ACTIVATE THE FOLLOWING LINES
      PARAMETER (IE77=20000)
      DIMENSION Y(IE77)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      DATA ONE /1./
C
      REWIND 1
      READ(1,END=10) ITIME,INAME,ILEV,ISIZE
      IF(ISIZE.EQ.0) GOTO 12
      REWIND 1
C
C     * MAIN LOOP
C
      ICNT=0
      WRITE(6,6015) 
  100 CONTINUE
         READ(1,END=9920) ITIME,INAME,ILEV,NWDS
         ICNT=ICNT+1
         IF(NWDS.GT.IE)THEN
            WRITE(6,6020) ICNT
            GOTO 9910
         ENDIF
         READ(1,END=13) (Y(I),I=1,NWDS) 
         YBAR=0.0
         HELP = 0.
         NPO = 0
         NUNDEF = 0
         DO 150 I=1,NWDS
            IF(Y(I) .LT. UNDEFE) THEN
             YBAR=YBAR+Y(I)
             HELP = HELP +Y(I)*Y(I)
             NPO = NPO + 1
            ELSE
             NUNDEF = NUNDEF + 1
            ENDIF
  150    CONTINUE
         IF(NPO.NE.0) THEN
           YBAR=YBAR/NPO
           HELP = HELP/NPO - YBAR*YBAR
           YSD = 0. 
           IF(HELP.GT.0.) YSD = SQRT(HELP)
         ELSE
           YBAR = UNDEF
           YSD  = UNDEF
         ENDIF
         WRITE(6,6040)ICNT,ITIME,INAME,ILEV,NWDS,YBAR,YSD,NUNDEF 
         GOTO 100
C
C     * ABNORMAL TERMINATION
C
   10 WRITE(6,6010) ' INPUT FILE IS EMPTY '
         GOTO 9910
   12 WRITE(6,6010) ' FIELD LENGTH ZERO '
         GOTO 9910
   13 WRITE(6,6030) ICNT
         GOTO 9910
 9910 CONTINUE
C
C        * CAUSE AN ABORT
C
         PRINT *,'*GGSTAT*: FATAL END '
         CALL MINONE(ONE,XNEGO)
         Z=SQRT(XNEGO)
         STOP
C
C     * NORMAL TERMINATION
C
 9920 CONTINUE
         STOP
C
C
 6010 FORMAT( A, I5)
 6015 FORMAT('0  RECORD   IDAT     IVAR    LEVEL   LEN ',
     1       '        MEAN     STD DEV  N*UNDEF'/
     2       ' ----------------------------------------',
     3       '---------------------------------')
 6020 FORMAT('0GGSTAT - RECORD ',I10,' IS TOO LONG.')
 6030 FORMAT('0GGSTAT -- UNEXPECTED EOF READING RECORD ',I10)
 6040 FORMAT(' ',I7,1X,4I8,2E12.4,I8) 
      END 
      SUBROUTINE MINONE(XIN,XOUT)
        XOUT=-XIN
      RETURN
      END
EOF
f77 ggstat.f -o ggstat.x
ggstat.x $1
rm fort.1
rm ggstat.x ggstat.f
exit