<html xmlns="http://www.w3.org/1999/xhtml">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=utf-8">

<title></title>

</head>

<body>

<div name="messageBodySection">

<div dir="auto">Hello helpful MITgcm community,<br>

<br>

I’m using the DDrake coupled configuration with a cube-sphere (C24) grid with 15 vertical levels (more details in Ferreira et al. 2010, <a href="https://doi.org/10.1175/2009JCLI3197.1" target="_blank">https://doi.org/10.1175/2009JCLI3197.1</a>), and I am trying

 to close the vertical heat balance in the model ocean.<br>

<br>

I have computed vertical profiles of the advective vertical heat flux (diagnostic ADVr_TH) and the total diffusive vertical heat flux (diagnostics DFrE_TH and DFrI_TH). These advective and diffusive fluxes almost exactly balance throughout the water column

 (see attached figure), consistent with the traditional ‘upwelling-diffusion’ vertical balance of Munk ’66 and others (I’ve plotted it such that a positive flux is downwards i.e. a warming effect).<br>

<br>

However, there is a strong departure from this balance in the top ~130 m of the ocean, indicating strong net upward cooling. Apparently penetrating shortwave radiation in the model could balance this cooling, but I am not sure how to properly confirm this.

 Are there existing model diagnostics to plot the shortwave penetration? Or could I estimate it from other outputs? <br>

<br>

Or is there something else that could provide a warming effect near the ocean surface? Or should I expect this balance to even hold at all?<br>

<br>

Thanks and all the best,<br>

<br>

Peter</div>

</div>

<div name="messageSignatureSection"><br>

<div class="matchFont">--

<div>Peter Shatwell</div>

<div>PhD student at MPECDT</div>

<div>Space and Atmospheric Physics</div>

<div>Imperial College London</div>

<div>Huxley 714</div>

</div>

</div>

</body>

</html>