[MITgcm-support] Heat Budget in tutorial (Global Ocean Simulation at 4 deg)

[Jean-Michel Campin]

Hi,

Just a detail: with Linear FreeSurf, there is a small loss of heat/salt/tracer
 at the surface (d/dt(Eta) inconsistent with fixed surf grid cell thickness):
a) if linFSConserveTr is set to true (file "data"), then the model 
  correct for this (see S/R calc_wsurf_tr.F).
b) if linFSConserveTr=False (=the default), need to account for this term
 in the budget, which can be obtain from the 1rst level of WTHMASS / WSLTMASS 
 diagnostics (diagnostics ADVr_TH / ADVr_SLT is not the right one for 
 this purpose, it's just zero at k=1).

Cheers,
Jean-Michel

On Fri, May 20, 2011 at 06:56:05AM -0700, Menemenlis, Dimitris (3248) wrote:
> Martin, you are probably aware of this but for support list benefit, JM wrote a very nice summary of surface diagnostic subtleties, which is available here: MITgcm/doc/diags_changes.txt
> 
> TFLUX is the correct diagnostic for temperature budget closure:
> 
> a) if useRealFreshWaterFlux=F or (nonlinFreeSurf=0 & usingZCoords=T)
>  TFLUX = oceQnet  + TRELAX + oceFreez
> 
> b) if useRealFreshWaterFlux=T & (nonlinFreeSurf>0 or usingPCoords=T),
>  TFLUX = oceQnet  + TRELAX + PmEpR*temp_EvPrRn*Cp + oceFreez
> 
> Dimitris Menemenlis
> 
> On May 20, 2011, at 5:39 AM, Martin Losch wrote:
> 
> > In order to get these diagnostics, you need to save (in data.diagnostics)
> > something like oceQnet (surface heat flux, but without relaxing contribution, as far as I remember), or TFLUX (which contains everything)
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support