1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
|
program open_close_read_write
implicit none
character(len=20) :: msg = 'Some error occured'
real, dimension(100) :: x, y
real, dimension(100) :: p, q
integer, target :: i, u = 1
integer, pointer :: u_ptr;
u_ptr => u
do i = 1, 100
x(i) = i * 0.1
y(i) = sin(x(i)) * (1 - cos(x(i)/3.0))
end do
open(unit=u_ptr, file='data.dat', status='replace')
do i = 1, 100
write(u, '(10F8.2)') x(i), y(i)
end do
close(1, err=999, iomsg=msg, iostat=u)
999 open (2, file='data.dat', status='old')
do i = 1, 100
read(2, *) p(i), q(i)
end do
close(2)
do i = 1, 100
write(*, *) p(i), q(i)
end do
end program
|
