!> clips.F
!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!> clip(fhf,f0) 
!!
      subroutine clip(fhf,f0)
!
! ... Local variables ..................................................
!
#include "param.inc"
#include "mdims.inc"
!
! ... Parameter variables ..............................................
!
      dimension fhf(LIp1,LJ)
!> Description of fhf
      dimension f0(LIp4,LJ)
!> Description of f0
!
! ... Begin ............................................................
!
#ifdef DEBUG_MODE_ON
      write(*,*) "DEBUG:  in clips.F::clip(fhf,f0)"
#endif
        do 200 j=1,MJ
        do 200 i=1,MIp1
  200   fhf(i,j) = amax1( amin1(f0(i+2,j),f0(i+1,j)),
     $          amin1( fhf(i,j), amax1(f0(i+2,j),f0(i+1,j))))
C
!
! ... End subroutine clip(fhf,f0) ......................................
!
      return
      end
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
!> CLIP2(X,X0)
      SUBROUTINE CLIP2(X,X0)
!
! ... Local variables ..................................................
!
#include "param.inc"
#include "mdims.inc"
!
! ... Parameter variables ..............................................
!
      DIMENSION X(lIP1,lJP1)
!> Description of X
      DIMENSION X0(lIP4,lJP4)
!> Description of X0
!
! ... Begin ............................................................
!
#ifdef DEBUG_MODE_ON
      write(*,*) "DEBUG:  in clips.F::CLIP2(X,X0)"
#endif
        DO 200 J=1,MJP1
        DO 200 I=1,MIP1
  200   X(I,J) = AMIN1(
     $           AMAX1(X0(I+1,J+1),X0(I+2,J+1),X0(I+1,J+2),X0(I+2,J+2)),
     $          AMAX1( X(I,J),
     $          AMIN1(X0(I+1,J+1),X0(I+2,J+1),X0(I+1,J+2),X0(I+2,J+2))))
C
!
! ... End subroutine CLIP2(X,X0) .......................................
!
        RETURN
        END
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
!> Iclip(fhf,f0)
!!
       subroutine Iclip(fhf,f0)
!
! ... Local variables ..................................................
!
#include "param.inc"
#include "mdims.inc"
!
! ... Parameter variables ..............................................
!
      dimension fhf(lIp1,lJP1)
!> description of fhf
      dimension f0(lIp4,lJP1)
!> description of f0
!
! ... Begin ............................................................
!
#ifdef DEBUG_MODE_ON
      write(*,*) "DEBUG:  in clips.F::Iclip(fhf,f0)"
#endif
        do 200 j=1,MJP1
        do 200 i=1,MIp1
  200   fhf(i,j) = amax1( amin1(f0(i+2,j),f0(i+1,j)),
     $          amin1( fhf(i,j), amax1(f0(i+2,j),f0(i+1,j))))
!
! ... End subroutine Iclip(fhf,f0) .....................................
!
        return
        end
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!> jclip(fhf,f0)
!!
      subroutine jclip(fhf,f0)
!
! ... Local variables ..................................................
!
#include "param.inc"
#include "mdims.inc"
!
! ... Parameter variables ..............................................
!
      dimension fhf(lIP1,lJp1)
!> description of fhf
      dimension f0(lIP1,lJp4)
!> description of f0
!
! ... Begin ............................................................
!
#ifdef DEBUG_MODE_ON
      write(*,*) "DEBUG:  in clips.F::jclip(fhf,f0)"
#endif
        do 200 j=1,MJp1
        do 200 i=1,MIP1
  200   fhf(i,j) = amax1( amin1(f0(i,j+2),f0(i,j+1)),
     $          amin1( fhf(i,j), amax1(f0(i,j+2),f0(i,j+1))))
!
! ... End subroutine jclip(fhf,f0) .....................................
!
        return
        end
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!