#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 4 ]
then
  echo " "
  echo "               *****     SELMON     *****      "
  echo " "
  echo "         The program selects a month or season"
  echo " "
  echo " Syntax: selmon <file1> <file2> <month_first> <month_last>"
  echo " "
  exit 0
else
 echo " "
 echo "               *****     SELMON     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "First month="$3
 echo "Last month="$4
 echo " "
 cp $1 fort.1
fi
cat > selmon.f <<EOF
      PROGRAM PSELMON
      IMPLICIT REAL (A-H,O-Z)
      
      PRINT *,'0*SELMON(INOPUT,OUTPUT,TAPE1,TAPE2)*'
      PRINT *,'  '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL SELMON (NDIM)       
      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
      RETURN    
 1000 FORMAT (1X,' TAPE=',I2,',  NDIM=',I5,',  NSET=',I5,/)
C
      END
C*******************************************************************
      SUBROUTINE SELMON(IE)
      IMPLICIT REAL (A-H,O-Z)
C
C   SELMON = select data from specified months 
C
C   TAPE1:   INPUT-FILE, daily data
C
C   TAPE2:   OUTPUT-FILE selected data
C
C
C   HVS 26/7/90 
C   MODIFIED: 4.6.91 H.BORGERT  CRAY2 VARIABLE FIELDS
C   modification 10/11/91 Heiko Borgert: now running on SUN CRAY CONVEX
C
      PARAMETER (UNDEF=0.9E+10,UNDEFE=.9*UNDEF)
      INTEGER ID(3),IT,MFIRST,MLAST
      LOGICAL LOG
C IF YOU DONT USE THIS PROGRAM ON A CRAY-2 DELETE THE FOLLOWING LINE
C      REAL X(IE)
C      POINTER (NPX,X)
C      ALLOCATE (NPX,IE)
C AND ACTIVATE THE NEXT TWO LINES
      PARAMETER (IE77=20000)
      REAL X(IE77)
C NOW YOU HAVE DONE ALL CHANGES FOR FORTRAN 77
      LOG = .TRUE.
      NALL = 0
      NW = 0
C
      PRINT *, '*SELMON(TAPE1,TAPE2)*'
      PRINT *,' '
      PRINT *,' ENTER THE SELECTED INTERVAL: MFIRST MLAST'
      READ(5,*) MFIRST, MLAST
      IF(MFIRST.GT.MLAST) LOG=.FALSE.
      PRINT *, ' SELECTED MONTHS: ', MFIRST, ' TO ', MLAST
      REWIND 1
      READ(1,END=1000) IT, ID 
      IF(ID(3).GT.IE)                             GOTO 1001 
      PRINT *, 'first record: ',ID
C
      REWIND 1
    1 READ(1,END=800) IT, ID
      READ(1) (X(I),I=1,ID(3))
      NALL = NALL + 1
      CALL SPLIDAT (IT,NY,NM,ND)
      IF(      LOG .AND. (NM.LT.MFIRST  .OR. NM.GT.MLAST)) GOTO 1
      IF(.NOT. LOG .AND. (NM.LT.MFIRST .AND. NM.GT.MLAST)) GOTO 1
      NW = NW + 1
      WRITE(2) IT,ID
      WRITE(2) (X(I),I=1,ID(3))
      GOTO 1
C
C  END-OF-FILE FOUND
  800 PRINT 9900, ' NUMBER OF RECORDS FOUND : ',NALL
      PRINT 9900, ' NUMBER OF RECORDS WRITTEN : ',NW
C
      PRINT 9900, ' REGULAR END '
      GOTO 1999
C
C   IRREGULAR ENDING -----------------------------------------------
 1000 PRINT 9910, ' TAPE1 EMPTY '
                                                           GOTO 1999
 1001 PRINT 9910, ' TOO MANY DATA ON TAPE1: ', IE, IT, ID
                                                           GOTO 1999
 1999 STOP 
C
 9900 FORMAT (1X, A, 8I10)
 9910 FORMAT (/' PROGRAM SELMON **** FATAL END :',/,
     &        1X,A,/,1X,4I10/1X,4I10)
      END 
C***************************************************
      SUBROUTINE SPLIDAT (IT, IYY, MM, IDD)
      IMPLICIT REAL (A-H,O-Z)
C
C  PURPOSE: SPLIT T21-FORMAT DATE IN PIECES
C  HVS  24/7/1989
      IYY = IT/10000 + 0.1
      II  = IT - IYY*10000
      MM  = II/100 + 0.1
      IDD = II - MM*100
      RETURN
      END 
EOF
f77 selmon.f -o selmon.x
selmon.x $1 $2 << M
$3
$4
M
cp fort.2 $2
rm fort.1 fort.2
rm selmon.x selmon.f
exit