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

Last change on this file since 315 was 12, checked in by lemeur, 12 years ago

parmela pspa initial

File size: 1.2 KB
Line 
1      subroutine out1(ne,nbuffer,outcor,nsize)
2c---print coordinates at end of time step (ne=0),
3c---or at end of element ne.
4c------------------------------------------------------------------------
5      save
6c
7      include 'param_sz.h'
8      include 'coordcom.h'
9      include 'ncordscom.h'
10      include 'outcom.h'
11      include 'syscom.h'
12      include 'tstepcom.h'
13      include 'ucom.h'
14c
15      dimension outcor(8,nsize)
16c--------------------------------------------------------------------------
17c*
18      if(ngood.le.0)return
19      if(ne.le.0)go to 30
20c---end of element
21      write(nout1,10) nbuf,ne,zloc(ne)
22   10 format(' nbuf=',i4,', ne=',i4,', z=',f6.1/
23     1 7x,'x',6x,'xp',7x,'y',6x,'yp',5x,'phi',7x,'w'/)
24      write(nout1,20) ((outcor(i,j),i=1,8),j=1,nbuffer)
25   20 format('+',2(0pf8.2,3pf8.4),0pf8.1,0pf8.3,i5,/)
26      return
27c---end of time step
28   30 continue
29      write(nnout,40) ngood,wt
30   40 format(' ngood=',i4,', omegat=',f8.1,/,
31     1 7x,'x',5x,'bgx',7x,'y',5x,'bgy',7x,'z',5x,'bgz   ne',/)
32      write(nnout,50) ((cord(i,j),i=1,7),j=1,ngood)
33   50 format('+',4f8.4,f8.2,f8.2,f5.0/)
34      return
35      end
36c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
Note: See TracBrowser for help on using the repository browser.