fixing unitialized values in refine_fn and get_element_order
energy computation and time stepping