/[MITgcm]/MITgcm/pkg/bulk_force/BULKF_ICE_CONSTANTS.h
ViewVC logotype

Contents of /MITgcm/pkg/bulk_force/BULKF_ICE_CONSTANTS.h

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

Revision 1.1 - (show annotations) (download)
Thu Nov 21 19:11:42 2002 UTC (21 years, 5 months ago) by cheisey
Branch: MAIN
CVS Tags: checkpoint47b_post, checkpoint51f_post, checkpoint48i_post, checkpoint51l_post, checkpoint47a_post, checkpoint51k_post, checkpoint48d_post, checkpoint50b_post, checkpoint47i_post, checkpoint51o_post, checkpoint48g_post, branchpoint-genmake2, checkpoint50c_pre, checkpoint50, checkpoint51j_post, branch-exfmods-tag, checkpoint51q_post, checkpoint47e_post, checkpoint50f_post, checkpoint50a_post, checkpoint48e_post, checkpoint47c_post, checkpoint50f_pre, checkpoint52a_pre, checkpoint48b_post, checkpoint47j_post, checkpoint47d_pre, checkpoint50d_pre, checkpoint47h_post, checkpoint51d_post, checkpoint51, checkpoint51r_post, checkpoint48c_pre, checkpoint52, checkpoint52b_pre, checkpoint48c_post, checkpoint50d_post, checkpoint51o_pre, checkpoint47f_post, checkpoint51t_post, checkpoint51b_pre, checkpoint52a_post, checkpoint51i_post, checkpoint50e_post, checkpoint47g_post, checkpoint50h_post, checkpoint50c_post, checkpoint51a_post, checkpoint51n_pre, checkpoint47d_post, checkpoint50e_pre, checkpoint50b_pre, checkpoint48d_pre, checkpoint50i_post, checkpoint51p_post, checkpoint51n_post, checkpoint51e_post, checkpoint51b_post, checkpoint48a_post, checkpoint51h_pre, checkpoint48f_post, checkpoint51i_pre, checkpoint51l_pre, checkpoint50g_post, checkpoint51u_post, checkpoint51c_post, checkpoint51g_post, checkpoint51m_post, ecco_c52_e35, checkpoint48, checkpoint49, checkpoint51f_pre, checkpoint51s_post, checkpoint48h_post
Branch point for: checkpoint51n_branch, branch-nonh, tg2-branch, branch-genmake2, branch-exfmods-curt
File MIME type: text/plain 
Two packages:  bulk_force (Bulk forcing)
and therm_seaice (thermodynamic_seaice) - adopted from LANL CICE.v2.0.2
Earlier integration from Stephaine Dutkiewicz
and Patrick Heimbach.

Two ifdef statements for compile time,
ALLOW_THERM_SEAICE and ALLOW_BULK_FORCE

Two switches in data.pkg to turn on at run-time:

cat data.pkg
# Packages
 &PACKAGES
 useBulkForce=.TRUE.,
 useThermSeaIce=.TRUE.,
 &

WARNING:  useSEAICE and useThermSEAICE are mutually exclusive.

The bulk package requires an additional parameter file
with two namelists, data.ice and data.blk.

c ADAPTED FROM:
c LANL CICE.v2.0.2
c-----------------------------------------------------------------------
c.. thermodynamics (vertical physics) based on M. Winton 3-layer model
c.. See Bitz, C. M. and W. H. Lipscomb, 1999:  "An energy-conserving
c..       thermodynamic sea ice model for climate study."  J. Geophys.
c..       Res., 104, 15669 - 15677.
c..     Winton, M., 1999:  "A reformulated three-layer sea ice model."
c..       Submitted to J. Atmos. Ocean. Technol.

c.. authors Elizabeth C. Hunke and William Lipscomb
c..         Fluid Dynamics Group, Los Alamos National Laboratory
c-----------------------------------------------------------------------
1 c----------------------------------------------------------------------------
2 c.. Common blocks for almost everything that the sea ice model passes around.
3 c----------------------------------------------------------------------------
4 c.. number layers of ice
5 integer nlyr ! maximum number of ice layers
6 _RL rnlyr ! real value of nlyr
7 parameter (nlyr = 2)
8 parameter (rnlyr = 2.)
9
10 c.. densities
11 _RL rhoa ! density of air (kg/m^3)
12 _RL rhos ! density of snow (kg/m^3)
13 _RL rhoi ! density of ice (kg/m^3)
14 _RL rhosw ! density of seawater (kg/m^3)
15 _RL rhoiw ! ocean-ice density difference (kg/m^3)
16
17 parameter (rhoa =1.3)
18 parameter (rhos = 330.)
19 parameter (rhoi = 900.)
20 parameter (rhosw = 1026.)
21 parameter (rhoiw = rhosw - rhoi)
22
23 c.. specific heats
24 _RL cpair ! specific heat of air (J/kg/K)
25 _RL cpice ! specific heat of fresh ice (J/kg/K)
26 _RL cpwater ! specific heat of water (J/kg/K)
27 _RL cpwv ! specific heat of water vapour (J/kg/KC)
28 _RL cpvir
29
30 parameter(cpair = 1004.)
31 parameter(cpice = 2106.)
32 parameter (cpwater = 4218.)
33 parameter (cpwv = 1.81e3)
34 parameter (cpvir = cpwv/cpair -1.)
35
36 c .. thermal conductivity. QQ check units
37 _RL kice ! thermal conductivity of pure ice (W/m/K)
38 _RL ksnow ! thermal conductivity of snow (W/m/K)
39
40 parameter (kice = 2.03) !QQ originally 2.03)
41 parameter (ksnow = 0.30) !QQ originally 0.30)
42
43 c .. heat transfer coefficient
44 _RL transcoef ! transfer coef between ice and water (unitless)
45
46 parameter (transcoef=0.006) !QQ originally 0.006
47
48 c .. latent heat
49 _RL Lvap ! latent heat of vaporizn at 0 C (J/kg)
50 _RL Lfresh ! latent heat of melting of pure ice (J/kg)
51
52 parameter(Lvap = 2.5e+6)
53 parameter (Lfresh = 3.34e5)
54
55 c.. Enthalpy
56 _RL qsnow ! snow enthalpy (J/kg)
57 parameter (qsnow = Lfresh)
58
59 c .. Albedo
60 _RL albsnodry ! albedo of dry snow (Tsfc < 0)
61 _RL albsnowet ! albedo of melting snow (Tsfc = 0)
62 _RL albicemax ! max albedo of bare ice
63 _RL albicemin ! minimum ice albedo (thin melting ice)
64 _RL halb ! melt rate parameter for albedo QQ units
65
66 parameter (albsnodry = 0.85)
67 parameter (albsnowet = 0.75)
68 parameter (albicemax = 0.65)
69 parameter (albicemin = 0.20)
70 parameter (halb = .5)
71
72 c.. Solar parameters
73 _RL i0 ! fraction of penetrating solar rad
74 _RL ksolar ! bulk solar abs coeff of sea ice (m-1)
75
76 parameter (i0 = 0.3) !QQQQQQ was 0.3
77 parameter (ksolar = 1.5) !QQQQQQ was1.5
78
79
80 c .. Salinity
81 _RL saltice ! salinity of ice (o/oo)
82 _RL S_winton ! winton salinity of ice (o/oo)
83 _RL mu_Tf ! Tf:brine salinity ratio (C/ppt)
84
85 parameter (saltice = 4.)
86 parameter (S_winton = 1.)
87 parameter (mu_Tf = 0.054)
88
89 c .. melting
90 _RL Tf0kel ! Freezing temp of fresh ice in Kelvin = 273.15
91 _RL Tmlt1 ! melting temp; depends on S (C)
92
93 parameter (Tf0kel = 273.15)
94 parameter (Tmlt1=-mu_Tf*S_winton)
95
96 c .. wind drag
97 c cdrag_[n] - n = 1,2,3 coefficients used to evaluate
98 c drag coefficient
99 _RL cdrag_1, cdrag_2, cdrag_3
100 parameter ( cdrag_1 = 0.0027000,
101 & cdrag_2 = 0.0001420,
102 & cdrag_3 = 0.0000764)
103
104 c .. constants
105 _RL stefan ! Stefan-Boltzmann constant (W/m^2 K^4)
106 _RL xkar ! Von Karman constant QQ units?
107 _RL Rvap ! gas constant for H2O vapor (J/kg/K)
108
109 parameter (stefan = 5.67e-8)
110 parameter (xkar = 0.4 )
111 parameter(Rvap = 461.)
112
113 c.. Miscellaneous
114 _RL p0 ! surface pressure (mb)
115
116 parameter(p0 = 1013.)
117
118 c.. Combinations used for efficiency
119 _RL rhoi_Lfresh ! rhoi*Lfresh
120 _RL Qcoef ! another constant for latent heat flux
121 _RL lvrrv ! (Lvap + Lfresh) / Rvap
122
123 parameter (rhoi_Lfresh = rhoi*Lfresh)
124 parameter (Qcoef = 6.11*0.622/p0)
125 parameter (lvrrv = (Lvap + Lfresh) / Rvap)
126
127 c .. Min/Max
128 _RL himin ! minimum thickness for ice (m)
129 _RL Terrmax ! temperature convergence precision (C)
130
131 parameter (himin = 0.01) !was 0.01?
132 parameter (Terrmax = 5.0e-1) !was 5.0e-3????
133
134 c .. for bulk formula
135 _RL humid_fac ! const. for the evaluation of the virtual temp.
136 _RL saltsat ! reduction of sat. vapor pressure over salt water
137 _RL gamma_blk ! adiabatic lapse rate
138
139
140 parameter (humid_fac = 0.606)
141 parameter (saltsat = 0.980)
142 parameter (gamma_blk = 0.010)
143
144
145
146 c .. for Vince bulk formula QQ check units
147 _RL Lvap_ice
148 _RL Rgas
149 _RL Sha
150
151 parameter (Lvap_ice = 2.83e6) ! latent heat from sublimation
152 parameter (Rgas = 287.e0) ! gas constant for dry air
153 parameter (Sha = Rgas/.286) !
154
155 c .. emissivities QQQ find
156 _RL atm_emissivity
157 _RL ocean_emissivity
158 _RL snow_emissivity
159 _RL ice_emissivity
160
161 cQQQQQQ find real values
162 parameter(atm_emissivity=.90d0) !QQ
163 parameter(ocean_emissivity=.985d0) !QQ
164 parameter(snow_emissivity=.98d0) !QQ
165 parameter(ice_emissivity=.98d0) !QQ
  ViewVC Help
Powered by ViewVC 1.1.22