1 |
subroutine instore( n, fc, gnorm0, is, m, jmin, jmax, |
2 |
& cold, ifail ) |
3 |
c |
4 |
c arguments |
5 |
c |
6 |
integer n, is, m, jmin, jmax, ifail |
7 |
double precision fc, gnorm0 |
8 |
logical cold |
9 |
c---- |
10 |
c |
11 |
common /xxstore/ itape, ntape, isize |
12 |
integer itape, ntape, isize |
13 |
|
14 |
ifail = 0 |
15 |
itape = 91 |
16 |
ntape = 92 |
17 |
isize = is |
18 |
cold = .true. |
19 |
|
20 |
cc one record = two arrays of real*(isize) |
21 |
|
22 |
open( itape |
23 |
$ , file = 'OPWARMI' |
24 |
$ , status = 'unknown' |
25 |
$ , form = 'formatted' |
26 |
$ , access = 'sequential' |
27 |
$ ) |
28 |
|
29 |
read( itape, *, end=800, err=900 ) n, fc, gnorm0, |
30 |
$ isize, m, jmin, jmax |
31 |
cold = .false. |
32 |
|
33 |
800 continue |
34 |
close(itape) |
35 |
|
36 |
print* |
37 |
write(*,'(a,i10,a,i2)') |
38 |
$ ' opening direct access with recl = ',n,' * ',isize |
39 |
print* |
40 |
|
41 |
open( ntape |
42 |
$ , file = 'OPWARMD' |
43 |
$ , status = 'unknown' |
44 |
$ , form = 'unformatted' |
45 |
$ , access = 'direct' |
46 |
$ , recl = n * isize |
47 |
$ ) |
48 |
|
49 |
return |
50 |
|
51 |
900 continue |
52 |
close(itape) |
53 |
ifail = 1 |
54 |
|
55 |
return |
56 |
end |