--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/gcc/testsuite/gfortran.dg/pr71642.f90	Fri Oct 27 22:46:09 2017 +0900
@@ -0,0 +1,27 @@
+! PR debug/71642
+! { dg-do compile }
+! { dg-options "-g" }
+
+MODULE gauss_colloc
+  INTEGER, PARAMETER :: dp=8
+CONTAINS
+SUBROUTINE collocGauss(h,h_inv,grid,poly,alphai,posi,max_r2,&
+        periodic,gdim,local_bounds,local_shift,poly_shift,scale,lgrid,error)
+    REAL(dp), DIMENSION(0:, 0:, 0:), &
+      INTENT(inout)                          :: grid
+    INTEGER, INTENT(inout), OPTIONAL :: lgrid
+    CONTAINS
+    SUBROUTINE kloop6
+    IF (kJump/=1 .AND. (ikstart+kmax-kstart>=ndim(2)+l_shift(2) .OR.&
+        ikstart2+kmin-kstart2<=l_ub(2)-ndim(2))) THEN
+        DO
+            DO k=kstart2,kend2,-1
+                IF ( PRESENT ( lgrid ) ) THEN
+                  grid(ik,ij,ii) = grid(ik,ij,ii) + p_v*res_k
+                END IF
+            END DO
+        END DO
+    END IF
+    END SUBROUTINE
+END SUBROUTINE
+END MODULE gauss_colloc