/[MITgcm]/MITgcm_contrib/torge/itd/code/SEAICE_PARAMS.h
ViewVC logotype

Annotation of /MITgcm_contrib/torge/itd/code/SEAICE_PARAMS.h

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


Revision 1.4 - (hide annotations) (download)
Mon Oct 22 19:26:42 2012 UTC (12 years, 9 months ago) by heimbach
Branch: MAIN
Changes since 1.3: +9 -11 lines
File MIME type: text/plain
Separate common blocks for ITD vs non-ITD

1 heimbach 1.4 C $Header: /u/gcmpack/MITgcm_contrib/torge/itd/code/SEAICE_PARAMS.h,v 1.3 2012/10/22 19:21:34 torge Exp $
2 dimitri 1.1 C $Name: $
3    
4     C *==========================================================*
5     C | SEAICE_PARAMS.h
6     C | o Basic parameter header for sea ice model.
7     C *==========================================================*
8    
9     C-- COMMON /SEAICE_PARM_L/ Logical parameters of sea ice model.
10     C
11     C SEAICEwriteState :: If true, write sea ice state to file;
12     C default is false.
13     C SEAICEuseDYNAMICS :: If false, do not use dynamics;
14     C default is to use dynamics.
15     C SEAICEuseTEM :: to use truncated ellipse method (see Geiger et al.
16     C 1998) set this parameter to true, default is false
17     C SEAICEuseEVP :: If false, use Zhangs LSR solver for VP equations
18     C if true use elastic viscous plastic solver
19     C SEAICEuseFREEDRIFT :: If True use free drift velocity instead of EVP
20     C or LSR
21 torge 1.3 C SEAICEuseJFNK :: If true, use Jacobi-free Newton-Krylov solver
22     C instead of LSR (default: false)
23     C SEAICEuseTILT :: If true then include surface tilt term in dynamics
24 dimitri 1.1 C SEAICEheatConsFix :: If true then fix ocn<->seaice advective heat flux.
25     C SEAICEuseEVPpickup :: Set to false in order to start EVP solver with
26     C non-EVP pickup files. Default is true.
27     C Applied only if SEAICEuseEVP=.TRUE.
28     C SEAICEuseFluxForm :: use flux form for advection and diffusion
29     C of seaice
30     C SEAICEuseMetricTerms :: use metric terms for dynamics solver
31     C (default = .true. )
32     C SEAICEuseFlooding :: turn on scheme to convert submerged snow into ice
33     C SEAICEadvHeff :: turn on advection of effective thickness
34     C (default = .true.)
35     C SEAICEadvArea :: turn on advection of fraction area
36     C (default = .true.)
37     C SEAICEadvSnow :: turn on advection of snow (does not work with
38     C non-default Leap-frog scheme for advection)
39     C SEAICEadvSalt :: turn on advection of salt (does not work with
40     C non-default Leap-frog scheme for advection)
41     C useHB87stressCoupling :: use an intergral over ice and ocean surface
42     C layer to define surface stresses on ocean
43     C following Hibler and Bryan (1987, JPO)
44     C usePW79thermodynamics :: use "0-layer" thermodynamics as described in
45     C Parkinson and Washington (1979) and Hibler (1979)
46     C useMaykutSatVapPoly :: use Maykut Polynomial for saturation vapor pressure
47     C instead of extended temp-range exponential law; def=F.
48     C SEAICE_mcPheeStepFunc :: use step function (not linear tapering) in
49     C ocean-ice turbulent flux
50 torge 1.3 C SEAICE_useMultDimSnow :: use same fixed pdf for snow as for MULITCATEGORY ICE
51 dimitri 1.1 C SEAICE_doOpenWaterGrowth :: use open water heat flux directly to grow ice
52     C (when false cool ocean, and grow later if needed)
53     C SEAICE_doOpenWaterMelt :: use open water heat flux directly to melt ice
54     C (when false warm ocean, and melt later if needed)
55     C SEAICE_salinityTracer :: use SItracer to exchange and trace ocean
56     C salt in ice
57     C SEAICE_age Tracer :: use SItracer to trace the age of ice
58     C SEAICErestoreUnderIce :: restore surface T/S also underneath ice
59     C ( default is false )
60     C SEAICE_no_slip :: apply no slip boundary conditions to seaice velocity
61     C SEAICE_clipVelocities :: clip velocities to +/- 40cm/s
62     C SEAICE_maskRHS :: mask the RHS of the solver where there is no ice
63     C SEAICE_tave_mdsio :: write TimeAverage output using MDSIO
64     C SEAICE_dump_mdsio :: write snap-shot output using MDSIO
65     C SEAICE_mon_stdio :: write monitor to std-outp
66     C SEAICE_tave_mnc :: write TimeAverage output using MNC
67     C SEAICE_dump_mnc :: write snap-shot output using MNC
68     C SEAICE_mon_mnc :: write monitor to netcdf file
69     LOGICAL
70     & SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseEVP,
71 torge 1.3 & SEAICEuseFREEDRIFT, SEAICEuseTILT, SEAICEuseTEM,
72     & SEAICEuseJFNK,
73 dimitri 1.1 & SEAICEheatConsFix,
74     & SEAICEuseMetricTerms,
75     & SEAICEuseEVPpickup, SEAICEuseFlooding,
76     & SEAICEadvHeff, SEAICEadvArea,
77     & SEAICEadvSnow, SEAICEadvSalt,
78     & SEAICEuseFluxForm, useHB87stressCoupling,
79     & usePW79thermodynamics, useMaykutSatVapPoly,
80 torge 1.3 & SEAICE_mcPheeStepFunc, SEAICE_useMultDimSnow,
81 dimitri 1.1 & SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt,
82     & SEAICE_salinityTracer, SEAICE_ageTracer,
83     & SEAICErestoreUnderIce,
84     & SEAICE_no_slip, SEAICE_clipVelocities, SEAICE_maskRHS,
85     & SEAICE_tave_mdsio, SEAICE_dump_mdsio, SEAICE_mon_stdio,
86     & SEAICE_tave_mnc, SEAICE_dump_mnc, SEAICE_mon_mnc
87     COMMON /SEAICE_PARM_L/
88     & SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseEVP,
89 torge 1.3 & SEAICEuseFREEDRIFT, SEAICEuseTILT, SEAICEuseTEM,
90     & SEAICEuseJFNK,
91 dimitri 1.1 & SEAICEheatConsFix,
92     & SEAICEuseMetricTerms,
93     & SEAICEuseEVPpickup, SEAICEuseFlooding,
94     & SEAICEadvHeff, SEAICEadvArea,
95     & SEAICEadvSnow, SEAICEadvSalt,
96     & SEAICEuseFluxForm, useHB87stressCoupling,
97     & usePW79thermodynamics, useMaykutSatVapPoly,
98 torge 1.3 & SEAICE_mcPheeStepFunc, SEAICE_useMultDimSnow,
99 dimitri 1.1 & SEAICE_doOpenWaterGrowth, SEAICE_doOpenWaterMelt,
100     & SEAICE_salinityTracer, SEAICE_ageTracer,
101     & SEAICErestoreUnderIce,
102     & SEAICE_no_slip, SEAICE_clipVelocities, SEAICE_maskRHS,
103     & SEAICE_tave_mdsio, SEAICE_dump_mdsio, SEAICE_mon_stdio,
104     & SEAICE_tave_mnc, SEAICE_dump_mnc, SEAICE_mon_mnc
105    
106     C-- COMMON /SEAICE_PARM_I/ Integer valued parameters of sea ice model.
107     C IMAX_TICE :: number of iterations for ice surface temp (default=10)
108     C postSolvTempIter :: select flux calculation after surf. temp solver iteration
109     C 0 = none, i.e., from last iter ; 2 = full non-lin form
110     C 1 = use linearized approx (consistent with tsurf finding)
111     C SOLV_MAX_ITERS :: maximum number of allowed LSR-solver iterations
112     C SOLV_NCHECK :: iteration interval for solver convergence test
113     C NPSEUDOTIMESTEPS :: number of extra pseudo time steps (>= 2)
114 torge 1.3 C SEAICEnewtonIterMax :: maximum number of allowed Newton iterations
115     C in JFNK-solver
116     C SEAICEkrylovIterMax :: maximum number of allowed Krylov iterations
117     C in JFNK-solver
118 dimitri 1.1 C LSR_mixIniGuess :: control mixing of free-drift sol. into LSR initial guess
119     C :: =0 : no mix ; =2,4 : mix with (1/err)^2,4 factor
120 torge 1.3 C SEAICEpresPow0 :: HEFF exponent for ice strength below SEAICEpresH0
121     C SEAICEpresPow1 :: HEFF exponent for ice strength above SEAICEpresH0
122 dimitri 1.1 C SEAICEadvScheme :: sets the advection scheme for thickness and area
123     C SEAICEadvSchArea :: sets the advection scheme for area
124     C SEAICEadvSchHeff :: sets the advection scheme for effective thickness
125     C (=volume), snow thickness, and salt if available
126     C SEAICEadvSchSnow :: sets the advection scheme for snow on sea-ice
127     C SEAICEadvSchSalt :: sets the advection scheme for sea ice salinity
128     C SEAICEadvSchSnow :: sets the advection scheme for snow on sea-ice
129     C SEAICE_areaLossFormula :: selects formula for ice cover loss from melt
130     C :: 1=from all but only melt conributions by ATM and OCN
131     C :: 2=from net melt-growth>0 by ATM and OCN
132     C :: 3=from predicted melt by ATM
133     C SEAICE_areaGainFormula :: selects formula for ice cover gain from open water growth
134     C :: 1=from growth by ATM
135     C :: 2=from predicted growth by ATM
136     C SEAICE_multDim :: number of ice categories
137     C SEAICE_debugPointI :: I,J index for seaice-specific debuggin
138     C SEAICE_debugPointJ
139     C
140     INTEGER IMAX_TICE, postSolvTempIter
141     INTEGER SOLV_MAX_ITERS, SOLV_NCHECK
142     INTEGER NPSEUDOTIMESTEPS
143     INTEGER LSR_mixIniGuess
144 torge 1.3 INTEGER SEAICEnewtonIterMax, SEAICEkrylovIterMax
145 dimitri 1.1 INTEGER SEAICEadvScheme
146     INTEGER SEAICEadvSchArea
147     INTEGER SEAICEadvSchHeff
148     INTEGER SEAICEadvSchSnow
149     INTEGER SEAICEadvSchSalt
150     INTEGER SEAICEadjMODE
151     INTEGER SEAICE_areaLossFormula
152     INTEGER SEAICE_areaGainFormula
153     INTEGER SEAICE_multDim
154     INTEGER SEAICE_debugPointI
155     INTEGER SEAICE_debugPointJ
156 torge 1.3 INTEGER SEAICEpresPow0, SEAICEpresPow1
157 dimitri 1.1 COMMON /SEAICE_PARM_I/
158     & IMAX_TICE, postSolvTempIter,
159     & SOLV_MAX_ITERS, SOLV_NCHECK,
160     & NPSEUDOTIMESTEPS,
161     & LSR_mixIniGuess,
162 torge 1.3 & SEAICEnewtonIterMax, SEAICEkrylovIterMax,
163     & SEAICEpresPow0, SEAICEpresPow1,
164 dimitri 1.1 & SEAICEadvScheme,
165     & SEAICEadvSchArea,
166     & SEAICEadvSchHeff,
167     & SEAICEadvSchSnow,
168     & SEAICEadvSchSalt,
169     & SEAICEadjMODE,
170     & SEAICE_areaLossFormula,
171     & SEAICE_areaGainFormula,
172     & SEAICE_multDim,
173     & SEAICE_debugPointI,
174     & SEAICE_debugPointJ
175    
176     C-- COMMON /SEAICE_PARM_C/ Character valued sea ice model parameters.
177     C AreaFile :: File containing initial sea-ice concentration
178     C HsnowFile :: File containing initial snow thickness
179     C HsaltFile :: File containing initial sea ice salt content
180     C HeffFile :: File containing initial sea-ice thickness
181     C uIceFile :: File containing initial sea-ice U comp. velocity
182     C vIceFile :: File containing initial sea-ice V comp. velocity
183     C !!! NOTE !!! Initial sea-ice thickness can also be set using
184     C SEAICE_initialHEFF below. But a constant initial condition
185     C can mean large artificial fluxes of heat and freshwater in
186     C the surface layer during the first model time step.
187     C
188     CHARACTER*(MAX_LEN_FNAM) AreaFile
189     CHARACTER*(MAX_LEN_FNAM) HsnowFile
190     CHARACTER*(MAX_LEN_FNAM) HsaltFile
191     CHARACTER*(MAX_LEN_FNAM) HeffFile
192     CHARACTER*(MAX_LEN_FNAM) uIceFile
193     CHARACTER*(MAX_LEN_FNAM) vIceFile
194     COMMON /SEAICE_PARM_C/
195     & AreaFile, HsnowFile, HsaltFile, HeffFile,
196     & uIceFile, vIceFile
197    
198     C-- COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model.
199     C SEAICE_deltaTtherm :: Seaice timestep for thermodynamic equations (s)
200     C SEAICE_deltaTdyn :: Seaice timestep for dynamic solver (s)
201     C SEAICE_deltaTevp :: Seaice timestep for EVP solver (s)
202     C SEAICE_elasticParm :: parameter that sets relaxation timescale
203     C tau = SEAICE_elasticParm * SEAICE_deltaTdyn
204     C SEAICE_evpTauRelax :: relaxation timescale tau (s)
205     C SEAICE_evpDampC :: evp damping constant (Hunke,JCP,2001) (kg/m^2)
206 torge 1.3 C JFNKgamma_nonlin :: non-linear tolerance parameter for JFNK solver
207     C JFNKgamma_lin_min/max :: tolerance parameters for linear JFNK solver
208     C JFNKres_t :: tolerance parameter for FGMRES residual
209 dimitri 1.1 C SEAICE_zetaMaxFac :: factor determining the maximum viscosity (s)
210     C (default = 5.e+12/2.e4 = 2.5e8)
211     C SEAICE_zetaMin :: lower bound for viscosity (default = 0) (N s/m^2)
212 torge 1.3 C SEAICEpresH0 :: HEFF threshold for ice strength (m)
213 dimitri 1.1 C SEAICE_monFreq :: SEAICE monitor frequency. (s)
214     C SEAICE_dumpFreq :: SEAICE dump frequency. (s)
215     C SEAICE_taveFreq :: SEAICE time-averaging frequency. (s)
216     C SEAICE_initialHEFF :: initial sea-ice thickness (m)
217     C SEAICE_rhoAir :: density of air (kg/m^3)
218     C SEAICE_rhoIce :: density of sea ice (kg/m^3)
219     C SEAICE_rhoSnow :: density of snow (kg/m^3)
220     C ICE2WATR :: ratio of sea ice density to water density
221     C OCEAN_drag :: air-ocean drag coefficient
222     C SEAICE_cpAir :: specific heat of air (J/kg/K)
223     C
224     C SEAICE_drag :: air-ice drag coefficient
225     C SEAICE_waterDrag :: water-ice drag coefficient * water density
226     C SEAICE_dryIceAlb :: winter albedo
227     C SEAICE_wetIceAlb :: summer albedo
228     C SEAICE_drySnowAlb :: dry snow albedo
229     C SEAICE_wetSnowAlb :: wet snow albedo
230     C HO :: AKA "lead closing parameter", demarcation thickness
231     C between thin and thick ice. Alternatively, HO (in
232     C meters) can be interpreted as the thickness of ice
233     C formed in open water.
234     C HO is a key ice-growth parameter that determines
235     C the partition between vertical and lateral growth.
236     C The default is 0.5m, increasing this value leads
237     C slower formation of a closed ice cover and thus to
238     C more ice (and thicker) ice, decreasing to faster
239     C formation of a closed ice cover (leads are closing
240     C faster) and thus less (thinner) ice.
241     C
242     C SEAICE_drag_south :: Southern Ocean SEAICE_drag
243     C SEAICE_waterDrag_south :: Southern Ocean SEAICE_waterDrag
244     C SEAICE_dryIceAlb_south :: Southern Ocean SEAICE_dryIceAlb
245     C SEAICE_wetIceAlb_south :: Southern Ocean SEAICE_wetIceAlb
246     C SEAICE_drySnowAlb_south :: Southern Ocean SEAICE_drySnowAlb
247     C SEAICE_wetSnowAlb_south :: Southern Ocean SEAICE_wetSnowAlb
248     C HO_south :: Southern Ocean HO
249     C
250     C SEAICE_wetAlbTemp :: Temp (deg.C) above which wet-albedo values are used
251     C SEAICE_waterAlbedo :: water albedo
252     C SEAICE_strength :: sea-ice strength Pstar
253     C SEAICE_eccen :: sea-ice eccentricity of the elliptical yield curve
254     C SEAICE_lhFusion :: latent heat of fusion for ice and snow (J/kg)
255     C SEAICE_lhEvap :: latent heat of evaporation for water (J/kg)
256     C SEAICE_dalton :: Dalton number (= sensible heat transfer coefficient)
257     C SEAICE_iceConduct :: sea-ice conductivity
258     C SEAICE_snowConduct :: snow conductivity
259     C SEAICE_emissivity :: longwave ocean-surface emissivity (-)
260     C SEAICE_ice_emiss :: longwave ice-surface emissivity (-)
261     C SEAICE_snow_emiss :: longwave snow-surface emissivity (-)
262     C SEAICE_boltzmann :: Stefan-Boltzman constant (not a run time parameter)
263     C SEAICE_snowThick :: cutoff snow thickness (for snow-albedo)
264     C SEAICE_shortwave :: ice penetration shortwave radiation factor
265     C SEAICE_saltFrac :: salinity of newly formed seaice defined as a
266     C fraction of the ocean surface salinity at the time
267     C of freezing
268     C SEAICE_salt0 :: prescribed salinity of seaice (in g/kg).
269     C facOpenGrow :: 0./1. version of logical SEAICE_doOpenWaterGrowth
270     C facOpenMelt :: 0./1. version of logical SEAICE_doOpenWaterMelt
271     C SEAICE_mcPheePiston:: ocean-ice turbulent flux "piston velocity" (m/s)
272     C that sets melt efficiency.
273     C SEAICE_mcPheeTaper :: tapering down of turbulent flux term with ice
274     C concentration. The 100% cover turb. flux is
275     C multiplied by 1.-SEAICE_mcPheeTaper
276     C SEAICE_frazilFrac :: Fraction of surface level negative heat content
277     C anomalies (relative to the local freezing point)
278     C may contribute as frazil over one time step.
279     C SEAICE_tempFrz0 :: sea water freezing point is
280     C SEAICE_dTempFrz_dS :: tempFrz = SEAICE_tempFrz0 + salt*SEAICE_dTempFrz_dS
281     C SEAICEstressFactor :: factor by which ice affects wind stress (default=1)
282     C LSR_ERROR :: sets accuracy of LSR solver
283     C DIFF1 :: parameter used in advect.F
284     C SEAICE_area_max :: usually set to 1. Seeting areaMax below 1 specifies
285     C the minimun amount of leads (1-areaMax) in the
286     C ice pack.
287     C SEAICE_area_floor :: usually set to 1x10^-5. Specifies a minimun
288     C ice fraction in the ice pack.
289     C SEAICE_area_reg :: usually set to 1x10^-5. Specifies a minimun
290     C ice fraction for the purposes of regularization
291     C SEAICE_hice_reg :: usually set to 5 cm. Specifies a minimun
292     C ice thickness for the purposes of regularization
293     C SEAICEdiffKhArea :: sets the diffusivity for area (m^2/s)
294     C SEAICEdiffKhHeff :: sets the diffusivity for effective thickness (m^2/s)
295     C SEAICEdiffKhSnow :: sets the diffusivity for snow on sea-ice (m^2/s)
296     C SEAICEdiffKhSalt :: sets the diffusivity for sea ice salinity (m^2/s)
297     C SEAICE_airTurnAngle :: turning angles of air-ice interfacial stress
298     C SEAICE_waterTurnAngle :: and ice-water interfacial stress (in degrees)
299 dimitri 1.2 C SEAICE_tauAreaObsRelax :: Timescale of relaxation to observed
300     C sea ice concentration (s), default=unset
301 dimitri 1.1 C
302     _RL SEAICE_deltaTtherm, SEAICE_deltaTdyn, SEAICE_deltaTevp
303     _RL SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq
304     _RL SEAICE_initialHEFF
305     _RL SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR
306     _RL SEAICE_cpAir
307     _RL SEAICE_drag, SEAICE_waterDrag, SEAICE_dryIceAlb
308     _RL SEAICE_wetIceAlb, SEAICE_drySnowAlb, SEAICE_wetSnowAlb, HO
309     _RL SEAICE_drag_south, SEAICE_waterDrag_south
310     _RL SEAICE_dryIceAlb_south, SEAICE_wetIceAlb_south
311     _RL SEAICE_drySnowAlb_south, SEAICE_wetSnowAlb_south, HO_south
312     _RL SEAICE_wetAlbTemp, SEAICE_waterAlbedo
313     _RL SEAICE_strength, SEAICE_eccen
314     _RL SEAICE_lhFusion, SEAICE_lhEvap
315     _RL SEAICE_dalton
316     _RL SEAICE_iceConduct, SEAICE_snowConduct
317     _RL SEAICE_emissivity, SEAICE_ice_emiss, SEAICE_snow_emiss
318     _RL SEAICE_boltzmann
319     _RL SEAICE_snowThick, SEAICE_shortwave
320     _RL SEAICE_saltFrac, SEAICE_salt0, SEAICEstressFactor
321     _RL SEAICE_mcPheeTaper, SEAICE_mcPheePiston
322     _RL SEAICE_frazilFrac, SEAICE_availHeatFrac
323     _RL facOpenGrow, facOpenMelt
324     _RL SEAICE_tempFrz0, SEAICE_dTempFrz_dS
325     _RL OCEAN_drag, LSR_ERROR, DIFF1
326 torge 1.3 _RL JFNKgamma_nonlin, JFNKres_t
327     _RL JFNKgamma_lin_min, JFNKgamma_lin_max
328 dimitri 1.1 _RL SEAICE_area_reg, SEAICE_hice_reg
329     _RL SEAICE_area_floor, SEAICE_area_max
330     _RL SEAICE_airTurnAngle, SEAICE_waterTurnAngle
331     _RL SEAICE_elasticParm, SEAICE_evpTauRelax
332     _RL SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac
333 torge 1.3 _RL SEAICEpresH0
334 dimitri 1.1 _RL SEAICEdiffKhArea, SEAICEdiffKhHeff, SEAICEdiffKhSnow
335     _RL SEAICEdiffKhSalt
336 dimitri 1.2 _RL SEAICE_tauAreaObsRelax
337 dimitri 1.1
338     COMMON /SEAICE_PARM_RL/
339     & SEAICE_deltaTtherm, SEAICE_deltaTdyn,
340     & SEAICE_deltaTevp, SEAICE_elasticParm, SEAICE_evpTauRelax,
341     & SEAICE_evpDampC, SEAICE_zetaMin, SEAICE_zetaMaxFac,
342 torge 1.3 & SEAICEpresH0,
343 dimitri 1.1 & SEAICE_monFreq, SEAICE_dumpFreq, SEAICE_taveFreq,
344     & SEAICE_initialHEFF,
345     & SEAICE_rhoAir, SEAICE_rhoIce, SEAICE_rhoSnow, ICE2WATR,
346     & SEAICE_drag, SEAICE_waterDrag, SEAICE_dryIceAlb,
347     & SEAICE_wetIceAlb, SEAICE_drySnowAlb, SEAICE_wetSnowAlb, HO,
348     & SEAICE_drag_south, SEAICE_waterDrag_south,
349     & SEAICE_dryIceAlb_south, SEAICE_wetIceAlb_south,
350     & SEAICE_drySnowAlb_south, SEAICE_wetSnowAlb_south, HO_south,
351     & SEAICE_wetAlbTemp, SEAICE_waterAlbedo,
352     & SEAICE_strength, SEAICE_eccen,
353     & SEAICE_lhFusion, SEAICE_lhEvap,
354     & SEAICE_dalton, SEAICE_cpAir,
355     & SEAICE_iceConduct, SEAICE_snowConduct,
356     & SEAICE_emissivity, SEAICE_ice_emiss, SEAICE_snow_emiss,
357     & SEAICE_boltzmann,
358     & SEAICE_snowThick, SEAICE_shortwave,
359     & SEAICE_saltFrac, SEAICE_salt0, SEAICEstressFactor,
360     & SEAICE_mcPheeTaper, SEAICE_mcPheePiston,
361     & SEAICE_frazilFrac, SEAICE_availHeatFrac,
362     & facOpenGrow, facOpenMelt,
363     & SEAICE_tempFrz0, SEAICE_dTempFrz_dS,
364     & OCEAN_drag, LSR_ERROR, DIFF1,
365 torge 1.3 & JFNKgamma_nonlin, JFNKres_t,
366     & JFNKgamma_lin_min, JFNKgamma_lin_max,
367 dimitri 1.1 & SEAICE_area_reg, SEAICE_hice_reg,
368     & SEAICE_area_floor, SEAICE_area_max,
369     & SEAICEdiffKhArea, SEAICEdiffKhHeff, SEAICEdiffKhSnow,
370 dimitri 1.2 & SEAICEdiffKhSalt, SEAICE_tauAreaObsRelax,
371 dimitri 1.1 & SEAICE_airTurnAngle, SEAICE_waterTurnAngle
372    
373     C-- COMMON /SEAICE_BOUND_RL/ Various bounding values
374     C MAX_HEFF :: maximum ice thickness (m)
375     C MIN_ATEMP :: minimum air temperature (deg C)
376     C MIN_LWDOWN :: minimum downward longwave (W/m^2)
377     C MIN_TICE :: minimum ice temperature (deg C)
378     C SEAICE_EPS :: small number used to reduce derivative singularities
379     C SEAICE_EPS_SQ :: small number square
380     C
381     _RL MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MIN_TICE
382     _RL SEAICE_EPS, SEAICE_EPS_SQ
383     COMMON /SEAICE_BOUND_RL/
384     & MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MIN_TICE,
385 torge 1.3 & SEAICE_EPS, SEAICE_EPS_SQ
386 heimbach 1.4
387     #ifdef SEAICE_ITD
388 dimitri 1.2 CToM<<<
389 heimbach 1.4 C Hlimit :: ice thickness category limits (m), array of size nITD+1
390     C Hlimit_c1,_c2,_c3 :: coefficients set in seaice_readparams.F to calculate Hlimit in seaice_init_fixed.F
391     _RL Hlimit(0:nITD)
392     _RL Hlimit_c1, Hlimit_c2, Hlimit_c3
393     COMMON /SEAICE_BOUND_ITD_RL/
394 dimitri 1.2 & Hlimit,
395     & Hlimit_c1,Hlimit_c2,Hlimit_c3
396     C>>>ToM
397 heimbach 1.4 #endif /* SEAICE_ITD */
398 dimitri 1.1
399     C-- Constants used by sea-ice model
400     _RL ZERO , ONE , TWO
401     PARAMETER ( ZERO = 0.0 _d 0, ONE = 1.0 _d 0, TWO = 2.0 _d 0 )
402     _RL QUART , HALF
403     PARAMETER ( QUART = 0.25 _d 0, HALF = 0.5 _d 0 )
404     _RL siEps
405     PARAMETER ( siEps = 1. _d -5 )
406     INTEGER MPSEUDOTIMESTEPS
407     PARAMETER (MPSEUDOTIMESTEPS=2)
408    
409     C-- Constants needed by McPhee formulas for turbulent ocean fluxes :
410     C Stanton number (dimensionless), typical friction velocity
411     C beneath sea ice (m/s), and tapering factor (dimensionless)
412     _RL STANTON_NUMBER, USTAR_BASE, MCPHEE_TAPER_FAC
413     PARAMETER ( MCPHEE_TAPER_FAC = 12.5 _d 0 , STANTON_NUMBER =
414     & 0.0056 _d 0, USTAR_BASE = 0.0125 _d 0 )
415    
416     C-- identifiers for advected properties
417     INTEGER GAD_HEFF,GAD_AREA,GAD_QICE1,GAD_QICE2,GAD_SNOW
418     INTEGER GAD_SALT,GAD_SITR
419     PARAMETER ( GAD_HEFF = 1,
420     & GAD_AREA = 2,
421     & GAD_SNOW = 3,
422     & GAD_SALT = 4,
423     & GAD_QICE1 = 5,
424     & GAD_QICE2 = 6,
425     & GAD_SITR = 7)
426    
427     CEH3 ;;; Local Variables: ***
428     CEH3 ;;; mode:fortran ***
429     CEH3 ;;; End: ***

  ViewVC Help
Powered by ViewVC 1.1.22