# Ocean Modeling Discussion

ROMS/TOMS

Search for:
 It is currently Mon Dec 11, 2017 11:19 am

 Page 1 of 1 [ 1 post ]
 All times are UTC
Author Message
 Post subject: MPDATA AdjointPosted: Wed Oct 09, 2013 7:00 pm

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:
Quote:
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?

Thanks,
Austin

Top

 Display posts from previous: All posts1 day7 days2 weeks1 month3 months6 months1 year Sort by AuthorPost timeSubject AscendingDescending
 Page 1 of 1 [ 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 forumYou cannot reply to topics in this forumYou cannot edit your posts in this forumYou cannot delete your posts in this forumYou cannot post attachments in this forum

Search for:
 Jump to:  Select a forum ------------------ News, Events & Job Opportunities    Meetings/Workshops    Job Opportunities    Ocean News ROMS/TOMS    ROMS Adjoint    ROMS Benchmarks    ROMS Bugs    ROMS Discussion    ROMS Documentation    ROMS Ecosystem    ROMS FAQ    ROMS Ice    ROMS Information    ROMS Installation    ROMS Known Problems    ROMS Messages    ROMS Problems    ROMS Releases    ROMS Results    ROMS Sediment    ROMS Source    ROMS Tips    ROMS Tools and Techniques    ROMS Trivia    ROMS Usage    ROMS Webinar    ROMS Wish List ROMS/TOMS Applications    User Applications    Adriatic Sea Ocean Modeling    Ocean Modeling FAQ