#! /bin/sh
#
#    Last change: Mihai Dima 2001
#
if [ $# -lt 3 ]
then
  echo " "
  echo "            *****     BOXAVG     *****      "
  echo " "
  echo " The program averages over areas of given dimensions"
  echo " "
  echo "Syntax: boxavg <file1> <file2> <int_1> <int_2> <int_3>"
  echo " "
  echo "    <int_1> - lines number in input file "
  echo "    <int_2> - lines number of the averaging box"
  echo "    <int_3> - columns number of the averaging box"
  echo " "
  exit 0
else
 echo " "
 echo "            *****     BOXAVG     *****      "
 echo " "
 echo "file1="$1
 echo "file2="$2
 echo "Lines number in the input file       "$3
 echo "Lines number of the averaging box    "$4
 echo "Columns number of the averaging box  "$5
 echo " "
 cp $1 fort.1
fi
cat > boxavg.f <<EOF
      PROGRAM PSELBOX
      IMPLICIT REAL (A-H,O-Z)
      PRINT *,'0**** BOXAVG (INPUT,OUTPUT,TAPE1,TAPE2)****'
      PRINT *,'   '
      NCHAN=1
      CALL PARLES (NCHAN,NDIM,NSET)
      CALL BOXAVG (NDIM)     
      STOP 
      END
C******************************************************************
      SUBROUTINE PARLES (NCHAN,NDIM,NSET)
      IMPLICIT REAL (A-H,O-Z)
C
      INTEGER I1,I2,I3

      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) 
      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 BOXAVG (IE)
      IMPLICIT REAL (A-H,O-Z)
C
C****  *BOXAVG* - AVERAGE SPATIAL FIELD OVER COARSER GRID.
C
C      PARAMETER (TMISS=0.9E+10,UNDEFE=.9*TMISS)
      INTEGER ROWLEN,COLLEN
      REAL REMAIN,TOT
      INTEGER XPOS,YPOS,XIND,NCOUNT
	
C IF YOU USE THIS PROGRAM NOT ON A CRAY-2 DELETE THE
C FOLLOWING LINES
C      PARAMETER(IE77=9E12)
C      DIMENSION X(IE),Y(IE)
C      POINTER (NPX,X),(NPY,Y)
C      ALLOCATE (NPX,IE)
C      ALLOCATE (NPY,IE)
C AND ACTIVATE THE FOLLOWING LINES
      PARAMETER(IE77=20000,IE78=20000,TMISS=0.9E+10,UNDEFE=0.9*TMISS)
      DIMENSION X(IE77),Y(IE78)
C NOW YOU HAVE DONE ALL CHANGES FOR STANDARD FORTRAN.
C X IS FOR DATA READ IN FROM INPUT FILE
C Y IS FOR AVERAGED DATA WRITTEN TO OUTPUT FILE
C 
C
C* GET THE BOXAVG CARD.
C
C      PRINT *,'*BOXAVG* ENTER: NROWS, BOXROWS, BOXCOLS'
C      PRINT *,'  FOR EXAMPLE:   90     2     5  IF COADS 2X2 TO 4X10'
C      READ(5, *,END=12) COLLEN,NROWS,NCOLS
      READ(*,*)COLLEN,NROWS,NCOLS
C
C LEN    = TOTAL LENGTH OF RECORD
C ROWLEN = NR DATA IN A ROW
C (ROWLEN MUST GO INTO LEN)
C COLLEN = LEN / ROWLEN
C NROWS  = NR ROWS TO BE AVERAGED
C (NROWS MUST GO INTO COLLEN)
C NCOLS  = NR COLS TO BE AVERAGED
C (NCOLS MUST GO INTO ROWLEN)
C NEWLEN = TOTAL LENGTH OF NEW RECORD (=LEN/(NROWS*NCOLS))
C
C  CHECK THE INPUT FILE.
C
      REWIND 1
      READ(1,END=10) ITIME,INAME,ILEV,LEN
      IF(LEN.GT.IE) GOTO 11
      IE=LEN
      REWIND 1
      WRITE(6,6008) ' FIRST RECORD:',ITIME,INAME,ILEV,LEN
C
C *******************************************************************
C *** LOOP TO READ EACH RECORD ***
C ********************************
C ICNT IS THE RECORD COUNTER VARIABLE
	ICNT=1
 980	CONTINUE
C
C READ A RECORD
  100 READ(1,END=9920) ITIME,INAME,ILEV,LEN
C                                 NORMAL BREAK OUT OF LOOP HERE -->
      IF(LEN.GT.IE) GOTO 14
      READ(1,END=15) (X(I),I=1,IE)
C PRINT SAMPLE RECORD
C	IF (ICNT.EQ.2.OR.ICNT.EQ.7) PRINT 107,' RECORD NR ',ICNT,
C     * ' 1ST DATA 5X2 BLOCKS ARE: ',( (X(I+J*NCOLS),I=1,5),
C     * (X(I+J*NCOLS),I=ROWLEN+1,ROWLEN+5), J=0,ROWLEN/NCOLS-1 )
C
C CHECK THAT VARIOUS RATIOS ARE INTEGERS
	REMAIN = MOD(LEN,COLLEN)
	IF (REMAIN.NE.0.) THEN
          WRITE(6,6010) 
     *          'NUMBER OF ROWS IS NOT EXACT FRACTION OF TOTAL DATA'
     *          ,COLLEN,LEN
		GOTO 9910
	ENDIF
	ROWLEN=LEN/COLLEN
	REMAIN = MOD(COLLEN,NROWS)
	IF (REMAIN.NE.0.) THEN
	WRITE(6,6010)
     * 'NR ROWS TO AVERAGE IS NOT EXACT FRACTN OF TOTAL ROWS'
     * ,NROWS,COLLEN
		GOTO 9910
	ENDIF
	REMAIN = MOD(ROWLEN,NCOLS)
	IF (REMAIN.NE.0.) THEN
	WRITE (6,6010)
     * 'NR COLS TO AVERAGE IS NOT EXACT FRACTN OF TOTAL COLS'
     * ,NCOLS,ROWLEN
		GOTO 9910
	ENDIF
 99	CONTINUE
	NEWLEN = LEN / (NROWS*NCOLS)
C
C WRITE HEADER OF NEW RECORD
        WRITE(2) ITIME,INAME,ILEV,NEWLEN
C
C ***************************************************
C --- LOOP TO DO AVERAGING ---
C ----------------------------
C SETUP
C XPOSA IS CURRENT START OF A BLOCK OF ROWS:
	XPOSA=1
C YPOS IS CURRENT Y INDEX POSITION:
	YPOS=1
C XPOS IS CURRENT X INDEX POSITION:
	XPOS=XPOSA
C XIND IS CURRENT X-VALUE BEING USED
C TOT IS CURRENT TOTAL OF VALUES BEING SUMMED
C NCOUNT IS NR OF GOOD VALUES SO FAR SUMMED
C
C *************************
C LOOP TO DO EACH NEW VALUE
 960	CONTINUE
	TOT=0.
	NCOUNT=0
	DO 220 JROW=1,NROWS
	XIND = XPOS + (JROW-1)*ROWLEN
	DO 210 JCOL=1,NCOLS
	IF (X(XIND).LT.UNDEFE) THEN
		TOT=TOT+X(XIND)
		NCOUNT=NCOUNT+1
	ENDIF
	XIND = XIND+1
 210	CONTINUE
 220	CONTINUE
C
	IF (NCOUNT.EQ.0) THEN
		Y(YPOS)=TMISS
	ELSE
		Y(YPOS)=TOT/NCOUNT
	ENDIF
C
	XPOS = XPOS+NCOLS
	YPOS=YPOS+1
C IF SAME ROW --
	IF (XPOS.LT.XPOSA+ROWLEN) GOTO 960
C IF NEW ROW --
	XPOSA = XPOSA+(ROWLEN*NROWS)
	XPOS  = XPOSA
C IF STILL WITHIN SET --
	IF (XPOSA.LT.LEN) GOTO 960
C   BREAK OUT OF LOOP HERE WHEN LEN IS REACHED -->
C
C END OF LOOP TO DO EACH NEW VALUE
C ********************************
C
C ---------------------------
C END OF LOOP TO DO AVERAGING
C *****************************************************
C
C WRITE NEW DATA
	WRITE(2) (Y(J),J=1,NEWLEN)
C PRINT SAMPLE RECORD
C	IF (ICNT.EQ.2.OR.ICNT.EQ.7) THEN
C	PRINT 107,' RECORD NR ',ICNT,
C     * ' NEW DATA 1ST 5 ITEMS ARE: ',(Y(I),I=1,5)
 107	FORMAT(A,I4,A/(5F8.2/5F8.2/))
C	ENDIF
C
      ICNT=ICNT+1
	GOTO 980
C *******************************************
C ***** END OF LOOP TO READ EACH RECORD *****
C *******************************************************************
C
C FATAL ERRORS   
   10	WRITE(6,6010) ' INPUT FILE IS EMPTY '
         GOTO 9910
   11 WRITE(6,6010) ' INPUT GRIDS TOO LARGE '
         GOTO 9910
   12 WRITE(6,6010) ' PARAMETER CARD TOO SHORT OR MISSING '
         GOTO 9910
   14 WRITE(6,6010) ' INVALID RECORD LENGTH:', LEN
         GOTO 9910
   15 WRITE(6,6010) ' UNEXPECTED *EOF* DURING READ '
         GOTO 9910
C
C*  ABNORMAL TERMINATION.
 9910 CONTINUE
         PRINT *, '*BOXAVG*: FATAL END '
         STOP
C
C*  NORMAL TERMINATION.
C
 9920 WRITE(6,6008) ' LAST RECORD READ:',ITIME,INAME,ILEV,LEN
      WRITE(6,6060) ICNT-1
	PRINT 6008,' LAST RECORD WRITTEN',ITIME,INAME,ILEV,NEWLEN
      PRINT *, '*BOXAVG*: NORMAL END OF PROGRAM'
      stop
 6008 FORMAT('*BOXAVG*: ', A,/,(6I10))
 6010 FORMAT('*******************************************************'/
     * '*BOXAVG*: ', A,/,(8X,I10,26X,I10/
     * '************************************************************')) 
 6011 FORMAT(' *BOXAVG*: ', A)
 6060 FORMAT(' READ ',I10,' RECORDS.')
      END
EOF
f77 boxavg.f -o boxavg.x
boxavg.x $1 $2 $3 $4 $5<< M
$3
$4
$5
M
cp fort.2 $2
rm fort.1 fort.2
rm boxavg.x boxavg.f
exit