source: PSPA/parmelaPSPA/trunk/file.f @ 12

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

parmela pspa initial

File size: 4.4 KB
Line 
1c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
2c     system dependant                                                 *
3c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
4      subroutine ouvrir
5c
6      include 'ucom.h'
7      include 'debug.h'
8      include 'flagcom.h'
9c
10c------------------------- open --------------------------------------
11c Partie transferee dans le main 14/09/00
12c      open(unit=nin,file='parmin',status='old',
13c     *     access='sequential',form='formatted')
14c      idebug=0
15c      if(idebug.eq.0) then
16c         open(unit=nnout,status='scratch',
17c     *        access='sequential',form='formatted')
18c      else
19c         open(unit=nnout,file='parmint',status='new',
20c     *        access='sequential',form='formatted')
21c      endif
22c-----------------------------------------------------------------------
23c open modifie le 14/11/2007 pour utilisation sous Linux
24c declaration des gros fichiers sur le disque tmp de la machine et non du
25c cluster
26c
27      open(unit=nav,access='sequential',
28     *     form='formatted',status='scratch')
29      open(unit=ndes1,file='parmdesz',form='formatted',
30     *     status='unknown')
31      open(unit=ndes2,file='/tmp/parmdest',form='unformatted',
32     *     status='unknown')
33      open(unit=nemit,access='sequential',
34     *     form='formatted',status='scratch')
35      open(unit=ndiag,form='formatted',status='scratch')
36      open(unit=nsemit,file='/tmp/parmemit2',form='unformatted',
37     *     status='unknown')
38      open(unit=nsnap,file='/tmp/parmsnap',form='unformatted',
39     *     status='unknown')
40c--------------------------------------------------------------------------
41      return
42      end
43c*********************************************************************
44c     vax version
45      subroutine orest1(i)
46c
47      include 'ucom.h'
48c
49      character*1 kfic1
50      character*36 bmnom1
51      write(kfic1,8337)i
52 8337 format(i1)
53      bmnom1='savecor'//kfic1
54      open (unit=nsav,file=bmnom1,access='sequential',
55     *      form='unformatted',status='old')
56      return
57      end
58c*********************************************************************
59c     vax version
60      subroutine orest2
61c
62      include 'ucom.h'
63c
64      open (unit=nsav,file='/tmp/savecor',access='sequential',
65     *      form='unformatted',status='old')
66      return
67      end
68c*********************************************************************
69c     vax version
70      subroutine osave1(i)
71c
72      include 'ucom.h'
73c
74      character*1 kfic2
75      character*36 bmnom2
76      write(kfic2,8338)i
77 8338 format(i1)
78      bmnom2='savecor'//kfic2
79      open (unit=nsav,file='/tmp/bmnom2',access='sequential',
80     *      form='unformatted',status='new')
81      return
82      end
83c*********************************************************************
84c     vax version
85      subroutine inscgrid(fexist) !02/92
86      logical fexist
87      inquire(file='scgrid',exist=fexist)
88      return
89      end
90c*********************************************************************
91c     vax version
92c     file back pour le back bombardment
93      subroutine opeback(iope)
94c
95      include 'ucom.h'
96c
97      if(iope.eq.1) then
98      open(unit=nback,file='parmback',status='new',form='unformatted')
99      open(unit=nback1,file='parmback1',status='new',form='unformatted')
100      else
101      open(unit=nback1,file='parmback1',status='new',form='unformatted')
102      endif
103      return
104      end
105c*********************************************************************
106c     vax version
107      subroutine appendparm
108c
109      include 'ucom.h'
110c
111      character*1 k(80)
112      integer bmout
113c---------------------------------------------------------------------
114      bmout=90
115      open(unit=bmout,file='parmout',access='sequential',
116     *     status='unknown',form='formatted')
117      rewind(nnout)
118      rewind(nemit)
119      rewind(nav)
120      rewind(ndiag)
121 100  format(80a1)
122 10   continue
123      read(nnout,100,end=20) (k(i),i=1,80)
124      write(bmout,100) (k(i),i=1,80)
125      go to 10
126 20   continue
127      read(nemit,100,end=30) (k(i),i=1,80)
128      write(bmout,100) (k(i),i=1,80)
129      go to 20
130 30   continue
131      read(nav,100,end=40) (k(i),i=1,80)
132      write(bmout,100) (k(i),i=1,80)
133      go to 30
134 40   continue
135      read(ndiag,100,end=50) (k(i),i=1,80)
136      write(bmout,100) (k(i),i=1,80)     
137      go to 40
138 50   continue
139c-----------------------------------------------------------------------
140      end
141c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*
142
Note: See TracBrowser for help on using the repository browser.