!!
SUBROUTINE ana_tobc (ng, tile, model)
!
!! git $Id$
!! svn $Id: ana_tobc.h 1151 2023-02-09 03:08:53Z arango $
!!======================================================================
!! Copyright (c) 2002-2023 The ROMS/TOMS Group !
!! Licensed under a MIT/X style license !
!! See License_ROMS.md !
!=======================================================================
! !
! This routine sets tracer-type variables open boundary conditions !
! using analytical expressions. !
! !
!=======================================================================
!
USE mod_param
USE mod_boundary
USE mod_grid
USE mod_ncparam
USE mod_ocean
USE mod_stepping
!
! Imported variable declarations.
!
integer, intent(in) :: ng, tile, model
!
! Local variable declarations.
!
character (len=*), parameter :: MyFile = &
& __FILE__
!
#include "tile.h"
!
CALL ana_tobc_tile (ng, tile, model, &
& LBi, UBi, LBj, UBj, &
& IminS, ImaxS, JminS, JmaxS, &
& nstp(ng), &
& GRID(ng) % z_r, &
& OCEAN(ng) % t)
!
! Set analytical header file name used.
!
#ifdef DISTRIBUTE
IF (Lanafile) THEN
#else
IF (Lanafile.and.(tile.eq.0)) THEN
#endif
ANANAME(34)=MyFile
END IF
!
RETURN
END SUBROUTINE ana_tobc
!
!***********************************************************************
SUBROUTINE ana_tobc_tile (ng, tile, model, &
& LBi, UBi, LBj, UBj, &
& IminS, ImaxS, JminS, JmaxS, &
& nstp, &
& z_r, t)
!***********************************************************************
!
USE mod_param
USE mod_scalars
USE mod_boundary
USE mod_ncparam
USE mod_ocean
#ifdef SEDIMENT
USE mod_sediment
#endif
!
! Imported variable declarations.
!
integer, intent(in) :: ng, tile, model
integer, intent(in) :: LBi, UBi, LBj, UBj
integer, intent(in) :: IminS, ImaxS, JminS, JmaxS
integer, intent(in) :: nstp
#ifdef ASSUMED_SHAPE
real(r8), intent(in) :: z_r(LBi:,LBj:,:)
real(r8), intent(in) :: t(LBi:,LBj:,:,:,:)
#else
real(r8), intent(in) :: z_r(LBi:UBi,LBj:UBj,N(ng))
real(r8), intent(in) :: t(LBi:UBi,LBj:UBj,N(ng),3,NT(ng))
#endif
!
! Local variable declarations.
!
integer :: i, ised, itrc, j, k
!
real(r8) :: cff
#include "set_bounds.h"
!
!-----------------------------------------------------------------------
! Tracers open boundary conditions.
!-----------------------------------------------------------------------
!
#ifdef ESTUARY_TEST
IF (ANY(LBC(ieast,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Eastern_Edge(tile)) THEN
DO k=1,N(ng)
DO j=JstrT,JendT
BOUNDARY(ng)%t_east(j,k,itemp)=T0(ng)
# ifdef SALINITY
BOUNDARY(ng)%t_east(j,k,isalt)=0.0_r8
# endif
# ifdef SEDIMENT
DO ised=1,NST
BOUNDARY(ng)%t_east(j,k,idsed(ised))=0.0_r8
END DO
# endif
END DO
END DO
END IF
IF (ANY(LBC(iwest,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Western_Edge(tile)) THEN
DO k=1,N(ng)
DO j=JstrT,JendT
BOUNDARY(ng)%t_west(j,k,itemp)=T0(ng)
# ifdef SALINITY
BOUNDARY(ng)%t_west(j,k,isalt)=30.0_r8
# endif
# ifdef SEDIMENT
DO ised=1,NST
BOUNDARY(ng)%t_west(j,k,idsed(ised))=0.0_r8
END DO
# endif
END DO
END DO
END IF
#elif defined NJ_BIGHT
IF (ANY(LBC(ieast,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Eastern_Edge(tile)) THEN
DO k=1,N(ng)
DO j=JstrT,JendT
IF (z_r(Iend+1,j,k).ge.-15.0_r8) THEN
BOUNDARY(ng)%t_east(j,k,itemp)=2.04926425772840E+01_r8- &
& z_r(Iend+1,j,k)* &
& (2.64085084879392E-01_r8+ &
& z_r(Iend+1,j,k)* &
& (2.75112532853521E-01_r8+ &
& z_r(Iend+1,j,k)* &
& (9.20748976164887E-02_r8+ &
& z_r(Iend+1,j,k)* &
& (1.44907572574284E-02_r8+ &
& z_r(Iend+1,j,k)* &
& (1.07821568591208E-03_r8+ &
& z_r(Iend+1,j,k)* &
& (3.24031805390397E-05_r8+ &
& 1.26282685769027E-07_r8*
& z_r(Iend+1,j,k)))))))
# ifdef SALINITY
BOUNDARY(ng)%t_east(j,k,isalt)=3.06648914919313E+01_r8- &
& z_r(Iend+1,j,k)* &
& (1.47672526294673E-01_r8+ &
& z_r(Iend+1,j,k)* &
& (1.12645576031340E-01_r8+ &
& z_r(Iend+1,j,k)* &
& (3.90092328187102E-02_r8+ &
& z_r(Iend+1,j,k)* &
& (6.93901493744710E-03_r8+ &
& z_r(Iend+1,j,k)* &
& (6.60443669679294E-04_r8+ &
& z_r(Iend+1,j,k)* &
& (3.19179236195422E-05_r8+ &
& 6.17735263440932E-07_r8*
& z_r(Iend+1,j,k)))))))
# endif
ELSE
cff=TANH(1.1_r8*z_r(Iend+1,j,k)+15.9_r8)
t_east(j,k,itemp)=14.6_r8+6.70_r8*cff
# ifdef SALINITY
t_east(j,k,isalt)=31.3_r8-0.55_r8*cff
# endif
END IF
END DO
END DO
END IF
IF (ANY(LBC(isouth,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Southern_Edge(tile)) THEN
DO k=1,N(ng)
DO i=IstrT,IendT
IF (z_r(i,Jstr-1,k).ge.-15.0_r8) THEN
BOUNDARY(ng)%t_south(i,k,itemp)=2.04926425772840E+01_r8- &
& z_r(i,Jstr-1,k)* &
& (2.64085084879392E-01_r8+ &
& z_r(i,Jstr-1,k)* &
& (2.75112532853521E-01_r8+&
& z_r(i,Jstr-1,k)* &
& (9.20748976164887E-02_r8+&
& z_r(i,Jstr-1,k)* &
& (1.44907572574284E-02_r8+&
& z_r(i,Jstr-1,k)* &
& (1.07821568591208E-03_r8+&
& z_r(i,Jstr-1,k)* &
& (3.24031805390397E-05_r8+&
& 1.26282685769027E-07_r8*
& z_r(i,Jstr-1,k)))))))
# ifdef SALINITY
BOUNDARY(ng)%t_south(i,k,isalt)=3.06648914919313E+01_r8- &
& z_r(i,Jstr-1,k)* &
& (1.47672526294673E-01_r8+ &
& z_r(i,Jstr-1,k)* &
& (1.12645576031340E-01_r8+&
& z_r(i,Jstr-1,k)* &
& (3.90092328187102E-02_r8+&
& z_r(i,Jstr-1,k)* &
& (6.93901493744710E-03_r8+&
& z_r(i,Jstr-1,k)* &
& (6.60443669679294E-04_r8+&
& z_r(i,Jstr-1,k)* &
& (3.19179236195422E-05_r8+&
& 6.17735263440932E-07_r8*
& z_r(i,Jstr-1,k)))))))
# endif
ELSE
cff=TANH(1.1_r8*depth+15.9_r8)
BOUNDARY(ng)%t_south(i,k,itemp)=14.6_r8+6.70_r8*cff
# ifdef SALINITY
BOUNDARY(ng)%t_south(i,k,isalt)=31.3_r8-0.55_r8*cff
# endif
END IF
END DO
END DO
END IF
#elif defined SED_TEST1
IF (ANY(LBC(ieast,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Eastern_Edge(tile)) THEN
DO k=1,N(ng)
DO j=JstrT,JendT
BOUNDARY(ng)%t_east(j,k,itemp)=20.0_r8
# ifdef SALINITY
BOUNDARY(ng)%t_east(j,k,isalt)=0.0_r8
# endif
END DO
END DO
END IF
#else
IF (ANY(LBC(ieast,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Eastern_Edge(tile)) THEN
DO itrc=1,NT(ng)
DO k=1,N(ng)
DO j=JstrT,JendT
BOUNDARY(ng)%t_east(j,k,itrc)=0.0_r8
END DO
END DO
END DO
END IF
IF (ANY(LBC(iwest,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Western_Edge(tile)) THEN
DO itrc=1,NT(ng)
DO k=1,N(ng)
DO j=JstrT,JendT
BOUNDARY(ng)%t_west(j,k,itrc)=0.0_r8
END DO
END DO
END DO
END IF
IF (ANY(LBC(isouth,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Southern_Edge(tile)) THEN
DO itrc=1,NT(ng)
DO k=1,N(ng)
DO i=IstrT,IendT
BOUNDARY(ng)%t_south(i,k,itrc)=0.0_r8
END DO
END DO
END DO
END IF
IF (ANY(LBC(inorth,isTvar(:),ng)%acquire).and. &
& DOMAIN(ng)%Northern_Edge(tile)) THEN
DO itrc=1,NT(ng)
DO k=1,N(ng)
DO i=IstrT,IendT
BOUNDARY(ng)%t_north(i,k,itrc)=0.0_r8
END DO
END DO
END DO
END IF
#endif
!
RETURN
END SUBROUTINE ana_tobc_tile