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

[Martin Losch]

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