anchor.diff:      if(nSplit > 5) call chr_cast(lnSplit(6), iParams(2), cErr)
   
     case("GAUSS")
  -    if(nSplit /= 6 .and. nSplit /= 7 .and. nSplit /= 9) then
  -      write(lerr,"(a,i0)") "DIST> ERROR FILL format GAUSS expected 3, 4 or 6 arguments, got ",nSplit-3
  -      write(lerr,"(a)")    "DIST>       FILL column GAUSS scale sigma mu [cut] [first last]"
  +    if(nSplit /= 5 .and. nSplit /= 6 .and. nSplit /= 8) then
  +      write(lerr,"(a,i0)") "DIST> ERROR FILL format GAUSS expected 2, 3 or 5 arguments, got ",nSplit-3
  +      write(lerr,"(a)")    "DIST>       FILL column GAUSS sigma mu [cut] [first last]"
         iErr = .true.
         return
       end if
  -    
  -    allocate(fParams(4))
  +
  +    allocate(fParams(3))
       allocate(iParams(2))
  -    
  +
       fillMethod = dist_fillGAUSS
       fParams(1) = one
  -    fParams(2) = one
  +    fParams(2) = zero
       fParams(3) = zero
  -    fParams(4) = zero
  -    iParams(1) =  1
  -    iParams(2) = -1
  -
  -    if(nSplit > 3) call chr_cast(lnSplit(4), fParams(1), cErr)
  -    if(nSplit > 4) call chr_cast(lnSplit(5), fParams(2), cErr)
  -    if(nSplit > 5) call chr_cast(lnSplit(6), fParams(3), cErr)
  -    if(nSplit > 6) call chr_cast(lnSplit(7), fParams(4), cErr)
  -    if(nSplit > 7) call chr_cast(lnSplit(8), iParams(1), cErr)
  -    if(nSplit > 8) call chr_cast(lnSplit(9), iParams(2), cErr)
  -
  -  case("UNIFORM")
  -    if(nSplit /= 6 .and. nSplit /= 8) then
  -      write(lerr,"(a,i0)") "DIST> ERROR FILL format UNIFORM expected 3 or 5 arguments, got ",nSplit-3
  -      write(lerr,"(a)")    "DIST>       FILL column UNIFORM scale lower upper [first last]"
  -      iErr = .true.
  -      return
  -    end if
  -
  -    allocate(fParams(4))
  -    allocate(iParams(2))
  -
  -    fillMethod = dist_fillUNIFORM
  -    fParams(1) = one
  -    fParams(2) = zero
  -    fParams(3) = one
       iParams(1) =  1
       iParams(2) = -1
   
  
positive.code: end subroutine dist_saveParticle
  
  ! ================================================================================================ !
  !  Post-Processing of Particle Arrays
  !  V.K. Berglyd Olsen, BE-ABP-HSS
  !  Created: 2019-07-05
  
negative.code:       iErr = .true.
        return
      end if
  
      allocate(fParams(1))
      allocate(iParams(2))
  
      fillMethod = dist_fillFLOAT
      fParams(1) = zero
      iParams(1) =  1
  