mod_boundschk Module Reference

Functions/Subroutines
subroutine	setup_boundschk
subroutine	boundschk

Variables
type(infovar), dimension(:), allocatable	vlist

Function/Subroutine Documentation

boundschk()

subroutine mod_boundschk::boundschk

(

)

Definition at line 146 of file mod_boundschk.f90.

    use all_vars
    use mod_utils, only : pstop
    use mod_wd
    use mod_ncdio, only : archive
    use mod_par,   only : egid,ngid
    implicit none
    type(infovar), allocatable :: vlist_global(:)
    integer  :: icnt,pt,lay,i,j,k,nviolations,printproc,ierr
    integer  :: iiint 
    real(sp) :: val
    integer, allocatable :: violations(:),tmp(:)
    character(len=20) :: vname

    
    iiint=iint
    if(.not.boundschk_on) return
    if(mod(iiint,chk_interval)/= 0) return
!    if(.not.boundschk_on .or. (mod(iint,chk_interval)/= 0)) return
    allocate(tmp(nprocs))
    tmp = 0

    !------------------------------------------------------
    !Check if variables are in user-defined bounds
    !------------------------------------------------------
    icnt = 0 

    ! vert-averaged x-velocity
    val = maxval(abs(ua(1:n)))
    if(val > veloc_mag_max)then
      icnt = icnt + 1
      pt = 1
      vname = 'vert-averaged u'
      ualoop: do i=1,n
        if(abs(ua(i)) > veloc_mag_max)then 
          pt = i
          val = abs(ua(i))
          exit ualoop
        endif
      end do ualoop 
      call set_infovar(vlist(icnt),vname,2,val,veloc_mag_max,xc(pt)+vxmin,yc(pt)+vymin,myid,pt,1)
    endif

    ! vert-averaged y-velocity
    val = maxval(abs(va(1:n)))
    if(val > veloc_mag_max)then      
      icnt = icnt + 1
      pt = 1
      vname = 'vert-averaged v'
      valoop: do i=1,n
        if(abs(va(i)) > veloc_mag_max)then
          pt = i
          val = abs(va(i))
          exit valoop
        endif
      end do valoop
      call set_infovar(vlist(icnt),vname,2,val,veloc_mag_max,xc(pt)+vxmin,yc(pt)+vymin,myid,pt,1)
    endif

    ! x-velocity
    val = maxval(abs(u(1:n,1:kbm1)))
    if(val > veloc_mag_max)then
      icnt = icnt + 1
      pt = 1
      vname = 'u'
      uloop: do k=1,kbm1
      do i=1,n
        if(abs(u(i,k)) > veloc_mag_max)then
          pt = i
          lay = k
          val = abs(u(i,k))
          exit uloop
        endif
      end do 
      end do uloop
      call set_infovar(vlist(icnt),vname,2,val,veloc_mag_max,xc(pt)+vxmin,yc(pt)+vymin,myid,pt,lay)
    endif

    !hydrographic vars
    if(.not. barotropic)then

    ! salinity - max 
    val = maxval(s1(1:m,1:kbm1))
    if(val > salt_max)then
      icnt = icnt + 1
      pt = 1
      vname = 'salinity'
      smaxloop: do k=1,kbm1
      do i=1,m
        if(s1(i,k) > salt_max)then
          pt = i
          lay = k
          val = s1(i,k)
          exit smaxloop
        endif
      end do
      end do smaxloop
      call set_infovar(vlist(icnt),vname,1,val,salt_max,vx(pt)+vxmin,vy(pt)+vymin,myid,pt,lay)
    endif

    ! salinity - min 
    val = minval(s1(1:m,1:kbm1))
    if(val < salt_min)then
      icnt = icnt + 1
      pt = 1
      vname = 'salinity'
      sminloop: do k=1,kbm1
      do i=1,m
        if(s1(i,k) < salt_min)then
          pt = i
          lay = k
          val = s1(i,k)
          exit sminloop
        endif
      end do
      end do sminloop
      call set_infovar(vlist(icnt),vname,1,val,salt_min,vx(pt)+vxmin,vy(pt)+vymin,myid,pt,lay)
    endif

    ! temp - max 
    val = maxval(t1(1:m,1:kbm1))
    if(val > temp_max)then
      icnt = icnt + 1
      pt = 1
      vname = 'temperature'
      tmaxloop: do k=1,kbm1
      do i=1,m
        if(t1(i,k) > temp_max)then
          pt = i
          lay = k
          val = t1(i,k)
          exit tmaxloop
        endif
      end do
      end do tmaxloop
      call set_infovar(vlist(icnt),vname,1,val,temp_max,vx(pt)+vxmin,vy(pt)+vymin,myid,pt,lay)
    endif

    ! temperature - min 
    val = minval(t1(1:m,1:kbm1))
    if(val < temp_min)then
      icnt = icnt + 1
      pt = 1
      vname = 'temperature'
      tminloop: do k=1,kbm1
      do i=1,m
        if(t1(i,k) < temp_min)then
          pt = i
          lay = k
          val = t1(i,k)
          exit tminloop
        endif
      end do
      end do tminloop
      call set_infovar(vlist(icnt),vname,1,val,temp_min,vx(pt)+vxmin,vy(pt)+vymin,myid,pt,lay)
    endif


    endif !.not. barotropic





    !-----------------------------------------------------------
    !Collect number of violations from each proc into vector tmp
    !-----------------------------------------------------------
    nviolations = icnt
    tmp(1) = icnt 


    !check total violations and return if none found
    if(sum(tmp) == 0)then
      deallocate(tmp)
      return 
    endif

    !look at violations, just pick lowest procid with violation to print violation
    printproc = 1
    do i=1,nprocs
      if(tmp(i) > 0)then
        printproc = i
        exit
      endif
    end do

    !write violations to screen
    if(printproc==myid)then
      write(ipt,*)'WARNING: Variable(s) have exceeded user-defined thresholds'
      do i=1,nviolations
        call print_infovar(vlist(i),ipt)
      end do 
      write(ipt,*)'ARCHIVING FRAME AND HALTING'
    endif

    !dump frame to netcdf file
    force_archive = .true.
    call archive

    !shutdown
    call pstop

    deallocate(tmp)    

Here is the call graph for this function:

Here is the caller graph for this function:

setup_boundschk()

subroutine mod_boundschk::setup_boundschk

(

)

Definition at line 118 of file mod_boundschk.f90.

    use all_vars
    integer :: nvars

    !report the bounds
    if(boundschk_on)then
      write(ipt,*)'bounds checking  :  on'
      write(ipt,*)'checking interval: ',chk_interval
      write(ipt,*)'veloc_mag_max    : ',veloc_mag_max
      write(ipt,*)'zeta_mag_max     : ',zeta_mag_max
      write(ipt,*)'temp_max         : ',temp_max
      write(ipt,*)'temp_min         : ',temp_min
      write(ipt,*)'salt_max         : ',salt_max
      write(ipt,*)'salt_min         : ',salt_min

      !allocate vlist
      nvars = 7
      allocate(vlist(nvars))

    else
      allocate(vlist(1))
    endif

Here is the caller graph for this function:

Variable Documentation

vlist

type(infovar), dimension(:), allocatable mod_boundschk::vlist

Definition at line 112 of file mod_boundschk.f90.

  type(infovar), allocatable :: vlist(:)