--- MITgcm_contrib/dfer/matlab_stuff/calcHeatTransDirect.m 2018/03/07 22:01:18 1.2 +++ MITgcm_contrib/dfer/matlab_stuff/calcHeatTransDirect.m 2018/04/19 00:03:59 1.3 @@ -39,7 +39,11 @@ % nFld = 5 = Residual [5=2-4] % Ocn only: % If GM advective form: -% nFld = 6 = HT by horizontal diffusion +% nFld = 6 = total (GM+Redi) eddy transport: K*dT/dy +% WARNING: This assumes K_gm=K_iso to get +% the cancellation of the off-diagonal term +% (K_iso-K_gm)*S_y*dT/dz +% nFld = 7 = GM (advective) eddy transport % If Skew flux form: % nFld = 6 = total (GM+Redi) eddy transport % nFld = 7 = GM (advective) eddy transport