Fail to active T and S boundary conditons

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
cloudy4568
Posts: 39
Joined: Thu Oct 17, 2019 8:56 pm
Location: Clarkson University

Fail to active T and S boundary conditons

#1 Unread post by cloudy4568 »

Hi all,

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
Data of boundary conditions are well prepared.

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
Do you have any idea about this? Thank you.

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: Fail to active T and S boundary conditons

#2 Unread post by kate »

You only have two tracers, right?

Code: Select all

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)
You need to change this to:

Code: Select all

LBC(isTvar) ==   Clo     Clo     RadNud  Clo    \    ! temperature
                    Clo     Clo     RadNud  Clo 

cloudy4568
Posts: 39
Joined: Thu Oct 17, 2019 8:56 pm
Location: Clarkson University

Re: Fail to active T and S boundary conditons

#3 Unread post by cloudy4568 »

That's the reason! Thank you, Kate!

Post Reply