SPEED
Functions/Subroutines
MAKE_SDOF_OUTPUT_FILES.f90 File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine make_sdof_output_files
 Creates files to store oscillator motion.
 

Function/Subroutine Documentation

◆ make_sdof_output_files()

subroutine make_sdof_output_files

Creates files to store oscillator motion.

Author
Aline Herlin
Date
November, 2020
Version
1.0

Definition at line 24 of file MAKE_SDOF_OUTPUT_FILES.f90.

25
26 use speed_sci
27 use speed_timeloop
28 use speed_par
29
30 implicit none
31 character*100 :: file
32 integer*4 :: temp, sfs, ii, SDOFmon
33
34 sdofmon = 701 + mpi_id
35
36 if (n_bld.gt.0) then
37
38 sfs = 0
39 do ii = 1, n_bld
40 temp = sys(ii)%SFS
41 if(temp.eq.1) sfs = temp
42 enddo
43
44 if (sdofout(1).eq.1) then ! displacement
45
46 if(sfs.eq.0) then
47
48 !!! structure relative displacement
49
50 file="SDOF000000.DX"
51
52 if(len_trim(monitor_file) .ne. 70) then
53 sdofdisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
54 else
55 sdofdisplx = file
56 endif
57
58 open(sdofmon,file=sdofdisplx,status='replace')
59 close(sdofmon)
60
61 file="SDOF000000.DY"
62
63 if(len_trim(monitor_file) .ne. 70) then
64 sdofdisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
65 else
66 sdofdisply = file
67 endif
68
69 open(sdofmon,file=sdofdisply,status='replace')
70 close(sdofmon)
71
72 file="SDOF000000.DZ"
73
74 if(len_trim(monitor_file) .ne. 70) then
75 sdofdisplz = monitor_file(1:len_trim(monitor_file)) // '/' // file
76 else
77 sdofdisplz = file
78 endif
79
80 open(sdofmon,file=sdofdisplz,status='replace')
81 close(sdofmon)
82
83 !!! soil displacement
84
85 file="SDOF000000.GDX"
86
87 if(len_trim(monitor_file) .ne. 70) then
88 sdofgrdisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
89 else
90 sdofgrdisplx = file
91 endif
92
93 open(sdofmon,file=sdofgrdisplx,status='replace')
94 close(sdofmon)
95
96 file="SDOF000000.GDY"
97
98 if(len_trim(monitor_file) .ne. 70) then
99 sdofgrdisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
100 else
101 sdofgrdisply = file
102 endif
103
104 open(sdofmon,file=sdofgrdisply,status='replace')
105 close(sdofmon)
106
107 file="SDOF000000.GDZ"
108
109 if(len_trim(monitor_file) .ne. 70) then
110 sdofgrdisplz = monitor_file(1:len_trim(monitor_file)) // '/' // file
111 else
112 sdofgrdisplz = file
113 endif
114
115 open(sdofmon,file=sdofgrdisplz,status='replace')
116 close(sdofmon)
117
118 elseif(sfs.eq.1) then
119
120 !!! structure relative displacement
121
122 file="STR000000.DX"
123
124 if(len_trim(monitor_file) .ne. 70) then
125 strdisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
126 else
127 strdisplx = file
128 endif
129
130 open(sdofmon,file=strdisplx,status='replace')
131 close(sdofmon)
132
133 file="STR000000.DY"
134
135 if(len_trim(monitor_file) .ne. 70) then
136 strdisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
137 else
138 strdisply = file
139 endif
140
141 open(sdofmon,file=strdisply,status='replace')
142 close(sdofmon)
143
144 !!! soil displacement
145
146 file="GRD000000.DX"
147
148 if(len_trim(monitor_file) .ne. 70) then
149 grddisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
150 else
151 grddisplx = file
152 endif
153
154 open(sdofmon,file=grddisplx,status='replace')
155 close(sdofmon)
156
157 file="GRD000000.DY"
158
159 if(len_trim(monitor_file) .ne. 70) then
160 grddisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
161 else
162 grddisply = file
163 endif
164
165 open(sdofmon,file=grddisply,status='replace')
166 close(sdofmon)
167
168 file="GRD000000.DZ"
169
170 if(len_trim(monitor_file) .ne. 70) then
171 grddisplz = monitor_file(1:len_trim(monitor_file)) // '/' // file
172 else
173 grddisplz = file
174 endif
175
176 open(sdofmon,file=grddisplz,status='replace')
177 close(sdofmon)
178
179 !!! foundation displacement
180
181 file="FND000000.DX"
182
183 if(len_trim(monitor_file) .ne. 70) then
184 fnddisplx = monitor_file(1:len_trim(monitor_file)) // '/' // file
185 else
186 fnddisplx = file
187 endif
188
189 open(sdofmon,file=fnddisplx,status='replace')
190 close(sdofmon)
191
192 file="FND000000.DY"
193
194 if(len_trim(monitor_file) .ne. 70) then
195 fnddisply = monitor_file(1:len_trim(monitor_file)) // '/' // file
196 else
197 fnddisply = file
198 endif
199
200 open(sdofmon,file=fnddisply,status='replace')
201 close(sdofmon)
202
203 file="FND000000.DRX"
204
205 if(len_trim(monitor_file) .ne. 70) then
206 fnddisplrx = monitor_file(1:len_trim(monitor_file)) // '/' // file
207 else
208 fnddisplrx = file
209 endif
210
211 open(sdofmon,file=fnddisplrx,status='replace')
212 close(sdofmon)
213
214 file="FND000000.DRY"
215
216 if(len_trim(monitor_file) .ne. 70) then
217 fnddisplry = monitor_file(1:len_trim(monitor_file)) // '/' // file
218 else
219 fnddisplry = file
220 endif
221
222 open(sdofmon,file=fnddisplry,status='replace')
223 close(sdofmon)
224
225 file="FND000000.DZX"
226
227 if(len_trim(monitor_file) .ne. 70) then
228 fnddisplzx = monitor_file(1:len_trim(monitor_file)) // '/' // file
229 else
230 fnddisplzx = file
231 endif
232
233 open(sdofmon,file=fnddisplzx,status='replace')
234 close(sdofmon)
235
236 file="FND000000.DZY"
237
238 if(len_trim(monitor_file) .ne. 70) then
239 fnddisplzy = monitor_file(1:len_trim(monitor_file)) // '/' // file
240 else
241 fnddisplzy = file
242 endif
243
244 open(sdofmon,file=fnddisplzy,status='replace')
245 close(sdofmon)
246 endif
247 endif
248
249 if (sdofout(2).eq.1) then ! acceleration
250
251 if(sfs.eq.0) then
252
253 !!! structure total acceleration
254
255 file="SDOF000000.AX"
256
257 if(len_trim(monitor_file) .ne. 70) then
258 sdofaccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
259 else
260 sdofaccx = file
261 endif
262
263 open(sdofmon,file=sdofaccx,status='replace')
264 close(sdofmon)
265
266 file="SDOF000000.AY"
267
268 if(len_trim(monitor_file) .ne. 70) then
269 sdofaccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
270 else
271 sdofaccy = file
272 endif
273
274 open(sdofmon,file=sdofaccy,status='replace')
275 close(sdofmon)
276
277 file="SDOF000000.AZ"
278
279 if(len_trim(monitor_file) .ne. 70) then
280 sdofaccz = monitor_file(1:len_trim(monitor_file)) // '/' // file
281 else
282 sdofaccz = file
283 endif
284
285 open(sdofmon,file=sdofaccz,status='replace')
286 close(sdofmon)
287
288 !!! soil acceleration
289
290 file="SDOF000000.GAX"
291
292 if(len_trim(monitor_file) .ne. 70) then
293 sdofgraccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
294 else
295 sdofgraccx = file
296 endif
297
298 open(sdofmon,file=sdofgraccx,status='replace')
299 close(sdofmon)
300
301 file="SDOF000000.GAY"
302
303 if(len_trim(monitor_file) .ne. 70) then
304 sdofgraccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
305 else
306 sdofgraccy = file
307 endif
308
309 open(sdofmon,file=sdofgraccy,status='replace')
310 close(sdofmon)
311
312 file="SDOF000000.GAZ"
313
314 if(len_trim(monitor_file) .ne. 70) then
315 sdofgraccz = monitor_file(1:len_trim(monitor_file)) // '/' // file
316 else
317 sdofgraccz = file
318 endif
319
320 open(sdofmon,file=sdofgraccz,status='replace')
321 close(sdofmon)
322 elseif(sfs.eq.1) then
323
324 !!! structure acceleration
325
326 file="STR000000.AX"
327
328 if(len_trim(monitor_file) .ne. 70) then
329 straccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
330 else
331 straccx = file
332 endif
333
334 open(sdofmon,file=straccx,status='replace')
335 close(sdofmon)
336
337 file="STR000000.AY"
338
339 if(len_trim(monitor_file) .ne. 70) then
340 straccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
341 else
342 straccy = file
343 endif
344
345 open(sdofmon,file=straccy,status='replace')
346 close(sdofmon)
347
348 !!! soil acceleration
349
350 file="GRD000000.AX"
351
352 if(len_trim(monitor_file) .ne. 70) then
353 grdaccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
354 else
355 grdaccx = file
356 endif
357
358 open(sdofmon,file=grdaccx,status='replace')
359 close(sdofmon)
360
361 file="GRD000000.AY"
362
363 if(len_trim(monitor_file) .ne. 70) then
364 grdaccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
365 else
366 grdaccy = file
367 endif
368
369 open(sdofmon,file=grdaccy,status='replace')
370 close(sdofmon)
371
372 file="GRD000000.AZ"
373
374 if(len_trim(monitor_file) .ne. 70) then
375 grdaccz = monitor_file(1:len_trim(monitor_file)) // '/' // file
376 else
377 grdaccz = file
378 endif
379
380 open(sdofmon,file=grdaccz,status='replace')
381 close(sdofmon)
382
383 !!! foundation acceleration
384
385 file="FND000000.AX"
386
387 if(len_trim(monitor_file) .ne. 70) then
388 fndaccx = monitor_file(1:len_trim(monitor_file)) // '/' // file
389 else
390 fndaccx = file
391 endif
392
393 open(sdofmon,file=fndaccx,status='replace')
394 close(sdofmon)
395
396 file="FND000000.AY"
397
398 if(len_trim(monitor_file) .ne. 70) then
399 fndaccy = monitor_file(1:len_trim(monitor_file)) // '/' // file
400 else
401 fndaccy = file
402 endif
403
404 open(sdofmon,file=fndaccy,status='replace')
405 close(sdofmon)
406
407 file="FND000000.ARX"
408
409 if(len_trim(monitor_file) .ne. 70) then
410 fndaccrx = monitor_file(1:len_trim(monitor_file)) // '/' // file
411 else
412 fndaccrx = file
413 endif
414
415 open(sdofmon,file=fndaccrx,status='replace')
416 close(sdofmon)
417
418 file="FND000000.ARY"
419
420 if(len_trim(monitor_file) .ne. 70) then
421 fndaccry = monitor_file(1:len_trim(monitor_file)) // '/' // file
422 else
423 fndaccry = file
424 endif
425
426 open(sdofmon,file=fndaccry,status='replace')
427 close(sdofmon)
428
429 file="FND000000.AZX"
430
431 if(len_trim(monitor_file) .ne. 70) then
432 fndacczx = monitor_file(1:len_trim(monitor_file)) // '/' // file
433 else
434 fndacczx = file
435 endif
436
437 open(sdofmon,file=fndacczx,status='replace')
438 close(sdofmon)
439
440 file="FND000000.AZY"
441
442 if(len_trim(monitor_file) .ne. 70) then
443 fndacczy = monitor_file(1:len_trim(monitor_file)) // '/' // file
444 else
445 fndacczy = file
446 endif
447
448 open(sdofmon,file=fndacczy,status='replace')
449 close(sdofmon)
450 endif
451 endif
452
453 if (sdofout(3).eq.1) then ! interaction force
454
455 if (sfs.eq.0) then
456
457 !!! structure shear force
458
459 file="SDOF000000.FX"
460
461 if(len_trim(monitor_file) .ne. 70) then
462 sdoffx = monitor_file(1:len_trim(monitor_file)) // '/' // file
463 else
464 sdoffx = file
465 endif
466
467 open(sdofmon,file=sdoffx,status='replace')
468 close(sdofmon)
469
470 file="SDOF000000.FY"
471
472 if(len_trim(monitor_file) .ne. 70) then
473 sdoffy = monitor_file(1:len_trim(monitor_file)) // '/' // file
474 else
475 sdoffy = file
476 endif
477
478 open(sdofmon,file=sdoffy,status='replace')
479 close(sdofmon)
480
481 file="SDOF000000.FZ"
482
483 if(len_trim(monitor_file) .ne. 70) then
484 sdoffz = monitor_file(1:len_trim(monitor_file)) // '/' // file
485 else
486 sdoffz = file
487 endif
488
489 open(sdofmon,file=sdoffz,status='replace')
490 close(sdofmon)
491 elseif(sfs.eq.1) then
492
493 !!! structure shear force
494
495 file="STR000000.FX"
496
497 if(len_trim(monitor_file) .ne. 70) then
498 strfx = monitor_file(1:len_trim(monitor_file)) // '/' // file
499 else
500 strfx = file
501 endif
502
503 open(sdofmon,file=strfx,status='replace')
504 close(sdofmon)
505
506 file="STR000000.FY"
507
508 if(len_trim(monitor_file) .ne. 70) then
509 strfy = monitor_file(1:len_trim(monitor_file)) // '/' // file
510 else
511 strfy = file
512 endif
513
514 open(sdofmon,file=strfy,status='replace')
515 close(sdofmon)
516
517 !!! foundation shear force
518
519 file="FND000000.FX"
520
521 if(len_trim(monitor_file) .ne. 70) then
522 fndfx = monitor_file(1:len_trim(monitor_file)) // '/' // file
523 else
524 fndfx = file
525 endif
526
527 open(sdofmon,file=fndfx,status='replace')
528 close(sdofmon)
529
530 file="FND000000.FY"
531
532 if(len_trim(monitor_file) .ne. 70) then
533 fndfy = monitor_file(1:len_trim(monitor_file)) // '/' // file
534 else
535 fndfy = file
536 endif
537
538 open(sdofmon,file=fndfy,status='replace')
539 close(sdofmon)
540
541 !!! interaction force
542
543 file="INT000000.FX"
544
545 if(len_trim(monitor_file) .ne. 70) then
546 intfx = monitor_file(1:len_trim(monitor_file)) // '/' // file
547 else
548 intfx = file
549 endif
550
551 open(sdofmon,file=intfx,status='replace')
552 close(sdofmon)
553
554 file="INT000000.FY"
555
556 if(len_trim(monitor_file) .ne. 70) then
557 intfy = monitor_file(1:len_trim(monitor_file)) // '/' // file
558 else
559 intfy = file
560 endif
561
562 open(sdofmon,file=intfy,status='replace')
563 close(sdofmon)
564
565 file="INT000000.FZ"
566
567 if(len_trim(monitor_file) .ne. 70) then
568 intfz = monitor_file(1:len_trim(monitor_file)) // '/' // file
569 else
570 intfz = file
571 endif
572
573 open(sdofmon,file=intfz,status='replace')
574 close(sdofmon)
575 endif
576 endif
577 endif
578
579 return
speed_par
Contains SPEED PARAMETERS used in (SPEED, READ_INPUT_FILES, MAKE_PARTION_AND_MPI_FILES,...
Definition MODULES.f90:207
speed_par::monitor_file
character *70 monitor_file
Definition MODULES.f90:243
speed_par::mpi_id
integer *4 mpi_id
Definition MODULES.f90:308
speed_sci
Contains parameters for MDOF.
Definition MODULES.f90:725
speed_sci::grddisplx
character *100 grddisplx
Definition MODULES.f90:794
speed_sci::sdofgraccz
character *100 sdofgraccz
SDOF base ground acceleration.
Definition MODULES.f90:790
speed_sci::fnddisplzy
character *100 fnddisplzy
4DOF foundation displacement
Definition MODULES.f90:795
speed_sci::grdaccz
character *100 grdaccz
4DOF ground acceleration
Definition MODULES.f90:797
speed_sci::fnddisplx
character *100 fnddisplx
Definition MODULES.f90:795
speed_sci::sdofgraccx
character *100 sdofgraccx
Definition MODULES.f90:790
speed_sci::sdofout
integer *4, dimension(3) sdofout
displ, acc, f_react
Definition MODULES.f90:776
speed_sci::fndacczy
character *100 fndacczy
4DOF foundation acceleration
Definition MODULES.f90:798
speed_sci::strfy
character *100 strfy
Definition MODULES.f90:799
speed_sci::sdofaccz
character *100 sdofaccz
SDOF total acceleration.
Definition MODULES.f90:789
speed_sci::straccx
character *100 straccx
Definition MODULES.f90:796
speed_sci::fndfy
character *100 fndfy
4DOF superstructure and foundation shear force
Definition MODULES.f90:799
speed_sci::fndaccy
character *100 fndaccy
Definition MODULES.f90:798
speed_sci::grdaccy
character *100 grdaccy
Definition MODULES.f90:797
speed_sci::sdofgrdisplx
character *100 sdofgrdisplx
Definition MODULES.f90:788
speed_sci::sdofdisplz
character *100 sdofdisplz
SDOF displacement.
Definition MODULES.f90:787
speed_sci::sdofgrdisply
character *100 sdofgrdisply
Definition MODULES.f90:788
speed_sci::n_bld
integer *4 n_bld
Definition MODULES.f90:773
speed_sci::grdaccx
character *100 grdaccx
Definition MODULES.f90:797
speed_sci::fnddisplrx
character *100 fnddisplrx
Definition MODULES.f90:795
speed_sci::straccy
character *100 straccy
4DOF structure acceleration
Definition MODULES.f90:796
speed_sci::sdofgrdisplz
character *100 sdofgrdisplz
SDOF base displacement.
Definition MODULES.f90:788
speed_sci::strdisply
character *100 strdisply
4DOF structure displacement
Definition MODULES.f90:793
speed_sci::sdofgraccy
character *100 sdofgraccy
Definition MODULES.f90:790
speed_sci::strdisplx
character *100 strdisplx
Definition MODULES.f90:793
speed_sci::fndfx
character *100 fndfx
Definition MODULES.f90:799
speed_sci::sdoffy
character *100 sdoffy
Definition MODULES.f90:791
speed_sci::sdofdisply
character *100 sdofdisply
Definition MODULES.f90:787
speed_sci::sdofaccy
character *100 sdofaccy
Definition MODULES.f90:789
speed_sci::sdoffx
character *100 sdoffx
Definition MODULES.f90:791
speed_sci::fnddisplry
character *100 fnddisplry
Definition MODULES.f90:795
speed_sci::sdoffz
character *100 sdoffz
SDOF reaction force.
Definition MODULES.f90:791
speed_sci::sdofaccx
character *100 sdofaccx
Definition MODULES.f90:789
speed_sci::fnddisplzx
character *100 fnddisplzx
Definition MODULES.f90:795
speed_sci::fndaccx
character *100 fndaccx
Definition MODULES.f90:798
speed_sci::strfx
character *100 strfx
Definition MODULES.f90:799
speed_sci::grddisplz
character *100 grddisplz
4DOF ground displacement
Definition MODULES.f90:794
speed_sci::intfy
character *100 intfy
Definition MODULES.f90:800
speed_sci::intfx
character *100 intfx
Definition MODULES.f90:800
speed_sci::fndaccry
character *100 fndaccry
Definition MODULES.f90:798
speed_sci::fnddisply
character *100 fnddisply
Definition MODULES.f90:795
speed_sci::fndaccrx
character *100 fndaccrx
Definition MODULES.f90:798
speed_sci::fndacczx
character *100 fndacczx
Definition MODULES.f90:798
speed_sci::sys
type(system), dimension(:), allocatable sys
SDOF system.
Definition MODULES.f90:771
speed_sci::sdofdisplx
character *100 sdofdisplx
Definition MODULES.f90:787
speed_sci::intfz
character *100 intfz
4DOF interaction forces
Definition MODULES.f90:800
speed_sci::grddisply
character *100 grddisply
Definition MODULES.f90:794
speed_timeloop
Contains a subset of SPEED paramters (used in TIME_LOOP)
Definition MODULES.f90:519

References speed_sci::fndaccrx, speed_sci::fndaccry, speed_sci::fndaccx, speed_sci::fndaccy, speed_sci::fndacczx, speed_sci::fndacczy, speed_sci::fnddisplrx, speed_sci::fnddisplry, speed_sci::fnddisplx, speed_sci::fnddisply, speed_sci::fnddisplzx, speed_sci::fnddisplzy, speed_sci::fndfx, speed_sci::fndfy, speed_sci::grdaccx, speed_sci::grdaccy, speed_sci::grdaccz, speed_sci::grddisplx, speed_sci::grddisply, speed_sci::grddisplz, speed_sci::intfx, speed_sci::intfy, speed_sci::intfz, speed_par::monitor_file, speed_par::mpi_id, speed_sci::n_bld, speed_sci::sdofaccx, speed_sci::sdofaccy, speed_sci::sdofaccz, speed_sci::sdofdisplx, speed_sci::sdofdisply, speed_sci::sdofdisplz, speed_sci::sdoffx, speed_sci::sdoffy, speed_sci::sdoffz, speed_sci::sdofgraccx, speed_sci::sdofgraccy, speed_sci::sdofgraccz, speed_sci::sdofgrdisplx, speed_sci::sdofgrdisply, speed_sci::sdofgrdisplz, speed_sci::sdofout, speed_sci::straccx, speed_sci::straccy, speed_sci::strdisplx, speed_sci::strdisply, speed_sci::strfx, speed_sci::strfy, and speed_sci::sys.