endif
sigma_odl(k,ii)=1.0d0/(sigma_odl(k,ii)*sigma_odl(k,ii))
else
- ii=ii+1
- l_homo(k,ii)=.false.
+c ii=ii+1
+c l_homo(k,ii)=.false.
endif
enddo
enddo
read(ientin,*) (inclust(k,i),k=1,ninclust(i))
read(ientin,*) (iresclust(k,i),k=1,nresclust(i))
enddo
+ close(ientin)
c
c Loop over clusters
c