[MITgcm-support] GM-Redi question

[Julian Mak]

I think it's "GM_isopycK" which is the Redi part? By default if that is not specified I think it takes on the value of GM (although I don't know if this would be "GM_background_K" or whatever the kGM value that the Visbeck scheme spits out in this case).


A data.gmredi I used is given below; ignore the MMB stuff which is our GEOMETRIC implementation where we hijacked the existing Visbeck scheme in a rather disgusting fashion (I haven't tried very hard tidying it up...)


Best wishes,
Julian


%%%%%%%%%%%%%%%%%%%%


&GM_PARM01
  GM_AdvForm         = .FALSE.
  GM_background_K    = 0.,
  GM_isopycK         = 200.,
  GM_taper_scheme    = 'gkw91',
  GM_maxSlope        = 5.e-3,
  GM_Kmin_horiz      = 0.,
  GM_Scrit           = 4.e-3,
  GM_Sd              = 1.e-3,
  GM_Visbeck_alpha   = 0.016,
  GM_Visbeck_minval_K= 50.,
  GM_Visbeck_maxval_K= 1.e+4,
# MMB params
  MMB_lambda         = 1.e-7,
  MMB_alpha          = 0.04,
  MMB_minval_K       = 50.,
  MMB_maxval_K       = 1.e+4,
#  MMB_minval_K       = 1155,
#  MMB_maxval_K       = 1155,
  energy_init        = 1.e-2,
  energy_kappa       = 0,
  energy_local       = .false.,


%%%%%%%%%%%%%%%%%%%%


On 20/02/2019 10:34, Naughten, Kaitlin A. wrote:

Hello everyone,


I am a bit confused about the GM-Redi scheme, and how to turn on the "GM" part (bolus or skew fluxes) versus just the "Redi" part (tracer diffusion along isentropes). An old message from Martin (http://mailman.mitgcm.org/pipermail/mitgcm-support/2016-March/010355.html) said that he only used the Redi part, but I'm not sure what options control this.


I switched on all the GM-Redi diagnostics in my configuration and discovered that many variables are zero everywhere (GM_PsiX, GM_PsiY, GM_ubT, GM_vbT, GM_hTrsl, GM_baseS, GM_rLamb, plus two components of the Redi tensor which I think are supposed to be zero: GM_Kuz and GM_Kvz). Does this mean that I only have the Redi part switched on? And if so, why is the GM part not on?


Here are the active lines of my data.gmredi:


 &GM_PARM01
 GM_background_K=10.,
 GM_taper_scheme='ldd97',
 GM_Kmin_horiz=50.,
 GM_Visbeck_alpha=0.07,
 GM_Visbeck_length=50.E3,
 GM_Visbeck_depth=2000.,
 GM_Visbeck_maxVal_K=300.,


And here are the active lines of my GM_REDI_OPTIONS.h:


#ifndef GMREDI_OPTIONS_H

#define GMREDI_OPTIONS_H
#include "PACKAGES_CONFIG.h"
#include "CPP_OPTIONS.h"
#ifdef ALLOW_GMREDI
#define GM_VISBECK_VARIABLE_K
#undef OLD_VISBECK_CALC
#undef GM_K3D
#undef GM_K3D_PASSIVE
#define GM_NON_UNITY_DIAGONAL
#define GM_EXTRA_DIAGONAL
#define GM_BOLUS_ADVEC
#undef GM_BOLUS_BVP
#endif /* ALLOW_GMREDI */
#endif /* GMREDI_OPTIONS_H */

This is probably a stupid question but hopefully that means it has a quick and easy answer!

Many thanks,
Kaitlin


Dr Kaitlin Naughten
British Antarctic Survey


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 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. UK Research and Innovation does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses.
Opinions, conclusions or other information in this message and attachments that are not related directly to UK Research and Innovation business are solely those of the author and do not represent the views of UK Research and Innovation.




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


--
Julian Mak,
Department of Physics,
University of Oxford,
Clarendon Laboratory, Parks Road
Oxford, OX1 3PU
e-mail: julian.mak at physics.ox.ac.uk<mailto:julian.mak at physics.ox.ac.uk>
www: https://sites.google.com/site/julianclmak/home
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20190220/585a3a3b/attachment-0001.html>