[MITgcm-support] Spinup in a 2-D channel surface waves?

Dmitri Leonov dleonov at u.washington.edu
Thu Oct 5 16:56:38 EDT 2006


Does your topography go all the way to the surface? (Otherwise your 
channel will be periodic in the cross-channel direction).
That's the only hint I can think of right now.

Dmitri

Jody Klymak wrote:
>
> Hello,
>
> I want to perform spinup of a 2-D channel to along-channel winds for a 
> class I'm teaching.  While I'm sure the MITgcm is overkill for this, 
> I'm at least vaguely familiar with it.
>
> I have the winds blowing, the channel seems to behave, and an Ekman 
> layer seems to form.  However, the surface pressure is not behaving 
> and I'm finding the docs on the surface pressure solver opaque enough 
> that I thought I'd ask for help here.
>
> My vertical resolution near the surface is 0.75m, for 40 bins.  The 
> total water depth is as deep as 200 m, with sloping sides in a 20 
> km-wide channel.  dx = 200 m.
>
> What happens is I blow a south wind, the Ekman layers pile water up in 
> the west and move it away from the East. However, the depth anomalies 
> do not seem to get large enough, and the ones that do form
> do not propagate more than a couple of grid cells away at anything 
> like the surface wave speed.  For instance, the cross-channel Ekman 
> flux reaches 1 m^2/s, which should be changing the surface tilt by 
> something like 36 cm over an hour. However, after 7 h or so (in sim 
> time), I still have 10 cm-high anomalies in the grids adjacent to the 
> edge of the channel, and nothing more than 3 grid cells away.
>
> Any hints what I've done wrong or suggestions for getting the 
> barotropic pressure gradient to set up properly across the channel 
> would be very helpful!
>
> Thanks a lot,  Jody
>
> In "data" I hope the relevant parameters are:
>
> viscAz=1.E-3,
> viscAh=1.E-3,
> no_slip_sides=.FALSE.,
> no_slip_bottom=.TRUE.,
> viscA4=0.E12,
>
> gBaro=9.81,
> rigidLid=.FALSE.,
> implicitFreeSurface=.TRUE.,
>
> nonHydrostatic=.FALSE.,
> hFacMin=0.5,
> implicSurfPress=0.5,
> implicDiv2DFlow=0.5,
> nonlinFreeSurf=3,
> hFacInf=0.1,
> hFacSup=2.0,
> exactConserv=.TRUE.,
>
>
>
>
> -- 
> Jody Klymak
> http://web.uvic.ca/~jklymak/
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support


-- 
------------------------------------------------
Dmitri A. Leonov
Research Associate, Postdoctoral
University of Washington, School of Oceanography





More information about the MITgcm-support mailing list