1! { dg-do compile } 2 subroutine ice(npoint, nterm, x, g) 3 implicit none 4 integer norder 5 parameter (norder=10) 6 integer j 7 integer k 8 integer ii 9 integer nterm 10 integer npoint 11 real b(norder) 12 real c(norder) 13 real d(norder) 14 real x(npoint) 15 real g(npoint) 16 real gg 17 real prev 18 real prev2 19 20 j = 1 21 100 continue 22 j = j+1 23 if (nterm == j) then 24 do ii=1,npoint 25 k = nterm 26 gg= d(k) 27 prev= 0.0 28 do k=k-1,1,-1 29 prev2= prev 30 prev= gg 31 gg = d(k)+(x(ii)-b(k))*prev-c(k+1)*prev2 32 enddo 33 g(ii) = gg 34 enddo 35 endif 36 go to 100 37 end 38