subroutine set_itv_161(ratio) use constraint_variables, only: fbeta_min real(kind(1.d0)) :: ratio fbeta_min = ratio end subroutine set_itv_161