Deallocate arrays
subroutine dealloc_arrays
implicit none
if(allocated(shat_arr)) deallocate(shat_arr)
if(allocated(beta_arr)) deallocate(beta_arr)
if(allocated(dbdrho_arr)) deallocate(dbdrho_arr)
if(allocated(stability)) deallocate(stability)
end subroutine dealloc_arrays