# Output internal sea ice temperatures from ACCESS-OM2 **URL:** \ **Category:** Sea Ice\ **Tags:** cosima\ **Created:** [20 October 2023 00:16 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531 "2023-10-20T00:16:07Z")\ **Posts on this page:** 11\ **Page:** 1
**Author:** ![kial](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/kial/32/246_2.png) [@kial](https://forum.access-hive.org.au/u/kial)\ **Post date:** [20 October 2023 00:16 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/1 "2023-10-20T00:16:07Z")
Howdy. I’m trying to output the internal sea ice temperatures from the standard config of ACCESS-OM2. The model is running and behaving normally, so I go into the cice\_in.nml file and modify: f\_tinz = ‘x’ to f\_tinz = ‘m’ Now the model isn’t able to get past a month (presumably where it is trying to write out this diagnostic). Has anyone been able to get internal sea ice temperatures out successfully? Thanks!! ---
**Author:** ![sofarrell](https://avatars.discourse-cdn.com/v4/letter/s/ea5d25/32.png) [@sofarrell](https://forum.access-hive.org.au/u/sofarrell)\ **Post date:** [20 October 2023 01:34 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/2 "2023-10-20T01:34:23Z")
Hi @kial we had them switched on in access-cm2 during testing but off in the main runs, they are a little trickier as they had an extra dimension to the output files, I am trying to remember how we got round that for all cases. I think it worked best with a new restart. ---
**Author:** ![kial](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/kial/32/246_2.png) [@kial](https://forum.access-hive.org.au/u/kial)\ **Post date:** [20 October 2023 01:43 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/3 "2023-10-20T01:43:47Z")
Thanks @sofarrell I’ve just seen that @aekiss had both Tinz and Sinz output for a couple of cycles with the IAF in ACCESS-OM2-01. I’m trawling through the cice\_in.nml files for other differences. Is it possible that they only work with ktherm=2 and NOT with ktherm=1? ---
**Author:** ![sofarrell](https://avatars.discourse-cdn.com/v4/letter/s/ea5d25/32.png) [@sofarrell](https://forum.access-hive.org.au/u/sofarrell)\ **Post date:** [20 October 2023 01:56 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/4 "2023-10-20T01:56:03Z")
We had it working for ktherm=1 in the access-cm2 runs ---
**Author:** ![aekiss](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aekiss/32/46_2.png) [@aekiss](https://forum.access-hive.org.au/u/aekiss)\ **Post date:** [20 October 2023 03:40 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/5 "2023-10-20T03:40:38Z")
Yes there’s Tinz, Sinz in the IAF ACCESS-OM2-01: [Data available: 0.1° 1958-2018 ACCESS-OM2 IAF runs – COSIMA](https://cosima.org.au/index.php/2020/07/29/data-available-0-1-1958-2018-access-om2-iaf-run/) - This used ktherm=2 (mushy ice). I thought Sinz is time-invariant with ktherm=1? Tinz, Sinz are much bigger files due to the extra dimension - perhaps there’s a IO bottleneck in your run? You may find this useful for comparing namelists: [GitHub - aekiss/nmltab: Python module and command-line tool to semantically tabulate, diff and superset Fortran namelist files.](https://github.com/aekiss/nmltab) ---
**Author:** ![sofarrell](https://avatars.discourse-cdn.com/v4/letter/s/ea5d25/32.png) [@sofarrell](https://forum.access-hive.org.au/u/sofarrell)\ **Post date:** [20 October 2023 04:38 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/6 "2023-10-20T04:38:52Z")
Yes @Kial would only need to look at Tinz as @aekiss says salinity profile doesn’t change in ktherm =1 option. ---
**Author:** ![kial](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/kial/32/246_2.png) [@kial](https://forum.access-hive.org.au/u/kial)\ **Post date:** [21 October 2023 01:54 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/7 "2023-10-21T01:54:35Z")
Still no luck. I’ve tried switching off all other CICE diagnostics, and can successfully toggle these on/off (including putting in additional 4D diags), but any run without f\_tinz = ‘x’ gets stuck at the point where it is trying to output these internal sea ice temperatures. The config is the latest 1deg\_jra55\_ryf if that helps… ---
**Author:** ![aekiss](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aekiss/32/46_2.png) [@aekiss](https://forum.access-hive.org.au/u/aekiss)\ **Post date:** [23 October 2023 04:55 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/8 "2023-10-23T04:55:05Z")
That’s odd - so even `f_tinz = ‘m’` fails? ---
**Author:** ![aekiss](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aekiss/32/46_2.png) [@aekiss](https://forum.access-hive.org.au/u/aekiss)\ **Post date:** [23 October 2023 04:57 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/9 "2023-10-23T04:57:55Z")
Note that the index order is incorrect with Tinz, Sinz: [diagnostic bug in CICE when outputting 4d fields · Issue #62 · COSIMA/cice5 · GitHub](https://github.com/COSIMA/cice5/issues/62#issuecomment-1097556730) ---
**Author:** ![aekiss](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/aekiss/32/46_2.png) [@aekiss](https://forum.access-hive.org.au/u/aekiss)\ **Post date:** [24 October 2023 02:43 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/10 "2023-10-24T02:43:39Z")
OK I can reproduce this problem. I’ve opened an issue: [Can't output Tinz · Issue #72 · COSIMA/cice5 · GitHub](https://github.com/COSIMA/cice5/issues/72) ---
**Author:** ![kial](https://sea2.discourse-cdn.com/flex020/user_avatar/forum.access-hive.org.au/kial/32/246_2.png) [@kial](https://forum.access-hive.org.au/u/kial)\ **Post date:** [24 October 2023 04:32 UTC](https://forum.access-hive.org.au/t/output-internal-sea-ice-temperatures-from-access-om2/1531/11 "2023-10-24T04:32:15Z")
Thank you @aekiss !!