# Ocean Modeling Discussion

ROMS/TOMS

Search for:
 It is currently Thu Apr 19, 2018 5:38 pm

 Page 1 of 1 [ 5 posts ]
 All times are UTC
Author Message
 Post subject: Bottom stress in ROMSPosted: Mon Nov 24, 2003 8:02 pm

Joined: Tue Nov 25, 2003 3:36 pm
Posts: 14
Location: CSIRO
Hernan,

I am a little confused with the bottom stress calculation in ROMS. There are 3 ways for this calculation using either linear drag coefficient, quadratic bottom drag coefficient, or logarithmic formulation via bottom roughtness. Only one method should be chosen. I find bottom stress is calculated in two places. One is in subroutine set_vbc.F, where bottom stress is proportional to the velocity at the lowest grid near the bottom. And antother is in subroutine step2d.F. Here bottom stress is proportional to the depth integrated velocity, which is true if only barotropic mode is simulated, and obviously it is different from that calculated in set_vbc.F in a 3D simulation. Bottom stress should be the same in both barotropic mode and baroclinic mode in a 3D simulation. Therefore it makes sense to me only in these two situations.

1) Only 2D not 3D simulation is done. In this case, set_vbc.F is not called. bottom stress is only calculated in step2d.F.

2) Choose logarithmic formulation for 3D simulation, then bottom stress is calculated only in set_vbc.F, since step2d.F does not compute bottom stress when linear or quadratic drag coefficient is zero.

If we simulate 3D circulation and choose linear or quadratic bottom drag coefficient instead of bottom roughness, we will run into this situation that bottom stress is different in barotropic and baroclinic mode.

I may misunderstand the code. If I am wrong, please correct me.

Liejun

Top

 Post subject: Posted: Mon Nov 24, 2003 8:05 pm

Joined: Wed Feb 26, 2003 4:41 pm
Posts: 1022
Location: IMCS, Rutgers University
Hi Liejun,

Yes, there are several ways to add bottom stress in ROMS:

Code:
1) Linear          RDRG > 0,  RDRG2 = 0,   Zob = 0
2) quadratic       RDRG = 0,  RDRG2 > 0,   Zob = 0
3) law of wall     RDRG = 0,  RDRG2 = 0,   Zob > 0
4) bottom boundary layer model

Notice that in schemes 1 to 3, only one value is greater than zero. This is a must in order to have the appropriate code executed.

The logarithmic law of the wall is only apply to 3D configurations. It does not make sense, I think, to apply this to shallow-water applications neither to barotropic components "ubar" and "vbar". It is only applied to bottom baroclinic velocities!!! In 3D applications, the logarithmic bottom stress or other BBL formulations enter in the 2D equations (step2d) in the vertically integrated residual terms "rufrc" and "rvfrc".

Do this makes sense?

Good luck,

Hernan G. Arango
arango@imcs.rutgers.edu

Top

 Post subject: Bottom VelocitiesPosted: Tue Jan 30, 2007 2:55 am

Joined: Tue Mar 09, 2004 3:21 pm
Posts: 15
Location: University of New South Wales, Australia
Hi,

In a 3D application with quadratic bottom friction what do the velocities at the lowest level in the model represent?

Are they a shear velocity (u*) or the velocity at the top of an assumed constant stress bottom boundary layer or something different?

cheers,
Deborah

Top

 Post subject: Posted: Tue Jan 30, 2007 8:36 pm

Joined: Tue Jul 01, 2003 4:12 am
Posts: 503
Location: NIWA
Quote:
If we simulate 3D circulation and choose linear or quadratic bottom drag coefficient instead of bottom roughness, we will run into this situation that bottom stress is different in barotropic and baroclinic mode.

I may misunderstand the code. If I am wrong, please correct me.

I think it does all work out correctly in the end, owing to the mysterious workings of the 2D/3D coupling. But I can't guarantee that. A while back I ran a test case with a doubly periodic domain initialized with a vertically uniform flow and confirmed that the deceleration of the flow matched the expected values. I could dig that up for you if you want to verify it.

Quote:
In a 3D application with quadratic bottom friction what do the velocities at the lowest level in the model represent?

Are they a shear velocity (u*) or the velocity at the top of an assumed constant stress bottom boundary layer or something different?

Velocity at top of an assumed boundary layer.

Top

 Post subject: Posted: Wed Jan 31, 2007 10:37 am

Joined: Wed Mar 26, 2003 4:14 pm
Posts: 11
Location: DOI/MMS
For a 3D application, the bottom stress is computed in set_vbc.F and is not computed in step2d.F at all.

From step2d.F:
Code:
# ifndef SOLVE3D
!
!-----------------------------------------------------------------------
!-----------------------------------------------------------------------
!

It enters the barotropic mode through the magic of 2D/3D coupling (as does surface momentum stress).

This is one of the things that took me quite a while to completely figure out when I was working on the diagnostics.

Top

 Display posts from previous: All posts1 day7 days2 weeks1 month3 months6 months1 year Sort by AuthorPost timeSubject AscendingDescending
 Page 1 of 1 [ 5 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 forumYou cannot reply to topics in this forumYou cannot edit your posts in this forumYou cannot delete your posts in this forumYou cannot post attachments in this forum

Search for:
 Jump to:  Select a forum ------------------ News, Events & Job Opportunities    Meetings/Workshops    Job Opportunities    Ocean News ROMS/TOMS    ROMS Adjoint    ROMS Benchmarks    ROMS Bugs    ROMS Discussion    ROMS Documentation    ROMS Ecosystem    ROMS FAQ    ROMS Ice    ROMS Information    ROMS Installation    ROMS Known Problems    ROMS Messages    ROMS Problems    ROMS Releases    ROMS Results    ROMS Sediment    ROMS Source    ROMS Tips    ROMS Tools and Techniques    ROMS Trivia    ROMS Usage    ROMS Webinar    ROMS Wish List ROMS/TOMS Applications    User Applications    Adriatic Sea Ocean Modeling    Ocean Modeling FAQ