# NOT RUN {
# Calculate normative temperature drop based on Minenergo-325 for pipe segment
pipeline <- list(
year = 1968,
laying = "channel",
d = 700,
l = 1000
)
operation_temperature <- c(130, 150) # [<U+00B0>C]
foo <- dropt(
temperature = operation_temperature,
flux = do.call(
m325nhl,
c(pipeline, temperature = list(operation_temperature))
)
)
foo
# [1] 1.366806 1.433840
# This is the same as using m325dropt:
bar <- m325dropt(temperature = operation_temperature,
year = 1968, laying = "channel", d = 700, len = 1000
)
bar
# [1] 1.366806 1.433840
# }
Run the code above in your browser using DataLab