/[MITgcm]/MITgcm/pkg/shelfice/SHELFICE.h
ViewVC logotype

Contents of /MITgcm/pkg/shelfice/SHELFICE.h

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


Revision 1.11 - (show annotations) (download)
Fri May 18 12:17:01 2012 UTC (13 years, 1 month ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64f, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63n, checkpoint63o, checkpoint64
Changes since 1.10: +8 -1 lines
File MIME type: text/plain
add parameterisation of advective-diffusive flux into the ice shelf
according to Holland and Jenkins (1999), eq.22-33, as suggested by Jan
De Ryd
Add some comments in SHELFICE.h about useGammaFrict

1 C $Header: /u/gcmpack/MITgcm/pkg/shelfice/SHELFICE.h,v 1.10 2012/01/06 13:45:59 heimbach Exp $
2 C $Name: $
3
4 #ifdef ALLOW_SHELFICE
5
6 CBOP
7 C !ROUTINE: SHELFICE.h
8
9 C !DESCRIPTION: \bv
10 C /==========================================================\
11 C | SHELFICE.h |
12 C | o Basic header thermodnynamic shelf ice package. |
13 C | Contains all SHELFICE field declarations. |
14 C \==========================================================/
15
16 C-----------------------------------------------------------------------
17 C
18 C-- Constants that can be set in data.shelfice
19 C SHELFICEtopoFile - File containing the topography of the
20 C shelfice draught (unit=m)
21 C SHELFICEloadAnomalyFile - name of shelfice load anomaly file
22 C SHELFICEDragLinear - linear drag at bottom shelfice (1/s)
23 C SHELFICEDragQuadratic - quadratic drag at bottom shelfice (1/m)
24 C SHELFICEheatTransCoeff - heat transfer coefficient that determines
25 C heat flux into shelfice (m/s)
26 C SHELFICEsaltTransCoeff - salinity transfer coefficient that determines
27 C salt flux into shelfice (m/s)
28 C SHELFICElatentHeat - latent heat of fusion (J/kg)
29 C useISOMIPTD - use simple ISOMIP thermodynamics
30 C SHELFICEconserve - use conservative form of H&O-thermodynamics
31 C following Jenkins et al. (2001, JPO)
32 C SHELFICEboundaryLayer - turn on vertical merging of cells to for a
33 C boundary layer of drF thickness
34 C SHELFICEadvDiffHeatFlux - use advective-diffusive heat flux into the ice shelf
35 C instead of diffusive heat flux (default), see Holland
36 C and Jenkins (1999), eq.21,22,26,31
37 C SHELFICEuseGammaFrict - use velocity dependent exchange coefficients, see Holland a
38 C and Jenkins (1999), eq.11-18
39 C no_slip_shelfice - set slip conditions for shelfice separately,
40 C (by default the same as no_slip_bottom)
41 C SHELFICEwriteState - enable output
42 C SHELFICE_dump_mnc - use netcdf for snapshot output
43 C SHELFICE_tave_mnc - use netcdf for time-averaged output
44 C SHELFICE_dumpFreq - analoguous to dumpFreq (= default)
45 C SHELFICE_taveFreq - analoguous to taveFreq (= default)
46 C
47 C-- Fields
48 C ktopC - index of the top "wet cell" (2D)
49 C R_shelfIce - shelfice topography [m]
50 C shelficeLoadAnomaly - pressure load anomaly of shelfice [Pa]
51 C shelficeHeatFlux - upward heat flux [W/m^2]
52 C shelficeFreshWaterFlux - upward fresh water flux (virt. salt flux) [m/s]
53 C shelficeForcingT - analogue of surfaceForcingT
54 C units are r_unit.Kelvin/s (=Kelvin.m/s if r=z)
55 C shelficeForcingS - analogue of surfaceForcingS
56 C units are r_unit.psu/s (=psu.m/s if r=z)
57 C-----------------------------------------------------------------------
58 C \ev
59 CEOP
60
61 COMMON /SHELFICE_PARMS_I/ kTopC
62 INTEGER kTopC (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
63
64 COMMON /SHELFICE_PARMS_R/
65 & SHELFICE_dumpFreq, SHELFICE_taveFreq,
66 & SHELFICEheatTransCoeff, SHELFICEsaltTransCoeff,
67 & rhoShelfice, SHELFICEkappa,
68 & SHELFICElatentHeat, recip_SHELFICElatentHeat,
69 & SHELFICEheatCapacity_Cp,
70 & SHELFICEthetaSurface,
71 & SHELFICEDragLinear, SHELFICEDragQuadratic,
72 & shiCdrag, shiZetaN, shiRc,
73 & shiPrandtl, shiSchmidt, shiKinVisc
74
75 _RL SHELFICE_dumpFreq, SHELFICE_taveFreq
76 _RL SHELFICEheatTransCoeff
77 _RL SHELFICEsaltTransCoeff
78 _RL SHELFICElatentHeat
79 _RL SHELFICEheatCapacity_Cp
80 _RL rhoShelfice
81 _RL SHELFICEkappa
82 _RL recip_SHELFICElatentHeat
83 _RL SHELFICEDragLinear
84 _RL SHELFICEDragQuadratic
85 _RL SHELFICEthetaSurface
86 _RL shiCdrag, shiZetaN, shiRc
87 _RL shiPrandtl, shiSchmidt, shiKinVisc
88
89
90 COMMON /SHELFICE_FIELDS_RL/
91 & shelficeForcingT, shelficeForcingS,
92 & shiTransCoeffT, shiTransCoeffS
93 _RL shelficeForcingT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
94 _RL shelficeForcingS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
95 _RL shiTransCoeffT (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
96 _RL shiTransCoeffS (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
97
98 COMMON /SHELFICE_FIELDS_RS/
99 & R_shelfIce,
100 & shelficeLoadAnomaly,
101 & shelficeHeatFlux,
102 & shelfIceFreshWaterFlux
103 _RS R_shelfIce (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
104 _RS shelficeLoadAnomaly (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
105 _RS shelficeHeatFlux (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
106 _RS shelficeFreshWaterFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
107
108 #ifdef ALLOW_SHIFWFLX_CONTROL
109 COMMON /SHELFICE_MASKS_CTRL/ maskSHI
110 _RS maskSHI (1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr,nSx,nSy)
111 #endif /* ALLOW_SHIFWFLX_CONTROL */
112
113 LOGICAL SHELFICEisOn
114 LOGICAL useISOMIPTD
115 LOGICAL SHELFICEconserve
116 LOGICAL SHELFICEboundaryLayer
117 LOGICAL no_slip_shelfice
118 LOGICAL SHELFICEwriteState
119 LOGICAL SHELFICE_dump_mdsio
120 LOGICAL SHELFICE_tave_mdsio
121 LOGICAL SHELFICE_dump_mnc
122 LOGICAL SHELFICE_tave_mnc
123 LOGICAL SHELFICEadvDiffHeatFlux
124 LOGICAL SHELFICEuseGammaFrict
125 COMMON /SHELFICE_PARMS_L/
126 & SHELFICEisOn,
127 & useISOMIPTD,
128 & SHELFICEconserve,
129 & SHELFICEboundaryLayer,
130 & no_slip_shelfice,
131 & SHELFICEwriteState,
132 & SHELFICE_dump_mdsio,
133 & SHELFICE_tave_mdsio,
134 & SHELFICE_dump_mnc,
135 & SHELFICE_tave_mnc,
136 & SHELFICEadvDiffHeatFlux,
137 & SHELFICEuseGammaFrict
138
139 CHARACTER*(MAX_LEN_FNAM) SHELFICEloadAnomalyFile
140 CHARACTER*(MAX_LEN_FNAM) SHELFICEtopoFile
141 COMMON /SHELFICE_PARM_C/
142 & SHELFICEloadAnomalyFile,
143 & SHELFICEtopoFile
144
145 #endif /* ALLOW_SHELFICE */

  ViewVC Help
Powered by ViewVC 1.1.22