!!
SUBROUTINE ana_wwave (ng, tile, model)
!
!! git $Id$
!! svn $Id: ana_wwave.h 1178 2023-07-11 17:50:57Z arango $
!!======================================================================
!! Copyright (c) 2002-2023 The ROMS/TOMS Group !
!! Licensed under a MIT/X style license !
!! See License_ROMS.md !
!=======================================================================
! !
! This subroutine sets wind induced wave amplitude, direction and !
! period to be used in the bottom boundary layer formulation. !
! !
!=======================================================================
!
USE mod_param
USE mod_forces
USE mod_grid
USE mod_ncparam
!
! Imported variable declarations.
!
integer, intent(in) :: ng, tile, model
!
! Local variable declarations.
!
character (len=*), parameter :: MyFile = &
& __FILE__
!
#include "tile.h"
!
CALL ana_wwave_tile (ng, tile, model, &
& LBi, UBi, LBj, UBj, &
& IminS, ImaxS, JminS, JmaxS, &
#ifdef WAVES_DIR
& FORCES(ng) % Dwave, &
#endif
#ifdef WAVES_HEIGHT
& FORCES(ng) % Hwave, &
#endif
#ifdef WAVES_LENGTH
& FORCES(ng) % Lwave, &
#endif
#ifdef WAVES_TOP_PERIOD
& FORCES(ng) % Pwave_top, &
#endif
#ifdef WAVES_BOT_PERIOD
& FORCES(ng) % Pwave_bot, &
#endif
#ifdef WAVES_UB
& FORCES(ng) % Uwave_rms, &
#endif
& GRID(ng) % angler, &
& GRID(ng) % h)
!
! Set analytical header file name used.
!
#ifdef DISTRIBUTE
IF (Lanafile) THEN
#else
IF (Lanafile.and.(tile.eq.0)) THEN
#endif
ANANAME(37)=MyFile
END IF
!
RETURN
END SUBROUTINE ana_wwave
!
!***********************************************************************
SUBROUTINE ana_wwave_tile (ng, tile, model, &
& LBi, UBi, LBj, UBj, &
& IminS, ImaxS, JminS, JmaxS, &
#ifdef WAVES_DIR
& Dwave, &
#endif
#ifdef WAVES_HEIGHT
& Hwave, &
#endif
#ifdef WAVES_LENGTH
& Lwave, &
#endif
#ifdef WAVES_TOP_PERIOD
& Pwave_top, &
#endif
#ifdef WAVES_BOT_PERIOD
& Pwave_bot, &
#endif
#ifdef WAVES_UB
& Uwave_rms, &
#endif
& angler, h)
!***********************************************************************
!
USE mod_param
USE mod_scalars
!
USE exchange_2d_mod, ONLY : exchange_r2d_tile
#ifdef DISTRIBUTE
USE mp_exchange_mod, ONLY : mp_exchange2d
#endif
!
! Imported variable declarations.
!
integer, intent(in) :: ng, tile, model
integer, intent(in) :: LBi, UBi, LBj, UBj
integer, intent(in) :: IminS, ImaxS, JminS, JmaxS
!
#ifdef ASSUMED_SHAPE
real(r8), intent(in) :: angler(LBi:,LBj:)
real(r8), intent(in) :: h(LBi:,LBj:)
# #ifdef WAVES_DIR
real(r8), intent(inout) :: Dwave(LBi:,LBj:)
# endif
# ifdef WAVES_HEIGHT
real(r8), intent(inout) :: Hwave(LBi:,LBj:)
# endif
# ifdef WAVES_LENGTH
real(r8), intent(inout) :: Lwave(LBi:,LBj:)
# endif
# ifdef WAVES_TOP_PERIOD
real(r8), intent(inout) :: Pwave_top(LBi:,LBj:)
# endif
# ifdef WAVES_BOT_PERIOD
real(r8), intent(inout) :: Pwave_bot(LBi:,LBj:)
# endif
# ifdef WAVES_UB
real(r8), intent(inout) :: Uwave_rms(LBi:,LBj:)
# endif
#else
real(r8), intent(in) :: angler(LBi:UBi,LBj:UBj)
real(r8), intent(in) :: h(LBi:UBi,LBj:UBj)
# #ifdef WAVES_DIR
real(r8), intent(inout) :: Dwave(LBi:UBi,LBj:UBj)
# endif
# ifdef WAVES_HEIGHT
real(r8), intent(inout) :: Hwave(LBi:UBi,LBj:UBj)
# endif
# ifdef WAVES_LENGTH
real(r8), intent(inout) :: Lwave(LBi:UBi,LBj:UBj)
# endif
# ifdef WAVES_TOP_PERIOD
real(r8), intent(inout) :: Pwave_top(LBi:UBi,LBj:UBj)
# endif
# ifdef WAVES_BOT_PERIOD
real(r8), intent(inout) :: Pwave_bot(LBi:UBi,LBj:UBj)
# endif
# ifdef WAVES_UB
real(r8), intent(inout) :: Uwave_rms(LBi:UBi,LBj:UBj)
# endif
#endif
!
! Local variable declarations.
!
integer :: i, j
!
real(r8) :: cff, wdir
#include "set_bounds.h"
!
!-----------------------------------------------------------------------
! Set wind induced wave amplitude (m), direction (radians) and
! period (s) at RHO-points.
!-----------------------------------------------------------------------
!
#if defined MY_APPLICATION
DO j=JstrT,JendT
DO i=IstrT,IendT
# ifdef WAVES_DIR
Dwave(i,j)=90.0_r8*deg2rad
# endif
# ifdef WAVES_HEIGHT
Hwave(i,j)=2.0_r8
# endif
# ifdef WAVES_LENGTH
Lwave(i,j)=20.0_r8
# endif
# ifdef WAVES_TOP_PERIOD
Pwave_top(i,j)=8.0_r8
# endif
# ifdef WAVES_BOT_PERIOD
Pwave_bot(i,j)=8.0_r8
# endif
# ifdef WAVES_UB
Uwave_rms(i,j)=0.2_r8
# endif
END DO
END DO
#else
ana_wwave: no values provided for Hwave, Dwave, Pwave, Lwave.
#endif
!
! Exchange boundary data.
!
#if defined WAVES_DIR
IF (EWperiodic(ng).or.NSperiodic(ng)) THEN
CALL exchange_r2d_tile (ng, tile, &
& LBi, UBi, LBj, UBj, &
& Dwave)
END IF
# ifdef DISTRIBUTE
CALL mp_exchange2d (ng, tile, model, 1, &
& LBi, UBi, LBj, UBj, &
& NghostPoints, &
& EWperiodic(ng), NSperiodic(ng), &
& Dwave)
# endif
#endif
#ifdef WAVES_HEIGHT
IF (EWperiodic(ng).or.NSperiodic(ng)) THEN
CALL exchange_r2d_tile (ng, tile, &
& LBi, UBi, LBj, UBj, &
& Hwave)
END IF
# ifdef DISTRIBUTE
CALL mp_exchange2d (ng, tile, model, 1, &
& LBi, UBi, LBj, UBj, &
& NghostPoints, &
& EWperiodic(ng), NSperiodic(ng), &
& Hwave)
# endif
#endif
#ifdef WAVES_LENGTH
IF (EWperiodic(ng).or.NSperiodic(ng)) THEN
CALL exchange_r2d_tile (ng, tile, &
& LBi, UBi, LBj, UBj, &
& Lwave)
END IF
# ifdef DISTRIBUTE
CALL mp_exchange2d (ng, tile, model, 1, &
& LBi, UBi, LBj, UBj, &
& NghostPoints, &
& EWperiodic(ng), NSperiodic(ng), &
& Lwave)
# endif
#endif
#ifdef WAVES_TOP_PERIOD
IF (EWperiodic(ng).or.NSperiodic(ng)) THEN
CALL exchange_r2d_tile (ng, tile, &
& LBi, UBi, LBj, UBj, &
& Pwave_top)
END IF
# ifdef DISTRIBUTE
CALL mp_exchange2d (ng, tile, model, 1, &
& LBi, UBi, LBj, UBj, &
& NghostPoints, &
& EWperiodic(ng), NSperiodic(ng), &
& Pwave_top)
# endif
#endif
#ifdef WAVES_BOT_PERIOD
IF (EWperiodic(ng).or.NSperiodic(ng)) THEN
CALL exchange_r2d_tile (ng, tile, &
& LBi, UBi, LBj, UBj, &
& Pwave_bot)
END IF
# ifdef DISTRIBUTE
CALL mp_exchange2d (ng, tile, model, 1, &
& LBi, UBi, LBj, UBj, &
& NghostPoints, &
& EWperiodic(ng), NSperiodic(ng), &
& Pwave_bot)
# endif
#endif
#ifdef WAVES_UB
IF (EWperiodic(ng).or.NSperiodic(ng)) THEN
CALL exchange_r2d_tile (ng, tile, &
& LBi, UBi, LBj, UBj, &
& Uwave_rms)
END IF
# ifdef DISTRIBUTE
CALL mp_exchange2d (ng, tile, model, 1, &
& LBi, UBi, LBj, UBj, &
& NghostPoints, &
& EWperiodic(ng), NSperiodic(ng), &
& Uwave_rms)
# endif
#endif
!
RETURN
END SUBROUTINE ana_wwave_tile