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++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*