using Glob
using DelimitedFiles
const lataus = -0.05
const purku = 0.25
const eps = 1.0
function error(file)
data_matrix = readdlm(file, ' ', Float64, comments=true)
Energy_in = data_matrix[end, 8]
E_solid = data_matrix[end, 6]
E_oil = data_matrix[end, 7]
return 100*(Energy_in - E_oil - E_solid) / Energy_in
end
function eurot(file)
data_matrix = readdlm(file, ' ', Float64, comments=true)
# P_in = data_matrix[end, 4]
dt = (data_matrix[2, 4] - data_matrix[1, 4])/3600.0
eurot = 0
for P_in in data_matrix[:,4]
if P_in > eps
eurot += P_in*lataus*dt
elseif P_in < eps
eurot -= P_in*purku*dt
end
end
return eurot
end
function main()
println("#############################################################3")
directory = length(ARGS) >= 1 ? ARGS[1] : ""
for file in glob("*_history.dat", directory)
println("\n$file")
println("virhe energiataseessa $(error(file)) %")
println("tuotto $(eurot(file))")
end
end
if abspath(PROGRAM_FILE) == @__FILE__
main()
end