[MITgcm-support] GMREDI without the small slope approximation
Martin Losch
Martin.Losch at awi.de
Tue Apr 25 04:44:40 EDT 2023
Hi Flynn,
I am not really qualified to answer this, so take my comments with a grain of salt.
As far as I can see, the missing tensor components (basically k12=Kuy and k21=Kvx) need to be defined and computed and you are right that these are small and straightforward changes to the code (gmredi_calc_tensor.F). But you will also need to code the action of these tensor components on the \nabla(u) in the right place (k12*dv/dy in xtranport and k21*du/dy in ytransport?, involves some averaging of the derivatives). I find this hard to figure out, because there are so many different options in pkg/gmredi, but it should be possible.
Martin
> On 21. Apr 2023, at 18:38, Flynn Ames <f.ames at pgr.reading.ac.uk> wrote:
>
> Dear MITgcm community,
>
> Hello - I'm wondering if anyone has ever tried implementing (or has implemented) GMREDI into the MITgcm without making the small slope approximation?
>
> This is something I would like to do as I'm modelling a deep, ice-covered ocean where the steep isopycnals that develop violate the small slope approximation over much of the domain.
>
> I may be being naive here, but looking at the MITgcm documentation, the terms that are neglected in the REDI projection tensor when making the small slope approximation appear to consist of quantities that are already computed in the model (Sx, Sy etc.).
> I'm hoping to ask: Does anyone know if it would just be a case of putting these extra terms into the REDI tensor in the model? or is there a lot that I'm missing here? - I'm guessing the latter :)
>
> Thank you for your time!
> Flynn
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org <mailto:MITgcm-support at mitgcm.org>
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230425/10fd5c4b/attachment.html>
More information about the MITgcm-support
mailing list