/[MITgcm]/MITgcm_contrib/verification_other/shelfice_remeshing/code/SHELFICE.h
ViewVC logotype

Diff of /MITgcm_contrib/verification_other/shelfice_remeshing/code/SHELFICE.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by dgoldberg, Fri Jan 22 18:09:42 2016 UTC revision 1.9 by jmc, Fri Dec 15 20:06:21 2017 UTC
# Line 18  C Line 18  C
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)
# Line 35  C     SHELFICEDynMassOnly      :: step i Line 32  C     SHELFICEDynMassOnly      :: step i
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
# Line 54  C     ---------------------------------- Line 48  C     ----------------------------------
48  C     SHELFICEuseGammaFrict    :: use velocity dependent exchange coefficients,  C     SHELFICEuseGammaFrict    :: use velocity dependent exchange coefficients,
49  C                                 see Holland and Jenkins (1999), eq.11-18,  C                                 see Holland and Jenkins (1999), eq.11-18,
50  C                                 with the following parameters (def: F):  C                                 with the following parameters (def: F):
51    C     SHELFICE_oldCalcUStar    :: use old uStar averaging expression
52  C     shiCdrag                 :: quadratic drag coefficient to compute uStar  C     shiCdrag                 :: quadratic drag coefficient to compute uStar
53  C                                 (def: 0.0015)  C                                 (def: 0.0015)
54  C     shiZetaN                 :: ??? (def: 0.052)  C     shiZetaN                 :: ??? (def: 0.052)
# Line 62  C     shiPrandtl, shiSchmidt   :: consta Line 57  C     shiPrandtl, shiSchmidt   :: consta
57  C                                 numbers used to compute gammaTurb  C                                 numbers used to compute gammaTurb
58  C     shiKinVisc               :: constant kinetic viscosity used to compute  C     shiKinVisc               :: constant kinetic viscosity used to compute
59  C                                 gammaTurb (def: 1.95e-5)  C                                 gammaTurb (def: 1.95e-5)
60  C     SHELFICERemeshFrequency  :: Frequency that size of etaN is checked to  C     SHELFICERemeshFrequency  :: Frequency that size of etaN is checked to
61  C                                 trigger remesh  C                                 trigger remesh
62  C     SHELFICESplitThreshold   :: Max size of etaN allowed before a remesh  C     SHELFICESplitThreshold   :: Max size of etaN allowed before a remesh
63  C     SHELFICEMergeThreshold   :: Min size of etaN allowed before a remesh  C     SHELFICEMergeThreshold   :: Min size of etaN allowed before a remesh
# Line 87  C Line 82  C
82  C--   Fields  C--   Fields
83  C     ktopC                  :: index of the top "wet cell" (2D)  C     ktopC                  :: index of the top "wet cell" (2D)
84  C     R_shelfIce             :: shelfice topography [m]  C     R_shelfIce             :: shelfice topography [m]
 C     R_Grounding            :: "grounded" ice topography [m]  
85  C     shelficeMassInit       :: ice-shelf mass (per unit area) (kg/m^2)  C     shelficeMassInit       :: ice-shelf mass (per unit area) (kg/m^2)
86  C     shelficeMass           :: ice-shelf mass (per unit area) (kg/m^2)  C     shelficeMass           :: ice-shelf mass (per unit area) (kg/m^2)
87  C     shelfIceMassDynTendency :: other mass balance tendency  (kg/m^2/s)  C     shelfIceMassDynTendency :: other mass balance tendency  (kg/m^2/s)
# Line 100  C     shelficeForcingT       :: analogue Line 94  C     shelficeForcingT       :: analogue
94  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)
95  C     shelficeForcingS       :: analogue of surfaceForcingS  C     shelficeForcingS       :: analogue of surfaceForcingS
96  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)
97    C     conserve_ssh           :: KS16. Use the obcs to conserve net open
98    C                               ocean eta to 0m
99  C-----------------------------------------------------------------------  C-----------------------------------------------------------------------
100  C \ev  C \ev
101  CEOP  CEOP
102    
103        COMMON /SHELFICE_PARMS_I/  kTopC,        COMMON /SHELFICE_PARMS_I/  kTopC,
104       &     SHELFICEselectDragQuadr       &     SHELFICEselectDragQuadr,
105         &     shelfice_etarestore_spongewidth
106        INTEGER kTopC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        INTEGER kTopC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
107        INTEGER SHELFICEselectDragQuadr        INTEGER SHELFICEselectDragQuadr
108          INTEGER shelfice_etarestore_spongewidth
109    
110        COMMON /SHELFICE_PARMS_R/        COMMON /SHELFICE_PARMS_R/
111       &     SHELFICE_dumpFreq, SHELFICE_taveFreq,       &     SHELFICE_dumpFreq, SHELFICE_taveFreq,
112       &     SHELFICEheatTransCoeff, SHELFICEsaltTransCoeff,       &     SHELFICEheatTransCoeff, SHELFICEsaltTransCoeff,
113       &     rhoShelfice, SHELFICEkappa,       &     rhoShelfice, SHELFICEkappa,
114       &     SHELFICElatentHeat,       &     SHELFICElatentHeat,
115       &     SHELFICEheatCapacity_Cp,       &     SHELFICEheatCapacity_Cp,
116       &     SHELFICEthetaSurface,       &     SHELFICEthetaSurface,
117       &     SHELFICEDragLinear, SHELFICEDragQuadratic,       &     SHELFICEDragLinear, SHELFICEDragQuadratic,
# Line 121  CEOP Line 119  CEOP
119       &     shiPrandtl, shiSchmidt, shiKinVisc,       &     shiPrandtl, shiSchmidt, shiKinVisc,
120       &     SHELFICERemeshFrequency,       &     SHELFICERemeshFrequency,
121       &     SHELFICESplitThreshold,       &     SHELFICESplitThreshold,
122       &     SHELFICEMergeThreshold       &     SHELFICEMergeThreshold,
123         &     shelficeEtaRelax
124        _RL SHELFICE_dumpFreq, SHELFICE_taveFreq        _RL SHELFICE_dumpFreq, SHELFICE_taveFreq
125        _RL SHELFICEheatTransCoeff        _RL SHELFICEheatTransCoeff
126        _RL SHELFICEsaltTransCoeff        _RL SHELFICEsaltTransCoeff
# Line 130  CEOP Line 129  CEOP
129        _RL rhoShelfice        _RL rhoShelfice
130        _RL SHELFICEkappa        _RL SHELFICEkappa
131        _RL SHELFICEDragLinear        _RL SHELFICEDragLinear
132        _RL SHELFICEDragQuadratic, SHELFICEMergeThreshold        _RL SHELFICEDragQuadratic
133          _RL SHELFICEMergeThreshold
134        _RL SHELFICEthetaSurface, SHELFICESplitThreshold        _RL SHELFICEthetaSurface, SHELFICESplitThreshold
135        _RL shiCdrag, shiZetaN, shiRc, SHELFICERemeshFrequency        _RL shiCdrag, shiZetaN, shiRc
136          _RL SHELFICERemeshFrequency
137        _RL shiPrandtl, shiSchmidt, shiKinVisc        _RL shiPrandtl, shiSchmidt, shiKinVisc
138        _RL SHELFICEGroundW, SHELFICEGroundC        _RL SHELFICEGroundW, SHELFICEGroundC, shelficeEtaRelax
139        COMMON /SHELFICE_FIELDS_RL/        COMMON /SHELFICE_FIELDS_RL/
140       &     shelficeMass, shelficeMassInit, shelficeGroundInit,       &     shelficeMass, shelficeMassInit,
141       &     shelficeLoadAnomaly,       &     shelficeLoadAnomaly,
142       &     shelficeForcingT, shelficeForcingS,       &     shelficeForcingT, shelficeForcingS,
143       &     shiTransCoeffT, shiTransCoeffS, EFFMASS       &     shiTransCoeffT, shiTransCoeffS, EFFMASS
144        _RL shelficeMass          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL shelficeMass          (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
145        _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)  
146        _RL shelficeLoadAnomaly   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL shelficeLoadAnomaly   (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
147        _RL shelficeForcingT      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL shelficeForcingT      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
148        _RL shelficeForcingS      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL shelficeForcingS      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
# Line 152  CEOP Line 152  CEOP
152        _RL SeaLevelRestore       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL SeaLevelRestore       (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
153    
154        COMMON /SHELFICE_FIELDS_RS/        COMMON /SHELFICE_FIELDS_RS/
155       &     R_shelfIce,       &     R_shelfIce, R_MWCT,
      &     R_Grounding,  
156       &     shelficeHeatFlux,       &     shelficeHeatFlux,
157       &     shelfIceFreshWaterFlux,       &     shelfIceFreshWaterFlux,
158       &     shelfIceMassDynTendency       &     shelfIceMassDynTendency
159        _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)
160        _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)
161        _RS shelficeHeatFlux      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS shelficeHeatFlux      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
162        _RS shelficeFreshWaterFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS shelficeFreshWaterFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
163        _RS        _RS
# Line 183  CEOP Line 182  CEOP
182        LOGICAL SHELFICE_tave_mnc        LOGICAL SHELFICE_tave_mnc
183        LOGICAL SHELFICEadvDiffHeatFlux        LOGICAL SHELFICEadvDiffHeatFlux
184        LOGICAL SHELFICEuseGammaFrict        LOGICAL SHELFICEuseGammaFrict
185          LOGICAL SHELFICE_oldCalcUStar
186        LOGICAL SHELFICEMassStepping        LOGICAL SHELFICEMassStepping
187        LOGICAL SHELFICEDynMassOnly        LOGICAL SHELFICEDynMassOnly
188          LOGICAL SHELFICEEtaSponge
189          LOGICAL SHELFICE_dig_ice
190          LOGICAL SHELFICE_massmin_truedens
191    C   KS16 put var here
192          LOGICAL conserve_ssh
193        COMMON /SHELFICE_PARMS_L/        COMMON /SHELFICE_PARMS_L/
194       &     SHELFICEisOn,       &     SHELFICEisOn,
195       &     useISOMIPTD,       &     useISOMIPTD,
# Line 200  CEOP Line 205  CEOP
205       &     SHELFICE_tave_mnc,       &     SHELFICE_tave_mnc,
206       &     SHELFICEadvDiffHeatFlux,       &     SHELFICEadvDiffHeatFlux,
207       &     SHELFICEuseGammaFrict,       &     SHELFICEuseGammaFrict,
208         &     SHELFICE_oldCalcUStar,
209       &     SHELFICEMassStepping,       &     SHELFICEMassStepping,
210       &     SHELFICEDynMassOnly       &     SHELFICEDynMassOnly,
211         &     SHELFICEEtaSponge,
212         &     SHELFICE_dig_ice,
213         &     SHELFICE_massmin_truedens,
214    C  KS16 and here;
215         &     conserve_ssh
216    
217        CHARACTER*(MAX_LEN_FNAM) SHELFICEloadAnomalyFile        CHARACTER*(MAX_LEN_FNAM) SHELFICEloadAnomalyFile
218        CHARACTER*(MAX_LEN_FNAM) SHELFICEmassFile        CHARACTER*(MAX_LEN_FNAM) SHELFICEmassFile
# Line 209  CEOP Line 220  CEOP
220        CHARACTER*(MAX_LEN_FNAM) SHELFICEtopoFile        CHARACTER*(MAX_LEN_FNAM) SHELFICEtopoFile
221        CHARACTER*(MAX_LEN_FNAM) SHELFICEMassDynTendFile        CHARACTER*(MAX_LEN_FNAM) SHELFICEMassDynTendFile
222        CHARACTER*(MAX_LEN_FNAM) SHELFICEGroundInitFile        CHARACTER*(MAX_LEN_FNAM) SHELFICEGroundInitFile
223        CHARACTER*(MAX_LEN_FNAM) SHELFICETransCoeffTFile        CHARACTER*(MAX_LEN_FNAM) SHELFICETransCoeffTFile
   
224    
225        COMMON /SHELFICE_PARM_C/        COMMON /SHELFICE_PARM_C/
226       &     SHELFICEloadAnomalyFile,       &     SHELFICEloadAnomalyFile,
# Line 219  CEOP Line 229  CEOP
229       &     SHELFICEGroundTopoFile,       &     SHELFICEGroundTopoFile,
230       &     SHELFICEMassDynTendFile,       &     SHELFICEMassDynTendFile,
231       &     SHELFICEGroundInitFile,       &     SHELFICEGroundInitFile,
232       &     SHELFICETransCoeffTFile       &     SHELFICETransCoeffTFile
233    
234  #endif /* ALLOW_SHELFICE */  #endif /* ALLOW_SHELFICE */
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22