[MITgcm-support] Re: pickup problem

chris hill cnh at mit.edu
Mon Apr 11 23:26:15 EDT 2005


Hi Tom,

  Have you tried the Smag scheme?
  there is an example in verification/MLAdjust/input.0.smag/data
  Baylor - what do you think?

Chris

Tom Haine wrote:
> Dear Alistair and Jean-Michel, 
> 
> I think I've found the source of the pickup problem with my version of
> the new horizontal viscosity: At pickup I removed the mean SSH from the
> SSH field (to avoid problems with very stiff barotropic modes when
> assimilating with an earlier version of this code - I should have cut
> this a while ago, but I forgot: sorry if I wasted your time, J-M.). 
> 
> Without this adjustment (i.e., the regular code) I can integrate past
> the pickup for several hundred more steps. Eventually, there's another
> divergence however:
> 
> WARNING: r*FacC < hFacInf at: i,j,bi,bj,Thid,Iter=  71  20   1   1   1
> 5977
> rStarFac,H,eta = -0.025522  5.000000E+00 -5.127608E+00
> STOP in CALC_R_STAR : too SMALL rStarFacC !
> 
> in STDERR.0032. Again, I guess this is the Leith scheme diverging
> despite the diagnostics about rStar (right?). I suppose my flow has
> become too energetic (I see more or less monotonic increase in KE prior
> to the crash). How should I control this?  Should I just increase
> viscC4Leith (currently set to 1.0)?
> 
> Thanks again for your help. I agree with Dimitris - this is a fabulous
> new scheme for the MITgcm.
> Tom.
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
> 




More information about the MITgcm-support mailing list