5#if defined NONLINEAR && defined SEDIMENT
44 integer,
intent(in) :: ng, tile
48 character (len=*),
parameter :: myfile = &
57 & lbi, ubi, lbj, ubj, &
58 & imins, imaxs, jmins, jmaxs, &
71 & LBi, UBi, LBj, UBj, &
72 & IminS, ImaxS, JminS, JmaxS, &
88 integer,
intent(in) :: ng, tile
89 integer,
intent(in) :: LBi, UBi, LBj, UBj
90 integer,
intent(in) :: IminS, ImaxS, JminS, JmaxS
91 integer,
intent(in) :: nstp, nnew
94 real(r8),
intent(inout) :: bed_frac(LBi:,LBj:,:,:)
95 real(r8),
intent(inout) :: bottom(LBi:,LBj:,:)
97 real(r8),
intent(inout) :: bed_frac(LBi:UBi,LBj:UBj,Nbed,NST)
98 real(r8),
intent(inout) :: bottom(LBi:UBi,LBj:UBj,MBOTP)
103 integer :: i, ised, j
105 real(r8) :: cff1, cff2, cff3, cff4
106 real(r8),
parameter :: eps = 1.0e-14_r8
108# include "set_bounds.h"
115 j_loop :
DO j=jstr,jend
122 cff1=cff1*
tau_ce(ised,ng)**bed_frac(i,j,1,ised)
123 cff2=cff2*
sd50(ised,ng)**bed_frac(i,j,1,ised)
124 cff3=cff3*(
wsed(ised,ng)+eps)**bed_frac(i,j,1,ised)
125 cff4=cff4*
srho(ised,ng)**bed_frac(i,j,1,ised)
127 bottom(i,j,
itauc)=cff1
129 bottom(i,j,
iwsed)=cff3
130 bottom(i,j,
idens)=max(cff4,1050.0_r8)
143 & lbi, ubi, lbj, ubj, 1, mbotp, &
147 & lbi, ubi, lbj, ubj, 1, mbotp, &
subroutine bc_r3d_tile(ng, tile, lbi, ubi, lbj, ubj, lbk, ubk, a)
logical, dimension(:), allocatable ewperiodic
logical, dimension(:), allocatable nsperiodic
real(r8), dimension(:), allocatable zob
type(t_sedbed), dimension(:), allocatable sedbed
real(r8), dimension(:,:), allocatable srho
real(r8), dimension(:,:), allocatable sd50
real(r8), dimension(:,:), allocatable wsed
real(r8), dimension(:,:), allocatable tau_ce
integer, dimension(:), allocatable nnew
integer, dimension(:), allocatable nstp
subroutine mp_exchange4d(ng, tile, model, nvar, lbi, ubi, lbj, ubj, lbk, ubk, lbt, ubt, nghost, ew_periodic, ns_periodic, a, b, c)
subroutine mp_exchange3d(ng, tile, model, nvar, lbi, ubi, lbj, ubj, lbk, ubk, nghost, ew_periodic, ns_periodic, a, b, c, d)
subroutine sed_surface_tile(ng, tile, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs, nstp, nnew, bed_frac, bottom)
subroutine, public sed_surface(ng, tile)
recursive subroutine wclock_off(ng, model, region, line, routine)
recursive subroutine wclock_on(ng, model, region, line, routine)