1 |
guez |
3 |
module in_out |
2 |
|
|
|
3 |
|
|
implicit none |
4 |
|
|
|
5 |
|
|
private prep_file, go_column, s_pr_mat, d_pr_mat |
6 |
|
|
interface pr_matrix |
7 |
|
|
module procedure s_pr_mat, d_pr_mat |
8 |
|
|
end interface |
9 |
|
|
|
10 |
|
|
contains |
11 |
|
|
|
12 |
|
|
!*********************************************************** |
13 |
|
|
|
14 |
|
|
integer function new_unit() |
15 |
|
|
|
16 |
|
|
logical opened, exist |
17 |
|
|
|
18 |
|
|
!------------------------------------------------------ |
19 |
|
|
|
20 |
|
|
new_unit = 0 |
21 |
|
|
do |
22 |
|
|
inquire(unit=new_unit, opened=opened, exist=exist) |
23 |
|
|
if (exist .and. .not. opened) exit |
24 |
|
|
new_unit = new_unit + 1 |
25 |
|
|
end do |
26 |
|
|
|
27 |
|
|
end function new_unit |
28 |
|
|
|
29 |
|
|
!*********************************************************** |
30 |
|
|
|
31 |
|
|
function csvread(file, first_r, first_c, last_r, last_c) |
32 |
|
|
|
33 |
|
|
! Reads comma-separated numeric values in a file. The |
34 |
|
|
! last column and/or last row parameters may be 0. This is |
35 |
|
|
! interpreted as "last in the file". |
36 |
|
|
|
37 |
|
|
real, pointer:: csvread(:,:) |
38 |
|
|
|
39 |
|
|
character(len=*), intent(in):: file |
40 |
|
|
integer, intent(in), optional:: first_r ! (first row to read) |
41 |
|
|
integer, intent(in), optional:: first_c ! (first column to read) |
42 |
|
|
integer, intent(in), optional:: last_r ! (last row to read) |
43 |
|
|
integer, intent(in), optional:: last_c ! (last column to read) |
44 |
|
|
|
45 |
|
|
! Variables local to the subprogram: |
46 |
|
|
integer i, unit |
47 |
|
|
integer f_r_loc ! (first row to read, local variable) |
48 |
|
|
integer f_c_loc ! (first column to read, local variable) |
49 |
|
|
integer l_r_loc ! (last row to read, local variable) |
50 |
|
|
integer l_c_loc ! (last column to read, local variable) |
51 |
|
|
|
52 |
|
|
!------------------------------------------------------ |
53 |
|
|
|
54 |
|
|
print *, 'Reading data from file "' // file // '"' |
55 |
|
|
unit = new_unit() |
56 |
|
|
open(unit, file=file, status='old', action='read', position='rewind') |
57 |
|
|
|
58 |
|
|
call prep_file(unit, first_r, first_c, last_r, last_c, f_r_loc, & |
59 |
|
|
f_c_loc, l_r_loc, l_c_loc) |
60 |
|
|
|
61 |
|
|
allocate(csvread(l_r_loc - f_r_loc + 1, l_c_loc - f_c_loc + 1)) |
62 |
|
|
|
63 |
|
|
do i = 1, l_r_loc - f_r_loc + 1 |
64 |
|
|
call go_column(unit, f_c_loc) |
65 |
|
|
read(unit, fmt=*) csvread(i, :) |
66 |
|
|
end do |
67 |
|
|
! (no implicit loop in read to allow partial reading of a line) |
68 |
|
|
|
69 |
|
|
close(unit) |
70 |
|
|
|
71 |
|
|
end function csvread |
72 |
|
|
|
73 |
|
|
!*********************************************************** |
74 |
|
|
|
75 |
|
|
function csvread_dp(file, first_r, first_c, last_r, last_c) |
76 |
|
|
|
77 |
|
|
! Reads comma-separated numeric values from a file, into a |
78 |
|
|
! double precision array. The last column and/or last row parameters may be |
79 |
|
|
! 0. This is interpreted as "last in the file". |
80 |
|
|
|
81 |
|
|
double precision, pointer:: csvread_dp(:,:) |
82 |
|
|
|
83 |
|
|
character(len=*), intent(in):: file |
84 |
|
|
integer, intent(in), optional:: first_r ! (first row to read) |
85 |
|
|
integer, intent(in), optional:: first_c ! (first column to read) |
86 |
|
|
integer, intent(in), optional:: last_r ! (last row to read) |
87 |
|
|
integer, intent(in), optional:: last_c ! (last column to read) |
88 |
|
|
|
89 |
|
|
! Variables local to the subprogram: |
90 |
|
|
integer i, unit |
91 |
|
|
integer f_r_loc ! (first row to read, local variable) |
92 |
|
|
integer f_c_loc ! (first column to read, local variable) |
93 |
|
|
integer l_r_loc ! (last row to read, local variable) |
94 |
|
|
integer l_c_loc ! (last column to read, local variable) |
95 |
|
|
|
96 |
|
|
!------------------------------------------------------ |
97 |
|
|
|
98 |
|
|
print *, 'Reading data from file "' // file // '"' |
99 |
|
|
unit = new_unit() |
100 |
|
|
open(unit, file=file, status='old', action='read', position='rewind') |
101 |
|
|
|
102 |
|
|
call prep_file(unit, first_r, first_c, last_r, last_c, f_r_loc, & |
103 |
|
|
f_c_loc, l_r_loc, l_c_loc) |
104 |
|
|
|
105 |
|
|
allocate(csvread_dp(l_r_loc - f_r_loc + 1, l_c_loc - f_c_loc + 1)) |
106 |
|
|
|
107 |
|
|
do i = 1, l_r_loc - f_r_loc + 1 |
108 |
|
|
call go_column(unit, f_c_loc) |
109 |
|
|
read(unit, fmt=*) csvread_dp(i, :) |
110 |
|
|
end do |
111 |
|
|
! (no implicit loop in read to allow partial reading of a line) |
112 |
|
|
|
113 |
|
|
close(unit) |
114 |
|
|
|
115 |
|
|
end function csvread_dp |
116 |
|
|
|
117 |
|
|
!*********************************************************** |
118 |
|
|
|
119 |
|
|
subroutine go_column(unit, column) |
120 |
|
|
|
121 |
|
|
! This subroutine is used by the various versions of "csvread". On |
122 |
|
|
! the current line of file, it advances to the input column. Columns are |
123 |
|
|
! assumend to be separated by commas. |
124 |
|
|
|
125 |
|
|
integer, intent(in):: unit ! logical unit for input file |
126 |
|
|
integer, intent(in):: column |
127 |
|
|
|
128 |
|
|
! Variables local to the subprogram: |
129 |
|
|
integer j |
130 |
|
|
character c |
131 |
|
|
|
132 |
|
|
!------------------------------------------------------ |
133 |
|
|
|
134 |
|
|
! Skip columns before "column" : |
135 |
|
|
j = 1 ! column index |
136 |
|
|
do while (j <= column - 1) |
137 |
|
|
read(unit, fmt='(a)', advance='no') c |
138 |
|
|
if (c == ',') j = j + 1 |
139 |
|
|
end do |
140 |
|
|
|
141 |
|
|
end subroutine go_column |
142 |
|
|
|
143 |
|
|
!*********************************************************** |
144 |
|
|
|
145 |
|
|
subroutine prep_file(unit, first_r, first_c, last_r, last_c, f_r_not_opt, & |
146 |
|
|
f_c_not_opt, l_r_not_opt, l_c_not_opt) |
147 |
|
|
|
148 |
|
|
! This subroutine is used by the various versions of "csvread". It |
149 |
|
|
! fills non-optional arguments: first and last row, first and last |
150 |
|
|
! column which will actually be read, taking information from the |
151 |
|
|
! file itself if necessary. It also positions the input file on the |
152 |
|
|
! first row to read. |
153 |
|
|
|
154 |
|
|
integer, intent(in):: unit ! logical unit for input file |
155 |
|
|
integer, intent(in), optional:: first_r ! (first row to read) |
156 |
|
|
integer, intent(in), optional:: first_c ! (first column to read) |
157 |
|
|
integer, intent(in), optional:: last_r ! (last row to read) |
158 |
|
|
integer, intent(in), optional:: last_c ! (last column to read) |
159 |
|
|
integer, intent(out):: f_r_not_opt ! (first row to read, not optional) |
160 |
|
|
integer, intent(out):: f_c_not_opt ! (first column to read, not optional) |
161 |
|
|
integer, intent(out):: l_r_not_opt ! (last row to read, not optional) |
162 |
|
|
integer, intent(out):: l_c_not_opt ! (last column to read, not optional) |
163 |
|
|
|
164 |
|
|
! Variables local to the subprogram: |
165 |
|
|
integer iostat, i |
166 |
|
|
character c |
167 |
|
|
|
168 |
|
|
!------------------------------------------------------ |
169 |
|
|
|
170 |
|
|
f_r_not_opt = opt_merge(first_r, 1) |
171 |
|
|
f_c_not_opt = opt_merge(first_c, 1) |
172 |
|
|
l_r_not_opt = opt_merge(last_r, 0) |
173 |
|
|
l_c_not_opt = opt_merge(last_c, 0) |
174 |
|
|
|
175 |
|
|
if (l_r_not_opt == 0) then |
176 |
|
|
! Count the number of lines in the file: |
177 |
|
|
i = 0 |
178 |
|
|
do |
179 |
|
|
read(unit, fmt=*, iostat=iostat) |
180 |
|
|
if (iostat /= 0) exit |
181 |
|
|
i = i + 1 |
182 |
|
|
end do |
183 |
|
|
l_r_not_opt = i |
184 |
|
|
if (l_r_not_opt == 0) stop 'Empty file.' |
185 |
|
|
|
186 |
|
|
rewind(unit) |
187 |
|
|
end if |
188 |
|
|
|
189 |
|
|
! Go to first row to read: |
190 |
|
|
do i = 1, f_r_not_opt - 1 |
191 |
|
|
read(unit, fmt=*) |
192 |
|
|
end do |
193 |
|
|
|
194 |
|
|
if (l_c_not_opt == 0) then |
195 |
|
|
! Count the number of values per line: |
196 |
|
|
i = 0 |
197 |
|
|
do |
198 |
|
|
read(unit, fmt='(a)', advance='no', iostat=iostat) c |
199 |
|
|
if (iostat /= 0) exit |
200 |
|
|
if (c == ',') i = i + 1 |
201 |
|
|
end do |
202 |
|
|
l_c_not_opt = i + 1 |
203 |
|
|
|
204 |
|
|
backspace(unit) |
205 |
|
|
end if |
206 |
|
|
|
207 |
|
|
print *, 'Reading column(s) ', f_c_not_opt, ':', l_c_not_opt, & |
208 |
|
|
', row(s) ', f_r_not_opt, ':', l_r_not_opt |
209 |
|
|
|
210 |
|
|
end subroutine prep_file |
211 |
|
|
|
212 |
|
|
!*********************************************************** |
213 |
|
|
|
214 |
|
|
integer function opt_merge(param, default) |
215 |
|
|
|
216 |
|
|
! Analogous to the intrinsic procedure "merge" : merges an |
217 |
|
|
! optional parameter and a default value depending on the |
218 |
|
|
! presence of the optional parameter. |
219 |
|
|
|
220 |
|
|
integer, intent(in), optional:: param |
221 |
|
|
integer, intent(in):: default |
222 |
|
|
|
223 |
|
|
!-------------- |
224 |
|
|
|
225 |
|
|
if (present(param)) then |
226 |
|
|
opt_merge = param |
227 |
|
|
else |
228 |
|
|
opt_merge = default |
229 |
|
|
end if |
230 |
|
|
|
231 |
|
|
end function opt_merge |
232 |
|
|
|
233 |
|
|
!*********************************************************** |
234 |
|
|
|
235 |
|
|
subroutine s_pr_mat(name, a) |
236 |
|
|
|
237 |
|
|
! This subroutine prints a rank 2 real matrix. |
238 |
|
|
|
239 |
|
|
character(len=*), intent(in):: name |
240 |
|
|
real, intent(in):: a(:,:) |
241 |
|
|
|
242 |
|
|
character(len=20) fmt |
243 |
|
|
integer n_lines, n_col, i |
244 |
|
|
|
245 |
|
|
!----------------- |
246 |
|
|
|
247 |
|
|
n_lines = size(a, 1) |
248 |
|
|
n_col = size(a, 2) |
249 |
|
|
if (n_lines <= 10 .and. n_col <= 5) then |
250 |
|
|
print *, name, ":" |
251 |
|
|
write(unit=fmt, fmt='("(1p, ", i0, "(g10.3: 1X))")') n_col |
252 |
|
|
do i = 1, n_lines |
253 |
|
|
print fmt, a(i, :) |
254 |
|
|
end do |
255 |
|
|
else |
256 |
|
|
print *, '"', name, '" is too big to print.' |
257 |
|
|
end if |
258 |
|
|
|
259 |
|
|
end subroutine s_pr_mat |
260 |
|
|
|
261 |
|
|
!*********************************************************** |
262 |
|
|
|
263 |
|
|
subroutine d_pr_mat(name, a) |
264 |
|
|
|
265 |
|
|
! This subroutine prints a rank 2 double precision matrix. |
266 |
|
|
|
267 |
|
|
character(len=*), intent(in):: name |
268 |
|
|
double precision, intent(in):: a(:,:) |
269 |
|
|
|
270 |
|
|
character(len=20) fmt |
271 |
|
|
integer n_lines, n_col, i |
272 |
|
|
|
273 |
|
|
!----------------- |
274 |
|
|
|
275 |
|
|
n_lines = size(a, 1) |
276 |
|
|
n_col = size(a, 2) |
277 |
|
|
if (n_lines <= 10 .and. n_col <= 5) then |
278 |
|
|
print *, name, ":" |
279 |
|
|
write(unit=fmt, fmt='("(1p, ", i0, "(g8.1: 1X))")') n_col |
280 |
|
|
do i = 1, n_lines |
281 |
|
|
print fmt, a(i, :) |
282 |
|
|
end do |
283 |
|
|
else |
284 |
|
|
print *, '"', name, '" is too big to print.' |
285 |
|
|
end if |
286 |
|
|
|
287 |
|
|
end subroutine d_pr_mat |
288 |
|
|
|
289 |
|
|
end module in_out |