<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body dir="auto">Hi Peter,<div><br></div><div>In case you use python, you may find this notebook useful</div><div><a href="https://github.com/miniufo/notebooks/blob/master/Budget%20in%20MITgcm%20Part%20I.ipynb">https://github.com/miniufo/notebooks/blob/master/Budget%20in%20MITgcm%20Part%20I.ipynb</a></div><div><br></div><div>This is a summary of my practice to close the heat budget, which is based on a previous notebook by Ryan Abernathy on Pangeo website.  Now the link seems broken.</div><div><br></div><div>The shortwave heating can be calculated offline if you have oceQsw output.  The code snippet<span style="font-size: 17px;"> </span><span style="color: blue; -webkit-text-size-adjust: 100%;"><span style="font-size: 17px;">__sw_fraction</span><span style="font-size: 14px;"> </span></span>is a copy from the model source code, and calculates the short wave penetration fraction to a specific depth.  Hope this helps.</div><div><br><div dir="ltr">Sent from my iPad</div><div dir="ltr"><br><blockquote type="cite">在 2020年8月25日,21:48,Shatwell, Peter A <peter.shatwell12@imperial.ac.uk> 写道:<br><br></blockquote></div><blockquote type="cite"><div dir="ltr">

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


<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>


<div><ctrl_climate_vertical_heat_balance.jpg></div><span>_______________________________________________</span><br><span>MITgcm-support mailing list</span><br><span>MITgcm-support@mitgcm.org</span><br><span>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</span><br></div></blockquote></div></body></html>