[MITgcm-cvs] MITgcm/pkg/seaice CVS Commit

[Martin Losch]

Update of /u/gcmpack/MITgcm/pkg/seaice
In directory forge:/tmp/cvs-serv26984/pkg/seaice

Modified Files:
	seaice_calc_lhs.F seaice_calc_rhs.F seaice_dynsolver.F 
	seaice_evp.F seaice_lsr.F seaice_preconditioner.F 
Log Message:
introduce new parameter SEAICEscaleSurfStress (default = .FALSE.)
  - if TRUE scale ice-ocean and ice-atmosphere stress acting on sea ice
    by AREA according to Connelly et al. (2004), JPO.
  - For EVP and the LSR solver, the implicit part of the drag term is 
    the only non-zero term in the denominator when there is no ice
    (which prevented running the model with zero drag). If
    SEAICEscaleSurfStress = .TRUE., this results in a division by zero
    (or zero main diagonals BU/BV) which needs to be caught. A practical
    consequence is that, for open water, the momentum equation reduce
    to trivially 0 = 0 + 0.