# Question about sea ice related heat fluxes in ACCESS-OM2

**URL:** <https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260>\
**Category:** Sea Ice\
**Created:** [13 April 2026 04:21 UTC](https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260 "2026-04-13T04:21:20Z")\
**Posts on this page:** 7\
**Page:** 1

<div class="post-metadata">

**Author:** ![rmholmes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/rmholmes/32/117_2.png) [@rmholmes](https://forum.access-hive.org.au/u/rmholmes)\
**Post date:** [13 April 2026 04:21 UTC](https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260/1 "2026-04-13T04:21:20Z")

</div>

Hi everyone,

I am trying to confirm that I’m not making a mistake when accounting for the sea ice related heat fluxes, `mh_flux` and `liceht`, in my OM2 mixed layer heat budgets. These diagnostics are somewhat difficult to track through the code since they have different names in CICE and in the coupler - so if someone who knows more about sea ice in OM2 than me can help I’d appreciate some guidance!

Specifically, I can see that the ice-\>ocean volume fluxes are included within [MOM5’s `pme` diagnostic](https://github.com/rmholmes/MOM5/blob/8b7ee7aba6f9cb5eb1b7b66a8ff783b0f15d8dc5/src/mom5/ocean_core/ocean_sbc.F90#L3642), and hence the heat flux into the ocean associated with these volume fluxes (assumed to enter the ocean at the surface ocean temperature) are included in the [`sfc_hflux_pme` heat budget term](https://github.com/rmholmes/MOM5/blob/8b7ee7aba6f9cb5eb1b7b66a8ff783b0f15d8dc5/src/mom5/ocean_core/ocean_sbc.F90#L4838). Hence, `mh_flux` and `liceht` should only contain heat fluxes that are not volume flux related (e.g. they should account only for differences between the melt water and surface ocean temperature, as well as other processes that I don’t understand). Can someone confirm whether that is the case or point me at somewhere which documents these terms from the CICE side in more detail?

Thanks!

---

<div class="post-metadata">

**Author:** ![sofarrell](https://avatars.discourse-cdn.com/v4/letter/s/ea5d25/32.png) [@sofarrell](https://forum.access-hive.org.au/u/sofarrell)\
**Post date:** [13 April 2026 07:30 UTC](https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260/2 "2026-04-13T07:30:35Z")

</div>

Hi @rmholmes I can probably make a quick attempt here with a bit more detail just looking at some basic issues. We can get into the CICE code if we need to.

The liceht code is for the access-cm2 model it has been set to zero for the access-om2 code so it shouldn’t be a factor in your calculations it represents the iceberg flux which we distributed away from the coast in the ACCESS-cm2 code and accounted for the latent heat of melting, in ACCESS-OM2 there is the calv flux term which represented solid runoff but was just added to the liquid runoff and no latent heat was accounted for.

The mh\_flux term can be traced back to the “iohtflx” term in the coupler and then to the “fhocn” term in the CICE code.

There are some old out of date comment statements referencing an fhnet field but that is no longer used in the coupler.

Hope this helps

Siobhan

---

<div class="post-metadata">

**Author:** ![rmholmes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/rmholmes/32/117_2.png) [@rmholmes](https://forum.access-hive.org.au/u/rmholmes)\
**Post date:** [13 April 2026 22:57 UTC](https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260/3 "2026-04-13T22:57:49Z")

</div>

Thanks Siobhan!

So tracking things further following your email and notes:

- In the `namcouple` `mh_flux` comes from `htflx_io` (and `wfiform` comes from `form_io`).
- Then in the [CICE coupling code](https://github.com/COSIMA/cice5/tree/master/drivers/auscom) `cpl_interface.f90` these come from `iohtflx` and `ioform`. These are described in `cpl_arrays_setup.f90`.
- In turn, in `cpl_forcing_handler.f90`, `iohtflx` and `ioform` come from `tiohtflx` and `tioform`, which themselves come from `fhocn` and `fresh` respectively.

`fhocn` is everywhere in the CICE code, so it’s hard to track and I’m not really sure how to proceed. However, from the [CICE model documentation](https://csdms.colorado.edu/w/images/CICE_documentation_and_software_user's_manual.pdf) we have “The portion of the melting potential actually used to melt ice is returned to the coupler in Fhocn” - which seems to support what I said above that it doesn’t include volume flux-related heat fluxes, only including “extra” heat required to melt sea ice. The [comment here](https://github.com/COSIMA/cice5/blob/edcfa6f9c76ed05b63196ce4b5355fa5a8f4fe3a/drivers/auscom/cpl_arrays_setup.F90#L46) also seems to support that. I guess maybe that’s good enough to answer my question…

---

<div class="post-metadata">

**Author:** ![anton](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/anton/32/1836_2.png) [@anton](https://forum.access-hive.org.au/u/anton)\
**Post date:** [13 April 2026 23:21 UTC](https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260/4 "2026-04-13T23:21:16Z")

</div>

I don’t think i can help further with the question, but the model document for CICE5 is here:

> **[cicedoc.pdf](https://github.com/ACCESS-NRI/cice5/blob/master/doc/cicedoc.pdf)**

It is for the standalone / CESM cice and doesn’t contain information re: the OASIS coupling for ACCESS.

---

<div class="post-metadata">

**Author:** ![sofarrell](https://avatars.discourse-cdn.com/v4/letter/s/ea5d25/32.png) [@sofarrell](https://forum.access-hive.org.au/u/sofarrell)\
**Post date:** [14 April 2026 03:19 UTC](https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260/5 "2026-04-14T03:19:57Z")

</div>

HI Ryan thats a pretty good description of how the variables change names through the namcouple, cpl\_arrays\_setup.F90, cpl\_forcing\_handler.F90 subroutines. As is pointed out in the comment statements in cpl\_arrays\_set\_up and in the CICE manual Fhocn is only the part of the heat flux that is used to melt the ice. The total flux is called ‘frzmlt’ in CICE and is not passed to the ocean model. I guess its partly confusing to the ocean community, its historic because the way frazil term of new ice formation is handled in CICE and added in later to the budget depending on which thermodynamic ice scheme is used.

---

<div class="post-metadata">

**Author:** ![rmholmes](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/rmholmes/32/117_2.png) [@rmholmes](https://forum.access-hive.org.au/u/rmholmes)\
**Post date:** [14 April 2026 04:37 UTC](https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260/6 "2026-04-14T04:37:48Z")

</div>

Thanks @sofarrell . Sorry I meant to say that the ocean also has the `frazil_3d` heat flux (although I guess that’s computed internally in MOM???), but that’s definitely not related to volume fluxes so I didn’t mention it.

---

<div class="post-metadata">

**Author:** ![sofarrell](https://avatars.discourse-cdn.com/v4/letter/s/ea5d25/32.png) [@sofarrell](https://forum.access-hive.org.au/u/sofarrell)\
**Post date:** [14 April 2026 09:42 UTC](https://forum.access-hive.org.au/t/question-about-sea-ice-related-heat-fluxes-in-access-om2/6260/7 "2026-04-14T09:42:37Z")

</div>

yes I guess you are using are already using the frazil\_3d heat in the calculation that you are trying to do but are missing the calculation from the congelation ice growth from the ice model side.
