anchor.diff:    if(spErr) goto 20
     if(nSplit == 0) goto 10
   
  -  if(isFirst) then
  -    if(nSplit < dist_nColumns) then
  -      dist_readCols = nSplit
  -      write(lout,"(2(a,i0))") "DIST> WARNING Distribution file has ",nSplit," columns, but format has ",dist_nColumns
  -    end if
  -    isFirst = .false.
  -  end if
  -  if(nSplit < dist_readCols) then
  +  if(nSplit /= dist_nColumns) then
       write(lerr,"(3(a,i0),a)") "DIST> ERROR Number of columns in file on line ",lineNo," is ",nSplit,&
  -      " but expected ",dist_readCols," columns"
  +      " but FORMAT defines ",dist_nColumns," columns"
       call prror
     end if
  -  do i=1,dist_readCols
  +  do i=1,nSplit
       call dist_saveParticle(j, i, lnSplit(i), cErr)
     end do
     if(cErr) goto 20
  +  dist_numPart = dist_numPart + 1
   
     goto 10
   
  
positive.code:   call f_close(fUnit)
    write(lout,"(a,i0,a)") "DIST> Read ",j," particles from file '"//trim(dist_distFile)//"'"
  
    if(j < napx) then
      write(lout,"(a,i0)") "DIST> WARNING Read a number of particles LOWER than requested: ",napx
      napx = j
    end if
  
  end subroutine dist_readDist
  
  ! ================================================================================================ !
  
negative.code: ! ================================================================================================ !
  subroutine dist_generateDist
  
    use crcoall
    use mod_common
    use mod_particles
  
    if(dist_nParticle > 0) then
      call dist_parseParticles
    end if
  
    if(dist_distFile /= " ") then
      call dist_readDist
    end if
  
    if(dist_numPart /= napx) then
      write(lerr,"(2(a,i0),a)") "DIST> ERROR Number of particles read or generated is ",dist_numPart,&
        " but the simulation setup requests ",napx," particles"
      call prror
    end if
  
    if(dist_numPart > 0) then
      call dist_postprParticles
      call dist_finaliseDist
      call part_applyClosedOrbit
  