[MITgcm-support] flux limiters and timestep advice...

Jean-Michel Campin jmc at ocean.mit.edu
Thu Apr 9 14:58:32 EDT 2009


Hi Jody,

>From your "data" file, it looks like you would need to turn on
staggerTimeStep (setting staggerTimeStep=.TRUE.,).

We added an error message a year ago (on Apr 08, 2008) for
this case:
> tail -2 STDERR.0000
> (PID.TID 0000.0001) ** WARNING ** GAD_CHECK: potentially unstable time-stepping (Internal Wave)
> (PID.TID 0000.0001) ** WARNING ** GAD_CHECK: need "staggerTimeStep=.TRUE." in "data", nml PARM01

Can you confirm that you see this error message in STDERR.##### 
or does your code is older than that ?

Thanks,
Jean-Michel

On Thu, Apr 09, 2009 at 11:27:49AM -0700, Klymak Jody wrote:
>
> Hi all,
>
> I'm running relatively high-resolution 2-D simulations of flow over  
> topography.  I was relatively happy using 2-nd order centered advection. 
> There was some gridscale noise, but the diffusivity tended to erase it 
> and the major features look fine.  Now I am being a little more fussy and 
> trying to balance energy with dissipation \epsilon = viscAz*(du/dz)^2 and 
> find that my dissipation is scale dependent such that fine-scale models 
> dissipate more than coarse scale.  My supposition is that this is due to 
> the extra gridscale noise.
>
> So I switched to a superbee flux limiting scheme (tempAdvScheme=77).  It 
> certainly removed much of the gridscale noise yet the resolved  
> dissipative features remain the same.
>
> However, after a while I start to get low-mode noise in the solution.   
> When I was using centered 2nd-order I could run 12.5 s timesteps w/o  
> difficulty.  If I run superbee I have to reduce this to 2 s, and still  
> about 20 h I start getting low-mode grid-scale noise.  The noise seems  
> largest near the edges of my topography, and extends through the whole  
> water column.
>
> Please see: http://web.uvic.ca/~jklymak/SuperbeeProblems.pdf
>
> FWIW, these are not huge signals, but completely destroy my energy  
> budgets.
>
> Perhaps my solution is simply to go down to 1-s timesteps (that will be 
> what I will do for now). However, I was curious if there were other  
> suggestions. Should I decrease just the barotropic time step?  Anything 
> else I've set up incorrectly in the attached?
>
> Thanks for any thoughts.
>
> Cheers,  Jody
>


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

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support




More information about the MITgcm-support mailing list