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