<div dir="ltr"><div>I think I answered a similar question from  you ? 6+ ? months ago.</div><div>---<br></div><div>Attached is a library of ocean related functions. <b><br></b></div><div><b><span style="color:rgb(255,0,0)">The library is not supported nor fully tested.</span></b></div><div>^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^<br></div><div>With a minor modification to the input temperature [degC], the function <b>heat_storage_ocn</b> will do what you request. <br></div><div>It calculates: Heat_Storage=ρCp∫Tdz (integral from 0 to z)</div><div><br></div><div><div>It can calculate TCHP [=ρCp∫(T−26)dz ]  by inputting  't-26' rather than 't'<br></div><div><br></div><div>   t = where(t.gt.26, t-26.0, 0.0)<br></div></div><div><br></div><div>Pay attention to units.</div><div><br></div><div>===</div><div>function <b>heat_storage_ocn</b>(t, rho, dz, opt)<br>;###################################<br>; Not fully tested: Use with caution<br>;###################################<br>;<br>; <a href="https://en.wikipedia.org/wiki/Ocean_heat_content">https://en.wikipedia.org/wiki/Ocean_heat_content</a><br>;<br>; Ocean heat storage is a non-linear quantity.<br>; Hence, it should be calculated using high-frequency data.<br>; However, if 'rho' is considered invariant over 'time', then<br>;          heat_storage(t(time))<br>;<br>; The density of seawater is about 1025 kg/m^3 <br>; The specific heat is about 3850 J/(kg C) <br>;<br>; t   - temperatue (C)        ; (time,lev,nlat,nlon)<br>; rho - density (kg/m3)       ; same shape as 't'<br>; dz  - layer thickness (m)   ; dz[lev] or (time,lev,nlat,nlon)<br>; opt - set to False ... not used</div><div>=======================================</div><div><br></div><div>Good luck<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Feb 28, 2020 at 5:19 AM Debashis Paul via ncl-talk <<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:monospace,monospace">Dear NCL folks,</div><div class="gmail_default" style="font-family:monospace,monospace">               Has anyone calculated TCHP, if done can you please how was it computed in ncl. TCHP is computed by the following given formula:</div><div class="gmail_default" style="font-family:monospace,monospace"><br></div><div class="gmail_default" style="font-family:monospace,monospace">TCHP=ρCp∫(T−26)dz (integral from 0 to z)<br></div><div class="gmail_default" style="font-family:monospace,monospace"><br></div><div class="gmail_default" style="font-family:monospace,monospace">Help will be much appreciated.</div><div class="gmail_default" style="font-family:monospace,monospace">With regards,</div><div class="gmail_default" style="font-family:monospace,monospace"><br></div><div class="gmail_default"><font face="monospace" color="#0000ff"><b><i>AJAY</i></b></font></div><div><div dir="ltr"><div dir="ltr"></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div>