The following code produces an ICE. This was using the binary download version. Douglas Cox program hp_pi implicit none integer, parameter :: realdp = selected_real_kind(p=8,r=30) real (kind=realdp) :: b,one,four one = 1.0_realdp four = 4.0_realdp b = atan(one)*four write(unit=*,fmt="(a,f30.25)") 'pi = ', b end program hp_pi