[MITgcm-support] Controlling the adjoint behaviour
Matthew Mazloff
mmazloff at ucsd.edu
Wed May 18 13:27:45 EDT 2011
Hi Chris,
You can use
adjDumpFreq
in PARM03 of your data file to write snapshots of the adjoint
sensitivities.
e.g.
# Time stepping parameters
&PARM03
nIter0 = 0,
nTimeSteps = 72,
deltaTmom = 1200.0,
deltaTtracer= 1200.0,
deltaTClock = 1200.0,
adjDumpFreq= 3600.0,
&
will dump sensitivities every hour
-Matt
On May 18, 2011, at 9:37 AM, Wilson, Christopher wrote:
> Hello all,
>
> I have set up a shallow water model and have hacked
> cost_atlantic_heat.F to set the cost function to be the maximum
> value of etaN in the domain.
>
> I would like to calculate the sensitivity of this cost function to
> bottom drag (BOTTOM_DRAG_CONTROL) at frequent intervals, but at
> present the adxx_* output contains only one time slice.
> I presume this is due to COST_ATLANTIC_HEAT only being called once,
> within cost_final.F.
>
> I suspect I need to insert the call to a subroutine evaluating the
> cost within cost_tile.F, but can't get this to work as expected:
> adxx_* still only produces one time slice.
> I'm generally mystified, so would be very grateful of any pointers.
>
> Thanks,
>
> Chris Wilson
>
>
>
>
>
> --
> This message (and any attachments) is for the recipient only. NERC
> is subject to the Freedom of Information Act 2000 and the contents
> of this email and any reply you make may be disclosed by NERC unless
> it is exempt from release under the Act. Any material supplied to
> NERC may be stored in an electronic records management system.
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list