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