Ocean Modeling Discussion

ROMS/TOMS

Search for:
It is currently Sat Jun 15, 2019 10:52 pm




Post new topic Reply to topic  [ 2 posts ] 

All times are UTC

Author Message
PostPosted: Fri Dec 05, 2003 6:13 am 
Offline
Site Admin
User avatar

Joined: Wed Feb 26, 2003 4:41 pm
Posts: 1076
Location: IMCS, Rutgers University
Dear all,

We found a bug in ROMS horizontal viscosity for momentum. This only affect applications with explicit horizontal viscosity. Actually, any value of viscosity chosen does not seem to do anything and the simulation remains inviscid.

This bug has been in ROMS for very long time. We did not noticed it before because by default ROMS uses 3rd-order upstream bias horizontal advection which includes implicit hyper-viscosity in its leading truncation term. In the past, all the horizontal dissipation was mainly specified for the tracer equations.

The problem was that the horizontal stress tensor terms were not included in the residual vertical-integrated terms rufrc and rvfrc used in the split-explicit coupling between baroclinic and barotropic momentum equations. Therefore, the barotropic equations remained horizontally inviscid and the net viscosity was canceled during the barotropic/baroclinic coupling. We corrected this bug and are now getting the correct damping of barotropic Rossby waves and viscous boundary layers.


The next release of ROMS will include the corrections for this bug. In the mean time, get the following tar file:

http://www.myroms.org/links/Hviscosity.tar.gz

and substitute the following files with the new ones:

Code:
    rhs3d.F
    uv3dmix2_s.F
    uv3dmix2_geo.F
    uv3dmix4_s.F
    uv3dmix4_geo.F
    MakeDepend


Notice that in rhs3d.F the CPP option DIAGNOSTIC has been renamed to TS_FIXED. The new version of the code has a new NetCDF output file for model diagnostics like term balances and budgets. I am using the CPP flag DIAGNOGTICS for this type of new output. The model diagnostics will be available in the future.

Many thanks to Genta Mizuta (WHOI) for pointing out this bug.

Good Luck,

Hernan G. Arango
arango@imcs.rutgers.edu


Top
 Profile  
Reply with quote  
 Post subject:
PostPosted: Thu Jan 08, 2004 11:24 pm 
Offline
User avatar

Joined: Tue Jul 01, 2003 4:12 am
Posts: 515
Location: NIWA
Hi Hernan

I added your files to my ROMS directory and tried to rebuild a ROMS job. The make process failed because a couple of files were missing: mod_eclight.f90 and mod_diags.f90.

Perhaps a new release is in order :)


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 5 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