18 |
C-- Constants that can be set in data.shelfice |
C-- Constants that can be set in data.shelfice |
19 |
C SHELFICEtopoFile :: File containing the topography of the |
C SHELFICEtopoFile :: File containing the topography of the |
20 |
C shelfice draught (unit=m) |
C shelfice draught (unit=m) |
|
C SHELFICEGroundTopoFile :: File containing the topography of "grounded" |
|
|
c ice (unit=m) |
|
21 |
C SHELFICEmassFile :: name of shelfice Mass file |
C SHELFICEmassFile :: name of shelfice Mass file |
|
C SHELFICEGroundInitFile :: name of GroundingMass file, equivalent to justC float ice at GroundTopo level |
|
22 |
C SHELFICEloadAnomalyFile :: name of shelfice load anomaly file |
C SHELFICEloadAnomalyFile :: name of shelfice load anomaly file |
23 |
C SHELFICEMassDynTendFile :: file name for other mass tendency |
C SHELFICEMassDynTendFile :: file name for other mass tendency |
24 |
C (e.g. dynamics) |
C (e.g. dynamics) |
32 |
C (not melting/freezing) def: F |
C (not melting/freezing) def: F |
33 |
C SHELFICEboundaryLayer :: turn on vertical merging of cells to for a |
C SHELFICEboundaryLayer :: turn on vertical merging of cells to for a |
34 |
C boundary layer of drF thickness, def: F |
C boundary layer of drF thickness, def: F |
|
C SHELFICEthickboundaryLayer |
|
|
C :: similar to boundary layer but allow for a thicker |
|
|
C layer, encompassing 3 cells in vertical, def: F |
|
35 |
C SHELFICErealFWflux :: ensure vert advective flux at bdry uses top cell |
C SHELFICErealFWflux :: ensure vert advective flux at bdry uses top cell |
36 |
C value rather than "boundary layer" value F |
C value rather than "boundary layer" value F |
37 |
C SHELFICEadvDiffHeatFlux :: use advective-diffusive heat flux into the |
C SHELFICEadvDiffHeatFlux :: use advective-diffusive heat flux into the |
81 |
C-- Fields |
C-- Fields |
82 |
C ktopC :: index of the top "wet cell" (2D) |
C ktopC :: index of the top "wet cell" (2D) |
83 |
C R_shelfIce :: shelfice topography [m] |
C R_shelfIce :: shelfice topography [m] |
|
C R_Grounding :: "grounded" ice topography [m] |
|
84 |
C shelficeMassInit :: ice-shelf mass (per unit area) (kg/m^2) |
C shelficeMassInit :: ice-shelf mass (per unit area) (kg/m^2) |
85 |
C shelficeMass :: ice-shelf mass (per unit area) (kg/m^2) |
C shelficeMass :: ice-shelf mass (per unit area) (kg/m^2) |
86 |
C shelfIceMassDynTendency :: other mass balance tendency (kg/m^2/s) |
C shelfIceMassDynTendency :: other mass balance tendency (kg/m^2/s) |
93 |
C units are r_unit.Kelvin/s (=Kelvin.m/s if r=z) |
C units are r_unit.Kelvin/s (=Kelvin.m/s if r=z) |
94 |
C shelficeForcingS :: analogue of surfaceForcingS |
C shelficeForcingS :: analogue of surfaceForcingS |
95 |
C units are r_unit.psu/s (=psu.m/s if r=z) |
C units are r_unit.psu/s (=psu.m/s if r=z) |
96 |
|
C conserve_ssh :: KS16. Use the obcs to conserve net open |
97 |
|
C ocean eta to 0m |
98 |
C----------------------------------------------------------------------- |
C----------------------------------------------------------------------- |
99 |
C \ev |
C \ev |
100 |
CEOP |
CEOP |
101 |
|
|
102 |
COMMON /SHELFICE_PARMS_I/ kTopC, |
COMMON /SHELFICE_PARMS_I/ kTopC, |
103 |
& SHELFICEselectDragQuadr |
& SHELFICEselectDragQuadr, |
104 |
|
& shelfice_etarestore_spongewidth |
105 |
INTEGER kTopC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
INTEGER kTopC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
106 |
INTEGER SHELFICEselectDragQuadr |
INTEGER SHELFICEselectDragQuadr |
107 |
|
INTEGER shelfice_etarestore_spongewidth |
108 |
|
|
109 |
COMMON /SHELFICE_PARMS_R/ |
COMMON /SHELFICE_PARMS_R/ |
110 |
& SHELFICE_dumpFreq, SHELFICE_taveFreq, |
& SHELFICE_dumpFreq, SHELFICE_taveFreq, |
119 |
& SHELFICERemeshFrequency, |
& SHELFICERemeshFrequency, |
120 |
& SHELFICESplitThreshold, |
& SHELFICESplitThreshold, |
121 |
& SHELFICEMergeThreshold, |
& SHELFICEMergeThreshold, |
122 |
|
& shelficeEtaRelax |
123 |
_RL SHELFICE_dumpFreq, SHELFICE_taveFreq |
_RL SHELFICE_dumpFreq, SHELFICE_taveFreq |
124 |
_RL SHELFICEheatTransCoeff |
_RL SHELFICEheatTransCoeff |
125 |
_RL SHELFICEsaltTransCoeff |
_RL SHELFICEsaltTransCoeff |
128 |
_RL rhoShelfice |
_RL rhoShelfice |
129 |
_RL SHELFICEkappa |
_RL SHELFICEkappa |
130 |
_RL SHELFICEDragLinear |
_RL SHELFICEDragLinear |
131 |
_RL SHELFICEDragQuadratic, SHELFICEMergeThreshold |
_RL SHELFICEDragQuadratic |
132 |
|
_RL SHELFICEMergeThreshold |
133 |
_RL SHELFICEthetaSurface, SHELFICESplitThreshold |
_RL SHELFICEthetaSurface, SHELFICESplitThreshold |
134 |
_RL shiCdrag, shiZetaN, shiRc, SHELFICERemeshFrequency |
_RL shiCdrag, shiZetaN, shiRc |
135 |
|
_RL SHELFICERemeshFrequency |
136 |
_RL shiPrandtl, shiSchmidt, shiKinVisc |
_RL shiPrandtl, shiSchmidt, shiKinVisc |
137 |
_RL SHELFICEGroundW, SHELFICEGroundC |
_RL SHELFICEGroundW, SHELFICEGroundC, shelficeEtaRelax |
138 |
COMMON /SHELFICE_FIELDS_RL/ |
COMMON /SHELFICE_FIELDS_RL/ |
139 |
& shelficeMass, shelficeMassInit, shelficeGroundInit, |
& shelficeMass, shelficeMassInit, |
140 |
& shelficeLoadAnomaly, |
& shelficeLoadAnomaly, |
141 |
& shelficeForcingT, shelficeForcingS, |
& shelficeForcingT, shelficeForcingS, |
142 |
& shiTransCoeffT, shiTransCoeffS, EFFMASS |
& shiTransCoeffT, shiTransCoeffS, EFFMASS |
143 |
_RL shelficeMass (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL shelficeMass (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
144 |
_RL shelficeMassInit (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL shelficeMassInit (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
_RL shelficeGroundInit (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
|
145 |
_RL shelficeLoadAnomaly (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL shelficeLoadAnomaly (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
146 |
_RL shelficeForcingT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL shelficeForcingT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
147 |
_RL shelficeForcingS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL shelficeForcingS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
151 |
_RL SeaLevelRestore (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RL SeaLevelRestore (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
152 |
|
|
153 |
COMMON /SHELFICE_FIELDS_RS/ |
COMMON /SHELFICE_FIELDS_RS/ |
154 |
& R_shelfIce, |
& R_shelfIce, R_MWCT, |
|
& R_Grounding, |
|
155 |
& shelficeHeatFlux, |
& shelficeHeatFlux, |
156 |
& shelfIceFreshWaterFlux, |
& shelfIceFreshWaterFlux, |
157 |
& shelfIceMassDynTendency, |
& shelfIceMassDynTendency |
158 |
_RS R_shelfIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS R_shelfIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
159 |
_RS R_Grounding (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS R_MWCT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
160 |
_RS shelficeHeatFlux (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS shelficeHeatFlux (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
161 |
_RS shelficeFreshWaterFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS shelficeFreshWaterFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
162 |
_RS |
_RS |
183 |
LOGICAL SHELFICEuseGammaFrict |
LOGICAL SHELFICEuseGammaFrict |
184 |
LOGICAL SHELFICEMassStepping |
LOGICAL SHELFICEMassStepping |
185 |
LOGICAL SHELFICEDynMassOnly |
LOGICAL SHELFICEDynMassOnly |
186 |
|
LOGICAL SHELFICEEtaSponge |
187 |
|
LOGICAL SHELFICE_dig_ice |
188 |
|
C KS16 put var here |
189 |
|
LOGICAL conserve_ssh |
190 |
COMMON /SHELFICE_PARMS_L/ |
COMMON /SHELFICE_PARMS_L/ |
191 |
& SHELFICEisOn, |
& SHELFICEisOn, |
192 |
& useISOMIPTD, |
& useISOMIPTD, |
203 |
& SHELFICEadvDiffHeatFlux, |
& SHELFICEadvDiffHeatFlux, |
204 |
& SHELFICEuseGammaFrict, |
& SHELFICEuseGammaFrict, |
205 |
& SHELFICEMassStepping, |
& SHELFICEMassStepping, |
206 |
& SHELFICEDynMassOnly |
& SHELFICEDynMassOnly, |
207 |
|
& SHELFICEEtaSponge, |
208 |
|
& SHELFICE_dig_ice, |
209 |
|
C KS16 and here; |
210 |
|
& conserve_ssh |
211 |
|
|
212 |
CHARACTER*(MAX_LEN_FNAM) SHELFICEloadAnomalyFile |
CHARACTER*(MAX_LEN_FNAM) SHELFICEloadAnomalyFile |
213 |
CHARACTER*(MAX_LEN_FNAM) SHELFICEmassFile |
CHARACTER*(MAX_LEN_FNAM) SHELFICEmassFile |