subroutine SwanTranspAc ( amat , rhs , leakcf, ac2 , ac1 , &
cgo , cax , cay , cad , cas , &
anybin, rdx , rdy , spcsig, spcdir, &
obredf, idcmin, idcmax, iscmin, iscmax, &
iddlow, iddtop, isslow, isstop, anyblk, &
trac0 , trac1 )
!
! --|-----------------------------------------------------------|--
! | Delft University of Technology |
! | Faculty of Civil Engineering and Geosciences |
! | Environmental Fluid Mechanics Section |
! | P.O. Box 5048, 2600 GA Delft, The Netherlands |
! | |
! | Programmer: Marcel Zijlema |
! --|-----------------------------------------------------------|--
!
!
! SWAN (Simulating WAves Nearshore); a third generation wave model
! Copyright (C) 1993-2017 Delft University of Technology
!
! This program is free software; you can redistribute it and/or
! modify it under the terms of the GNU General Public License as
! published by the Free Software Foundation; either version 2 of
! the License, or (at your option) any later version.
!
! This program is distributed in the hope that it will be useful,
! but WITHOUT ANY WARRANTY; without even the implied warranty of
! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
! GNU General Public License for more details.
!
! A copy of the GNU General Public License is available at
! http://www.gnu.org/copyleft/gpl.html#SEC3
! or by writing to the Free Software Foundation, Inc.,
! 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
!
!
! Authors
!
! 40.80: Marcel Zijlema
!
! Updates
!
! 40.80, August 2007: New subroutine
! 40.85, August 2008: add tranport terms for output purposes
! 41.00, February 2009: add GSE correction
! 41.07, July 2009: add explicit scheme for sigma space
!
! Purpose
!
! Computes the transport part of the action balance equation
!
! Modules used
!
use ocpcomm4
use swcomm2
use swcomm3
use swcomm4
use SwanGriddata
!
implicit none
!
! Argument variables
!
integer, intent(in) :: iddlow ! minimum direction bin that is propagated within a sweep
integer, intent(in) :: iddtop ! maximum direction bin that is propagated within a sweep
integer, intent(in) :: isslow ! minimum frequency that is propagated within a sweep
integer, intent(in) :: isstop ! maximum frequency that is propagated within a sweep
!
integer, dimension(MSC), intent(in) :: idcmax ! maximum frequency-dependent counter in directional space
integer, dimension(MSC), intent(in) :: idcmin ! minimum frequency-dependent counter in directional space
integer, dimension(MDC), intent(in) :: iscmax ! maximum direction-dependent counter in frequency space
integer, dimension(MDC), intent(in) :: iscmin ! minimum direction-dependent counter in frequency space
!
real, dimension(MDC,MSC,nverts), intent(in) :: ac1 ! action density at previous time level
real, dimension(MDC,MSC,nverts), intent(in) :: ac2 ! action density at current time level
real, dimension(MDC,MSC,5), intent(out) :: amat ! coefficient matrix of system of equations in (sigma,theta) space
! 1: correspond to point (l ,m )
! 2: correspond to point (l-1,m )
! 3: correspond to point (l+1,m )
! 4: correspond to point (l ,m-1)
! 5: correspond to point (l ,m+1)
real, dimension(MDC,MSC), intent(in) :: cad ! wave transport velocity in theta-direction
real, dimension(MDC,MSC), intent(in) :: cas ! wave transport velocity in sigma-direction
real, dimension(MDC,MSC,ICMAX), intent(in) :: cax ! wave transport velocity in x-direction
real, dimension(MDC,MSC,ICMAX), intent(in) :: cay ! wave transport velocity in y-direction
real, dimension(MSC,ICMAX), intent(in) :: cgo ! group velocity
real, dimension(MDC,MSC), intent(out) :: leakcf ! leak coefficient
real, dimension(MDC,MSC,2), intent(in) :: obredf ! action reduction coefficient based on transmission
real, dimension(2), intent(in) :: rdx ! first component of contravariant base vector rdx(b) = a^(b)_1
real, dimension(2), intent(in) :: rdy ! second component of contravariant base vector rdy(b) = a^(b)_2
real, dimension(MDC,MSC), intent(out) :: rhs ! right-hand side of system of equations in (sigma,theta) space
real, dimension(MDC,6), intent(in) :: spcdir ! (*,1): spectral direction bins (radians)
! (*,2): cosine of spectral directions
! (*,3): sine of spectral directions
! (*,4): cosine^2 of spectral directions
! (*,5): cosine*sine of spectral directions
! (*,6): sine^2 of spectral directions
real, dimension(MSC), intent(in) :: spcsig ! relative frequency bins
real, dimension(MDC,MSC,MTRNP), intent(out) :: trac0 ! explicit part of propagation in present vertex for output purposes
real, dimension(MDC,MSC,MTRNP), intent(out) :: trac1 ! implicit part of propagation in present vertex for output purposes
!
logical, dimension(MDC,MSC), intent(in) :: anybin ! true if bin is active in considered sweep
logical, dimension(MDC,MSC), intent(out) :: anyblk ! true if bin is blocked by a counter current based on a CFL criterion
!
! Local variables
!
integer, save :: ient = 0 ! number of entries in this subroutine
!
! Structure
!
! Description of the pseudo code
!
! Source text
!
if (ltrace) call strace (ient,'SwanTranspAc')
!
! set matrix and right-hand side to zero
!
amat = 0.
rhs = 0.
!
! set leak and transport coefficients to zero
!
leakcf = 0.
trac0 = 0.
trac1 = 0.
!
! compute transport in x-y space
!
!TIMG call SWTSTA(140)
call SwanTranspX ( amat , rhs , ac2 , ac1 , cax , cay , &
rdx , rdy , obredf, idcmin, idcmax, isslow, &
isstop , trac0, trac1 )
!
! add GSE correction, if appropriate
!
if ( WAVAGE > 0. ) call SwanGSECorr ( rhs, ac2, cgo, spcdir, idcmin, idcmax, isslow, isstop, trac0 )
!TIMG call SWTSTO(140)
!
! compute transport in theta space
!
!TIMG call SWTSTA(142)
if ( IREFR /= 0 ) then
!
call STRSD ( DDIR , idcmin , idcmax , cad , &
amat(1,1,4), amat(1,1,1), amat(1,1,5), rhs , &
ac2 , isstop , anybin , leakcf , &
trac0 , trac1 )
!
endif
!TIMG call SWTSTO(142)
!
! compute transport in sigma space
!
!TIMG call SWTSTA(141)
if ( (DYNDEP .OR. ICUR /= 0) .and. ITFRE /= 0 ) then
!
if ( int(PNUMS(8)) == 1 ) then
!
! implicit scheme
!
call STRSSI ( spcsig , cas , amat(1,1,2), amat(1,1,1), &
amat(1,1,3), anybin, rhs , ac2 , &
iscmin , iscmax, iddlow , iddtop , &
trac0 , trac1 )
!
elseif ( int(PNUMS(8)) == 2 ) then
!
! explicit scheme
!
call STRSSB ( iddlow, iddtop, idcmin, idcmax, isstop, &
cax , cay , cas , ac2 , spcsig, &
rhs , anyblk, rdx , rdy , trac0 )
!
endif
!
endif
!TIMG call SWTSTO(141)
!
end subroutine SwanTranspAc