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