#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 3 ]
then
  echo " "
  echo "            *****     ANOMAL     *****      "
  echo " "
  echo "       The program calculates anomalies (file3) "
  echo "from monthly data (file1) and long term means (file2)"
  echo " "
  echo " Syntax: anomal <file1> <file2> <file3>"
  echo " "
  exit 0
else
 echo " "
 echo "            *****     ANOMAL     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "file3="$3
 cp $1 fort.1
 cp $2 fort.2
fi
cat > anomal.f <<EOF
      PROGRAM Panomal
      IMPLICIT REAL (A-H,O-Z)
      LOGICAL ANOMA

      ANOMA=.TRUE.

      IF (ANOMA) THEN
        WRITE(6,'(A)') '0*ANOMAL(,,ACT,MEAN,ANOMALIES)'
      ELSE
        WRITE(6,'(A)') '0*RECON(,,ANOMALIES,MEAN,ACT.)'
      ENDIF
      PRINT *,'  '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      NCHAN=2
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL ANOMAL (NDIM,ANOMA)      
      STOP 
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
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 ANOMAL(IE,ANOMA)
      IMPLICIT REAL (A-H,O-Z)
C
      PARAMETER (UNDEF=0.9E+10, UNDEFE=0.9*UNDEF)
      DIMENSION ID(4)
      LOGICAL ANOMA
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 ACTIVATE THE
C FOLLOWING TWO LINES IN THIS PROGRAM !
      PARAMETER (IE77=20000) 
      DIMENSION X(IE77),XM(IE77,12)
C AND DELETE THE FOLLOWING TWO LINES
C      PARAMETER (IE77=65000) 
C      DIMENSION X(IE), XM(IE,12)
C      POINTER (NPX,X),(NPY,XM)
      NMO(II) = (II+.1 - (IFIX(II/10000.+0.1)*10000.))/100
C      ALLOCATE (NPX,IE)
C      ALLOCATE (NPY,IE*12)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      DATA ONE /1./
      DO 1 I = 1,IE
        DO 1 M = 1,12
    1   XM(I,M) = UNDEF
C
C*  CHECK IF FILES EMPTY AND WHETHER DIMENSIONS MATCH.
C
      REWIND 1
      REWIND 2
      READ(1,END=100) ID
      IE2 = ID(2)
      IE3 = ID(3)
      IE4 = ID(4)
      READ(2,END=200) ID
      IF(IE4.GT. IE) GOTO 101
      IF(ID(4).GT. IE) GOTO 201
      IF(IE2.NE.ID(2)) GOTO 300
      IF(IE3.NE.ID(3)) GOTO 301
      IF(IE4.NE.ID(4)) GOTO 302
C
C*  READ DATA FROM <TAPE2> UNTIL **EOF**
C   (LONG TERM MEANS).
C
      REWIND 2
   10 CONTINUE
      READ(2,END=20) ID
      READ(2) (X(I),I=1,ID(4))
C
C*  FIND NUMBER OF MONTH ID(1) = YYMMDD.
C
      MM = NMO(ID(1))
      IF(MM.LT.1 .OR. MM.GT.12) GOTO 102
      DO 11 I = 1,IE4
        IF(X(I) .LT. UNDEFE) XM(I,MM) = X(I)
   11   CONTINUE
      GOTO 10
C
C*  READ OF ACTUAL DATA FROM <TAPE1> UNTIL *EOF*.
C
   20 REWIND 1
   23 READ(1,END=9920) ID
      READ(1) (X(I),I=1,ID(4))
C
C*  FIND NUMBER OF MONTH ID(1) = YYMMDD.
C
      MM = NMO(ID(1))
      IF(MM.LT.1 .OR. MM.GT.12) GOTO 202
      IF (ANOMA) THEN
        DO 21 I = 1,IE4
          IF(X(I) .LT. UNDEFE .AND. XM(I,MM).LT.UNDEFE) THEN
            X(I) = X(I) - XM(I,MM)
          ELSE
            X(I) = UNDEF
          ENDIF
   21   CONTINUE
      ELSE
        DO 22 I = 1,IE4
          IF(X(I) .LT. UNDEFE .AND. XM(I,MM).LT.UNDEFE) THEN
            X(I) = X(I) + XM(I,MM)
          ELSE
            X(I) = UNDEF
          ENDIF
   22   CONTINUE
      ENDIF
      WRITE(3) ID
      WRITE(3) (X(I),I=1,ID(4))
      GOTO 23
C
C*  ABNORMAL TERMINATION.
C
  100 WRITE(6,6010) '<TAPE1>-FILE EMPTY'
         GOTO 9910
  101 WRITE(6,6010) '<TAPE1> DIMENSION TOO LARGE >',IE77
         GOTO 9910
  102 WRITE(6,6010) 'NUMBER OF MONTH WRONG', ID, MM
         GOTO 9910
  200 WRITE(6,6010) '<TAPE2> FILE EMPTY'
         GOTO 9910
  201 WRITE(6,6010) '<TAPE2> DIMENSION TOO LARGE >',IE77
         GOTO 9910
  202 WRITE(6,6010) 'NUMBER OF MONTH WRONG', ID, MM
         GOTO 9910
  300 WRITE(6,6010) '<TAPE2>/<TAPE1> INCOMPATBILE (ID(2))',IE2,ID(2)
         GOTO 9910
  301 WRITE(6,6010) '<TAPE2>/<TAPE1> INCOMPATBILE (ID(3))',IE3,ID(3)
         GOTO 9910
  302 WRITE(6,6010) '<TAPE2>/<TAPE1> INCOMPATBILE (ID(4))',IE4,ID(4)
         GOTO 9910
 9910 CONTINUE
C
C*  CAUSE AN ABORT.
C
      PRINT *,'*ANOMAL*: FATAL END '
      CALL MINONE(ONE,XNEGO)
      Z=SQRT(XNEGO)
      STOP 
C
C*  NORMAL TERMINATION.
C
 9920 PRINT *,'*ANOMAL*: O.K.'
      STOP
 6010 FORMAT(' *ANOMAL* :',1X,A,4I10,/,1X,4I10)
      END
      SUBROUTINE MINONE(XIN,XOUT)
        XOUT=-XIN
      RETURN
      END
EOF
f77 anomal.f -o anomal.x
anomal.x $1 $2 $3
cp fort.3 $3
rm fort.1 fort.2 fort.3
rm anomal.x anomal.f
exit