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