<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">your configuration is different but perhaps you can find some useful suggestions for closing the MITgcm heat budget here:<div class=""><a href="https://ecco-group.org/docs/v4r3_budgets_howto.pdf" class="">https://ecco-group.org/docs/v4r3_budgets_howto.pdf</a></div><div class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Aug 25, 2020, at 6:48 AM, Shatwell, Peter A <<a href="mailto:peter.shatwell12@imperial.ac.uk" class="">peter.shatwell12@imperial.ac.uk</a>> wrote:</div><br class="Apple-interchange-newline"><div class="">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" class="">
<title class=""></title>
<div class="">
<div name="messageBodySection" class="">
<div dir="auto" class="">Hello helpful MITgcm community,<br class="">
<br class="">
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://urldefense.us/v3/__https://doi.org/10.1175/2009JCLI3197.1__;!!PvBDto6Hs4WbVuu7!fkKI3OFZ0-ngQFY6Q7KpZNLELkMDsL8Eo9SVh0ggrz8G8hGTRJkg6lCiW9ItCq5KMmVmZqkujHI$" target="_blank" class="">https://doi.org/10.1175/2009JCLI3197.1</a>), and I am trying
to close the vertical heat balance in the model ocean.<br class="">
<br class="">
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 class="">
<br class="">
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 class="">
<br class="">
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 class="">
<br class="">
Thanks and all the best,<br class="">
<br class="">
Peter</div>
</div>
<div name="messageSignatureSection" class=""><br class="">
<div class="matchFont">--
<div class="">Peter Shatwell</div>
<div class="">PhD student at MPECDT</div>
<div class="">Space and Atmospheric Physics</div>
<div class="">Imperial College London</div>
<div class="">Huxley 714</div>
</div>
</div>
</div>
<span id="cid:B06FF55A51356C47BEF77538D01A4FEF@eurprd06.prod.outlook.com"><ctrl_climate_vertical_heat_balance.jpg></span>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PvBDto6Hs4WbVuu7!fkKI3OFZ0-ngQFY6Q7KpZNLELkMDsL8Eo9SVh0ggrz8G8hGTRJkg6lCiW9ItCq5KMmVmlr1_kpQ$ <br class=""></div></blockquote></div><br class=""></div></body></html>