[MITgcm-support] Lateral heat flux due to temperature restoring

[Jean-Michel Campin]

Hi Josephine,

I would like to know few other details about your set-up and configuration:
 1) are you using z* (i.e., select_rStar=1 or 2), and are you using implicitDiffusion=T ?
 2) and regarding this: 
> Temperature and salinity are restored at the open lateral boundaries.
   are you using pkg/rbcs or is it applied as part of the pkg/obcs "sponge" option ?
  In the first case (with pkg/rbcs), the tendency from the relaxation would be included
  inside the general (3-D) "forcing" tendency" diagnostics "gT_Forc ' and "gS_Forc ",
  but not if using pkg/obcs "sponge" option..

Cheers,
Jean-Michel

On Tue, Sep 17, 2024 at 04:14:38PM +0000, Josephine Anselin - BAS wrote:
> Hi everyone,
> 
> I am evaluating the heat budget in a 3D domain with nx=640 grid points in the x-direction; ny=1200 grid points in the y-direction and nz=150 grid points in the z direction. The domain includes open ocean, ice shelf and grounded ice cells. Temperature and salinity are restored at the open lateral boundaries.
> 
> The following parameter settings are used:
> implicitFreeSurface=.TRUE.
> nonlinFreeSurf=4.
> useRealFreshWaterFlux=.FALSE.
> 
> To begin with, I am trying to close the budget for the cells away from the ice shelf-ocean interface and away from the ocean surface. I was assuming that for these cells, the heat budget should be:
> TOTTEND/86400 = Adv_tend + Dif_tend,
> with Adv_tend and Dif_tend calculated as per page 3 in this document: https://mitgcm.org/download/daily_snapshot/MITgcm/doc/Heat_Salt_Budget_MITgcm.pdf
> (KPPg_TH and oceQsw are zero in the model configuration).
> 
> However, the heat budget does not close using this calculation.
> I would like to investigate the potential impact of the lateral heat flux due to temperature restoring to check if this might be causing the large residual.
> Is there a diagnostics analogous to TRELAX but for lateral heat flux instead of surface heat flux?
> If not, does anyone have any suggestions on how this could be calculated using built-in diagnostics?
> 
> Thank you
> Josephine
> 
> 
> This email and any attachments are intended solely for the use of the named recipients. If you are not the intended recipient you must not use, disclose, copy or distribute this email or any of its attachments and should notify the sender immediately and delete this email from your system. UK Research and Innovation (UKRI) has taken every reasonable precaution to minimise risk of this email or any attachments containing viruses or malware but the recipient should carry out its own virus and malware checks before opening the attachments. UKRI does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses.
> 

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support