[MITgcm-support] Adding Biogeochemistry (BGC) Cost/Objective Function

[Gael Forget]

Hi Yohei,

I always tend to to think that it’s better to use the adjoint controls as a way to diagnose adjoint sensitivities (adxx* derived by TAF) instead of the hand-coded ADJ files (which we never know to be correct). Just my 2 cents. Hope this helps, or you have already resolved the issue.

Cheers,
Gaël


> On Nov 20, 2024, at 6:28 AM, Yohei Takano - BAS <yokano at bas.ac.uk> wrote:
> 
> Dear all,
> 
>    I am trying to add new cost/objective function within pkg/ecco framework following the
> ptracer example from the codes. For the bgc objective function, I am thinking of introducing
> air-sea co2 flux, pco2 (DIC and BLING package), chlorophyll (BLING package)  for example and
> now working on air-sea co2 flux in DIC package as a starter.
> 
>    In the end, I have modified "cost_gencost_customize.F", "ecco_phys.F", and "ecco_check.F"
> codes to include  "m_DICCFLX" and "m_boxmean_DICCFLX" for gencost bar file and manage to
> compile and run the model. However, the adjoint part of the simulation does not generate any
> ADJ* file (adxx_* files are generated, I still need to check the output) and I am wondering if I need 
> to change other part of the code to include this such as codes in pkg/autodiff. Here I am thinking of 
> the default ADJ* files (such as taux, temperature etc. common for other adjoint simulations) 
> not something new. The costfunction0012 contains values so I think it is defining the objective function okay, 
> just the output behaviour is different from what I expected from other runs.
> 
>    In the end I understand I should separate bgc from physics (i.e. not to include in "ecco_phys.F") but
> just for a starter I am starting from modifying the existing codes. I am running the model for 2 years
> the adjDumpFreq    = 1209600.0, so I am expecting to have ADJ* output every two weeks.
> 
>    Is there any other part of the codes/pkg I need to update (or missing some settings)?
> Any comments, suggestions will be great.
> 
> Thank you,
> 
> Yohei
>    
> 
> 
> 
> This email and any attachments are intended solely for the use of the named recipients. If you are not the intended recipient you must not use, disclose, copy or distribute this email or any of its attachments and should notify the sender immediately and delete this email from your system. UK Research and Innovation (UKRI) has taken every reasonable precaution to minimise risk of this email or any attachments containing viruses or malware but the recipient should carry out its own virus and malware checks before opening the attachments. UKRI does not accept any liability for any losses or damages which the recipient may sustain due to presence of any viruses.
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org <mailto:MITgcm-support at mitgcm.org>
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20241202/5ec6f1be/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/pkcs7-signature
Size: 1859 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20241202/5ec6f1be/attachment.p7s>