#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 4 ]
then
  echo " "
  echo "             *****     SELVAR     *****      "
  echo " "
  echo "       The program selects a variable and level"
  echo " "
  echo "  Syntax: selvar <file1> <file2> <var_code> <var_level>"
  echo " "
  exit 0
else
 echo " "
 echo "             *****     SELVAR     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "Variable's code="$3
 echo "Variable's level="$4
 echo " "
 cp $1 fort.1
fi
cat > selvar.f <<EOF
      PROGRAM PSELVAR
      IMPLICIT REAL (A-H,O-Z)

      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL SELVAR (NDIM)
      STOP 
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
C
      INTEGER I1,I2,I3
      REAL F

      REWIND NCHAN
      READ(NCHAN,END=10) ITIME,INAME,ILEV,ISIZE
      REWIND NCHAN

      NSET=0
      NDIM=0  
    1 READ (NCHAN,END=100) I1,I2,I3,I4
      NDIM=MAX0(NDIM,I4)
      READ (NCHAN,END=101) F
      NSET=NSET+1
      GOTO1
      
  101 PRINT *,' MORE HEADER THEN DATA '     
  100 PRINT 1000,NCHAN,NDIM,NSET
      REWIND NCHAN
      RETURN    
   10 PRINT *,'INPUT FILE ON TAPE',NCHAN,' IS EMPTY'
      STOP 'MISTAKE'
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5)
      END

C*****************************************************************
      SUBROUTINE SELVAR (IE)
      IMPLICIT REAL (A-H,O-Z)
C
      PARAMETER (UNDEF=0.9E+10,UNDEFE=.9*UNDEF)
      LOGICAL TVAR,TLEV
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 DELETE THE
C FOLLOWING TWO LINES
C      PARAMETER(IE77=9E12)
C      DIMENSION X(IE), ID(4)
C      POINTER (NPX,X)
C      ALLOCATE (NPX,IE)
C AND ACTIVATE THE FOLLOWING TWO LINES
      PARAMETER(IE77=20000)
      DIMENSION X(IE77),ID(4)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
      DATA TVAR,TLEV/.FALSE.,.FALSE./
      DATA ONE/1.0/
      PRINT *,' ENTER CODE FOR VARIABLE AND LEVEL: IVAR LEVEL'
      READ(5,*,END=12) IVAR,LEVEL
      PRINT *,' VARIABLE:',IVAR,', LEVEL:', LEVEL
      IF (IVAR.EQ.-9999) TVAR=.TRUE.
      IF (LEVEL.EQ.-9999) TLEV=.TRUE.
C
C*  CHECK THE INPUT FILE.
C
      REWIND 1
      READ(1,END=10) ID
      IF(ID(4).GT.IE .OR. ID(4).GT.IE77) GOTO 11
      IEE=ID(4)
      REWIND 1
      WRITE(6,6010) ' FIRST RECORD:',ID
C
C*  MAIN LOOP.
C
      ICNT=0
      IOUT = 0
  100 READ(1,END=9920) ID
        IF(ID(4).GT.IE) GOTO 11

        IF(ID(2).EQ.IVAR .AND. TLEV) THEN
          READ(1,END=15) (X(I),I=1,ID(4))
          ICNT=ICNT+1
          WRITE(2) ID
          WRITE(2) (X(I),I=1,ID(4))
          IOUT = IOUT + 1
        ELSE IF(TVAR .AND. ID(3).EQ.LEVEL) THEN
          READ(1,END=15) (X(I),I=1,ID(4))
          ICNT=ICNT+1
          WRITE(2) ID
          WRITE(2) (X(I),I=1,ID(4))
          IOUT = IOUT + 1
        ELSE IF(ID(2).EQ.IVAR .AND. ID(3).EQ.LEVEL) THEN
          READ(1,END=15) (X(I),I=1,ID(4))
          ICNT=ICNT+1
          WRITE(2) ID
          WRITE(2) (X(I),I=1,ID(4))
          IOUT = IOUT + 1
        ELSE
          READ(1,END=15) 
          ICNT=ICNT+1
        ENDIF
        GOTO 100

C**** END OF MAIN LOOP

   10 WRITE(6,6010) ' INPUT FILE IS EMPTY '
         GOTO 9910
   11 WRITE(6,6010) ' INPUT GRID TOO LARGE ', ID(4)
         GOTO 9910
   12 WRITE(6,6010) ' PARAMETER CARD TOO SHORT OR MISSING '
         GOTO 9910
   15 WRITE(6,6010) ' UNEXPECTED *EOF* DURING READ '
         GOTO 9910
 9910 CONTINUE
C
C*  ABNORMAL TERMINATION.
C
         PRINT *,  '*SELVAR*: FATAL END '
         CALL MINONE(ONE,XNEGO)
         ZZ=SQRT(XNEGO)
         STOP
C
C*  NORMAL TERMINATION.
C
 9920 WRITE(6,6010) ' LAST RECORD:',ID
      WRITE(6,6060) ICNT, IOUT
      PRINT *,  '*SELVAR*: NORMAL '
          STOP
 6010 FORMAT(' *SELVAR*: ', A,/,(5I10))
 6060 FORMAT(' *SELVAR* -- READ ',I10,' RECORDS.',/,
     &       '          WRITTEN ',I10,' RECORDS')
      END
      SUBROUTINE MINONE(XIN,XOUT)
        XOUT=-XIN
      RETURN
      END
EOF
f77 selvar.f -o selvar.x
selvar.x $1 $2 << M
$3
$4
M
cp fort.2 $2
rm fort.1 fort.2
rm selvar.x selvar.f
exit