polaris.ocean.rpe.compute_rpe

polaris.ocean.rpe.compute_rpe(mesh_filename, initial_state_filename, output_filenames)[source]

Computes the reference (resting) potential energy for the whole domain

Parameters:

mesh_filename (str) – Name of the netCDF file containing the MPAS horizontal mesh variables
initial_state_filename (str) – Name of the netCDF file containing the initial state
output_filenames (list) – List of netCDF files containing output of forward steps

Returns:

rpe (numpy.ndarray) – the reference potential energy of size Time x len(output_files)