C $Header: /u/gcmpack/MITgcm/pkg/dic/dic_diagnostics_init.F,v 1.1 2010/05/14 21:52:33 jahn Exp $ C $Name: $ #include "GCHEM_OPTIONS.h" C !INTERFACE: ========================================================== SUBROUTINE DIC_DIAGNOSTICS_INIT(myThid ) C !DESCRIPTION: C define diagnostics for dic package C experiment C !USES: =============================================================== IMPLICIT NONE #include "SIZE.h" #include "EEPARAMS.h" C !INPUT PARAMETERS: =================================================== C myThid :: thread number INTEGER myThid CEOP #ifdef ALLOW_DIAGNOSTICS C !LOCAL VARIABLES: INTEGER diagNum CHARACTER*8 diagName CHARACTER*16 diagCode CHARACTER*16 diagUnits CHARACTER*(80) diagTitle C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| C Define diagnostics Names : diagName = 'DICBIOA ' diagTitle = 'Biological Productivity (mol/m3/s)' diagUnits = 'mol/m3/sec ' diagCode = 'SMRP MR ' CALL DIAGNOSTICS_ADDTOLIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) diagName = 'DICCARB ' diagTitle = 'Carbonate chg-biol prod and remin (mol eq/m3/s)' diagUnits = 'mol eq/m3/sec ' diagCode = 'SMRP MR ' CALL DIAGNOSTICS_ADDTOLIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) diagName = 'DICTFLX ' diagTitle = 'Tendency of DIC due to air-sea exch (mol/m3/s)' diagUnits = 'mol/m3/sec ' diagCode = 'SM P L1 ' CALL DIAGNOSTICS_ADDTOLIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) diagName = 'DICOFLX ' diagTitle = 'Tendency of O2 due to air-sea exch (mol/m3/s)' diagUnits = 'mol/m3/sec ' diagCode = 'SM P L1 ' CALL DIAGNOSTICS_ADDTOLIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) diagName = 'DICCFLX ' diagTitle = 'Flux of CO2 - air-sea exch (mol/m2/s)' diagUnits = 'mol/m2/sec ' diagCode = 'SM P L1 ' CALL DIAGNOSTICS_ADDTOLIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) diagName = 'DICPCO2 ' diagTitle = 'Partial Pressure of CO2 (atm)' diagUnits = 'atm ' diagCode = 'SM P M1 ' CALL DIAGNOSTICS_ADDTOLIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) diagName = 'DICPHAV ' diagTitle = 'pH (dimensionless)' diagUnits = 'dimensionless ' diagCode = 'SM P M1 ' CALL DIAGNOSTICS_ADDTOLIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, 0, myThid ) C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| #endif /* ALLOW_DIAGNOSTICS */ RETURN END