[MITgcm-support] Controlling the adjoint behaviour

[Wilson, Christopher]

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.