anchor.diff:    nzz(:)   = 0
     nucm(:)  = zero
   
  -! initialise particle counter
     jj = 0
  -
  -  open( unit=dist_readUnit, file=dist_readFile )
  -
  -! cycle on lines in file:
  +  ln = 0
  +  open(unit=dist_readUnit, file=dist_readFile)
  +
   1981 continue
  -  read(dist_readUnit,"(a)",end=1983,err=1982) tmp_line
  -  if( tmp_line(1:1).eq.comment_char ) goto 1981
  +  read(dist_readUnit,"(a)",end=1983,err=1982) inLine
  +  ln = ln+1
  +
  +  if(inLine(1:1) == "*") goto 1981
  +  if(inLine(1:1) == "#") goto 1981
  +  if(inLine(1:1) == "!") goto 1981
     jj = jj+1
   
  -  if( jj.gt.napx ) then
  -    write(lout,"(a,i0,a)") "DIST> Stopping reading file as ",napx," particles have been as requested in fort.3"
  +  if(jj > napx) then
  +    write(lout,"(a,i0,a)") "DIST> Stopping reading file as ",napx," particles have been read, as requested in fort.3"
       jj = napx
       goto 1983
     end if
   
  -  read(tmp_line, *, err=1982) id,gen,weight,xv(1,jj),xv(2,jj),z,yv(1,jj),yv(2,jj),zp,naa(jj),nzz(jj),nucm(jj),ejfv(jj),dt(jj)
  +  read(inLine, *, err=1982) id,gen,weight,xv(1,jj),xv(2,jj),z,yv(1,jj),yv(2,jj),zp,naa(jj),nzz(jj),nucm(jj),ejfv(jj),dt(jj)
   
     xv(1,jj)  = xv(1,jj)*c1e3
     xv(2,jj)  = xv(2,jj)*c1e3
  
positive.code: 
  ! error while parsing file:
  1982 continue
    write(lout,"(a)") "DIST> ERROR Reading particles from line: '"//trim(tmp_line)//"'"
    call prror(-1)
    return
  
  
negative.code:   write(lout,"(a,i0,a)") "DIST> Read ",jj," particles from file '"//trim(dist_readFile)//"'"
  
    if(jj < napx) then
      write(lout,"(a,i0)") "DIST> WARNING Read a number of particles LOWER than requested: ",napx
      napx = jj
    end if
  
  