[MITgcm-support] Understanding gT_Forc and gS_Forc diagnostics

[Dafydd Stephenson]

Hi all,
I’m hoping to get the external forcing impact on T,S,U, and V in terms of a temperature (K/s), salinity (g/kg/s) and velocity (m/s2) tendency. For the velocity fields I think this is fairly straightforward: Um_Ext and Vm_Ext seem to provide exactly what I’m looking for, and a crude rescaling of EXFtaux and EXFtauy [dividing by the product of a constant reference density (rho0) and upper layer thickness (10m)] compares well to Um_Ext and Vm_Ext.

I’m a bit stumped by the equivalent output for T and S, though. It seems from the manual that I want gT_Forc and gS_Forc, but these fields do not compare as well to -EXFqnet/(rho0*cp*10m) and (EXFempmr*S0)/(rho0*10m) at the surface, particularly in the polar oceans (which is perhaps expected). However, my gT_Forc and gS_Forc output is three-dimensional, which is not what I expected – the source code only seems to refer to, e.g., gtForc(i,j).

Are these fields what I’m looking for? Where does the extra dimension come from in the final output?

Thanks for your help!
Dafydd
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20210923/aefa7683/attachment.html>