Write system output files.
26
29
30 implicit none
31 integer*4 :: sfs, temp, i, SDOFmon, jdof
32 real*8 :: tt1tmp
33 character*32 :: fmt_spec
34
35 sdofmon=10*(mpi_id+1)+7
36
37 sfs = 0
40 if(temp.eq.1) sfs = temp
41 enddo
42
44
45 if (sfs.eq.0) then
46 open(sdofmon,file=
sdofdisplx,position=
'append')
47 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
50 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
51 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%tempU1(jdof,1), jdof=1,
sys(i)%NDOF)
52 else
53 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempU1(
sys(i)%NDOF,1)
54 endif
55 enddo
56
57 write(sdofmon,"(A1)") " "
58 close(sdofmon)
59
60 open(sdofmon,file=
sdofdisply,position=
'append')
61 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
64 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
65 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%tempU1(jdof,2), jdof=1,
sys(i)%NDOF)
66 else
67 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempU1(
sys(i)%NDOF,2)
68 endif
69 enddo
70
71 write(sdofmon,"(A1)") " "
72 close(sdofmon)
73
74 open(sdofmon,file=
sdofdisplz,position=
'append')
75 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
78 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
79 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%tempU1(jdof,3), jdof=1,
sys(i)%NDOF)
80 else
81 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempU1(
sys(i)%NDOF,3)
82 endif
83 enddo
84
85 write(sdofmon,"(A1)") " "
86 close(sdofmon)
87
88
89
91 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
93 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,1)
94 enddo
95
96 write(sdofmon,"(A1)") " "
97 close(sdofmon)
98
100 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
102 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,2)
103 enddo
104
105 write(sdofmon,"(A1)") " "
106 close(sdofmon)
107
109 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
111 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,3)
112 enddo
113
114 write(sdofmon,"(A1)") " "
115 close(sdofmon)
116 elseif(sfs.eq.1) then
117
118 open(sdofmon,file=
strdisplx,position=
'append')
119 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
121 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%u(1,1)
122 enddo
123
124 write(sdofmon,"(A1)") " "
125 close(sdofmon)
126
127 open(sdofmon,file=
strdisply,position=
'append')
128 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
130 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%u(1,2)
131 enddo
132
133 write(sdofmon,"(A1)") " "
134 close(sdofmon)
135
136 open(sdofmon,file=
grddisplx,position=
'append')
137 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
139 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,1)
140 enddo
141
142 write(sdofmon,"(A1)") " "
143 close(sdofmon)
144
145 open(sdofmon,file=
grddisply,position=
'append')
146 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
148 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,2)
149 enddo
150
151 write(sdofmon,"(A1)") " "
152 close(sdofmon)
153
154 open(sdofmon,file=
grddisplz,position=
'append')
155 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
157 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,3)
158 enddo
159
160 write(sdofmon,"(A1)") " "
161 close(sdofmon)
162
163 open(sdofmon,file=
fnddisplx,position=
'append')
164 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
166 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%u(2,1)
167 enddo
168
169 write(sdofmon,"(A1)") " "
170 close(sdofmon)
171
172 open(sdofmon,file=
fnddisply,position=
'append')
173 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
175 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%u(2,2)
176 enddo
177
178 write(sdofmon,"(A1)") " "
179 close(sdofmon)
180
181 open(sdofmon,file=
fnddisplrx,position=
'append')
182 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
184 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%u(3,1)
185 enddo
186
187 write(sdofmon,"(A1)") " "
188 close(sdofmon)
189
190 open(sdofmon,file=
fnddisplry,position=
'append')
191 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
193 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%u(3,2)
194 enddo
195
196 write(sdofmon,"(A1)") " "
197 close(sdofmon)
198
199 open(sdofmon,file=
fnddisplzx,position=
'append')
200 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
202 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%u(4,1)
203 enddo
204
205 write(sdofmon,"(A1)") " "
206 close(sdofmon)
207
208 open(sdofmon,file=
fnddisplzy,position=
'append')
209 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
211 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%u(4,2)
212 enddo
213
214 write(sdofmon,"(A1)") " "
215 close(sdofmon)
216 endif
217 endif
218
219
221
222 if(sfs.eq.0) then
223
224
225 open(sdofmon,file=
sdofaccx,position=
'append')
226 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
229 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
230 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%tempRA1(jdof,1), jdof=1,
sys(i)%NDOF)
231 else
232 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempRA1(
sys(i)%NDOF,1)
233 endif
234 enddo
235 write(sdofmon,"(A1)") " "
236 close(sdofmon)
237
238 open(sdofmon,file=
sdofaccy,position=
'append')
239 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
242 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
243 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%tempRA1(jdof,2), jdof=1,
sys(i)%NDOF)
244 else
245 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempRA1(
sys(i)%NDOF,2)
246 endif
247 enddo
248 write(sdofmon,"(A1)") " "
249 close(sdofmon)
250
251 open(sdofmon,file=
sdofaccz,position=
'append')
252 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
255 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
256 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%tempRA1(jdof,3), jdof=1,
sys(i)%NDOF)
257 else
258 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempRA1(
sys(i)%NDOF,3)
259 endif
260 enddo
261 write(sdofmon,"(A1)") " "
262 close(sdofmon)
263
264
265
266 open(sdofmon,file=
sdofgraccx,position=
'append')
267 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
269 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofag(i,1)
270 enddo
271 write(sdofmon,"(A1)") " "
272 close(sdofmon)
273
274 open(sdofmon,file=
sdofgraccy,position=
'append')
275 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
277 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofag(i,2)
278 enddo
279 write(sdofmon,"(A1)") " "
280 close(sdofmon)
281
282 open(sdofmon,file=
sdofgraccz,position=
'append')
283 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
285 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofag(i,3)
286 enddo
287 write(sdofmon,"(A1)") " "
288 close(sdofmon)
289 elseif(sfs.eq.1) then
290
291 open(sdofmon,file=
straccx,position=
'append')
292 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
294 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%a(1,1)
295 enddo
296
297 write(sdofmon,"(A1)") " "
298 close(sdofmon)
299
300 open(sdofmon,file=
straccy,position=
'append')
301 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
303 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%a(1,2)
304 enddo
305
306 write(sdofmon,"(A1)") " "
307 close(sdofmon)
308
309 open(sdofmon,file=
grdaccx,position=
'append')
310 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
312 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofag(i,1)
313 enddo
314
315 write(sdofmon,"(A1)") " "
316 close(sdofmon)
317
318 open(sdofmon,file=
grdaccy,position=
'append')
319 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
321 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofag(i,2)
322 enddo
323
324 write(sdofmon,"(A1)") " "
325 close(sdofmon)
326
327 open(sdofmon,file=
grdaccz,position=
'append')
328 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
330 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofag(i,3)
331 enddo
332
333 write(sdofmon,"(A1)") " "
334 close(sdofmon)
335
336 open(sdofmon,file=
fndaccx,position=
'append')
337 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
339 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%a(2,1)
340 enddo
341
342 write(sdofmon,"(A1)") " "
343 close(sdofmon)
344
345 open(sdofmon,file=
fndaccy,position=
'append')
346 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
348 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%a(2,2)
349 enddo
350
351 write(sdofmon,"(A1)") " "
352 close(sdofmon)
353
354 open(sdofmon,file=
fndaccrx,position=
'append')
355 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
357 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%a(3,1)
358 enddo
359
360 write(sdofmon,"(A1)") " "
361 close(sdofmon)
362
363 open(sdofmon,file=
fndaccry,position=
'append')
364 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
366 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%a(3,2)
367 enddo
368
369 write(sdofmon,"(A1)") " "
370 close(sdofmon)
371
372 open(sdofmon,file=
fndacczx,position=
'append')
373 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
375 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%a(4,1)
376 enddo
377
378 write(sdofmon,"(A1)") " "
379 close(sdofmon)
380
381 open(sdofmon,file=
fndacczy,position=
'append')
382 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
384 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%a(4,2)
385 enddo
386
387 write(sdofmon,"(A1)") " "
388 close(sdofmon)
389 endif
390 endif
391
393
394 if(sfs.eq.0) then
395 open(sdofmon,file=
sdoffx,position=
'append')
396 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
399 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
400 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%IntForce(jdof,1), jdof=1,
sys(i)%NDOF)
401 else
402 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%IntForce(
sys(i)%NDOF,1)
403 endif
404 end do
405 write(sdofmon,"(A1)") " "
406 close(sdofmon)
407
408 open(sdofmon,file=
sdoffy,position=
'append')
409 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
412 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
413 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%IntForce(jdof,2), jdof=1,
sys(i)%NDOF)
414 else
415 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%IntForce(
sys(i)%NDOF,2)
416 endif
417 end do
418 write(sdofmon,"(A1)") " "
419 close(sdofmon)
420
421 open(sdofmon,file=
sdoffz,position=
'append')
422 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
425 write(fmt_spec,
'(A1, I0, A6)')
'(',
sys(i)%NDOF,
'E16.7)'
426 write(sdofmon,fmt_spec,advance=
'NO') (
sys(i)%IntForce(jdof,3), jdof=1,
sys(i)%NDOF)
427 else
428 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%IntForce(
sys(i)%NDOF,3)
429 endif
430 end do
431 write(sdofmon,"(A1)") " "
432 close(sdofmon)
433 elseif(sfs.eq.1) then
434
435 open(sdofmon,file=
strfx,position=
'append')
436 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
438 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%fs(1)
439 enddo
440
441 write(sdofmon,"(A1)") " "
442 close(sdofmon)
443
444 open(sdofmon,file=
strfy,position=
'append')
445 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
447 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%fs(2)
448 enddo
449
450 write(sdofmon,"(A1)") " "
451 close(sdofmon)
452
453 open(sdofmon,file=
fndfx,position=
'append')
454 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
456 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%fb(1)
457 enddo
458
459 write(sdofmon,"(A1)") " "
460 close(sdofmon)
461
462 open(sdofmon,file=
fndfy,position=
'append')
463 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
465 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%fb(2)
466 enddo
467
468 write(sdofmon,"(A1)") " "
469 close(sdofmon)
470
471 open(sdofmon,file=
intfx,position=
'append')
472 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
474 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofforceinput(3*(i-1)+1)
475 end do
476 write(sdofmon,"(A1)") " "
477 close(sdofmon)
478
479 open(sdofmon,file=
intfy,position=
'append')
480 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
482 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofforceinput(3*(i-1)+2)
483 end do
484 write(sdofmon,"(A1)") " "
485 close(sdofmon)
486
487 open(sdofmon,file=
intfz,position=
'append')
488 write(sdofmon,"(E16.7)",advance='NO') tt1tmp
490 if(
sys(i)%SFS.eq.1)
write(sdofmon,
"(E16.7)",advance=
'NO')
sdofforceinput(3*(i-1)+3)
491 end do
492 write(sdofmon,"(A1)") " "
493 close(sdofmon)
494 endif
495 endif
496
497 return
Contains parameters for MDOF.
character *100 sdofgraccz
SDOF base ground acceleration.
character *100 fnddisplzy
4DOF foundation displacement
character *100 grdaccz
4DOF ground acceleration
character *100 sdofgraccx
integer *4, dimension(3) sdofout
displ, acc, f_react
character *100 fndacczy
4DOF foundation acceleration
character *100 sdofaccz
SDOF total acceleration.
character *100 fndfy
4DOF superstructure and foundation shear force
character *100 sdofgrdisplx
character *100 sdofdisplz
SDOF displacement.
character *100 sdofgrdisply
character *100 fnddisplrx
integer *4 flag_outatalldofs
character *100 straccy
4DOF structure acceleration
character *100 sdofgrdisplz
SDOF base displacement.
character *100 strdisply
4DOF structure displacement
character *100 sdofgraccy
character *100 sdofdisply
real *8, dimension(:), allocatable sdofforceinput
real *8, dimension(:,:), allocatable sdofgd
real *8, dimension(:,:), allocatable sdofag
character *100 fnddisplry
character *100 sdoffz
SDOF reaction force.
character *100 fnddisplzx
character *100 grddisplz
4DOF ground displacement
type(system), dimension(:), allocatable sys
SDOF system.
character *100 sdofdisplx
character *100 intfz
4DOF interaction forces
Contains a subset of SPEED paramters (used in TIME_LOOP)