Thanks Ryan,
I should have been more clear in exactly what I’m trying to accomplish.
I’m trying to compute the mixed layer heat budget point-wise offline from daily average output.
The output I have is MLD, total surface heat flux, u,v,w, and T. BRAN2020 is based on MOM5 on a 0.1x0.1 degree B-grid that differs slightly from the ACCESS-OM-01 grid.
I’ve made two serious attempts, both of which don’t really come close to closing:
- using xgcm to compute udT/dx and vdT/dy on T-points for each vertical layer, followed by vertical integration;
- carefull finite-volumn type calculation, where I considered every T-cell as a control volumn, mapped temperature and u/v to the relevant center cell edges (using xesmf), and computed the heat flux through each cell edge, ie:
advective term: 1/(cell area) * [vT (northern edge) * dx(northern_dege) - vT (southern_edge) * dx(southern_edge) + uT(eastern edge) * dy(eastern_edge)-uT(western edge)*dy_(western edge) ]
The grid terms (area, dx, dy) I get from the grid-spec file.
I then mask all points below the mixed layer, and integrate vertically.
In the region I working in (the EAC) I expect the advective terms to dominate.
The first implementation is very noisy, but at least dT/dt and the advective term are of the same order of magnitude. The second approach, even though I’ve tried to be very careful, gives me results that are an order of magnitude greater than dT/dt, so clearly I’ve made an error.
Any advice on how you’ve computed the advective terms?