sed_fluxes_mod Module Reference

Functions/Subroutines
subroutine, public	sed_fluxes (ng, tile)
subroutine	sed_fluxes_tile (ng, tile, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs, nstp, nnew, hz, rmask_wet, bustrc, bvstrc, bustrw, bvstrw, bustrcwmax, bvstrcwmax, bustr, bvstr, t, ero_flux, settling_flux, bed_thick, bed, bed_frac, bed_mass, bottom)

Function/Subroutine Documentation

sed_fluxes()

subroutine, public sed_fluxes_mod::sed_fluxes	(	integer, intent(in)	ng,
		integer, intent(in)	tile )

Definition at line 34 of file sed_fluxes.F.

!***********************************************************************
!
      USE mod_param
      USE mod_forces
      USE mod_grid
      USE mod_ocean
      USE mod_sedbed
      USE mod_stepping
      USE mod_bbl
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, tile
!
!  Local variable declarations.
!
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__
!
# include "tile.h"
!
# ifdef PROFILE
      CALL wclock_on (ng, inlm, 16, __line__, myfile)
# endif
      CALL sed_fluxes_tile (ng, tile,                                   &
     &                      lbi, ubi, lbj, ubj,                         &
     &                      imins, imaxs, jmins, jmaxs,                 &
     &                      nstp(ng), nnew(ng),                         &
     &                      grid(ng) % Hz,                              &
# ifdef WET_DRY
     &                      grid(ng) % rmask_wet,                       &
# endif
# ifdef BBL_MODEL
     &                      bbl(ng) % bustrc,                           &
     &                      bbl(ng) % bvstrc,                           &
     &                      bbl(ng) % bustrw,                           &
     &                      bbl(ng) % bvstrw,                           &
     &                      bbl(ng) % bustrcwmax,                       &
     &                      bbl(ng) % bvstrcwmax,                       &
# endif
     &                      forces(ng) % bustr,                         &
     &                      forces(ng) % bvstr,                         &
     &                      ocean(ng) % t,                              &
     &                      sedbed(ng) % ero_flux,                      &
     &                      sedbed(ng) % settling_flux,                 &
# if defined SED_MORPH
     &                      sedbed(ng) % bed_thick,                     &
# endif
     &                      sedbed(ng) % bed,                           &
     &                      sedbed(ng) % bed_frac,                      &
     &                      sedbed(ng) % bed_mass,                      &
     &                      sedbed(ng) % bottom)
# ifdef PROFILE
      CALL wclock_off (ng, inlm, 16, __line__, myfile)
# endif
!
      RETURN

References mod_bbl::bbl, mod_forces::forces, mod_grid::grid, mod_param::inlm, mod_stepping::nnew, mod_stepping::nstp, mod_ocean::ocean, sed_fluxes_tile(), mod_sedbed::sedbed, wclock_off(), and wclock_on().

Referenced by sediment_mod::sediment().

Here is the call graph for this function:

Here is the caller graph for this function:

sed_fluxes_tile()

subroutine sed_fluxes_mod::sed_fluxes_tile ( integer, intent(in) ng, integer, intent(in) tile, integer, intent(in) lbi, integer, intent(in) ubi, integer, intent(in) lbj, integer, intent(in) ubj, integer, intent(in) imins, integer, intent(in) imaxs, integer, intent(in) jmins, integer, intent(in) jmaxs, integer, intent(in) nstp, integer, intent(in) nnew, real(r8), dimension(lbi:,lbj:,:), intent(in) hz, real(r8), dimension(lbi:,lbj:), intent(in) rmask_wet, real(r8), dimension(lbi:,lbj:), intent(in) bustrc, real(r8), dimension(lbi:,lbj:), intent(in) bvstrc, real(r8), dimension(lbi:,lbj:), intent(in) bustrw, real(r8), dimension(lbi:,lbj:), intent(in) bvstrw, real(r8), dimension(lbi:,lbj:), intent(in) bustrcwmax, real(r8), dimension(lbi:,lbj:), intent(in) bvstrcwmax, real(r8), dimension(lbi:,lbj:), intent(in) bustr, real(r8), dimension(lbi:,lbj:), intent(in) bvstr, real(r8), dimension(lbi:,lbj:,:,:,:), intent(inout) t, real(r8), dimension(lbi:,lbj:,:), intent(inout) ero_flux, real(r8), dimension(lbi:,lbj:,:), intent(inout) settling_flux, real(r8), dimension(lbi:,lbj:,:), intent(inout) bed_thick, real(r8), dimension(lbi:,lbj:,:,:), intent(inout) bed, real(r8), dimension(lbi:,lbj:,:,:), intent(inout) bed_frac, real(r8), dimension(lbi:,lbj:,:,:,:), intent(inout) bed_mass, real(r8), dimension(lbi:,lbj:,:), intent(inout) bottom )

private

Definition at line 95 of file sed_fluxes.F.

!***********************************************************************
!
      USE mod_param
      USE mod_scalars
      USE mod_sediment
!
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, tile
      integer, intent(in) :: LBi, UBi, LBj, UBj
      integer, intent(in) :: IminS, ImaxS, JminS, JmaxS
      integer, intent(in) :: nstp, nnew
!
# ifdef ASSUMED_SHAPE
      real(r8), intent(in) :: Hz(LBi:,LBj:,:)
#  ifdef WET_DRY
      real(r8), intent(in) :: rmask_wet(LBi:,LBj:)
#  endif
#  ifdef BBL_MODEL
      real(r8), intent(in) :: bustrc(LBi:,LBj:)
      real(r8), intent(in) :: bvstrc(LBi:,LBj:)
      real(r8), intent(in) :: bustrw(LBi:,LBj:)
      real(r8), intent(in) :: bvstrw(LBi:,LBj:)
      real(r8), intent(in) :: bustrcwmax(LBi:,LBj:)
      real(r8), intent(in) :: bvstrcwmax(LBi:,LBj:)
#  endif
      real(r8), intent(in) :: bustr(LBi:,LBj:)
      real(r8), intent(in) :: bvstr(LBi:,LBj:)
#  if defined SED_MORPH
      real(r8), intent(inout):: bed_thick(LBi:,LBj:,:)
#  endif
      real(r8), intent(inout) :: t(LBi:,LBj:,:,:,:)
      real(r8), intent(inout) :: ero_flux(LBi:,LBj:,:)
      real(r8), intent(inout) :: settling_flux(LBi:,LBj:,:)
      real(r8), intent(inout) :: bed(LBi:,LBj:,:,:)
      real(r8), intent(inout) :: bed_frac(LBi:,LBj:,:,:)
      real(r8), intent(inout) :: bed_mass(LBi:,LBj:,:,:,:)
      real(r8), intent(inout) :: bottom(LBi:,LBj:,:)
# else
      real(r8), intent(in) :: Hz(LBi:UBi,LBj:UBj,N(ng))
#  ifdef WET_DRY
      real(r8), intent(in) :: rmask_wet(LBi:UBi,LBj:UBj)
#  endif
#  ifdef BBL_MODEL
      real(r8), intent(in) :: bustrc(LBi:UBi,LBj:UBj)
      real(r8), intent(in) :: bvstrc(LBi:UBi,LBj:UBj)
      real(r8), intent(in) :: bustrw(LBi:UBi,LBj:UBj)
      real(r8), intent(in) :: bvstrw(LBi:UBi,LBj:UBj)
      real(r8), intent(in) :: bustrcwmax(LBi:UBi,LBj:UBj)
      real(r8), intent(in) :: bvstrcwmax(LBi:UBi,LBj:UBj)
#  endif
      real(r8), intent(in) :: bustr(LBi:UBi,LBj:UBj)
      real(r8), intent(in) :: bvstr(LBi:UBi,LBj:UBj)
#  if defined SED_MORPH
      real(r8), intent(inout):: bed_thick(LBi:UBi,LBj:UBj,3)
#  endif
      real(r8), intent(inout) :: t(LBi:UBi,LBj:UBj,N(ng),3,NT(ng))
      real(r8), intent(inout) :: ero_flux(LBi:UBi,LBj:UBj,NST)
      real(r8), intent(inout) :: settling_flux(LBi:UBi,LBj:UBj,NST)
      real(r8), intent(inout) :: bed(LBi:UBi,LBj:UBj,Nbed,MBEDP)
      real(r8), intent(inout) :: bed_frac(LBi:UBi,LBj:UBj,Nbed,NST)
      real(r8), intent(inout) :: bed_mass(LBi:UBi,LBj:UBj,Nbed,1:2,NST)
      real(r8), intent(inout) :: bottom(LBi:UBi,LBj:UBj,MBOTP)
# endif
!
!  Local variable declarations.
!
      integer :: Ksed, i, indx, ised, j, k, ks
      integer :: bnew
 
      real(r8) :: cff, cff1, cff2, cff3, cff4
 
      real(r8), dimension(IminS:ImaxS,N(ng)) :: Hz_inv
 
      real(r8), dimension(IminS:ImaxS,JminS:JmaxS) :: tau_w
 
# include "set_bounds.h"
 
# ifdef BEDLOAD
      bnew=nnew
# else
      bnew=nstp
# endif
!
!-----------------------------------------------------------------------
!  Compute sediment deposition, resuspension, and erosion.
!-----------------------------------------------------------------------
!
# if defined BEDLOAD_MPM || defined SUSPLOAD
#  ifdef BBL_MODEL
      DO j=jstr-1,jend+1
        DO i=istr-1,iend+1
          tau_w(i,j)=sqrt(bustrcwmax(i,j)*bustrcwmax(i,j)+              &
     &                    bvstrcwmax(i,j)*bvstrcwmax(i,j))
#   ifdef WET_DRY
          tau_w(i,j)=tau_w(i,j)*rmask_wet(i,j)
#   endif
        END DO
      END DO
#  else
      DO j=jstrm1,jendp1
        DO i=istrm1,iendp1
          tau_w(i,j)=0.5_r8*sqrt((bustr(i,j)+bustr(i+1,j))*             &
     &                           (bustr(i,j)+bustr(i+1,j))+             &
     &                           (bvstr(i,j)+bvstr(i,j+1))*             &
     &                           (bvstr(i,j)+bvstr(i,j+1)))
#   ifdef WET_DRY
          tau_w(i,j)=tau_w(i,j)*rmask_wet(i,j)
#   endif
        END DO
      END DO
#  endif
# endif
!
!-----------------------------------------------------------------------
!  Sediment deposition and resuspension near the bottom.
!-----------------------------------------------------------------------
!
!  The deposition and resuspension of sediment on the bottom "bed"
!  is due to precepitation settling_flux, already computed, and the
!  resuspension (erosion, hence called ero_flux). The resuspension is
!  applied to the bottom-most grid box value qc(:,1) so the total mass
!  is conserved. Restrict "ero_flux" so that "bed" cannot go negative
!  after both fluxes are applied.
!
      j_loop : DO j=jstr,jend
        DO k=1,n(ng)
          DO i=istr,iend
            hz_inv(i,k)=1.0_r8/hz(i,j,k)
          END DO
        END DO
!
        sed_loop: DO ised=1,nst
          indx=idsed(ised)
          DO i=istr,iend
!
!  Calculate critical shear stress in Pa
!
# if defined COHESIVE_BED
            cff = rho0/bed(i,j,1,ibtcr)
# elif defined MIXED_BED
            cff = max(bottom(i,j,idprp)*bed(i,j,1,ibtcr)/rho0+          &
     &            (1.0_r8-bottom(i,j,idprp))*tau_ce(ised,ng),           &
     &            tau_ce(ised,ng))
            cff=1.0_r8/cff
# else
            cff=1.0_r8/tau_ce(ised,ng)
# endif
!
!  Compute erosion, ero_flux (kg/m2).
!
            cff1=(1.0_r8-bed(i,j,1,iporo))*bed_frac(i,j,1,ised)
            cff2=dt(ng)*erate(ised,ng)*cff1
            cff3=srho(ised,ng)*cff1
            cff4=bed_mass(i,j,1,bnew,ised)
            ero_flux(i,j,ised)=                                         &
     &               min(max(0.0_r8,cff2*(cff*tau_w(i,j)-1.0_r8)),      &
     &                   min(cff3*bottom(i,j,iactv),cff4)+              &
     &                   settling_flux(i,j,ised))
!
!  Update global tracer variables (m Tunits for nnew indx, Tuints for 3)
!  for erosive flux.
!
            t(i,j,1,nnew,indx)=t(i,j,1,nnew,indx)+ero_flux(i,j,ised)
          END DO
        END DO sed_loop
      END DO j_loop
!
      RETURN

References mod_scalars::dt, mod_sediment::erate, mod_sediment::iactv, mod_sediment::ibtcr, mod_sediment::idprp, mod_sediment::idsed, mod_sediment::iporo, mod_scalars::rho0, mod_sediment::srho, and mod_sediment::tau_ce.

Referenced by sed_fluxes().

Here is the caller graph for this function: