Upwelling Channel with 4 Open Boundaries

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
hpftcb
Posts: 30
Joined: Wed Apr 13, 2016 7:37 pm
Location: OceanPact

Upwelling Channel with 4 Open Boundaries

#1 Unread post by hpftcb »

Hi all,

I'D like to know if there is a way to set up an upwelling channel case with four open boundaries (on the contrary of Roms' upwelling test case) and zero inflow and outflow current.
I tried the following configuration, yet the wind generated current did not manage to radiate from the domain, it kept recirculating along the boundaries


LBC(isFsur) == Che Che Che Che ! free-surface
LBC(isUbar) == Shc Shc Shc Shc ! 2D U-momentum
LBC(isVbar) == Shc Shc Shc Shc ! 2D V-momentum
LBC(isUvel) == RadNud RadNud RadNud RadNud ! 3D U-momentum
LBC(isVvel) == RadNud RadNud RadNud RadNud ! 3D V-momentum
LBC(isMtke) == Gra Gra Gra Gra ! mixing TKE

LBC(isTvar) == RadNud RadNud RadNud RadNud ! temperature
RadNud RadNud RadNud RadNud ! salinity
RadNud RadNud RadNud RadNud ! dye_01

My domain is a 2.8 X 1.5 degree rectangle with 1.1 km resolution.

Regards,

Fernando T. C. Barreto

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

Re: Upwelling Channel with 4 Open Boundaries

#2 Unread post by kate »

What values are you using on the boundary? Is this a realistic domain or an idealized one? There's a reason the toy problems are all periodic channels or closed basins - except the Kelvin wave which looks awful if you run it longer than in its input file.

hpftcb
Posts: 30
Joined: Wed Apr 13, 2016 7:37 pm
Location: OceanPact

Re: Upwelling Channel with 4 Open Boundaries

#3 Unread post by hpftcb »

Hi, thanks for the reply

Yes,It's an idealized experiment (bathymetry is attached). It's a 3D case with 30 sigma layers (vt=2,vs=4), with salinity and temp field build from one WOA point (uniform horizontally), zero current velocity in all boundaries, and 5 m/s eastern wind input with ana_smflux.h. Boundaries conditions were build from initial condition file. Time step is 180.0d0.

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

Re: Upwelling Channel with 4 Open Boundaries

#4 Unread post by kate »

So your boundary conditions are level isopycnal surfaces. You blow a wind which is going to drive currents in balance with non-level isopycnals. If the boundary conditions aren't consistent with this, of course you'll get rim currents.

User avatar
wilkin
Posts: 872
Joined: Mon Apr 28, 2003 5:44 pm
Location: Rutgers University
Contact:

Re: Upwelling Channel with 4 Open Boundaries

#5 Unread post by wilkin »

If you are trying to set-up something that is effectively 1-dimensional in the vertical (you have no horizontal variation in density or wind forcing) then you might want a small double periodic domain. See the BIO_TOY and SED_TOY examples for guidance on how to do this.
John Wilkin: DMCS Rutgers University
71 Dudley Rd, New Brunswick, NJ 08901-8521, USA. ph: 609-630-0559 jwilkin@rutgers.edu

hpftcb
Posts: 30
Joined: Wed Apr 13, 2016 7:37 pm
Location: OceanPact

Re: Upwelling Channel with 4 Open Boundaries

#6 Unread post by hpftcb »

Thank you for the replies.
I'll also try Gradient BC (Gradient = 0)

regards

hpftcb
Posts: 30
Joined: Wed Apr 13, 2016 7:37 pm
Location: OceanPact

Re: Upwelling Channel with 4 Open Boundaries

#7 Unread post by hpftcb »

I tried with grad = 0 in all boundaries and the domain energy kept rising (velocity rising) until my case blew up. Does anyone know if what I did is right, if I can't use grad to termohaline or something like that.

Thank you in advance


LBC(isFsur) == Gra Gra Gra Gra ! free-surface
LBC(isUbar) == Gra Gra Gra Gra ! 2D U-momentum
LBC(isVbar) == Gra Gra Gra Gra ! 2D V-momentum
LBC(isUvel) == Gra Gra Gra Gra ! 3D U-momentum
LBC(isVvel) == Gra Gra Gra Gra ! 3D V-momentum
LBC(isMtke) == Gra Gra Gra Gra ! mixing TKE

LBC(isTvar) == Gra Gra Gra Gra ! temperature
Gra Gra Gra Gra ! salinity
Gra Gra Gra Gra ! dye_01

! Adjoint-based algorithms can have different lateral boundary
! conditions keywords.

ad_LBC(isFsur) == Gra Gra Gra Gra ! free-surface
ad_LBC(isUbar) == Gra Gra Gra Gra ! 2D U-momentum
ad_LBC(isVbar) == Gra Gra Gra Gra ! 2D V-momentum
ad_LBC(isUvel) == Gra Gra Gra Gra ! 3D U-momentum
ad_LBC(isVvel) == Gra Gra Gra Gra ! 3D V-momentum
ad_LBC(isMtke) == Gra Gra Gra Gra ! mixing TKE

ad_LBC(isTvar) == Gra Gra Gra Gra ! temperature
Gra Gra Gra Gra ! salinity
Gra Gra Gra Gra ! dye_01

Post Reply