[MITgcm-support] user-defined diagnostic filled with zeros

[Jean-Michel Campin]

Hi Alexandre,

There has been 3 reply to a previous, similar email sent on Aug 05:
http://mitgcm.org/pipermail/mitgcm-support/2015-August/010036.html

Could you try to answer these messages, and specially these 2:
http://mitgcm.org/pipermail/mitgcm-support/2015-August/010037.html
http://mitgcm.org/pipermail/mitgcm-support/2015-August/010049.html

Cheers,
Jean-Michel

On Mon, Aug 17, 2015 at 09:19:09AM +0200, Alexandre Pohl wrote:
> Dear all,
> 
> I am using a coupled configuration of the MITgcm with the pkg/DIC and pkg/diagnostics enabled. I would like to add a new diagnostic to the permanent set of diagnostics in the model: 'omegaC' (computed in 'dic/calcite_saturation.F').
> 
> I followed the instructions from the guidebook (http://mitgcm.org/public/r2_manual/latest/online_documents/node269.html <http://mitgcm.org/public/r2_manual/latest/online_documents/node269.html>):
> 
> 1] First, I define the new diagnostic in 'dic_diagnostics_init.F' :
> 
>      diagName  = 'OMEGAC  '
>      diagTitle = 'Calcite Omega   '
>      diagUnits = 'dimensionless   '
>      diagCode  = 'SMRP    MR      '
>      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
>     I       diagName, diagCode, diagUnits, diagTitle, 0, myThid )
> 
> 2] Then, I call the 'DIAGNOSTICS_FILL' routine in 'dic_biotic_forcing.F'.
> 
>        CALL DIAGNOSTICS_FILL(SURC   ,'DICTFLX ',0,1 ,2,bi,bj,myThid)    <==== one of pkg/DIC default diagnostics
>        CALL DIAGNOSTICS_FILL(omegaC ,'OMEGAC  ',0,Nr,1,bi,bj,myThid)    <==== my new diag
> #ifdef ALLOW_O2
>        CALL DIAGNOSTICS_FILL(SURO   ,'DICOFLX ',0,1 ,2,bi,bj,myThid)    <==== another pkg/DIC default diagnostic
> 
> Compilation succeeds, the model runs and diagnostics are successfully written in the model output mdsio files. The problem is that the newly added diagnostic, OMEGAC, is filled with zeros. I guess that it has something to do with either 'diagCode' from step 1 (which I set to  'SMRP    MR      ' like for DICBIOA, which seemed to be correct after Table 7.1: http://mitgcm.org/public/r2_manual/latest/online_documents/node268.html <http://mitgcm.org/public/r2_manual/latest/online_documents/node268.html>), or with the bibjFlg from step 2 (which I set to 1 because omegaC is not modified in 'dic_biotic_forcing.F', following instructions given at the top of the same page).
> 
> Could anyone tell me what I am doing wrong? 
> 
> Thanks!
> Alexandre
> 
> __
> 
> Alexandre POHL
> PhD Student
> Laboratoire des Sciences du Climat et de l???Environnement (LSCE)
> Centre d'Etude de Saclay
> CEA-Orme des Merisiers, bât. 701
> F-91191 GIF-SUR-YVETTE CEDEX
> FRANCE
> alexandre.pohl at lsce.ipsl.fr

> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support