[MITgcm-devel] problem with exch1_rl and seaice_read_pickup

[Martin Losch]

Hi there,

I find that a restart with pkg/seaice, and SEAICE_MULTICATEGORY undefined (i.e. SEAICE_multDim = 1) does not work; here's the error message:

STOP  S/R EXCH1_RL: myNz  .GT. MAX_NR_EXCH 

When I check the values: MAX_NR_EXCH = 4 and myNz = 7 (because in seaice_read_pickup, there's a
        CALL EXCH_3D_RL   ( TICES, MULTDIM, myThid )
with PARAMETER ( MULTDIM = 7 ).

Obviously I have Nr = 3, so that MAX_NR_EXCH = 4 in conflict with MULTDIM  = 7.
Having an ice model with only a 3layer ocean (or a 1layer ocean) is not very realistic, but for tests (of the seaice model) very convenient.

What's a good solution to this? exchange TICES only with SEAICE_multDim>1 (will not work either) or change the definition of MAX_NR_EXCH to MIN(Nr+1, MULTDIM) (will not work in a parameter statement, will it?)

Martin