[MITgcm-support] Hydrostatic as grid approaches isotropic?

Jody Klymak jklymak at uvic.ca
Tue Jul 20 12:27:52 EDT 2010


Hi again,

Any thoughts on this?  Does anyone have a working near-isotropic setup  
that they could share the data file for?

Thanks,  Jody

On Jul 12, 2010, at  15:25 PM, Jody Klymak wrote:

>
> Hi All,
>
> I have been banging my head against a wall with a simulation where  
> dz \approx 2 m and dx \approx 5 m (there is some telescoping away  
> from my topography and with depth).  I was running with  
> nonHydrostatic=.FALSE.  I expected mode-1 internal waves to zip  
> along at 1 m/s, but the fastest I was getting was 0.1 m/s.  I  
> thought it was my stratification etc etc.
>
> Then, or fun, I changed dx to 50 m.  Now I get 1 m/s mode-1 waves  
> and all is fine.  Is this a known feature of the code?  No doubt an  
> effect I should have learnt about somewhere...
>
> I include the pressure solving parts of my data file; maybe I'm  
> doing something wrong in here.  Note I use linear free surface  
> because I use Orlanski BCs. (ahem, though I really don't as I have a  
> big sponge in there.  I really only use the Orlanski code to set  
> OBNu, OBNv, etc, so if going to a non-linear surface works with OBCs  
> then I can make the effort).
>
> Thanks for any pointers....
>
> Cheers,  Jody
>
>
> gravity=9.81,
> gBaro=9.81,
> rigidLid=.FALSE.,
> implicitFreeSurface=.TRUE.,
> eosType='LINEAR',
> nonHydrostatic=.FALSE.,
> # minimum cell fraction.  This reduces steppiness..
> hFacMin=0.1,
> # implicSurfPress=0.5,
> # implicDiv2DFlow=0.5,
> # nonlinFreeSurf=3,
> # hFacInf=0.2,
> # hFacSup=1.8,
> #  exactConserv=.TRUE.,
>  implicitDiffusion=.TRUE.
>  implicitViscosity=.TRUE.
> # Superbee on:
> tempAdvScheme=77,
> staggerTimeStep=.TRUE.,
>
>
>
>
> --
> Jody Klymak
> http://web.uvic.ca/~jklymak/
>
>
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

--
Jody Klymak
http://web.uvic.ca/~jklymak/







More information about the MITgcm-support mailing list