[MITgcm-support] Adjoins advection schemes

Patrick Heimbach heimbach at MIT.EDU
Mon Jul 2 09:41:30 EDT 2012


our scheme of choice for the adjoint has become
3rd-order direct space and time DST3 in non-flux limited form, i.e. scheme 30.
DST3FL (33) causes problems for certain applications, and we haven't followed up on this
(the obvious issue being how the limiters are formulated).

SOM is very memory-consuming, so we've staid away from this one for the adjoint
(which amplifies memory-consumption).

OS7MP is also currently disabled but might be worth a try
(would need some tweak though before it can be attempted).


On Jul 2, 2012, at 3:58 AM, Martin Losch wrote:

> Christopher, 
> the easiest way to find out is to have a look at gad_advection or gad_calc_rhs (the calls to GAD_*_ADV_X/Y). You'll see that most schemes "work" (i.e., they are allowed to be used) but in particular OS7MP (and I think SOM) doesn't work. From my experience I recommend using linear advection schemes (which rules out all flux limited schemes) for more accurate gradients.  
> Martin
> On Jun 29, 2012, at 10:25 PM, Christopher L. P. Wolfe wrote:
>> Hi all,
>> Does anyone know if there exists a list of advection schemes which are known to "work" with adjoint? I'm particular interested in the higher order schemes like OS7MP and SOM and their flux limited forms. 
>> Thanks in advance,
>> Christopher
>> _______________________________________________
>> MITgcm-support mailing list
>> MITgcm-support at mitgcm.org
>> http://mitgcm.org/mailman/listinfo/mitgcm-support
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support

Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS 54-1420 | 77 Massachusetts Ave | Cambridge MA 02139 USA
FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach

-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1588 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120702/7ea37f18/attachment.p7s>

More information about the MITgcm-support mailing list