anchor.diff:  !--elliptic beam x>z
                 if(sigman(1,imbb(i)).gt.sigman(2,imbb(i))) then
                   if(nbeaux(imbb(i)).eq.1.or.nbeaux(imbb(i)).eq.3) then
  -                  write(lout,"(a)") "BEAMBEAM> ERROR At each interaction point the beam must be either "//&
  +                  write(lerr,"(a)") "BEAMBEAM> ERROR At each interaction point the beam must be either "//&
                       "round or elliptical for all particles"
                     call prror(-1)
                   else
  
positive.code: !--elliptic beam z>x
                if(sigman(1,imbb(i)).lt.sigman(2,imbb(i))) then
                  if(nbeaux(imbb(i)).eq.1.or.nbeaux(imbb(i)).eq.2) then
                    write(lout,"(a)") "BEAMBEAM> ERROR At each interaction point the beam must be either "//&
                      "round or elliptical for all particles"
                    call prror(-1)
                  else
  
negative.code: !--round beam
                if(sigman(1,imbb(i)).eq.sigman(2,imbb(i))) then
                  if(nbeaux(imbb(i)).eq.2.or.nbeaux(imbb(i)).eq.3) then
                    write(lerr,"(a)") "BEAMBEAM> ERROR At each interaction point the beam must be either "//&
                      "round or elliptical for all particles"
                    call prror(-1)
                  else
  