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