[MITgcm-support] Numerical diffusion

Martin Losch mlosch at awi-bremerhaven.de
Mon Oct 2 04:46:01 EDT 2006

Hi Lars,

now that Chris set me up I have to answer (o:

I don't know much about GOTM, yet, I have tried to use it to  
understand the KPP code in the MITgcm, but that did not work out so far.

In terms of channel experiments: I have played with re-entrant  
channels with wind forcing, buoyancy forcing or temperature/density  
restoring along one boundary. I'll happily make these configurations  
available, although, they are based on 3year-old code an may need  
some adjustments (there are some retired parameters, I guess).

I guess you can prescribe a constant pressure gradient in  
calc_grad_phi_surf.F, but I have never tried this. (I'm not sure if  
this is want you want.)


On 29 Sep 2006, at 10:51, Lars Inge Enstad wrote:

> Hello
> I am currently trying to implement an alternative vertical turbulence
> model in MITgcm (GOTM (General Ocean Turbulence Model Burchard et
> al.), has this been tried before?). My question is about the numerical
> diffusion in the momentum advection scheme in MITgcm. How does this
> depend on the spatial and temporal discretisation? What is the  
> momentum
> advection scheme with the lowest possible numerical diffusion in
> MITgcm?
> One of the test cases I want to verify the code with is channel
> flow. In previous models I have used for channel flow, I have used a
> constant pressure gradient to force the flow. Is there some easy way
> to set up MITgcm this way?
> Best Regards
> Lars Inge Enstad
> -- 
> Lars Inge Enstad
> Post Doc
> Computational Mathematics Unit
> Bergen Center for Computational Science
> UNIFOB/University of Bergen
> Thormøhlensgt. 55
> N-5008 Bergen, Norway
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

More information about the MITgcm-support mailing list