[MITgcm-support] Turbulent Ice-Ocean heat flux

[Martin Losch]

Alexander, which version of the code do you use (at the top of seaice_growth.F you'll find the Header-line, the most recent version is v 1.173)?

Martin
On Jul 31, 2012, at 6:49 PM, Alexander Slavin wrote:

> 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/m2: 
> 
> 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
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support