[MITgcm-support] New added Diagnostics-problem

kaveh Purkiani kavehpurkiani at googlemail.com
Sun Oct 25 09:36:08 EDT 2015


Hi all,

I followed the user manual to add new diagnostics to the code.

1- New diagnostic in mypackage_diagnostics_init.F is added like below:

      diagName  = 'new     '         (it is 8 characters long)
      diagTitle = 'New Concentration'
      diagUnits = 'mmol/m^3/d      '
      diagCode  = 'SMR     MR      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I     diagName, diagCode, diagUnits, diagTitle, 0,myThid )


2- At mypackage_forcing.F the diagnostics_fill call is added:

#ifdef ALLOW_DIAGNOSTICS
      CALL DIAGNOSTICS_FILL(new_pTracer,'new     ',0,Nr,1,bi,bj,myThid)
#endif /* ALLOW_DIAGNOSTICS */


3- DIAGNOSTICS_SIZE.h is considered at code folder and data.diagnostics has
been defined in input folder.
Model compiled successfully. I have this error massage at model run :

 (PID.TID 0001.0001) *** ERROR *** DIAGNOSTICS_SET_POINTERS: new      is
not a Diagnostic

When I check the available.diagnostic.log in the running folder there is no
new variable which I thought it has been provided through the implemented
code.

the diagnostics_set_ponters.F says :

If (ndCount.EQ.0) call this error.
ndCount is number of diasable diagnostics in fldsList which is unknown to
me.


Any clue to solve this is very much appreciated.


regards,

Kaveh.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20151025/2c19c33e/attachment.htm>


More information about the MITgcm-support mailing list