#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 6 ]
then
  echo " "
  echo "           *****     EDGE     *****      "
  echo " "
  echo "The program replaces values in the input file"
  echo " "
  echo "     ITYP set the substitution conditions:"
  echo " "
  echo "ITYP < 0 : all values smaller DG1 set to SPECIAL"
  echo "ITYP > 0 : all values greater DG1 set to SPECIAL"
  echo "ITYP = 0 : all values smaller DG1 and greater"
  echo "                              DG2 set to SPECIAL"
  echo " "
  echo "Syntax: edge <file1> <file2> <ITYP> <DG1> <DG2> <SPECIAL>"
  echo " "
  exit 0
else
 echo " "
 echo "           *****     EDGE     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "ITYPE="$3
 echo "DG1="$4
 echo "DG2="$5
 echo "NEW VALUE="$6
 cp $1 fort.1
fi
cat > edge.f <<EOF
      PROGRAM PEDGE 
      IMPLICIT REAL (A-H,O-Z)
      
      PRINT *,'0*EDGE(INPUT,OUTPUT,TAPE1,TAPE2)*'
      PRINT *,'  '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL EDGE (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=',I7,', NSET=',I7,/)
C
      END
C*******************************************************************
      SUBROUTINE EDGE (NDIM)
      IMPLICIT REAL (A-H,O-Z)
C
C****   *EDGE*  - values smaler than *DG* are set to 
C       *UVALUE* or input value (UNDEF)
C
      DIMENSION X(NDIM),ID(4)
C      POINTER (NX,X)
C      ALLOCATE(NX,NDIM)
      
      UNDEF=0.9E+10

C     PRINT *,'Enter value for ITYP and DG1, and If ITYP equal 0 DG2:'

      READ(5,*,ERR=10) ITYP,DG1,DG2
   10 PRINT *,' '
      PRINT *,'enter the value for SPECIAL, otherwise all values'
      PRINT *,'smaler/greater than *DGx* are set to ',UNDEF
      READ (5,*,ERR=11,END=11) UNDEF
      GOTO 12 
   11 UNDEF=0.9E+10
   12 PRINT *,' ITYP=',ITYP,' DG1=',DG1,' DG2=',DG2
      PRINT *,' new data values=',UNDEF
      REWIND 1

1     READ(1,END=9) ID
      READ(1) (X(I),I=1,ID(4))
      DO 2 I=1,ID(4)
      IF(ITYP.LT.0.AND.X(I).LT.DG1) X(I)=UNDEF
      IF(ITYP.GT.0.AND.X(I).GT.DG1) X(I)=UNDEF
      IF(ITYP.EQ.0.AND.(X(I).LT.DG1.OR.X(I).GT.DG2)) X(I)=UNDEF
  
2     CONTINUE
      WRITE(2) ID
      WRITE(2) (X(I),I=1,ID(4))
      GOTO 1
9     CONTINUE
      STOP
      END
EOF
f77 edge.f -o edge.x
edge.x $1 $2 << M
$3
$4
$5
$6
M
cp fort.2 $2
rm fort.1 fort.2
rm edge.x edge.f
exit