anchor.diff:      63.7e-9_fPrec, 166.0e-9_fPrec, 322.0e-9_fPrec, 727.0e-9_fPrec, 823.0e-9_fPrec, 78.0e-9_fPrec, 78.0e-9_fPrec, &
       87.1e-9_fPrec, 152.9e-9_fPrec, 424.0e-9_fPrec, 320.8e-9_fPrec, 682.2e-9_fPrec, zero, c1e10 ]
   
  -
  -  ! IN "+CD DBTRTHIN" and "+CD DBDATEN"
   !  real(kind=fPrec) remitx_dist,remity_dist,
   
     integer, private :: k
  
positive.code:   cdist_betaY    = mybetay
    cdist_emitX    = myemitx0_dist
    cdist_emitY    = myemity0_dist
    cdist_ampX     = mynex
    cdist_ampY     = myney
    cdist_smearX   = mdex
    cdist_smearY   = mdey
    cdist_spreadE  = enerror
    cdist_bunchLen = bunchlength
    cdist_fileName = filename_dis
    call cdist_makeDist(do_thisdis)
    ! if(radial) then
    !   call makedis_radial(myalphax, myalphay, mybetax, mybetay, myemitx0_dist, myemity0_dist, &
  
negative.code: 
      else if(mynex == zero .and. myney > zero) then ! halo in y
  20    continue
        myemity = myemity0*(myney+(rndm4()*cdist_smearY))**2
        ysigmay = sqrt(mybetay*myemity)
        xv2(j)  = ysigmay * sin_mb(twopi*rndm4())
        if(abs(xv2(j)).lt.cutoff) goto 20
  