[MITgcm-devel] seaice: variable TICE?

[Martin Losch]

Hi there,

can anyone think of a reason (or show me code) that makes the global variable “TICE” in pkg/seaice necessary? Grepping TICE gives:
>> grep -in ‘TICE(‘ pkg/seaice/*.F
cost_ice_test.F:183:     &                 (TICE(i,j,bi,bj) - 273.15 _d 0 + 1.96 _d 0 ) *
seaice_check_pickup.F:115:                TICES(i,j,k,bi,bj) = TICE(i,j,bi,bj)
seaice_diagnostics_state.F:72:            TICE(I,J,bi,bj) = 0. _d 0
seaice_diagnostics_state.F:80:             TICE(I,J,bi,bj) = TICE(I,J,bi,bj)
seaice_growth.F:739:CADJ STORE tice(:,:,bi,bj)
seaice_growth.F:759:           TICE(I,J,bi,bj) = ZERO
seaice_growth.F:779:           TICE(I,J,bi,bj) = ZERO
seaice_growth.F:842:           TICE(I,J,bi,bj) = TICE(I,J,bi,bj)
seaice_growth.F:845:           TICE(I,J,bi,bj) = TICE(I,J,bi,bj)
seaice_growth.F:1358:c          IF (TICE(I,J,bi,bj) .LT. TMIX) THEN
seaice_init_varia.F:143:          TICE(i,j,bi,bj)   = 0. _d 0
seaice_init_varia.F:231:          TICE(i,j,bi,bj)=273.0 _d 0
seaice_read_pickup.F:199:              TICES(i,j,k,bi,bj) = TICE(i,j,bi,bj)
seaice_reg_ridge.F:222:           TICE(I,J,bi,bj)=celsius2K
seaice_reg_ridge.F:343:           TICE(I,J,bi,bj)=tmpscal1itd(i,j)/tmpscal2itd(i,j)
seaice_reg_ridge.F:345:           TICE(I,J,bi,bj)=celsius2K

as far as I can see, TICE is not used or can be replaced by TICES(i,j,1,bi,bj) or local variables.
If you agree with me, I will remove the global variable, because it may be tricking TAF (and it is not necessary). In fact, removing TICE fixes the NaN problem discussed here: <http://mitgcm.org/pipermail/mitgcm-devel/2014-April/006062.html>

Martin