subroutine set_itv_10(ratio) use physics_variables, only: hfact implicit none real(kind(1.d0)) :: ratio hfact = ratio end subroutine set_itv_10