[MITgcm-support] Fwd: New added Diagnostics-problem

Jean-Michel Campin jmc at ocean.mit.edu
Wed Oct 28 15:46:29 EDT 2015


Hi Kaveh,

The error you are getting is consistent with the fact that 
the new diagnostics "new     " is not present in the log file:
 available_diagnostics.log
and I suspect that this piece of code you added in mypackage_diagnostics_init.F
is not run through.

I don't have better suggestion than to check basic things, such as:
1) is pkg/mypackage compiled ?
2) is this package turned on ? ( useMYPACKAGE=.TRUE. in data.pkg)
3) check that S/R MYPACKAGE_DIAGNOSTICS_INIT is called 
  by adding "debugMode=.TRUE." in eedata ; this print the sequence of
 subroutine calls in STDOUT files.
4) check that your changes are indeed present in the "small .f" file 
 version, i.e.: mypackage_diagnostics_init.f (in the build dir).

Cheers,
Jean-Michel

On Wed, Oct 28, 2015 at 04:12:36PM +0100, kaveh Purkiani wrote:
> 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.

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




More information about the MITgcm-support mailing list