# NOT RUN {
with(m325nhldata, {
## Linear extrapolation adopted in Minenergo's Order 325 using last two points:
temperature <- seq(0, 270, 10) # [<U+00B0>C]
flux <- m325nhl(1980, "underground", TRUE, 0, 73, temperature) # [kcal/m/h]
plot(temperature, flux, type = "b")
## Consider heat losses of fittings:
stopifnot(
## when beta becomes 1.15
all(
round(
m325nhl(1980, "underground", d = 73, temperature = 65,
beta = c(FALSE, TRUE)),
3
) == c(65.500, 75.325)
),
## when beta becomes 1.2
all(
round(
m325nhl(2000, "channel", d = 73, temperature = 65,
beta = c(FALSE, TRUE)),
3
) == c(17.533, 21.040)
)
)
})
# }
Run the code above in your browser using DataLab