source: PSPA/parmelaPSPA/trunk/out3.f @ 478

      subroutine out3(ne,nbuffer,outcor,nsize)
c        plot pout vs pin and wout vs pin
c-----------------------------------------------------------------------
      save
c
      include 'param_sz.h'
      include 'constcom.h'
      include 'coordcom.h'
      include 'outcom.h'
      include 'syscom.h'
      include 'ucom.h'
c
      common/com1/scale(10),ws,ps,zs,bgs,wts,ntape
      common/outbuf/tcor(8,imaa*imb)
      real x(2,imaa),y(2,imaa),x0,b0
      dimension outcor(8,nsize)
      equivalence (tcor(1,imaa+1),x(1,1))
      equivalence (tcor(1,imaa+imaa/2),y(1,1))
c--------------------------------------------------------------------------
c*
      if(ne.eq.0)return
      ntape=3
      x0=w0/erest
      b0=sqrt(x0*(x0+2.))/(1.+x0)
      scale(1)=optcon(2)*b0*wavel/360.
      scale(3)=scale(1)
      scale(2)=optcon(3)
      scale(4)=optcon(4)
      ps=pr(ne)
      ws=wr(ne)
      zs=zloc(ne)
      do 20 np=1,nbuffer
c           get index of particle
      ind=indx(outcor(7,np))
      x(1,np)=z0-cor(5,ind)
      x(2,np)=x(1,np)
      y(1,np)=outcor(5,np)-ps
   20 y(2,np)=outcor(6,np)-ws
      write(ntape,30)optcon(2),scale(2),zs,nbuffer,optcon(2),scale(4)
   30 format(/5x,' pout vs pin ',f4.0,' x ',f4.0,' zr=',f7.1,' nbuf=',
     1 i4, '  wout vs pin ',f4.0,' x ',f7.3)
      call ttyplot(x,y,nbuffer,scale,ntape)
      write(ntape,40)
   40 format(1x)
      return
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*