gfortran.dg/gamma_(1|2).f90 fail on Darwin (and probably on any system with real(16) instead of real(10)) because selected_real_kind(p=15,r=900) returns -2: maximum range is 291. Using selected_real_kind (precision (0.0_8) + 1) solves the problem. Dominique