[MITgcm-devel] pkg/seaice pseudo time stepping for LSOR

[Patrick Heimbach]

Hi Martin,

what are max. expected number of pseudo-timesteps.
As you probably guess, your changes are producing some
relatively ugly recomp. of seaice_lsr within adseaice_dynsolver.
For the storing we'll need a hard limit (to be used as
common block size).
What's a good number?

Cheers
-p.

On May 25, 2009, at 5:50 AM, Martin Losch wrote:

> Hi Patrick,
> I have checked in code to do pseudo time stepping for LSOR.  
> Basically, the two calls to seaice_lsr in seaice_dynsolver have  
> been replaced by a loop with a default of 2 iterations, and the  
> copy/interpolation between time levels has been moved to seaice_lsr.
>
> A potential problem for the adjoint are the store directives for  
> uicec and vicec in the 'predictor time step'. I have moved them to  
> seaice_lsr along with the rest of the code and I have deliberately  
> left out the "(icall-1)*ncklev_1 part of the key that is used for  
> uice/vice, because the storing is only applied for icall==1. I am  
> not sure if this will work properly and apologize beforehand for  
> all problems that this may cause.
>
> Martin
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel

---
Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS 54-1518 | 77 Massachusetts Ave | Cambridge MA 02139 USA
FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach