depth_check.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	depth_check

Function/Subroutine Documentation

depth_check()

subroutine depth_check

(

)

Definition at line 47 of file depth_check.f90.

  !==============================================================================|
  USE all_vars
  USE mod_utils
  USE mod_par  
  IMPLICIT NONE
  INTEGER, DIMENSION(NPROCS) :: SBUF,RBUF
  INTEGER  :: I,II,MLOC,IERR, RECV
  REAL(SP) :: DMIN
  !==============================================================================|

   IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "Start: depth_check"
!--Calculate Minimum Depth and Set Global Node Number if Min Depth < MIN_DEPTH
   sbuf = 0
   mloc = 0
   ierr = 0
   dmin = minval(d(1:m))
   ! NGID NOW EXISTS FOR BOTH SERIAL AND PARALLEL
   IF(dmin < min_depth) mloc = ngid(minloc(d(1:m),dim=1))

!--Reduce in Master Processor Array and Dump To Screen
   sbuf(myid) = mloc
   rbuf = sbuf


!--If Depth Condition is Violated Write Warning and Halt User
   IF(msr)THEN
     IF(sum(rbuf) /= 0) THEN
        DO i=1,nprocs
           ii = rbuf(i)
           IF(ii /= 0)THEN
              !         WRITE(*,*)'DEPTH IN NODE: ',II,' AT ',XG(II)+VXMIN,YG(II)+VYMIN, &
              !                        ' IS LESS THAN MIN_DEPTH'
              
              WRITE(*,*)'DEPTH IN NODE: ',ii,'; IS LESS THAN MIN_DEPTH'
              WRITE(*,*)'ADJUST BATHYMETRY AT THIS (THESE) LOCATION(S) OR'
              WRITE(*,*)'RECOMPILE FVCOM WITH FLOODING/DRYING FORMULATION'                 
              WRITE(*,*)'STOPPING....'
              CALL pstop
           END IF
        END DO
     END IF
  END IF

   IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "End: depth_check"

Here is the call graph for this function: