[MITgcm-support] Using the "gencost" and "genctrl" approach for adjoint sensitivities

[Dan Jones]

Hi Gael,

I would like to continue this discussion from a couple months ago.  I am
performing a sensitivity experiment using "gencost" in the recently
released ECCOv4 framework that you wrote about, with the following settings
in data.ecco:

    &ECCO_GENCOST_NML

    gencost_avgperiod(1)  = 'month',

    gencost_barfile(1) = 'm_boxmean_theta',

    gencost_errfile(1) = 'my_region_mask',

    gencost_name(1) = 'boxmean_theta',

    gencost_posproc(1,1)='boxmean',

    gencost_outputlevel(1)=1,

    mult_gencost(1) = 1.0,

   #

If, for example, I want to quantify the impact of a given independent
variable on the last month of the simulation, I need to provide a temporal
mask (i.e. 'my_region_maskT' that is "1" for the last month and "0" for all
other months, correct?

Best,
Dan

On Sat, Sep 17, 2016 at 5:00 PM, <mitgcm-support-request at mitgcm.org> wrote:

> Send MITgcm-support mailing list submissions to
>         mitgcm-support at mitgcm.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         http://mitgcm.org/mailman/listinfo/mitgcm-support
> or, via email, send a message with subject or body 'help' to
>         mitgcm-support-request at mitgcm.org
>
> You can reach the person managing the list at
>         mitgcm-support-owner at mitgcm.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of MITgcm-support digest..."
>
>
> Today's Topics:
>
>    1. Re: Using the "gencost" and "genctrl" approach for        adjoint
>       sensitivities (gael forget)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 16 Sep 2016 18:45:10 -0400
> From: gael forget <gforget at mit.edu>
> To: mitgcm-support at mitgcm.org
> Subject: Re: [MITgcm-support] Using the "gencost" and "genctrl"
>         approach for    adjoint sensitivities
> Message-ID: <E6C49C65-5F6D-4E79-8754-76E967CB57D7 at mit.edu>
> Content-Type: text/plain; charset="windows-1252"
>
> Hi Dan,
> I added documentation in chapter 8 of the MITgcm manual.
> It is still rather incomplete and rough around the edges
> but hopefully is already helpful.
> Cheers,
> Gael
>
> ps: kudos to C. Wolfe for providing most of the tables.
>
>
> On Sep 15, 2016, at 8:31 AM, gael forget <gforget at mit.edu> wrote:
>
> > Hi Dan,
> >
> > there is an example for what you are aiming to do with gencost in
> > http://mitgcm.org/viewvc/*checkout*/MITgcm/MITgcm_
> contrib/verification_other/global_oce_cs32/input_ad.sens/data.ecco
> > where the first cost term is for tracer, the second is for transport, and
> > data.profiles is empty (so it is effectively turned off ? but you can
> also do this via
> > packages.conf and data.pkg). The cost function location (space and time)
> is specified
> > at run-time via masks (south30_mask* and north10_mask*) that were
> created using
> > http://mitgcm.org/viewvc/*checkout*/MITgcm/MITgcm_
> contrib/verification_other/global_oce_cs32/input_ad.sens/some_mask.m
> >
> > The genctrl settings for this experiment are in
> > http://mitgcm.org/viewvc/*checkout*/MITgcm/MITgcm_
> contrib/verification_other/global_oce_cs32/input_ad/data.ctrl
> > and provide examples for forcing, initial condition, and internal
> parameter controls.
> >
> > I will add some documentation in the manual shortly about this.
> >
> > Cheers,
> > Gael
> >
> > On Sep 15, 2016, at 6:49 AM, Dan Jones <dcjones.work at gmail.com> wrote:
> >
> >> Hi all,
> >>
> >> I have some questions about using "gencost" and "genctrl".  To be
> concrete, I'll pick a specific example below.  In this example, I am *not*
> generating a state estimate - I am interested in examining the gradients
> themselves (i.e. \partial J / \partial x, where J is the cost function and
> x is some control variable).
> >>
> >> I would like to start with an up-to-date verification exercise
> (MITgcm_contrib/verification_other/global_oce_biogeo_bling).  I want to
> modify this setup such that:
> >>
> >> - The cost function is the value of one of the tracer fields, at a
> location specified at run-time, averaged over the last 'x' seconds
> >> - The cost function should *not* include any contributions from
> "profiles"; how do I turn off those parts entirely?  Does just dropping the
> package from the packages.conf list work?
> >> - The controls include physical variables (e.g. u, v, T, S) *and*
> tracer fields (e.g. DIC)
> >>
> >> How do I need to modify the headers/source code to specify this?  I am
> using TAF to generate the adjoint code.
> >>
> >> I have looked at the comments in the header files, but as this approach
> is new I am feeling a bit uncertain about it.  Any help/clarification is
> very much appreciated!
> >>
> >> Best,
> >> Dan
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >> _______________________________________________
> >> 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
>
> -------------- next part --------------
> An HTML attachment was scrubbed...
> URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/
> 20160916/7918153d/attachment-0001.htm>
> -------------- next part --------------
> A non-text attachment was scrubbed...
> Name: smime.p7s
> Type: application/pkcs7-signature
> Size: 1843 bytes
> Desc: not available
> URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/
> 20160916/7918153d/attachment-0001.p7s>
>
> ------------------------------
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>
>
> End of MITgcm-support Digest, Vol 159, Issue 21
> ***********************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20161104/6878b97e/attachment.htm>