[ncl-talk] computing tropical cyclone heat potential

[Dennis Shea]

I think I answered a similar question from  you ? 6+ ? months ago.
---
Attached is a library of ocean related functions.
*The library is not supported nor fully tested.*
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
With a minor modification to the input temperature [degC], the function
*heat_storage_ocn* will do what you request.
It calculates: Heat_Storage=ρCp∫Tdz (integral from 0 to z)

It can calculate TCHP [=ρCp∫(T−26)dz ]  by inputting  't-26' rather than 't'

   t = where(t.gt.26, t-26.0, 0.0)

Pay attention to units.

===
function *heat_storage_ocn*(t, rho, dz, opt)
;###################################
; Not fully tested: Use with caution
;###################################
;
; https://en.wikipedia.org/wiki/Ocean_heat_content
;
; Ocean heat storage is a non-linear quantity.
; Hence, it should be calculated using high-frequency data.
; However, if 'rho' is considered invariant over 'time', then
;          heat_storage(t(time))
;
; The density of seawater is about 1025 kg/m^3
; The specific heat is about 3850 J/(kg C)
;
; t   - temperatue (C)        ; (time,lev,nlat,nlon)
; rho - density (kg/m3)       ; same shape as 't'
; dz  - layer thickness (m)   ; dz[lev] or (time,lev,nlat,nlon)
; opt - set to False ... not used
=======================================

Good luck

On Fri, Feb 28, 2020 at 5:19 AM Debashis Paul via ncl-talk <
ncl-talk at ucar.edu> wrote:

> Dear NCL folks,
>                Has anyone calculated TCHP, if done can you please how was
> it computed in ncl. TCHP is computed by the following given formula:
>
> TCHP=ρCp∫(T−26)dz (integral from 0 to z)
>
> Help will be much appreciated.
> With regards,
>
> *AJAY*
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200228/0bf6fedc/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ocean_lib.ncl
Type: application/octet-stream
Size: 24904 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200228/0bf6fedc/attachment.obj>