<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">Wow, thank you so much for pointing me to this link, Yu Kun! And thanks also to Martin Losch for pointing me to swfrac.F in the source code. <br>
<br>
I’ve now calculated the shortwave penetration offline, and the balance is looking a lot better! (figure attached) There’s still a small positive bump, but perhaps this will be closed by other terms as outlined in Yu Kun’s python notebook below.<br>
<br>
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 name="messageReplySection">On 25 Aug 2020, 18:13 +0100, Yu-Kun Qian <qianyk@mail3.sysu.edu.cn>, wrote:<br>
<blockquote type="cite" style="border-left-color: grey; border-left-width: thin; border-left-style: solid; margin: 5px 5px;padding-left: 10px;">
<table class="MsoNormalTable" style="mso-cellspacing: 0cm; mso-yfti-tbllook: 1184; mso-padding-alt: 0cm 0cm 0cm 0cm" cellspacing="0" cellpadding="0" border="1">
<tbody>
<tr style="mso-yfti-irow: 0; mso-yfti-firstrow: yes; mso-yfti-lastrow: yes">
<td style="BORDER-TOP: #f0f0f0; BORDER-RIGHT: #f0f0f0; BACKGROUND: #ffffcc; BORDER-BOTTOM: #f0f0f0; PADDING-BOTTOM: 0cm; PADDING-TOP: 0cm; PADDING-LEFT: 0cm; BORDER-LEFT: #f0f0f0; PADDING-RIGHT: 0cm" wrap="">
<p class="MsoNormal" style="MARGIN: 0cm 0cm 0pt"><span style="FONT-SIZE: 12pt; FONT-FAMILY: "Arial",sans-serif"><font size="2">This email from qianyk@mail3.sysu.edu.cn originates from outside Imperial. Do not click on links and attachments unless you recognise
the sender. If you trust the sender, add them to your <a href="https://spam.ic.ac.uk/SpamConsole/Senders.aspx">
safe senders list</a> to disable email stamping for this address.</font></span></p>
</td>
</tr>
</tbody>
</table>
<p></p>
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">
<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>
</blockquote>
</div>
</body>
</html>