diff --git a/kf.for b/kf.for
index 971e2a834ae9642c039ca08a54930e11c69e7b71..e7760c066d546c53462ba21d6ee1b6844f5bfcbd 100644
--- a/kf.for
+++ b/kf.for
@@ -58,7 +58,7 @@ C LOCALS
      2 RG(nx,ny),HPV(nx,ny)
 
       LIKE(d(1)+1:nobs) = 0.D0
-      DO 1000 imain = d(1)+1,nobs
+      DO imain = d(1)+1,nobs
 	 iny = IYK(imain,ny+1)
 
 C ------------------------------------
@@ -171,7 +171,7 @@ c       CALL F03ABF(V(1:iny,1:iny),iny,iny,DETV,COM(1:iny,1),IFAIL)
 	  XT(imain,1:nx)      = X1(1:nx)
 	  PT(imain,1:nx,1:nx) = P1(1:nx,1:nx)
 
-1000	 ENDIF
-
-	RETURN
-	END
+      END IF
+      END DO
+      RETURN
+      END SUBROUTINE