C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/dgoldberg/streamice/streamice_init_diagnostics.F,v 1.2 2012/04/20 14:09:13 dgoldberg Exp $
C $Name:  $

#include "STREAMICE_OPTIONS.h"

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
CBOP 0
C !ROUTINE: STREAMICE_INIT_FIXED

C !INTERFACE:
      SUBROUTINE STREAMICE_INIT_DIAGNOSTICS( myThid )

C     !DESCRIPTION:
C     Initialize STREAMICE variables that are kept fixed during the run.

C     !USES:
      IMPLICIT NONE
#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "STREAMICE.h"
#include "STREAMICE_CG.h"
#include "STREAMICE_BDRY.h"
#include "GRID.h"

      INTEGER myThid

#ifdef ALLOW_STREAMICE
#ifdef ALLOW_DIAGNOSTICS

C     myThid ::  my Thread Id number

      INTEGER       diagNum
      INTEGER       diagMate
      CHARACTER*8   diagName
      CHARACTER*16  diagCode
      CHARACTER*16  diagUnits
      CHARACTER*(80) diagTitle

      diagName  = 'SI_Uvel '
      diagTitle = 'Ice stream x-velocity'
      diagUnits = 'm/a        '
      diagCode  = 'UZ      L1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'SI_Vvel '
      diagTitle = 'Ice stream y-velocity'
      diagUnits = 'm/a        '
      diagCode  = 'VZ      L1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'SI_Thick'
      diagTitle = 'Ice stream thickness'
      diagUnits = 'm          '
      diagCode  = 'SM      L1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'SI_area '
      diagTitle = 'Ice stream cell area coverage'
      diagUnits = 'm^2        '
      diagCode  = 'SM      L1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'SI_float'
      diagTitle = 'Ice stream grounding ind'
      diagUnits = 'none       '
      diagCode  = 'SM      L1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )

      diagName  = 'SI_hmask'
      diagTitle = 'Ice stream thickness mask'
      diagUnits = 'none       '
      diagCode  = 'SM      L1      '
      CALL DIAGNOSTICS_ADDTOLIST( diagNum,
     I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )


      PRINT *, "ADDED DIAGS TO LIST"


#endif
#endif
      RETURN
      END