[MITgcm-support] Lock exchange problem with the MITgcm

Flavien Gouillon gouillon at coaps.fsu.edu
Thu Oct 29 15:03:51 EDT 2009

*Dear all,*

*I am starting with the MITgcm but have experience with ROMS and HYCOM. 
I am running a similar simulation of the well known lock exchange 
problem by Haidvogel and Beckmann (and more recently Burchard and 
Rennau, 2009). In a closed flat 2D (x-z) basin (20m depth), I have 2 
initial densities, 5kg/m3 difference (salinity is passive), separated by 
a vertical wall that I removed at t=0. I turn off all explicit mixing 
and diffusion, lateral and bottom free-slip. I have the latest version 
of the MITgcm and run in a hydrostatic fashion. Resolution is 500m (cst) 
in dx (128pts total points) and 40 levels (dz=0.5m, cst). It seems to 
work well but I have few questions since I am quite unfamiliar with this 

* *- *It seems that, for any advection scheme I choose, the MITgcm is 
quite “dispersive” (by dispersive I mean that oscillations appear at the 
interface of the 2 densities) if I do not prescribe high enough 
viscosity (1e-3 for viscA[hz] and 1e4 for viscA4). I did not notice this 
using ROMS. Is there any documentation on this that I could read or 
anybody has some info to give me? I have attached 2 figures (run with 
the 1st order upwind) for the case of no viscosity (all viscosity terms 
at 0) and viscosity at 1e-3 as well as my data file.*

* *- *The other question is: when running with the 3^rd order upwind 
advection scheme, my density profile has, in some places, a persistent 
non-monotic (lighter to denser) feature (unstable conditions, denser 
water above a lighter one). I thought that running with a hydrostatic 
model would prevent such thing and will not allow any ‘overturning’ (by 
quickly mixing any unstable profile) or will blow-up. Is the MITgcm 
capable of handling such thing even running with the hydrostatic option 
and without any mixing?*

* *- *If anybody tried to reproduce the lock exchange problem from 
Haidvogel and Beckmann (or any other lock exchange prob) I would greatly 
appreciate any advices/info/etc...!*

*I hope I was clear enough, please don't hesitate to ask me if I need to 
provide more info, Thanks in advance for all your help, *

*Flavien *

Flavien Gouillon
GRA, Physical Oceanography

-------------- next part --------------
A non-text attachment was scrubbed...
Name: MIT_no_viscosity.jpg
Type: image/jpeg
Size: 122369 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20091029/3eacd956/attachment-0002.jpg>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: MIT_viscosity.jpg
Type: image/jpeg
Size: 85730 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20091029/3eacd956/attachment-0003.jpg>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: data
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20091029/3eacd956/attachment-0001.el>

More information about the MITgcm-support mailing list