[MITgcm-support] Diagnosing internal ice stress with the seaice pkg

[Leng Hengling]

Dear MITgcm community,

I am studying sea ice-ocean interactions with the seaice pkg. Note that the internal ice stress (divergence) is calculated in seaice_calc_stressdiv.F but not defined as a diagnostic. I have tried adding it to the list of diagnostics, via calling DIAGNOSTICS_ADDTOLIST in seaice_diagnostics_init.F as:
...
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
      diagName  = 'SIfx    '
      diagTitle = 'SEAICE zonal internal stress'
      diagUnits = 'N/m^2           '
      diagCode  = 'UU      M1      '
      diagMate  = diagNum + 2
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
      diagName  = 'SIfy    '
      diagTitle = 'SEAICE meridional internal stress'
      diagUnits = 'N/m^2           '
      diagCode  = 'VV      M1      '
      diagMate  = diagNum
...
and filling the fields in seaice_diagnostics_state.F:
...
       CALL DIAGNOSTICS_FILL(
     I stressDivergenceX, 'SIfx    ',0,1, 0,1,1,myThid)
       CALL DIAGNOSTICS_FILL(
     I stressDivergenceY, 'SIfy    ',0,1, 0,1,1,myThid)
...
The model ran smoothly for 100 days without any errors. However, the ouput internal ice stress (both zonal and meridional components) is zero everywhere during the model period. Other output variables (sea ice concentration, thickness, velocity, etc) look fine. I would appreciate any insights or suggestions on what might be causing the internal ice stress to output as zero.

Best regards,

Hengling Leng
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20241010/8ef9e302/attachment.html>