source: PSPA/parmelaPSPA/trunk/sortie.f @ 315

      subroutine sortie(np,ngood,wtp)
c-----------------------------------------------------------------------
c   The subroutine is for the graphic at LAL 
c-----------------------------------------------------------------------
c
      include 'param_sz.h'
      include 'constcom.h'
      include 'coordcom.h'
      include 'flagcom.h'
      include 'misccom.h'
      include 'pcordcom.h'
      include 'syscom.h'
      include 'ucom.h'
c
      common/part/nparticle(imaa)
c--------------------------------------------------------------------------
c*
      ne=rne
      xx=x
      yy=y  
      xxp=bgx/bgz
      yyp=bgy/bgz
c---if there is a background bfield the beam is rotating.
c   therefor transform to the rotation frame of referance.
      if(ifld.ne.0)then
         if(ifoclal.eq.1) then
           cay=bfldlal(z)/(brhof*sqrt(gamma**2-1.))
         else
           cay=bfld(z,sqrt(x**2+y**2),brfld)/(brhof*sqrt(gamma**2-1.))
         endif
         xxp=bgx/bgz-.5*cay*y
         yyp=bgy/bgz+.5*cay*x
      endif
      phi=wtp    
      wz=(gamma-1.)*erest   
      if(ne.eq.0) go to 10
      write(ndes1,*) xx,xxp,bgx,yy,yyp,bgy,z,bgz,phi,wz,ne,np,
     1   ngood,nparticle(np),phizero(np),ksi1(np),ksi2(np),ksi3(np)
c     unit ndes1 end of element z=cte
 10   continue
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*