Ocean Modeling Discussion


Search for:
It is currently Thu Oct 19, 2017 11:06 am

Post new topic Reply to topic  [ 1 post ] 

All times are UTC

Author Message
 Post subject: MPDATA Adjoint
PostPosted: Wed Oct 09, 2013 7:00 pm 
User avatar

Joined: Wed Mar 11, 2009 4:24 pm
Posts: 25
Location: North Carolina State University
I've been doing some modeling of passive tracers, and I want to use the Adjoint to look at the backward advection of the passive tracers. Since my main interest lies in evaluating the tracer concentrations, I am using MPDATA for the forward integration to reduce overshoots. With other tracer advection schemes, I get values O(10^7-10^11) in both the positive and negative side for a tracer that should range from 0 to 1. MPDATA works well, with only small overshoots.

My problem then lies with the adjoint calculations, for which MPDATA routines are not yet coded. As I think about this more, I wonder if it is even possible to code an adjoint algorithm for a predictor-corrector tracer advection scheme like MPDATA? For example, how would we know the error of the first-pass a priori, if the adjoint requires calculating the anti-diffusive velocities first?

I noticed a track ticket posted a LONG time ago that referred to the difficulty of this calculation (https://www.myroms.org/projects/src/ticket/74), but haven't found any further discussion since. In this ticket, it is said:
The tangent linear and adjoint linearization of the MPDATA algorithm is not trivial. MPDATA, as implemented, is highly nonlinear and we don't know if the tangent linear approximation will hold. The adjoint is a sensitivity operator and it is likely that the TLM, RPM, and ADM will go unstable

Have there been any advances in ways to try to code the adjoint for MPDATA-type advection? It seems that the MPDATA code has been pretty extensively used since this 6-yr old post. If not, are there any thoughts on what the next best method would be to try to reduce the overshoots with the adjoint calculation?


Reply with quote  
Display posts from previous:  Sort by  
Post new topic Reply to topic  [ 1 post ] 

All times are UTC

Who is online

Users browsing this forum: No registered users and 1 guest

You cannot post new topics in this forum
You cannot reply to topics in this forum
You cannot edit your posts in this forum
You cannot delete your posts in this forum
You cannot post attachments in this forum

Search for:
Jump to:  
Powered by phpBB® Forum Software © phpBB Group