[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