31 integer*4 :: sfs, temp, i, SDOFmon, jdof
33 character*32 :: fmt_spec
35 sdofmon=10*(mpi_id+1)+7
40 if(temp.eq.1) sfs = temp
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)
53 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempU1(
sys(i)%NDOF,1)
57 write(sdofmon,
"(A1)")
" "
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)
67 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempU1(
sys(i)%NDOF,2)
71 write(sdofmon,
"(A1)")
" "
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)
81 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempU1(
sys(i)%NDOF,3)
85 write(sdofmon,
"(A1)")
" "
91 write(sdofmon,
"(E16.7)",advance=
'NO') tt1tmp
93 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,1)
96 write(sdofmon,
"(A1)")
" "
100 write(sdofmon,
"(E16.7)",advance=
'NO') tt1tmp
102 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,2)
105 write(sdofmon,
"(A1)")
" "
109 write(sdofmon,
"(E16.7)",advance=
'NO') tt1tmp
111 write(sdofmon,
"(E16.7)",advance=
'NO')
sdofgd(i,3)
114 write(sdofmon,
"(A1)")
" "
116 elseif(sfs.eq.1)
then
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)
124 write(sdofmon,
"(A1)")
" "
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)
133 write(sdofmon,
"(A1)")
" "
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)
142 write(sdofmon,
"(A1)")
" "
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)
151 write(sdofmon,
"(A1)")
" "
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)
160 write(sdofmon,
"(A1)")
" "
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)
169 write(sdofmon,
"(A1)")
" "
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)
178 write(sdofmon,
"(A1)")
" "
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)
187 write(sdofmon,
"(A1)")
" "
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)
196 write(sdofmon,
"(A1)")
" "
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)
205 write(sdofmon,
"(A1)")
" "
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)
214 write(sdofmon,
"(A1)")
" "
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)
232 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempRA1(
sys(i)%NDOF,1)
235 write(sdofmon,
"(A1)")
" "
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)
245 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempRA1(
sys(i)%NDOF,2)
248 write(sdofmon,
"(A1)")
" "
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)
258 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%tempRA1(
sys(i)%NDOF,3)
261 write(sdofmon,
"(A1)")
" "
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)
271 write(sdofmon,
"(A1)")
" "
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)
279 write(sdofmon,
"(A1)")
" "
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)
287 write(sdofmon,
"(A1)")
" "
289 elseif(sfs.eq.1)
then
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)
297 write(sdofmon,
"(A1)")
" "
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)
306 write(sdofmon,
"(A1)")
" "
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)
315 write(sdofmon,
"(A1)")
" "
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)
324 write(sdofmon,
"(A1)")
" "
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)
333 write(sdofmon,
"(A1)")
" "
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)
342 write(sdofmon,
"(A1)")
" "
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)
351 write(sdofmon,
"(A1)")
" "
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)
360 write(sdofmon,
"(A1)")
" "
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)
369 write(sdofmon,
"(A1)")
" "
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)
378 write(sdofmon,
"(A1)")
" "
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)
387 write(sdofmon,
"(A1)")
" "
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)
402 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%IntForce(
sys(i)%NDOF,1)
405 write(sdofmon,
"(A1)")
" "
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)
415 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%IntForce(
sys(i)%NDOF,2)
418 write(sdofmon,
"(A1)")
" "
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)
428 write(sdofmon,
"(E16.7)",advance=
'NO')
sys(i)%IntForce(
sys(i)%NDOF,3)
431 write(sdofmon,
"(A1)")
" "
433 elseif(sfs.eq.1)
then
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)
441 write(sdofmon,
"(A1)")
" "
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)
450 write(sdofmon,
"(A1)")
" "
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)
459 write(sdofmon,
"(A1)")
" "
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)
468 write(sdofmon,
"(A1)")
" "
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)
476 write(sdofmon,
"(A1)")
" "
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)
484 write(sdofmon,
"(A1)")
" "
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)
492 write(sdofmon,
"(A1)")
" "