1 | subroutine out1(ne,nbuffer,outcor,nsize) |
---|
2 | c---print coordinates at end of time step (ne=0), |
---|
3 | c---or at end of element ne. |
---|
4 | c------------------------------------------------------------------------ |
---|
5 | save |
---|
6 | c |
---|
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' |
---|
14 | c |
---|
15 | dimension outcor(8,nsize) |
---|
16 | c-------------------------------------------------------------------------- |
---|
17 | c* |
---|
18 | if(ngood.le.0)return |
---|
19 | if(ne.le.0)go to 30 |
---|
20 | c---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 |
---|
27 | c---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 |
---|
36 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|