subroutine set_itv_123(ratio) use constraint_variables, only: foh_stress implicit none real(kind(1.d0)) :: ratio foh_stress = ratio end subroutine set_itv_123