/[MITgcm]/MITgcm/pkg/seaice/SEAICE_PARAMS.h
ViewVC logotype

Contents of /MITgcm/pkg/seaice/SEAICE_PARAMS.h

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

Revision 1.11 - (show annotations) (download)
Thu Oct 9 04:19:20 2003 UTC (20 years, 7 months ago) by edhill
Branch: MAIN
CVS Tags: checkpoint51k_post, checkpoint52l_pre, hrcube4, hrcube5, checkpoint52d_pre, checkpoint52j_pre, checkpoint51o_pre, checkpoint51l_post, checkpoint52l_post, checkpoint52k_post, checkpoint52, checkpoint52f_post, checkpoint51t_post, checkpoint51n_post, checkpoint52i_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post, checkpoint52e_pre, checkpoint52e_post, checkpoint51n_pre, checkpoint52b_pre, checkpoint51l_pre, checkpoint52m_post, checkpoint51q_post, checkpoint52b_post, checkpoint52c_post, checkpoint52f_pre, checkpoint51r_post, checkpoint51i_post, checkpoint52d_post, checkpoint52a_pre, checkpoint52i_post, checkpoint52h_pre, checkpoint52j_post, branch-netcdf, checkpoint51o_post, checkpoint52a_post, ecco_c52_e35, checkpoint51m_post, checkpoint51p_post, checkpoint51u_post
Branch point for: branch-nonh, tg2-branch, netcdf-sm0, checkpoint51n_branch
Changes since 1.10: +2 -1 lines
File MIME type: text/plain 
 o first check-in for the "branch-genmake2" merge
 o verification suite as run on shelley (gcc 3.2.2):

Wed Oct  8 23:42:29 EDT 2003
                T           S           U           V
G D M    c        m  s        m  s        m  s        m  s
E p a R  g  m  m  e  .  m  m  e  .  m  m  e  .  m  m  e  .
N n k u  2  i  a  a  d  i  a  a  d  i  a  a  d  i  a  a  d
2 d e n  d  n  x  n  .  n  x  n  .  n  x  n  .  n  x  n  .

OPTFILE=NONE

Y Y Y Y 13 16 16 16  0 16 16 16 16 16 16 16 16 13 12  0  0 pass  adjustment.128x64x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16  0  0 16 16  0  0 pass  adjustment.cs-32x32x1
Y Y Y Y 16 16 16 16  0 16 16 16 16 16 16 22  0 16 16 22  0 pass  adjust_nlfs.cs-32x32x1
Y Y Y Y -- 13 13 16 16 13 13 13 13 16 16 16 16 16 16 16 16 N/O   advect_cs
Y Y Y Y -- 22 16 16 16 16 16 16 13 16 16 16 16 16 16 16 16 N/O   advect_xy
Y Y Y Y -- 13 16 13 16 16 16 16 16 16 16 22 16 16 16 16 16 N/O   advect_xz
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  aim.5l_cs
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 16 16 16 16 13 16 pass  aim.5l_Equatorial_Channel
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 13 16 16 13 13 16 pass  aim.5l_LatLon
Y Y Y Y 13 16 16 16 16 16 16 16 16 16 13 12 13 13 16 13 16 pass  exp0
Y Y Y Y 14 16 16 16 16 16 16 16 22 16 16 16 13 16 16 22 16 pass  exp1
Y Y Y Y 13 13 16 13 16 16 16 16 16 13 13 16 16 13 13 13 13 pass  exp2
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  exp4
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 22 16 16 16 22 16 pass  exp5
Y Y Y Y 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 16 pass  front_relax
Y Y Y Y 14 16 16 13 13 16 16 13 13 16 13 13 16 12 13 13 16 pass  global_ocean.90x40x15
Y Y Y Y 10 16 16 13 13 16 13 16 16 13 13 13 13 16 16 13 16 FAIL  global_ocean.cs32x15
Y Y Y Y  6 11 12 13 13 12 13 16 13  9  9  9  9 10  9  9 11 FAIL  global_ocean_pressure
Y Y Y Y 14 16 16 13 16 16 16 13 13 13 13 13 16 12 16 13 16 pass  global_with_exf
Y Y Y Y 14 16 16 16 16 16 16 16 16 11 13 22 13 16 16  9 16 pass  hs94.128x64x5
Y Y Y Y 13 16 16 16 16 16 16 16 16 11 16 16 16 13 16 22 13 pass  hs94.1x64x5
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 16 13 13 16 16 22 13 pass  hs94.cs-32x32x5
Y Y Y Y 10 10 16 13 13 16 16 16 22 16 13 13 13 13 13 22 13 FAIL  ideal_2D_oce
Y Y Y Y  8 16 16 16 16 16 16 16 16 13 13  8 16 16 16 16 16 FAIL  internal_wave
Y Y Y Y 14 16 16 16 16 16 16 16 16 13 13 22 13 13 13 22 16 pass  inverted_barometer
Y Y Y Y 12 16 16 16 16 16 16 16 16 16 13 12 13 13 13 13 13 FAIL  lab_sea
Y Y Y Y 11 16 16 16 16 16 16 16 13 13 13 12 13 16 13 12 13 FAIL  natl_box
Y Y Y Y 16 16 16 16 16 16 16 16 22 16 16 16 16 16 16 16 16 pass  plume_on_slope
Y Y Y Y 13 16 16 16 16 13 16 16 16 16 16 16 16 13 16 16 16 pass  solid-body.cs-32x32x1
1 C $Header: /u/u3/gcmpack/MITgcm/pkg/seaice/SEAICE_PARAMS.h,v 1.10.2.1 2003/10/02 18:18:33 adcroft Exp $
2 C $Name: $
3
4 #ifdef ALLOW_SEAICE
5
6 C /==========================================================\
7 C | SEAICE_PARAMS.h |
8 C | o Basic parameter header for sea ice model. |
9 C \==========================================================/
10
11 C-- COMMON /SEAICE_PARM_L/ Logical parameters of sea ice model.
12 C
13 C SEAICEwriteState - If true, write sea ice state to file;
14 C default is false.
15 C
16 C SEAICEuseDYNAMICS - If false, do not use dynamics;
17 C default is to use dynamics.
18 C
19 C SEAICEuseADI - If true, use ADI rather than LSR solver;
20 C default is LSR solver.
21 C
22 LOGICAL SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseADI
23 COMMON /SEAICE_PARM_L/
24 & SEAICEwriteState, SEAICEuseDYNAMICS, SEAICEuseADI
25
26 C-- COMMON /SEAICE_PARM_I/ Integer valued parameters of sea ice model.
27 C LAD - time stepping used for sea-ice advection:
28 C 1 = LEAPFROG, 2 = BACKWARD EULER.
29 C IMAX_TICE - number of iterations for ice heat budget 10
30 C NPSEUDO - number of pseudo-timesteps used in dynsolver
31 C DO PSEUDO-TIMESTEPS TO OBTAIN AN ACCURATE VISCOUS-PLASTIC SOLUTION
32 C 10 PSEUDO-TIMESTEPS OR MORE ARE SUGGESTED FOR HIGH-RESOLUTION (~10KM)
33 C 1 PSEUDO-TIMESTEP CAN BE USED FOR LOW-RESOLUTION GLOBAL MODELING
34 C
35 INTEGER LAD, IMAX_TICE, NPSEUDO
36 COMMON /SEAICE_PARM_I/ LAD, IMAX_TICE, NPSEUDO
37
38 C-- COMMON /SEAICE_PARM_C/ Character valued sea ice model parameters.
39 C uwindFile - File containing uwind
40 C vwindFile - File containing vwind
41 C atempFile - File containing atemp
42 C aqhFile - File containing aqh
43 C lwdownFile - File containing lwdown
44 C swdownFile - File containing swdown
45 C precipFile - File containing precip
46 C evapFile - File containing evap
47 C runoffFile - File containing runoffF
48 C HeffFile - File containing initial sea-ice thickness
49 C !!! NOTE !!! Initial sea-ice thickness can also be set using
50 C SEAICE_initialHEFF below. But a constant initial condition
51 C can mean large artificial fluxes of heat and freshwater in
52 C the surface layer during the first model time step.
53 C
54 CHARACTER*(MAX_LEN_FNAM) uwindFile
55 CHARACTER*(MAX_LEN_FNAM) vwindFile
56 CHARACTER*(MAX_LEN_FNAM) atempFile
57 CHARACTER*(MAX_LEN_FNAM) aqhFile
58 CHARACTER*(MAX_LEN_FNAM) lwdownFile
59 CHARACTER*(MAX_LEN_FNAM) swdownFile
60 CHARACTER*(MAX_LEN_FNAM) precipFile
61 CHARACTER*(MAX_LEN_FNAM) evapFile
62 CHARACTER*(MAX_LEN_FNAM) runoffFile
63 CHARACTER*(MAX_LEN_FNAM) HeffFile
64 COMMON /SEAICE_PARM_C/ uwindFile, vwindFile, atempFile, aqhFile,
65 & lwdownFile, swdownFile, precipFile, evapFile, runoffFile,
66 & HeffFile
67
68 C-- COMMON /SEAICE_PARM_RL/ Real valued parameters of sea ice model.
69 C SEAICE_dumpFreq - SEAICE dump frequency. (s)
70 C SEAICE_taveFreq - SEAICE time-averaging frequency. (s)
71 C SEAICE_initialHEFF - initial sea-ice thickness (m)
72 C SEAICE_drag - air-ice drag coefficient
73 C OCEAN_drag - air-ocean drag coefficient
74 C SEAICE_waterDrag - water-ice drag coefficient * water density
75 C SEAICE_dryIceAlb - winter albedo
76 C SEAICE_wetIceAlb - summer albedo
77 C SEAICE_drySnowAlb - dry snow albedo
78 C SEAICE_wetSnowAlb - wet snow albedo
79 C SEAICE_waterAlbedo - water albedo
80 C SEAICE_strength - sea-ice strength
81 C SEAICE_sensHeat - sensible heat transfer coefficient
82 C SEAICE_latentWater - latent heat transfer coefficient for water
83 C SEAICE_latentIce - latent heat transfer coefficient for ice
84 C SEAICE_iceConduct - sea-ice conductivity
85 C SEAICE_snowConduct - snow conductivity
86 C SEAICE_emissivity - Stefan-Boltzman constant * emissivity
87 C SEAICE_snowThick - cutoff snow thickness
88 C SEAICE_shortwave - penetration shortwave radiation factor
89 C LSR_ERROR - sets accuracy of LSR solver
90 C DIFF1 - parameter used in advect.F
91 C A22 - parameter used in growth.F
92 C HO - demarcation thickness between thin and
93 C thick ice: HO is a key ice-growth parameter
94 C WindForcingStart - Time of first wind forcing record (s)
95 C WindForcingEnd - Time of last wind forcing record (s)
96 C WindForcingPeriod - Period between wind forcing records (s)
97 C FluxForcingStart - Time of first flux forcing record (s)
98 C FluxForcingEnd - Time of last flux forcing record (s)
99 C FluxForcingPeriod - Period between flux forcing records (s)
100 C SSTForcingStart - Time of first SST forcing record (s)
101 C SSTForcingEnd - Time of last SST forcing record (s)
102 C SSTForcingPeriod - Period between SST forcing records (s)
103 C SSSForcingStart - Time of first SSS forcing record (s)
104 C SSSForcingEnd - Time of last SSS forcing record (s)
105 C SSSForcingPeriod - Period between SSS forcing records (s)
106 C StartingYear - Starting year of integration
107 C EndingYear - Ending year of integration
108 C
109 _RL SEAICE_dumpFreq, SEAICE_taveFreq, SEAICE_initialHEFF
110 _RL SEAICE_drag, SEAICE_waterDrag, SEAICE_dryIceAlb
111 _RL SEAICE_wetIceAlb, SEAICE_drySnowAlb, SEAICE_wetSnowAlb
112 _RL SEAICE_waterAlbedo, SEAICE_strength
113 _RL SEAICE_sensHeat, SEAICE_latentWater, SEAICE_latentIce
114 _RL SEAICE_iceConduct, SEAICE_snowConduct, SEAICE_emissivity
115 _RL SEAICE_snowThick, SEAICE_shortwave, OCEAN_drag
116 _RL LSR_ERROR,DIFF1,A22,HO
117 _RL WindForcingStart, WindForcingEnd, WindForcingPeriod
118 _RL FluxForcingStart, FluxForcingEnd, FluxForcingPeriod
119 _RL SSTForcingStart, SSTForcingEnd, SSTForcingPeriod
120 _RL SSSForcingStart, SSSForcingEnd, SSSForcingPeriod
121 _RL StartingYear, EndingYear
122 COMMON /SEAICE_PARM_RL/
123 & SEAICE_dumpFreq, SEAICE_taveFreq, SEAICE_initialHEFF,
124 & SEAICE_drag, SEAICE_waterDrag, SEAICE_dryIceAlb,
125 & SEAICE_wetIceAlb, SEAICE_drySnowAlb, SEAICE_wetSnowAlb,
126 & SEAICE_waterAlbedo, SEAICE_strength,
127 & SEAICE_sensHeat, SEAICE_latentWater, SEAICE_latentIce,
128 & SEAICE_iceConduct, SEAICE_snowConduct, SEAICE_emissivity,
129 & SEAICE_snowThick, SEAICE_shortwave, OCEAN_drag,
130 & LSR_ERROR, DIFF1, A22, HO,
131 & WindForcingStart, WindForcingEnd, WindForcingPeriod,
132 & FluxForcingStart, FluxForcingEnd, FluxForcingPeriod,
133 & SSTForcingStart, SSTForcingEnd, SSTForcingPeriod,
134 & SSSForcingStart, SSSForcingEnd, SSSForcingPeriod,
135 & StartingYear, EndingYear
136
137 C-- COMMON /SEAICE_BOUND_RL/ Various bounding values
138 C MAX_HEFF - maximum ice thickness (m)
139 C MIN_ATEMP - minimum air temperature (deg C)
140 C MIN_LWDOWN - minimum downward longwave (W/m^2)
141 C MAX_TICE - maximum ice temperature (deg C)
142 C MIN_TICE - minimum ice temperature (deg C)
143 C SEAICE_EPS, SEAICE_EPS_SQ - used to reduce derivative singularities
144 C
145 _RL MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE
146 _RL SEAICE_EPS, SEAICE_EPS_SQ
147 COMMON /SEAICE_BOUND_RL/
148 & MAX_HEFF, MIN_ATEMP, MIN_LWDOWN, MAX_TICE, MIN_TICE,
149 & SEAICE_EPS, SEAICE_EPS_SQ
150
151 C-- Seaice time step
152 COMMON /SEAICE_R/ SEAICE_DT
153 _RL SEAICE_DT
154
155 C-- Constants used by sea-ice model
156 _RL ZERO , ONE , TWO
157 parameter ( ZERO = 0.0 _d 0, ONE = 1.0 _d 0, TWO = 2.0 _d 0 )
158 _RL QUART , HALF
159 parameter ( QUART = 0.25 _d 0, HALF = 0.5 _d 0 )
160
161 #endif /* ALLOW_SEAICE */
  ViewVC Help
Powered by ViewVC 1.1.22