Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Wed Jun 19, 2019 5:00 am




Post new topic Reply to topic  [ 7 posts ] 

All times are UTC

Author Message
PostPosted: Mon May 22, 2017 1:35 pm 
Offline

Joined: Wed Apr 13, 2016 7:37 pm
Posts: 25
Location: University of Sao Paulo (USP)
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


Top
 Profile  
Reply with quote  
PostPosted: Mon May 22, 2017 3:25 pm 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3601
Location: IMS/UAF, USA
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.


Top
 Profile  
Reply with quote  
PostPosted: Mon May 22, 2017 5:15 pm 
Offline

Joined: Wed Apr 13, 2016 7:37 pm
Posts: 25
Location: University of Sao Paulo (USP)
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.


Top
 Profile  
Reply with quote  
PostPosted: Mon May 22, 2017 5:33 pm 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3601
Location: IMS/UAF, USA
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.


Top
 Profile  
Reply with quote  
PostPosted: Tue May 23, 2017 11:01 am 
Offline
User avatar

Joined: Mon Apr 28, 2003 5:44 pm
Posts: 488
Location: Rutgers University
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


Top
 Profile  
Reply with quote  
PostPosted: Tue May 23, 2017 3:28 pm 
Offline

Joined: Wed Apr 13, 2016 7:37 pm
Posts: 25
Location: University of Sao Paulo (USP)
Thank you for the replies.
I'll also try Gradient BC (Gradient = 0)

regards


Top
 Profile  
Reply with quote  
PostPosted: Wed May 24, 2017 1:58 pm 
Offline

Joined: Wed Apr 13, 2016 7:37 pm
Posts: 25
Location: University of Sao Paulo (USP)
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


Top
 Profile  
Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 7 posts ] 

All times are UTC


Who is online

Users browsing this forum: No registered users and 2 guests


You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group