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