1 | subroutine out3(ne,nbuffer,outcor,nsize) |
---|
2 | c plot pout vs pin and wout vs pin |
---|
3 | c----------------------------------------------------------------------- |
---|
4 | save |
---|
5 | c |
---|
6 | include 'param_sz.h' |
---|
7 | include 'constcom.h' |
---|
8 | include 'coordcom.h' |
---|
9 | include 'outcom.h' |
---|
10 | include 'syscom.h' |
---|
11 | include 'ucom.h' |
---|
12 | c |
---|
13 | common/com1/scale(10),ws,ps,zs,bgs,wts,ntape |
---|
14 | common/outbuf/tcor(8,imaa*imb) |
---|
15 | real x(2,imaa),y(2,imaa),x0,b0 |
---|
16 | dimension outcor(8,nsize) |
---|
17 | equivalence (tcor(1,imaa+1),x(1,1)) |
---|
18 | equivalence (tcor(1,imaa+imaa/2),y(1,1)) |
---|
19 | c-------------------------------------------------------------------------- |
---|
20 | c* |
---|
21 | if(ne.eq.0)return |
---|
22 | ntape=3 |
---|
23 | x0=w0/erest |
---|
24 | b0=sqrt(x0*(x0+2.))/(1.+x0) |
---|
25 | scale(1)=optcon(2)*b0*wavel/360. |
---|
26 | scale(3)=scale(1) |
---|
27 | scale(2)=optcon(3) |
---|
28 | scale(4)=optcon(4) |
---|
29 | ps=pr(ne) |
---|
30 | ws=wr(ne) |
---|
31 | zs=zloc(ne) |
---|
32 | do 20 np=1,nbuffer |
---|
33 | c get index of particle |
---|
34 | ind=indx(outcor(7,np)) |
---|
35 | x(1,np)=z0-cor(5,ind) |
---|
36 | x(2,np)=x(1,np) |
---|
37 | y(1,np)=outcor(5,np)-ps |
---|
38 | 20 y(2,np)=outcor(6,np)-ws |
---|
39 | write(ntape,30)optcon(2),scale(2),zs,nbuffer,optcon(2),scale(4) |
---|
40 | 30 format(/5x,' pout vs pin ',f4.0,' x ',f4.0,' zr=',f7.1,' nbuf=', |
---|
41 | 1 i4, ' wout vs pin ',f4.0,' x ',f7.3) |
---|
42 | call ttyplot(x,y,nbuffer,scale,ntape) |
---|
43 | write(ntape,40) |
---|
44 | 40 format(1x) |
---|
45 | return |
---|
46 | end |
---|
47 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|