[MITgcm-support] adjoint run

[Patrick Heimbach]

Hi Mirek,

yes, there's a template for this in place that you should be able to adapt without too much pain:
At each timestep, the model calls a generic subroutine cost_tile
(see pkg/cost/ )
This routine 
1. calls individual cost routines, depending on your choice (CPP option),
2. evaluates this cost function term as long as
      IF ( myTime .GT. (EndTime - lastinterval) ) THEN

I.e. the cost function is evaluated over the last "lastinterval" seconds.
Default is 2592000 seconds (= 30 days).
You can easily change this in data.cost by setting "lastinterval" 
to something different (e.g., 31104000 = 3600*24*360 = one 360-day year).

Also, in cost_tile.F add a call to whichever CO2-related cost routine you have,
and make sure to sum them up in cost_final.F (also in pkg/cost/ ).
A good generic example is the cost function used for the shelfice package
(grep for ALLOW_COST_SHELFICE).

Cheers
-Patrick

On Mar 8, 2013, at 8:03 AM, Miroslaw Andrejczuk <Andrejczuk at atm.ox.ac.uk> wrote:

> 
> Dear Support,
> 
> I was wondering if you could help me to set up an adjoint simulation. I need to calculate the sensitivity for the integral over the whole globe of pCO2. However, I want to keep this function (that is the value of the integral) constant for the duration of the integration (~10 years).
> 
> I would appreciate if you could either help mi to set up such simulation or provide me with the contact to person who could do that.
> 
> Best Regards,
> Mirek Andrejczuk
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support


---
Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS 54-1420 | 77 Massachusetts Ave | Cambridge MA 02139 USA
FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach

-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1588 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130310/3f36f79f/attachment.p7s>