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