Pascal/Turbo Pascal - Matrices

Hola que tal?
Estoy haciendo una práctica, en la cual , te dan una matriz, y tienes k hacer el perfil por columnas es decir, es el primer elemento, de la fila k no es nulo:
0 0 2 0 0 sería: la posición 3
1 0 0 5 0 posición 1
0 5 0 0 6 2
5 0 0 0 0 1
(3,1,2,1)

Yo lo tengo hecho, xo me da un fallo y en la tercera fila cuenta lo k kiere.Si alguien me puede ayudar se lo agradezco. El código es este:

PROGRAM practica2(FDATOS,INPUT,OUTPUT);
CONST
MAX= 6;
TYPE
TMATRIZ = ARRAY[1..MAX,1..MAX] OF INTEGER;
VAR
FDATOS: TEXT;
MATRIZ: TMATRIZ;
I,J:INTEGER;

PROCEDURE PERFIL_COLUMNAS(VAR F:TEXT;MAT:TMATRIZ;LONG:INTEGER);
VAR
I,J,K:INTEGER;
SIGUIENTE:BOOLEAN;
BEGIN
I:=1;
WHILE (I<=LONG) DO
BEGIN
SIGUIENTE:=FALSE;
J:=1;
WHILE (SIGUIENTE=FALSE) AND (J<=LONG) DO
BEGIN
READ(F,MAT[I,J]);
IF (MAT[I,J] = 0) THEN
SIGUIENTE:= FALSE
ELSE
BEGIN
SIGUIENTE:=TRUE;
WRITE(J);
END;
J:=J+1
END;
I:=I+1;
end;
end;

BEGIN
ASSIGN(FDATOS,'F:\datosp2.txt');
RESET(FDATOS);
WRITELN;
PERFIL_COLUMNAS(FDATOS,MATRIZ,MAX);
END.