[MITgcm-devel] 3 time levels in seaice

[Martin Losch]

Jean-Michel, Patrick

I have started to get rid of the 3 time levels in pkg/seaice.
So far I am still using the old 3d fields in the code, but only  
level1. This works fine; the next step will be to replace uice(:,:, 
3,:,:) etc with uice(:,:,:,:), but that will require changing the  
pickup files. I am not sure how to do that in order to be backward  
compatible. Is there a way to tell from the meta files? Should I  
introduce a flag useSeaicePickupPreCkpt?? and write code for these  
different cases?

Further, the biggest problem was dealing with the old advect.F and  
diffus.F routines. The most important application for these is the  
adjoint (because the flux-limited schemes are not yet differentiated).  
I had to change them quite a bit (among that removing the leap-frog  
code that hopefully nobody uses), so i am afraid that I might have  
broken the adjoint.

Should I make a tag, check in the stuff that I have now, and let  
Patrick fix the adjoint before we worry about the pickup files?

Martin