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

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

parmela pspa initial

File size: 1.7 KB
Line 
1      subroutine out5(ne,nbuffer,outcor,nsize)
2c---write coordinates on tape (tape2 or tape4)
3c-----------------------------------------------------------------------
4      save
5c
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'
16c
17      common/chars/title
18      common/errors/eraln(5,0:lmx)
19      common/tapes/ifl(2      )
20      dimension outcor(8,nsize)
21c--------------------------------------------------------------------------
22c*
23      ntape=2
24      if(ifl(ntape).ne.0)go to 10
25c---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'
38c--- 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
54c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
Note: See TracBrowser for help on using the repository browser.