[MITgcm-devel] seaice free-drift

[Jean-Michel Campin]

Hi,

I saw what I think is a bug in seaice_freedrift.F (discovered around 
Oct 27 and checked the code again two weeks ago), and so that 
everyone knows:
1) The wind-stress (already computed over seaice) is scaled a 2nd
   time by tmpscal1 = SEAICE_drag/OCEAN_drag
2) The Sea-level slope term is not included in the free-drift
   solution without any obvious reason. Using FORCEX0 & FORCEY0
   instead of TAUX,TAUY would fix this.

Cheers,
Jean-Michel