Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Sun Aug 18, 2019 9:32 am




Post new topic Reply to topic  [ 2 posts ] 

All times are UTC

Author Message
PostPosted: Fri Feb 19, 2016 2:27 pm 
Offline

Joined: Mon Jun 30, 2003 2:29 pm
Posts: 31
Location: CCPO/ODU, USA
I am trying to set up roms to represent alternating flow in a uniform depth channel (x-z). It is 40 km long and periodic in y. density is uniform. vertical viscosity is constant. The boundaries are open on east and west ends. I want the flow to have ubar = uo sin(2 pi t/period) for different choices of uo.

I have tried using ana_m2obc to set boundary values of ubar to the desired flow. I have also set ana_fsobc to have the desired changing elevation at either end. I finally tried to set ana_m3obc to the steady vertical frictional structure. OBC choices are

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

Because of the friction, the analytical solution is complicated (although possible). Do I really need to give zeta, ubar and u(z) at either boundary from an analytical solution?

Is there a way to simply set ubar at the ends and have the model calculate zeta and u(z)? I thought about using ana_m2clima and nudging, but it is unclear how to have a time dependent climatology using ana_m2clima.

Has anyone done such a model? I thought that this would be easy to set up. I am clearly not understanding something about ROMS. Hints and suggestions are appreciated.


Top
 Profile  
Reply with quote  
PostPosted: Fri Feb 19, 2016 7:55 pm 
Offline
User avatar

Joined: Wed Jul 02, 2003 5:29 pm
Posts: 3641
Location: IMS/UAF, USA
It's tricky, for sure. Since you have it periodic in y, you can't support a north-south surface slope in balance with the ubar. This used to be easy with a walled channel and the streamfunction solver.

If you want to apply the ana_m3obc, you will have to pick a different LBC for M3, say RadNud.

Why wouldn't the ana_m2clima be the same as your ana_m2obc? Are you trying to propagate something from one end to the other or is it all simultaneous?


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

All times are UTC


Who is online

Users browsing this forum: No registered users and 1 guest


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