pkg/streamice/streamice_init_diagnostics.F

5ca83cd8f7 Dani*0001 #include "STREAMICE_OPTIONS.h"
                0002
                0003 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0004 CBOP 0
                0005 C !ROUTINE: STREAMICE_INIT_FIXED
                0006
                0007 C !INTERFACE:
                0008       SUBROUTINE STREAMICE_INIT_DIAGNOSTICS( myThid )
                0009
                0010 C     !DESCRIPTION:
                0011 C     Initialize STREAMICE variables that are kept fixed during the run.
                0012
                0013 C     !USES:
                0014       IMPLICIT NONE
                0015 #include "EEPARAMS.h"
                0016 #include "SIZE.h"
                0017 #include "PARAMS.h"
                0018 #include "STREAMICE.h"
                0019 #include "STREAMICE_CG.h"
                0020 #include "STREAMICE_BDRY.h"
                0021 #include "GRID.h"
                0022
07e785229e dngo*0023 C     myThid ::  my Thread Id number
5ca83cd8f7 Dani*0024       INTEGER myThid
                0025
                0026 #ifdef ALLOW_STREAMICE
                0027 #ifdef ALLOW_DIAGNOSTICS
                0028
07e785229e dngo*0029       CHARACTER*(MAX_LEN_MBUF) msgBuf
5ca83cd8f7 Dani*0030       INTEGER       diagNum
07e785229e dngo*0031 c     INTEGER       diagMate
5ca83cd8f7 Dani*0032       CHARACTER*8   diagName
                0033       CHARACTER*16  diagCode
                0034       CHARACTER*16  diagUnits
                0035       CHARACTER*(80) diagTitle
                0036
                0037       diagName  = 'SI_Uvel '
                0038       diagTitle = 'Ice stream x-velocity'
                0039       diagUnits = 'm/a        '
                0040       diagCode  = 'UZ      L1      '
                0041       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0042      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0043
                0044       diagName  = 'SI_Vvel '
                0045       diagTitle = 'Ice stream y-velocity'
                0046       diagUnits = 'm/a        '
                0047       diagCode  = 'VZ      L1      '
                0048       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0049      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0050
                0051       diagName  = 'SI_Thick'
                0052       diagTitle = 'Ice stream thickness'
                0053       diagUnits = 'm          '
                0054       diagCode  = 'SM      L1      '
                0055       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0056      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0057
                0058       diagName  = 'SI_area '
                0059       diagTitle = 'Ice stream cell area coverage'
                0060       diagUnits = 'm^2        '
                0061       diagCode  = 'SM      L1      '
                0062       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0063      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0064
                0065       diagName  = 'SI_float'
                0066       diagTitle = 'Ice stream grounding ind'
                0067       diagUnits = 'none       '
                0068       diagCode  = 'SM      L1      '
                0069       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0070      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0071
                0072       diagName  = 'SI_hmask'
                0073       diagTitle = 'Ice stream thickness mask'
                0074       diagUnits = 'none       '
                0075       diagCode  = 'SM      L1      '
                0076       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0077      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0078
                0079       diagName  = 'SI_usurf'
                0080       diagTitle = 'Ice stream surface x-vel'
                0081       diagUnits = 'none       '
                0082       diagCode  = 'SM      L1      '
                0083       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0084      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0085
                0086       diagName  = 'SI_vsurf'
                0087       diagTitle = 'Ice stream surface y-vel'
                0088       diagUnits = 'none       '
                0089       diagCode  = 'SM      L1      '
                0090       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0091      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0092
7924606b84 Dani*0093       diagName  = 'SI_ubase'
                0094       diagTitle = 'Ice stream basal x-vel'
                0095       diagUnits = 'none       '
                0096       diagCode  = 'SM      L1      '
                0097       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0098      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0099
                0100       diagName  = 'SI_vbase'
                0101       diagTitle = 'Ice stream basal y-vel'
                0102       diagUnits = 'none       '
                0103       diagCode  = 'SM      L1      '
                0104       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0105      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0106
                0107       diagName  = 'SI_taubx'
                0108       diagTitle = 'Ice stream basal x-stress'
                0109       diagUnits = 'none       '
                0110       diagCode  = 'SM      L1      '
                0111       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0112      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0113
                0114       diagName  = 'SI_tauby'
                0115       diagTitle = 'Ice stream basal y-stress'
                0116       diagUnits = 'none       '
                0117       diagCode  = 'SM      L1      '
                0118       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0119      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0120
5ca83cd8f7 Dani*0121       diagName  = 'SI_selev'
                0122       diagTitle = 'Ice stream surface elev'
                0123       diagUnits = 'none       '
                0124       diagCode  = 'SM      L1      '
                0125       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0126      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0127
                0128       diagName  = 'SI2Dtrac'
                0129       diagTitle = 'Ice stream 2d tracer'
                0130       diagUnits = 'none       '
                0131       diagCode  = 'SM      L1      '
                0132       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0133      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0134
07e785229e dngo*0135       diagName  = 'BDOT_SI '
                0136       diagTitle = 'Ice stream 2d tracer'
                0137       diagUnits = 'none       '
                0138       diagCode  = 'SM      L1      '
                0139       CALL DIAGNOSTICS_ADDTOLIST( diagNum,
                0140      I      diagName, diagCode, diagUnits, diagTitle, 0, myThid )
                0141
5ca83cd8f7 Dani*0142       WRITE(msgBuf,'(A)') 'ADDED DIAGS TO LIST'
                0143        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
                0144      &                     SQUEEZE_RIGHT , 1)
                0145
07e785229e dngo*0146 #endif /* ALLOW_DIAGNOSTICS */
                0147 #endif /* ALLOW_STREAMICE */
5ca83cd8f7 Dani*0148       RETURN
                0149       END