[MITgcm-support] Turbulent Ice-Ocean heat flux

[Alexander Slavin]

Hello MITgcm users,

I have a question concerning turbulent Ice-Ocean heat flux.
Do MITgcm have some diagnostics available to output it?

The analytical formula based on bulk formulation can be written as:
Focn = Csens * Cp * Rho * (Ui-Uw)^2 * (Ti -Tw)

In pkg/seaice/seaice_growth.F I have found this variable for turbulent
ocean-ice heat flux converted to meters of potential ice melt:

a_QbyOCN(i,j) = -tmpscal1 * (HeatCapacity_Cp*rhoConst/QI)
&                                 * (theta(I,J,kSurface,bi,bj)-TBC)

After that I can divide it by convertQ2HI to obtain heat flux in W/m^2:

TurbIOCN(I,J) = a_QbyOCN(I,J) / convertQ2HI

Is it correct? If yes, why was the dependence on the relative ocean-ice
velocity (or the Ustar) dropped? I see that Ustar was calculated for the
McPhee option, but this part of the code does not seem to be operational.

Any help that you may provide regarding this issue would be welcome.

Thank you,
Alexander Slavin,
McGill University, Canada

This electronic communication is governed by the terms and conditions at
http://www.mun.ca/cc/policies/electronic_communications_disclaimer_2012.php
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120731/c3cba64b/attachment.htm>