1 |
cheisey |
1.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 |