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

      subroutine out2(ne,nbuffer,outcor,nsize)
c         output subroutine number 2
c         plot 2 profiles of all particles on the tty
c         plot their initial coordinates if indicated
c         plot the intermediate particle coordinates when indicated
c         output parameters
c             1=two initial graph types (defined on first output call)
c             2=two profile types       (defined on first output call)
c             3=two current graph types
c             4=unused
c         the tens digit of each parameter specifies the left hand
c         graph type while the units digit specifies the right hand grap
c         the initial and current graph type codes are%
c              0=no graph
c              1=x-xprime space
c              2=y-yprime space
c              3=phase-energy space (or z-zprime, if ne not zero)
c              4=x-y space
c          the profile type codes are%
c              0=  no profile
c              1-6=x,xp,y,yp,phi (or z), w (or zp)
c          optcon card parameters are graph scales for
c            x, xp, y, yp, phi, w, z, zp
c
c-----------------------------------------------------------------------
      save
c
      include 'param_sz.h'
      include 'constcom.h'
      include 'coordcom.h'
      include 'misccom.h'
      include 'ncordscom.h'
      include 'outcom.h'
      include 'tstepcom.h'
      include 'syscom.h'
      include 'var_char.h'
      include 'ucom.h'
c
      common/chars/title
      common/com1/scale(10),ws,ps,zs,bgs,wts,ntape
      common/outbuf/tcor(8,imaa*imb)
      real x(2      ,imaa),y(2,imaa)
      dimension outcor(8,nsize)
      equivalence (tcor(1,(imb-1)*imaa+1),x(1,1))
      equivalence (tcor(1,(imb-1)*imaa+imaa/2),y(1,1))
c--------------------------------------------------------------------------
c*
      data iflagg/0/,i7,i0,i1,i6/7,0,1,6/
c--------when iflagg = 0, title is printed, and nstart instead of ncell
c------printed on initial graph. iflagg becomes 1 when current graph is
c------called, and ncell is printed.  iflag controls print of heading on
c------profile plots.
      do 10 n=1,10
      scale(n)=optcon(n+4)
   10 continue
      ntape=nout2
      if(ne.ne.0)ntape=nout2
      if (iflagg.gt.0) go to 30
c         plot initial graphs
      write(nout2,*) title
      ws=w0
      ps=wt0
      zs=z0
      bgs=sqrt(ws/erest*(2.+ws/erest))
      wts=wt
      option=optcon(2)
      if (option.eq.0.0) go to 20
      call ggplot (i0,option,cor,i6,ngood,i1)
   20 continue
      iflagg=1
c         plot profiles
   30 continue
      option=optcon(3)
      opt=optcon(4)
      zs=cord(5,1)
      bgs=cord(6,1)
      wts=wt
      if(ne.eq.0)go to 35
      ws=wr(ne)
      ps=pr(ne)
      zs=zloc(ne)
   35 continue
      if(option.eq.0.)go to 40
      if(ne.eq.0)call pplot(ne,option,cord,8,ngood)
      if(ne.ne.0)call pplot(ne,option,outcor,8,nbuffer)
   40 continue
c         plot current particle coordinates
      if (opt.le.0.)go to 50
      if (ne.eq.0) call ggplot(ne,opt,cord,i7,ngood,i7)
      if (ne.ne.0) call ggplot(ne,opt,outcor,8,nbuffer,8)
   50 continue
      write(nnout,55)
   55 format(' returning from out2')
      return
c
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*