Re: [ferret_users] Calculate OHC time series of water column

[Sudheer Joseph <sjo@xxxxxxxxxxxxx>]

Hi Abhishek,

If you are looking monthly integrated ohc ( sum of daily ohcs) that cannot be done from monthly data.( it is similar to expecting vertically integrated ohc when one has only one layer's temperature)  If you use monthly data it can give only monthly mean ohc or 15 day snapshot based on how monthly data was formed from daily data.  
Sudheer
Sent from my iPhone

On 27-Jul-2017, at 6:48 PM, saurabh rathore <rohitsrb2020@xxxxxxxxx> wrote:

Dear Abhishek sir 

It is always good to use potential temperature for ohc calculation. In which rho.cp.T is the ocean heat J/m^3 and it's integration in z will give the OHC J/m^2. I think you did the right thing except the use of potential temperature. If it is a 4d data in x,y,z,t then ohc[x=@din,y=@din,z=@din] will give you the time series of OHC in joules. 

I hope this will help. Please let us know if it solved. 

Regards Saurabh 

On Jul 27, 2017 10:31 PM, "'Abhisek Chatterjee' via _OAR PMEL Ferret Users" <ferret_users@xxxxxxxx> wrote:

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