polaris.ocean.conservation.compute_total_energy

polaris.ocean.conservation.compute_total_energy(ds_mesh, ds, model)[source]

Compute the total heat content in an ocean model output file

Parameters:

ds_mesh (xarray.Dataset) – The mesh dataset, containing areaCell
ds (xarray.Dataset) – The dataset containing temperature and a layer thickness for a single time slice

Returns:

total_energy (float) – The total heat content in J over the whole domain