pkg/steep_icecavity/STIC.h

00f81e6785 Ou W*0001 CBOP
                0002 C     !ROUTINE: STIC.h
                0003 C     !INTERFACE:
                0004 C     #include "STIC.h"
                0005
                0006 C !DESCRIPTION:
                0007 C     *================================================================*
                0008 C     | STIC.h
                0009 C     | o Basic header thermodnynamic steep ice cavity package.
                0010 C     |   Contains all STEEP_ICECAVITY field declarations.
                0011 C     | o this package can only be used as an add-on to pkg/shelfice
                0012 C     |   because it uses most of the shelfice-infrastructure
                0013 C     *================================================================*
                0014
                0015 C-----------------------------------------------------------------------
                0016 C--   Constants that can be set in data.shelfice
                0017 C     icfThetaHorizDiffusionLength ::
                0018 C     icfThetaInterior  ::
                0019 C     STICdepthFile     :: name of icefront depth file (m)
                0020 C                          2D file containing depth of the ice front
                0021 C                          at each model grid cell
                0022 C     STIClengthFile    :: name of icefront length file (m/m^2)
                0023 C                          2D file containing the ratio of the horizontal
                0024 C                          length of the ice front in each model grid cell
                0025 C                          divided by the grid cell area
                0026 C-----------------------------------------------------------------------
                0027 CEOP
                0028
                0029       COMMON /STIC_PARMS_R/
                0030      &     icfThetaHorizDiffusionLength,
                0031      &     icfThetaInterior
                0032       _RL icfThetaHorizDiffusionLength
                0033       _RL icfThetaInterior
                0034
                0035       CHARACTER*(MAX_LEN_FNAM) STIClengthFile
                0036       CHARACTER*(MAX_LEN_FNAM) STICdepthFile
                0037       COMMON /STIC_PARM_C/
                0038      &     STIClengthFile,
                0039      &     STICdepthFile
                0040
                0041 C     ow - 06/29/2018
                0042 C     maskSHI of pkg/shelfice is not consistent with the spirit of gencost.
                0043 C       Use the following masks below instead.
                0044 C     mask2dSHIICF: 2d shelf-ice & ice-front mask:
                0045 C       1 for having shelf-ice and/or ice-front at one or more vertical levels
                0046 C       and 0 otherwise.
                0047 C     mask3dSHIICF: 3d shelf-ice & ice-front mask.
                0048 C     mask2dSHI: 2d shelf-ice mask
                0049 C     mask3dSHI: 3d shelf-ice mask
                0050 C     mask2dICF: 2d ice-front mask: 1 for having ice-front at one or more
                0051 C                vertical levels.
                0052 C     mask3dICF: 3d ice-front mask
                0053       COMMON /STIC_MASKS/ mask2dSHIICF, mask3dSHIICF,
                0054      &        mask2dSHI, mask3dSHI, mask2dICF, mask3dICF
                0055       _RS mask2dSHIICF(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0056       _RS mask3dSHIICF(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0057       _RS mask2dSHI   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0058       _RS mask3dSHI   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0059       _RS mask2dICF   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0060       _RS mask3dICF   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0061
                0062 C     INTEGER arrays
                0063 C     CURI_ARR       :: i-index for neighboring ice-front points
                0064 C     CURJ_ARR       :: j-index for neighboring ice-front points
                0065 C     sticfWidth_arr :: ice-front width in meters
                0066 C     kIcf           :: index of the bottommost ice front cell (2D)
                0067       COMMON /STIC_ICEFRONT_I/ kIcf, CURI_ARR, CURJ_ARR
                0068       INTEGER kIcf    (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0069       INTEGER CURI_ARR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy,4)
                0070       INTEGER CURJ_ARR(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy,4)
                0071
                0072       COMMON /STIC_FIELDS_RL/
                0073      &     stic_gT, stic_gS,
                0074      &     sticfWidth_arr
                0075       _RL stic_gT   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0076       _RL stic_gS   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0077       _RL sticfWidth_arr(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy,4)
                0078
                0079 #ifdef ALLOW_SHITRANSCOEFF_3D
                0080        COMMON /SHELFICE_GAMMA3D_RL/
                0081      &     shiTransCoeffT3D, shiTransCoeffS3D
                0082        _RL shiTransCoeffT3D(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0083        _RL shiTransCoeffS3D(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0084 #endif
                0085
                0086       COMMON /STIC_FIELDS_RS/
                0087      &     R_stic, sticfLength,
                0088      &     icfHeatFlux, icfFreshWaterFlux,
                0089      &     sticfHeatFlux, sticfFreshWaterFlux
                0090       _RS R_stic             (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0091       _RS sticfLength        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0092       _RL icfHeatFlux        (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0093       _RL icfFreshWaterFlux  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
                0094       _RL sticfHeatFlux      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0095       _RL sticfFreshWaterFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)