[12] | 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++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|