I have a confused question. The boundary conditions for T and S are set (RadNud) in .in file:
Code: Select all
! W S E N
! e o a o
! s u s r
! t t t t
! h h
!
! 1 2 3 4
LBC(isFsur) == Clo Clo Che Clo ! free-surface
LBC(isUbar) == Clo Clo Shc Clo ! 2D U-momentum
LBC(isVbar) == Clo Clo Shc Clo ! 2D V-momentum
LBC(isUvel) == Clo Clo RadNud Clo ! 3D U-momentum
LBC(isVvel) == Clo Clo RadNud Clo ! 3D V-momentum
LBC(isMtke) == Clo Clo Clo Clo ! mixing TKE
LBC(isTvar) == Clo Clo RadNud Clo \ ! temperature
Clo Clo RadNud Clo \ ! salinity
Clo Clo Clo Clo \ ! inert(1)
Clo Clo Clo Clo \ ! inert(2)
Clo Clo Clo Clo \ ! inert(3)
Clo Clo Clo Clo \ ! inert(4)
Clo Clo Clo Clo \ ! inert(5)
Clo Clo Clo Clo ! inert(6)
! Ice boundary conditions
LBC(isAice) == Clo Clo Clo Clo ! ice concentration
LBC(isHice) == Clo Clo Clo Clo ! ice thickness
LBC(isHsno) == Clo Clo Clo Clo ! snow thickness
LBC(isTice) == Clo Clo Clo Clo ! ice temperature
LBC(isApond)== Clo Clo Clo Clo ! surface water
LBC(isHpond)== Clo Clo Clo Clo ! surface water
LBC(isSig11)== Clo Clo Clo Clo ! sigma-11
LBC(isSig12)== Clo Clo Clo Clo ! sigma-12
LBC(isSig22)== Clo Clo Clo Clo ! sigma-22
LBC(isUice) == Gra Gra Gra Gra ! ice U-momentum
LBC(isVice) == Gra Gra Gra Gra ! ice V-momentum
! Adjoint-based algorithms can have different lateral boundary
! conditions keywords.
ad_LBC(isFsur) == Per Clo Per Clo ! free-surface
ad_LBC(isUbar) == Per Clo Per Clo ! 2D U-momentum
ad_LBC(isVbar) == Per Clo Per Clo ! 2D U-momentum
ad_LBC(isUvel) == Per Clo Per Clo ! 3D U-momentum
ad_LBC(isVvel) == Per Clo Per Clo ! 3D V-momentum
ad_LBC(isMtke) == Per Clo Per Clo ! mixing TKE
ad_LBC(isTvar) == Per Clo Per Clo \ ! temperature
Per Clo Per Clo ! salinity
! Set lateral open boundary edge volume conservation switch for
! nonlinear model and adjoint-based algorithms. Usually activated
! with radiation boundary conditions to enforce global mass
! conservation, except if tidal forcing is enabled. [1:Ngrids].
VolCons(west) == F ! western boundary
VolCons(east) == F ! eastern boundary
VolCons(south) == F ! southern boundary
VolCons(north) == F ! northern boundary
ad_VolCons(west) == F ! western boundary
ad_VolCons(east) == F ! eastern boundary
ad_VolCons(south) == F ! southern boundary
ad_VolCons(north) == F ! northern boundary
But in the log.file, 'temp' and 'salt' boundary conditions are closed:
Code: Select all
Lateral Boundary Conditions: NLM
============================
Variable Grid West Edge South Edge East Edge North Edge
--------- ---- ---------- ---------- ---------- ----------
zeta 1 Closed Closed Chapman Exp Closed
ubar 1 Closed Closed Shchepetkin Closed
vbar 1 Closed Closed Shchepetkin Closed
u 1 Closed Closed Rad + Nud Closed
v 1 Closed Closed Rad + Nud Closed
temp 1 Closed Closed Closed Closed
salt 1 Closed Closed Closed Closed
tke 1 Closed Closed Closed Closed
uice 1 Gradient Gradient Gradient Gradient
vice 1 Gradient Gradient Gradient Gradient
aice 1 Closed Closed Closed Closed
hice 1 Closed Closed Closed Closed
tisrf 1 Closed Closed Closed Closed
snow_thick 1 Closed Closed Closed Closed
apond 1 Closed Closed Closed Closed
hpond 1 Closed Closed Closed Closed
sig11 1 Closed Closed Closed Closed
sig12 1 Closed Closed Closed Closed
sig22 1 Closed Closed Closed Closed
