[ferret_users] Calculate OHC time series of water column

["'Abhisek Chatterjee' via _OAR PMEL Ferret Users" <ferret_users@xxxxxxxx>]

Hi,

I am trying to calculate time series of OHC(t) of the entire water column from monthly mean temperature (z,t) simulated by the model. I calculated OHC from temperature by the following formula:

let cp= 3850 ! J/kg C

let rho=rho_un(temp,salt,z[gz=temp])

let ohc=rho*cp*temp

let ohc_intz= ohc[z=@din]

But, I think the derived OHC represents a snapshot of a particular month and not the total heat content of a month. Thus, my derived series likely represent only a series of snapshots OHC and not the total OHC.

It will be a great help if someone can help me to create monthly integrated time series of OHC from monthly mean temperature? 

Thanks in advance,

Regards,

Abhisek