CalculiX  2.13
A Free Software Three-Dimensional Structural Finite Element Program
Functions/Subroutines
printout.f File Reference

Go to the source code of this file.

Functions/Subroutines

subroutine printout (set, nset, istartset, iendset, ialset, nprint, prlab, prset, v, t1, fn, ipkon, lakon, stx, eei, xstate, ener, mi, nstate_, ithermal, co, kon, qfx, ttime, trab, inotr, ntrans, orab, ielorien, norien, nk, ne, inum, filab, vold, ikin, ielmat, thicke, eme, islavsurf, mortar, time, ielprop, prop, veold, orname, nelemload, nload, sideload, xload)
 

Function/Subroutine Documentation

◆ printout()

subroutine printout ( character*81, dimension(*)  set,
integer  nset,
integer, dimension(*)  istartset,
integer, dimension(*)  iendset,
integer, dimension(*)  ialset,
integer  nprint,
character*6, dimension(*)  prlab,
character*81, dimension(*)  prset,
real*8, dimension(0:mi(2),*)  v,
real*8, dimension(*)  t1,
real*8, dimension(0:mi(2),*)  fn,
integer, dimension(*)  ipkon,
character*8, dimension(*)  lakon,
real*8, dimension(6,mi(1),*)  stx,
real*8, dimension(6,mi(1),*)  eei,
real*8, dimension(nstate_,mi(1),*)  xstate,
real*8, dimension(mi(1),*)  ener,
integer, dimension(*)  mi,
integer  nstate_,
integer  ithermal,
real*8, dimension(3,*)  co,
integer, dimension(*)  kon,
real*8, dimension(3,mi(1),*)  qfx,
real*8  ttime,
real*8, dimension(7,*)  trab,
integer, dimension(2,*)  inotr,
integer  ntrans,
real*8, dimension(7,*)  orab,
integer, dimension(mi(3),*)  ielorien,
integer  norien,
integer  nk,
integer  ne,
integer, dimension(*)  inum,
character*87, dimension(*)  filab,
real*8, dimension(0:mi(2),*)  vold,
integer  ikin,
integer, dimension(mi(3),*)  ielmat,
real*8, dimension(mi(3),*)  thicke,
real*8, dimension(6,mi(1),*)  eme,
integer, dimension(2,*)  islavsurf,
integer  mortar,
real*8  time,
integer, dimension(*)  ielprop,
real*8, dimension(*)  prop,
real*8, dimension(0:mi(2),*)  veold,
character*80, dimension(*)  orname,
integer, dimension(2,*)  nelemload,
integer  nload,
character*20, dimension(*)  sideload,
real*8, dimension(2,*)  xload 
)
25 !
26 ! stores results in the .dat file
27 !
28  implicit none
29 !
30  logical force
31 !
32  character*1 cflag
33  character*6 prlab(*)
34  character*8 lakon(*)
35  character*20 sideload(*)
36  character*80 noset,elset,orname(*)
37  character*81 set(*),prset(*)
38  character*87 filab(*)
39 !
40  integer nset,istartset(*),iendset(*),ialset(*),nprint,ipkon(*),
41  & mi(*),nstate_,ii,jj,iset,l,limit,node,ipos,ithermal,ielem,
42  & nelem,kon(*),inotr(2,*),ntrans,ielorien(mi(3),*),norien,nk,ne,
43  & inum(*),nfield,ikin,nodes,ne0,nope,mt,ielmat(mi(3),*),iface,
44  & jfaces,mortar,islavsurf(2,*),ielprop(*),nload,
45  & nelemload(2,*)
46 !
47  real*8 v(0:mi(2),*),t1(*),fn(0:mi(2),*),stx(6,mi(1),*),bhetot,
48  & eei(6,mi(1),*),xstate(nstate_,mi(1),*),ener(mi(1),*),energytot,
49  & volumetot,co(3,*),qfx(3,mi(1),*),rftot(0:3),ttime,time,
50  & trab(7,*),orab(7,*),vold(0:mi(2),*),enerkintot,thicke(mi(3),*),
51  & eme(6,mi(1),*),prop(*),veold(0:mi(2),*),xload(2,*)
52 !
53  mt=mi(2)+1
54 !
55 ! interpolation in the original nodes of 1d and 2d elements
56 !
57  do ii=1,nprint
58  if((prlab(ii)(1:4).eq.'U ').or.
59  & ((prlab(ii)(1:4).eq.'NT ').and.(ithermal.gt.1))) then
60  if(filab(1)(5:5).ne.' ') then
61  nfield=mt
62  cflag=' '
63  force=.false.
64  call map3dto1d2d(v,ipkon,inum,kon,lakon,nfield,nk,
65  & ne,cflag,co,vold,force,mi)
66  endif
67  exit
68  endif
69  enddo
70  do ii=1,nprint
71  if((prlab(ii)(1:4).eq.'NT ').and.(ithermal.le.1)) then
72  if(filab(2)(5:5).ne.' ') then
73  nfield=1
74  cflag=' '
75  force=.false.
76  call map3dto1d2d(t1,ipkon,inum,kon,lakon,nfield,nk,
77  & ne,cflag,co,vold,force,mi)
78  endif
79  exit
80  endif
81  enddo
82  do ii=1,nprint
83  if(prlab(ii)(1:2).eq.'RF') then
84  if(filab(1)(5:5).ne.' ') then
85  nfield=mt
86  cflag=' '
87  force=.true.
88  call map3dto1d2d(fn,ipkon,inum,kon,lakon,nfield,nk,
89  & ne,cflag,co,vold,force,mi)
90  endif
91  exit
92  endif
93  enddo
94 !
95  do ii=1,nprint
96 !
97 ! nodal values
98 !
99  if((prlab(ii)(1:4).eq.'U ').or.(prlab(ii)(1:4).eq.'NT ').or.
100  & (prlab(ii)(1:4).eq.'RF ').or.(prlab(ii)(1:4).eq.'RFL ').or.
101  & (prlab(ii)(1:4).eq.'PS ').or.(prlab(ii)(1:4).eq.'PN ').or.
102  & (prlab(ii)(1:4).eq.'MF ').or.(prlab(ii)(1:4).eq.'V ').or.
103  & (prlab(ii)(1:4).eq.'TS '))
104  & then
105 !
106  ipos=index(prset(ii),' ')
107  noset=' '
108  noset(1:ipos-1)=prset(ii)(1:ipos-1)
109 !
110 ! printing the header
111 !
112  if(prlab(ii)(1:4).eq.'U ') then
113  write(5,*)
114  if(mi(2).eq.3) then
115  write(5,100) noset(1:ipos-2),ttime+time
116  else
117  write(5,133) noset(1:ipos-2),ttime+time
118  endif
119  100 format(' displacements (vx,vy,vz) for set ',a,
120  & ' and time ',e14.7)
121  133 format(' displacements (v(i),i=1..ndof) for set ',a,
122  & ' and time ',e14.7)
123  write(5,*)
124  elseif((prlab(ii)(1:4).eq.'NT ').or.
125  & (prlab(ii)(1:4).eq.'TS ')) then
126  write(5,*)
127  write(5,101) noset(1:ipos-2),ttime+time
128  101 format(' temperatures for set ',a,' and time ',e14.7)
129  write(5,*)
130  elseif((prlab(ii)(1:5).eq.'RF ').or.
131  & (prlab(ii)(1:5).eq.'RF T')) then
132  write(5,*)
133  write(5,102) noset(1:ipos-2),ttime+time
134  102 format(' forces (fx,fy,fz) for set ',a,
135  & ' and time ',e14.7)
136  write(5,*)
137  elseif((prlab(ii)(1:5).eq.'RFL ').or.
138  & (prlab(ii)(1:5).eq.'RFL T')) then
139  write(5,*)
140  write(5,103) noset(1:ipos-2),ttime+time
141  103 format(' heat generation for set ',a,' and time ',e14.7)
142  write(5,*)
143  elseif(prlab(ii)(1:4).eq.'PS ') then
144  write(5,*)
145  write(5,115) noset(1:ipos-2),ttime+time
146  115 format(' static pressures for set ',a,' and time ',e14.7)
147  write(5,*)
148  elseif(prlab(ii)(1:4).eq.'PN ') then
149  write(5,*)
150  write(5,117) noset(1:ipos-2),ttime+time
151  117 format(' network pressures (total pressure for gases, sta
152  &tic pressure for liquids and fluid depth for channels) for set ',
153  &a,' and time ',e14.7)
154  write(5,*)
155  elseif(prlab(ii)(1:4).eq.'MF ') then
156  write(5,*)
157  write(5,118) noset(1:ipos-2),ttime+time
158  118 format(' mass flows for set ',a,' and time ',e14.7)
159  write(5,*)
160  elseif(prlab(ii)(1:4).eq.'V ') then
161  write(5,*)
162  write(5,119) noset(1:ipos-2),ttime+time
163  119 format(' velocities (vx,vy,vz) for set ',a,
164  & ' and time ',e14.7)
165  write(5,*)
166  endif
167 !
168 ! printing the data
169 !
170  do iset=1,nset
171  if(set(iset).eq.prset(ii)) exit
172  enddo
173  do jj=0,3
174  rftot(jj)=0.d0
175  enddo
176  do jj=istartset(iset),iendset(iset)
177  if(ialset(jj).lt.0) cycle
178  if(jj.eq.iendset(iset)) then
179  node=ialset(jj)
180  call printoutnode(prlab,v,t1,fn,ithermal,ii,node,
181  & rftot,trab,inotr,ntrans,co,mi,veold)
182  elseif(ialset(jj+1).gt.0) then
183  node=ialset(jj)
184  call printoutnode(prlab,v,t1,fn,ithermal,ii,node,
185  & rftot,trab,inotr,ntrans,co,mi,veold)
186  else
187  do node=ialset(jj-1)-ialset(jj+1),ialset(jj),
188  & -ialset(jj+1)
189  call printoutnode(prlab,v,t1,fn,ithermal,ii,node,
190  & rftot,trab,inotr,ntrans,co,mi,veold)
191  enddo
192  endif
193  enddo
194 !
195 ! writing total values to file
196 !
197  if((prlab(ii)(1:5).eq.'RF O').or.
198  & (prlab(ii)(1:5).eq.'RF T')) then
199  write(5,*)
200  write(5,104) noset(1:ipos-2),ttime+time
201  104 format(' total force (fx,fy,fz) for set ',a,
202  & ' and time ',e14.7)
203  write(5,*)
204  write(5,'(6x,1p,3(1x,e13.6))') rftot(1),rftot(2),rftot(3)
205  elseif((prlab(ii)(1:5).eq.'RFL O').or.
206  & (prlab(ii)(1:5).eq.'RFL T')) then
207  write(5,*)
208  write(5,105)noset(1:ipos-2),ttime+time
209  105 format(' total heat generation for set ',a,
210  & ' and time ',e14.7)
211  write(5,*)
212  write(5,'(6x,1p,1x,e13.6)') rftot(0)
213  endif
214 !
215 ! integration point values
216 !
217  elseif((prlab(ii)(1:4).eq.'S ').or.
218  & (prlab(ii)(1:4).eq.'E ').or.
219  & (prlab(ii)(1:4).eq.'ME ').or.
220  & (prlab(ii)(1:4).eq.'PEEQ').or.
221  & (prlab(ii)(1:4).eq.'ENER').or.
222  & (prlab(ii)(1:4).eq.'SDV ').or.
223  & (prlab(ii)(1:4).eq.'HFL ')) then
224 !
225  ipos=index(prset(ii),' ')
226  elset=' '
227  elset(1:ipos-1)=prset(ii)(1:ipos-1)
228 !
229  limit=1
230 !
231  do l=1,limit
232 !
233 ! printing the header
234 !
235  if(prlab(ii)(1:4).eq.'S ') then
236  write(5,*)
237  write(5,106) elset(1:ipos-2),ttime+time
238  106 format(' stresses (elem, integ.pnt.,sxx,syy,szz,sxy,sx
239  &z,syz) for set ',a,' and time ',e14.7)
240  write(5,*)
241  elseif(prlab(ii)(1:4).eq.'E ') then
242  write(5,*)
243  write(5,107) elset(1:ipos-2),ttime+time
244  107 format(' strains (elem, integ.pnt.,exx,eyy,ezz,exy,exz
245  &,eyz) forset ',a,' and time ',e14.7)
246  write(5,*)
247  elseif(prlab(ii)(1:4).eq.'PEEQ') then
248  write(5,*)
249  write(5,108) elset(1:ipos-2),ttime+time
250  108 format(' equivalent plastic strain (elem, integ.pnt.,p
251  &e)for set ',a,' and time ',e14.7)
252  write(5,*)
253  elseif(prlab(ii)(1:4).eq.'ENER') then
254  write(5,*)
255  write(5,109) elset(1:ipos-2),ttime+time
256  109 format(' internal energy density (elem, integ.pnt.,energy) for
257  &set ',a,' and time ',e14.7)
258  write(5,*)
259  elseif(prlab(ii)(1:4).eq.'SDV ') then
260  write(5,*)
261  write(5,111) elset(1:ipos-2),ttime+time
262  111 format
263  & (' internal state variables (elem, integ.pnt.,values) f
264  &or set ',a,' and time ',e14.7)
265  write(5,*)
266  elseif(prlab(ii)(1:4).eq.'HFL ') then
267  write(5,*)
268  write(5,112) elset(1:ipos-2),ttime+time
269  112 format(' heat flux (elem, integ.pnt.,qx,qy,qz) for set
270  & ',a,' and time ',e14.7)
271  write(5,*)
272  elseif(prlab(ii)(1:4).eq.'ME ') then
273  write(5,*)
274  write(5,130) elset(1:ipos-2),ttime+time
275  130 format(' mechanical strains (elem, integ.pnt.,exx,eyy,
276  &ezz,exy,exz,eyz) forset ',a,' and time ',e14.7)
277  write(5,*)
278  endif
279 !
280 ! printing the data
281 !
282  do iset=1,nset
283  if(set(iset).eq.prset(ii)) exit
284  enddo
285  do jj=istartset(iset),iendset(iset)
286  if(ialset(jj).lt.0) cycle
287  if(jj.eq.iendset(iset)) then
288  nelem=ialset(jj)
289  call printoutint(prlab,ipkon,lakon,stx,eei,xstate,
290  & ener,mi(1),nstate_,ii,nelem,qfx,
291  & orab,ielorien,norien,co,kon,ielmat,thicke,eme,
292  & ielprop,prop,nelem,ithermal,orname)
293  elseif(ialset(jj+1).gt.0) then
294  nelem=ialset(jj)
295  call printoutint(prlab,ipkon,lakon,stx,eei,xstate,
296  & ener,mi(1),nstate_,ii,nelem,qfx,orab,
297  & ielorien,norien,co,kon,ielmat,thicke,eme,
298  & ielprop,prop,nelem,ithermal,orname)
299  else
300  do nelem=ialset(jj-1)-ialset(jj+1),ialset(jj),
301  & -ialset(jj+1)
302  call printoutint(prlab,ipkon,lakon,stx,eei,
303  & xstate,ener,mi(1),nstate_,ii,nelem,
304  & qfx,orab,ielorien,norien,co,kon,ielmat,
305  & thicke,eme,ielprop,prop,nelem,ithermal,
306  & orname)
307  enddo
308  endif
309  enddo
310 !
311  enddo
312 !
313 ! whole element values
314 !
315  elseif((prlab(ii)(1:4).eq.'ELSE').or.
316  & (prlab(ii)(1:4).eq.'ELKE').or.
317  & (prlab(ii)(1:4).eq.'EVOL').or.
318  & (prlab(ii)(1:4).eq.'EBHE').or.
319  & (prlab(ii)(1:4).eq.'CSTR').or.
320  & (prlab(ii)(1:4).eq.'CDIS').or.
321  & (prlab(ii)(1:4).eq.'CNUM').or.
322  & (prlab(ii)(1:4).eq.'CELS')) then
323 !
324  ipos=index(prset(ii),' ')
325  elset=' '
326  elset(1:ipos-1)=prset(ii)(1:ipos-1)
327 !
328 ! printing the header
329 !
330  if((prlab(ii)(1:5).eq.'ELSE ').or.
331  & (prlab(ii)(1:5).eq.'ELSET')) then
332  write(5,*)
333  write(5,113) elset(1:ipos-2),ttime+time
334  113 format(' internal energy (element, energy) for set ',a,
335  & ' and time ',e14.7)
336  write(5,*)
337  elseif((prlab(ii)(1:5).eq.'ELKE ').or.
338  & (prlab(ii)(1:5).eq.'ELKET')) then
339  write(5,*)
340  write(5,110) elset(1:ipos-2),ttime+time
341  110 format(' kinetic energy (elem, energy) for set '
342  & ,a,' and time ',e14.7)
343  write(5,*)
344  elseif((prlab(ii)(1:5).eq.'EVOL ').or.
345  & (prlab(ii)(1:5).eq.'EVOLT')) then
346  write(5,*)
347  write(5,114) elset(1:ipos-2),ttime+time
348  114 format(' volume (element, volume) for set ',a,
349  & ' and time ',e14.7)
350  write(5,*)
351  elseif((prlab(ii)(1:5).eq.'EBHE ').or.
352  & (prlab(ii)(1:5).eq.'EBHET')) then
353  write(5,*)
354  write(5,131) elset(1:ipos-2),ttime+time
355  131 format(' body heating (element, volume) for set ',a,
356  & ' and time ',e14.7)
357  write(5,*)
358  elseif((prlab(ii)(1:5).eq.'CSTR ').or.
359  & (prlab(ii)(1:5).eq.'CSTRT')) then
360  write(5,*)
361  if(mortar.eq.0) then
362  write(5,122) ttime+time
363  122 format(' contact stress (slave node,press,'
364  & 'tang1,tang2) for all contact elements and time',
365  & e14.7)
366  elseif(mortar.eq.1) then
367  write(5,126) ttime+time
368  126 format(' contact stress (slave element+face,press,'
369  & 'tang1,tang2) for all contact elements and time',
370  & e14.7)
371  endif
372  write(5,*)
373  elseif((prlab(ii)(1:5).eq.'CDIS ').or.
374  & (prlab(ii)(1:5).eq.'CDIST')) then
375  write(5,*)
376  if(mortar.eq.0) then
377  write(5,123) ttime+time
378  123 format(' relative contact displacement (slave node,'
379  & 'normal,tang1,tang2) for all contact elements and '
380  & 'time',e14.7)
381  elseif(mortar.eq.1) then
382  write(5,127) ttime+time
383  127 format(
384  & ' relative contact displacement (slave element+face,'
385  & 'normal,tang1,tang2) for all contact elements and '
386  & 'time',e14.7)
387  endif
388  write(5,*)
389  elseif((prlab(ii)(1:5).eq.'CELS ').or.
390  & (prlab(ii)(1:5).eq.'CELST')) then
391  write(5,*)
392  if(mortar.eq.0) then
393  write(5,124) ttime+time
394  124 format(' contact print energy (slave node,energy) for'
395  & 'all contact elements and time',e14.7)
396  elseif(mortar.eq.1) then
397  write(5,128) ttime+time
398  128 format(
399  & ' contact print energy (slave element+face,energy)'
400  & 'for all contact elements and time',e14.7)
401  endif
402  endif
403 !
404 ! printing the data
405 !
406 
407  volumetot=0.d0
408  bhetot=0.d0
409  energytot=0.d0
410  enerkintot=0.d0
411 
412  if ((prlab(ii)(1:4).eq.'CSTR').or.
413  & (prlab(ii)(1:4).eq.'CDIS').or.
414  & (prlab(ii)(1:4).eq.'CNUM').or.
415  & (prlab(ii)(1:4).eq.'CELS')) then
416 !
417 ! ne0 is the number of the first contact element
418 !
419  do jj=ne,1,-1
420  if((lakon(jj)(2:2).ne.'S').or.
421  & (lakon(jj)(7:7).ne.'C')) then
422  ne0=jj+1
423  exit
424  endif
425  enddo
426 !
427  if(prlab(ii)(1:4).ne.'CNUM') then
428  if(mortar.eq.0) then
429  do nelem=ne0,ne
430  read(lakon(nelem)(8:8),'(i1)') nope
431  nope=nope+1
432  nodes=kon(ipkon(nelem)+nope)
433  call printoutelem(prlab,ipkon,lakon,kon,co,
434  & ener,mi(1),ii,nelem,energytot,volumetot,
435  & enerkintot,ne,stx,nodes,thicke,ielmat,
436  & ielem,iface,mortar,ielprop,prop,
437  & sideload,nload,nelemload,xload,bhetot)
438  enddo
439  elseif(mortar.eq.1) then
440  do nelem=ne0,ne
441  jfaces=
442  & islavsurf(1,kon(ipkon(nelem)+kon(ipkon(nelem))+2))
443  ielem=int(jfaces/10.d0)
444  iface=jfaces-10*ielem
445  call printoutelem(prlab,ipkon,lakon,kon,co,
446  & ener,mi(1),ii,nelem,energytot,volumetot,
447  & enerkintot,ne,stx,nodes,thicke,ielmat,
448  & ielem,iface,mortar,ielprop,prop,
449  & sideload,nload,nelemload,xload,bhetot)
450  enddo
451  endif
452  endif
453  else
454  do iset=1,nset
455  if(set(iset).eq.prset(ii)) exit
456  enddo
457  do jj=istartset(iset),iendset(iset)
458  if(ialset(jj).lt.0) cycle
459  if(jj.eq.iendset(iset)) then
460  nelem=ialset(jj)
461  call printoutelem(prlab,ipkon,lakon,kon,co,
462  & ener,mi(1),ii,nelem,energytot,volumetot,
463  & enerkintot,ne,stx,nodes,thicke,ielmat,
464  & ielem,iface,mortar,ielprop,prop,
465  & sideload,nload,nelemload,xload,bhetot)
466  elseif(ialset(jj+1).gt.0) then
467  nelem=ialset(jj)
468  call printoutelem(prlab,ipkon,lakon,kon,co,
469  & ener,mi(1),ii,nelem,energytot,volumetot,
470  & enerkintot,ne,stx,nodes,thicke,ielmat,
471  & ielem,iface,mortar,ielprop,prop,
472  & sideload,nload,nelemload,xload,bhetot)
473  else
474  do nelem=ialset(jj-1)-ialset(jj+1),ialset(jj),
475  & -ialset(jj+1)
476  call printoutelem(prlab,ipkon,lakon,kon,co,
477  & ener,mi(1),ii,nelem,energytot,volumetot,
478  & enerkintot,ne,stx,nodes,thicke,ielmat,
479  & ielem,iface,mortar,ielprop,prop,
480  & sideload,nload,nelemload,xload,bhetot)
481  enddo
482  endif
483  enddo
484  endif
485 !
486 ! writing total values to file
487 !
488  if((prlab(ii)(1:5).eq.'ELSEO').or.
489  & (prlab(ii)(1:5).eq.'ELSET')) then
490  write(5,*)
491  write(5,116) elset(1:ipos-2),ttime+time
492  116 format(' total internal energy for set ',a,' and time ',
493  & e14.7)
494  write(5,*)
495  write(5,'(6x,1p,1x,e13.6)') energytot
496  elseif((prlab(ii)(1:5).eq.'ELKEO').or.
497  & (prlab(ii)(1:5).eq.'ELKET')) then
498  write(5,*)
499  write(5,120) elset(1:ipos-2),ttime+time
500  120 format(' total kinetic energy for set ',a,' and time ',
501  & e14.7)
502  write(5,*)
503  write(5,'(6x,1p,1x,e13.6)') enerkintot
504  elseif((prlab(ii)(1:5).eq.'EVOLO').or.
505  & (prlab(ii)(1:5).eq.'EVOLT')) then
506  write(5,*)
507  write(5,121) elset(1:ipos-2),ttime+time
508  121 format(' total volume for set ',a,' and time ',e14.7)
509  write(5,*)
510  write(5,'(6x,1p,1x,e13.6)') volumetot
511  elseif((prlab(ii)(1:5).eq.'EBHEO').or.
512  & (prlab(ii)(1:5).eq.'EBHET')) then
513  write(5,*)
514  write(5,132) elset(1:ipos-2),ttime+time
515  132 format(' total body heating for set ',a,' and time ',
516  & e14.7)
517  write(5,*)
518  write(5,'(6x,1p,1x,e13.6)') bhetot
519  elseif((prlab(ii)(1:5).eq.'CELSO').or.
520  & (prlab(ii)(1:5).eq.'CELST')) then
521  write(5,*)
522  write(5,125) ttime+time
523  125 format(' total contact spring energy for time ',e14.7)
524  write(5,*)
525  write(5,'(6x,1p,1x,e13.6)') energytot
526  elseif(prlab(ii)(1:4).eq.'CNUM') then
527  write(5,*)
528  write(5,129) ttime+time
529  129 format
530  & (' total number of contact elements for time ',e14.7)
531  write(5,*)
532  write(5,'(6x,1p,1x,i10)') ne-ne0+1
533 !
534  endif
535  endif
536  enddo
537 !
538  return
printoutelem
subroutine printoutelem(prlab, ipkon, lakon, kon, co, ener, mi, ii, nelem, energytot, volumetot, enerkintot, ne, stx, nodes, thicke, ielmat, ielem, iface, mortar, ielprop, prop, sideload, nload, nelemload, xload, bhetot)
Definition: printoutelem.f:23
printoutint
subroutine printoutint(prlab, ipkon, lakon, stx, eei, xstate, ener, mi, nstate_, ii, nelem, qfx, orab, ielorien, norien, co, konf, ielmat, thicke, eme, ielprop, prop, nelel, ithermal, orname)
Definition: printoutint.f:22
printoutnode
subroutine printoutnode(prlab, v, t1, fn, ithermal, ii, node, rftot, trab, inotr, ntrans, co, mi, veold)
Definition: printoutnode.f:21
map3dto1d2d
subroutine map3dto1d2d(yn, ipkon, inum, kon, lakon, nfield, nk, ne, cflag, co, vold, force, mi)
Definition: map3dto1d2d.f:21
nodes
subroutine nodes(inpc, textpart, co, nk, nk_, set, istartset, iendset, ialset, nset, nset_, nalset, nalset_, istep, istat, n, iline, ipol, inl, ipoinp, inp, ipoinpc)
Definition: nodes.f:22
Hosted by OpenAircraft.com, (Michigan UAV, LLC)