parameter(m=25,m1=25) INTEGER IFL3(m,m1),IBD(m,m1) INTEGER J25,I25 OPEN(1,file='C:\Documents and Settings\IFL.txt') DO 10 I=1,25 DO 20 J=1,25 READ(1,*)IFL(I,J) 20 CONTINUE 10 CONTINUE CLOSE(1) C J25=25 I25=25 DO 222 J = 1,J25 DO 222 I = 1,I25 IBD (I,J) = 0 222 CONTINUE C DO 230 J = 1,J25 DO 231 I = 1,I25-1 IF(IFL(I+1,J)-IFL(I,J).EQ.-1) IBD(I,J) = 1 231 CONTINUE DO 232 I = I25,2,-1 IF(IFL(I,J)-IFL(I-1,J).EQ.1) IBD(I,J) = 1
232 CONTINUE 230 CONTINUE DO 250 I = 1,I25 DO 251 J=1,J25-1 IF(IFL(I,J)-IFL(I,J+1).EQ.-1) IBD(I,J) =1 251 CONTINUE DO 252 J = J25,2,-1 IF(IFL(I,J-1)-IFL(I,J).EQ.1) IBD(I,J)=1 252 CONTINUE 250 CONTINUE OPEN(1,file='IBD.dat') DO 253 I=1,I25 DO 254 J=1,J25 WRITE(1,*) IBD(I,J) 254 CONTINUE 253 CONTINUE CLOSE(1) STOP END