0.40
10
TEMPs
0.35
5
TEMP
0.30
0
TRHO II
TRHO DEL
0.25
5
0.20
10
0
10
20
30
40
50
60
70
80
90
Time (h)
Figure 17. Behavior of TRHO_DEL and TRHO_II. TRHO_II steadily in-
creases for the entire simulation period. TRHO_DEL steadily increases except
when TEMPs drops below zero from 22 to 30 hours and when TEMPs drops
below zero again from 87 to 90 hours. This figure shows that TRHO_DEL is
more sensitive than TRHO_II to changes in TEMPs, because TRHO_DEL sig-
nificantly increases twice when TEMPs drops below zero while TRHO_II
steadily increases despite the value of TEMPs.
Table 5. Conservation of mass within SWC.
Hours
SWC
MELT I
MELT II
REFREZ
SNOWE
SUBLM
TEMPs
0
46.8
0
0
0
0
0.1
0.4
1
46.7
0
0
0
0
0.1
2.3
2
46.6
0
0
0
0
0.1
1.3
3
46.5
0
0
0
0
0.1
0.6
4
46.4
0
0
0
0
0.1
0.4
5
46.3
0
0
0
0
0.1
0.4
6
46.2
0
0
0
0
0.1
0.2
7
46.1
0
0
0
0
0.1
0.1
8
46
0.04
0
0
0
0
0.2
9
45.96
0.21
0
0
0
0
1
10
45.75
0.31
0
0
0
0
1.5
11
45.44
0.38
0
0
0
0
1.8
12
45.06
0.33
0
0
0
0
1.6
13
44.73
0.21
0
0
0
0
1
14
44.52
0.15
0
0
0
0
0.7
15
44.37
0.02
0
0
0
0
0.1
16
44.35
0
0
0.04
0
0.1
0.2
17
44.3
0
0
0.04
0
0.1
0.2
18
44.24
0
0
0.06
0
0.1
0.3
19
44.2
0
0
0.06
0
0.1
0.3
20
44.16
0
0
0.08
0
0.1
0.4
demonstrates that conservation of mass occurs with SWC for every hour, because the value of SWC
minus MELT_I, minus MELT_II, minus SUBLM, plus REFREZ, plus SNOWE equals the value of
SWC for the next hour (eq 2). SUBLM was set to 0.1 to help demonstrate conservation of mass.
Table 6 shows that mass is conserved within LWC. For any hour, LWC plus MELT, minus
EXCESS_LWC, minus REFREZ, plus RAIN_IN_I equal LWC for the following hour. According to
Table 6 the value of LWC at 16 hours should be 5.15, but is 5.16 due to round off error.
Table 7 shows that mass is conserved in SDEP. SNOW_IN was set to ten from five to 15 hours to
better demonstrate conservation of mass. For any hour, SDEP minus COMPACTION, minus
SNOWMELT, plus SNOW_IN equals SDEP for the following hour.
25