1 | subroutine out5(ne,nbuffer,outcor,nsize) |
---|
2 | c---write coordinates on tape (tape2 or tape4) |
---|
3 | c----------------------------------------------------------------------- |
---|
4 | save |
---|
5 | c |
---|
6 | include 'param_sz.h' |
---|
7 | include 'constcom.h' |
---|
8 | include 'coordcom.h' |
---|
9 | include 'misccom.h' |
---|
10 | include 'ncordscom.h' |
---|
11 | include 'outcom.h' |
---|
12 | include 'tstepcom.h' |
---|
13 | include 'syscom.h' |
---|
14 | include 'var_char.h' |
---|
15 | include 'ucom.h' |
---|
16 | c |
---|
17 | common/chars/title |
---|
18 | common/errors/eraln(5,0:lmx) |
---|
19 | common/tapes/ifl(2 ) |
---|
20 | dimension outcor(8,nsize) |
---|
21 | c-------------------------------------------------------------------------- |
---|
22 | c* |
---|
23 | ntape=2 |
---|
24 | if(ifl(ntape).ne.0)go to 10 |
---|
25 | c---write basic information on tape |
---|
26 | open(unit=ntape,file='tape2',form='unformatted', |
---|
27 | * access='sequential',status='unknown') |
---|
28 | ifl(ntape)=1 |
---|
29 | write(ntape,err=15) irun,title,ac,bc,cc,nsize,ngood |
---|
30 | write(ntape,err=15) ((cor(i,j),i=1,6),j=1,nsize) |
---|
31 | 10 continue |
---|
32 | if(ne.eq.0)go to 16 |
---|
33 | write(ntape,err=15) nbuffer,ne,pr(ne),wr(ne),outcor |
---|
34 | return |
---|
35 | 15 continue |
---|
36 | call appendparm |
---|
37 | stop ' Abnormal stop error writing to tape2' |
---|
38 | c--- output cordinates at end of time step on unit 3 |
---|
39 | 16 if(ifl(1).eq.0) then |
---|
40 | ntape4=4 |
---|
41 | open(unit=ntape4,file='tape4',form='unformatted', |
---|
42 | * access='sequential',status='unknown') |
---|
43 | ifl(1)=1 |
---|
44 | nstep=1 |
---|
45 | endif |
---|
46 | write(ntape4,err=17)nstep,ngood, |
---|
47 | * ((cord(i,j),i=1,5,2),j=1,ngood) |
---|
48 | nstep=nstep+1 |
---|
49 | return |
---|
50 | 17 continue |
---|
51 | call appendparm |
---|
52 | stop ' Abnormal stop error writing to tape4' |
---|
53 | end |
---|
54 | c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++* |
---|