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