add if statment to init_histogram
update init_histogram with else statement
add gnuplot_head
move floating point operation outside of loop
move vector_type outside of loop
move xyz and 0.1_sp back inside loops. move if statements outside of loops
add do concurrent