# [ncl-talk] computing tropical cyclone heat potential

Debashis Paul debgeo668 at gmail.com
Fri Feb 28 23:38:05 MST 2020

Thank you Dennis Sir, will go through it.
With regards,
DP
********************************
*Debashis Paul*

*Research Scholar,NIT-RKL*
*Ph:+916002713628*
*E-mail: Debgeo668 at gmail.com <Debgeo668 at gmail.com>*

On Fri, Feb 28, 2020 at 11:59 PM Dennis Shea <shea at ucar.edu> wrote:

> 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/20200229/60d061f5/attachment.html>