_RL lhSublim
lhSublim = SEAICE_lhEvap + SEAICE_lhFusion
latentHeatFluxMaxMult(I,J,IT) = lhSublim*recip_deltaTtherm *
latentHeatFluxMax(I,J) = lhSublim * recip_deltaTtherm *
DIAGarrayA(I,J) = a_FWbySublim(I,J) * lhSublim
D1I=SEAICE_dalton*lhSublim*SEAICE_rhoAir
FWsublim(I,J) = F_lh(I,J)/lhSublim