OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
genoutp.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| genoutp ../engine/source/output/sty/genoutp.F
25!||--- called by ------------------------------------------------------
26!|| imp_buck ../engine/source/implicit/imp_buck.F
27!|| sortie_main ../engine/source/output/sortie_main.F
28!||--- calls -----------------------------------------------------
29!|| count_arsz_cs ../engine/source/output/sty/outp_c_s.F
30!|| count_arsz_ct ../engine/source/output/sty/outp_c_t.F
31!|| count_arsz_rs ../engine/source/output/sty/outp_r_s.F
32!|| count_arsz_rt ../engine/source/output/sty/outp_r_t.F
33!|| count_arsz_sps ../engine/source/output/sty/outp_sp_s.F
34!|| count_arsz_spt ../engine/source/output/sty/outp_sp_t.F
35!|| count_arsz_sptt ../engine/source/output/sty/outp_sp_t.F
36!|| count_arsz_ss ../engine/source/output/sty/outp_s_s.F
37!|| count_arsz_st ../engine/source/output/sty/outp_s_t.F
38!|| filelen_sys ../common_source/tools/input_output/write_routines.c
39!|| outp_arsz_cs ../engine/source/mpi/interfaces/spmd_outp.F
40!|| outp_arsz_ct ../engine/source/mpi/interfaces/spmd_outp.F
41!|| outp_arsz_rs ../engine/source/mpi/interfaces/spmd_outp.F
42!|| outp_arsz_rt ../engine/source/mpi/interfaces/spmd_outp.F
43!|| outp_arsz_sps ../engine/source/mpi/interfaces/spmd_outp.F
44!|| outp_arsz_spt ../engine/source/mpi/interfaces/spmd_outp.F
45!|| outp_arsz_sptt ../engine/source/mpi/interfaces/spmd_outp.F
46!|| outp_arsz_ss ../engine/source/mpi/interfaces/spmd_outp.F
47!|| outp_arsz_st ../engine/source/mpi/interfaces/spmd_outp.F
48!|| outp_c_s ../engine/source/output/sty/outp_c_s.F
49!|| outp_c_t ../engine/source/output/sty/outp_c_t.F
50!|| outp_c_tf ../engine/source/output/sty/outp_c_t.F
51!|| outp_mt ../engine/source/output/sty/outp_mt.F
52!|| outp_n_v ../engine/source/output/sty/outp_n_v.F
53!|| outp_n_v2 ../engine/source/output/sty/outp_n_v.F
54!|| outp_n_vc ../engine/source/output/sty/outp_n_v.F
55!|| outp_no ../engine/source/output/sty/outp_no.F
56!|| outp_r_s ../engine/source/output/sty/outp_r_s.F
57!|| outp_r_t ../engine/source/output/sty/outp_r_t.F
58!|| outp_s_s ../engine/source/output/sty/outp_s_s.F
59!|| outp_s_t ../engine/source/output/sty/outp_s_t.F
60!|| outp_s_tt ../engine/source/output/sty/outp_s_t.F
61!|| outp_sp_s ../engine/source/output/sty/outp_sp_s.F
62!|| outp_sp_t ../engine/source/output/sty/outp_sp_t.F
63!|| outp_sp_tt ../engine/source/output/sty/outp_sp_t.F
64!|| spmd_outpitab ../engine/source/mpi/interfaces/spmd_outp.F
65!||--- uses -----------------------------------------------------
66!|| drape_mod ../engine/share/modules/drape_mod.F
67!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
68!|| element_mod ../common_source/modules/elements/element_mod.F90
69!|| inoutfile_mod ../common_source/modules/inoutfile_mod.F
70!|| output_mod ../common_source/modules/output/output_mod.F90
71!|| stack_mod ../engine/share/modules/stack_mod.F
72!||====================================================================
73 SUBROUTINE genoutp(X ,D ,V ,A ,
74 2 IXS ,IXQ ,IXC ,IXT ,IXP ,
75 3 IXR ,IXTG ,IPARG ,PM ,IGEO ,
76 4 MS ,CONT ,ITAB ,PARTSAV,FINT ,
77 5 FEXT ,TANI ,EANI ,IPART ,
78 6 VR ,ELBUF_TAB ,DD_IAD,WEIGHT,
79 7 IPM ,KXSP ,SPBUF ,NODGLOB,LENG ,
80 8 FOPT ,NOM_OPT ,NPBY ,FNCONT ,FTCONT,
81 9 GEO ,THKE ,STACK ,DRAPE_SH4N, DRAPE_SH3N,
82 A DRAPEG,OUTPUT )
83C-----------------------------------------------
84C M o d u l e s
85C-----------------------------------------------
86 USE elbufdef_mod
87 USE stack_mod
88 USE inoutfile_mod
89 USE drape_mod
90 USE output_mod , ONLY : output_
91 use element_mod , only : nixs,nixq,nixc,nixp,nixr,nixt,nixtg
92C-----------------------------------------------
93C I m p l i c i t T y p e s
94C-----------------------------------------------
95#include "implicit_f.inc"
96C-----------------------------------------------
97C C o m m o n B l o c k s
98C-----------------------------------------------
99#include "com01_c.inc"
100#include "com04_c.inc"
101!#include "com06_c.inc"
102#include "com08_c.inc"
103#include "com09_c.inc"
104#include "param_c.inc"
105#include "units_c.inc"
106#include "scr03_c.inc"
107#include "scr11_c.inc"
108#include "scr14_c.inc"
109#include "scr16_c.inc"
110#include "scr17_c.inc"
111#include "chara_c.inc"
112#include "task_c.inc"
113#include "sphcom.inc"
114#include "filescount_c.inc"
115C-----------------------------------------------
116C D u m m y A r g u m e n t s
117C-----------------------------------------------
118 INTEGER IPARG(*),IPART(LIPART1,*),
119 . IXS(NIXS,*),IXQ(NIXQ,*),IXC(NIXC,*),IXTG(NIXTG,*),
120 . IXT(NIXT,*),IXP(NIXP,*),IXR(NIXR,*),ITAB(*),
121 . WEIGHT(*), DD_IAD(*),IPM(*),IGEO(*),KXSP(NISP,*),
122 . LENG,NODGLOB(*), NOM_OPT(*), NPBY(*)
123 my_real
124 . X(*), D(*), V(*), VR(*), A(*),
125 . PM(NPROPM,*), MS(*),
126 . CONT(*), PARTSAV(*),FINT(*),FEXT(*),TANI(6,*),EANI(*),
127 . spbuf(*), fopt(*), fncont(3,*), ftcont(3,*),geo(*),thke(*)
128 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP) :: ELBUF_TAB
129 TYPE (STACK_PLY) :: STACK
130 TYPE(DRAPE_) :: DRAPE_SH4N(NUMELC_DRAPE), DRAPE_SH3N(NUMELTG_DRAPE)
131 TYPE(DRAPEG_) :: DRAPEG
132 TYPE(output_),INTENT(INOUT) :: OUTPUT
133C-----------------------------------------------
134C L o c a l V a r i a b l e s
135C-----------------------------------------------
136 CHARACTER CHOUTP*4,FILNAM*100,T10*10,MES*40,VAR_CHAR*11,CHAR_JJ*2,
137 . titre*40
138 INTEGER FILEN,I,INNODA,J,JJ
139 INTEGER SIZLOC,SIZP0
140 INTEGER I161,I16A,I16B,I16C,I16D,I16E,I16F,I16G,I16H,I16I
141 INTEGER SIZWR_CT(3),SIZP0_CT(3),SIZLOC_CT(3),
142 . SIZWR_CS(2),SIZP0_CS(2),SIZLOC_CS(2),
143 . SIZWR_SS(2),SIZP0_SS(2),SIZLOC_SS(2),
144 . SIZWR_ST(3),SIZP0_ST(3),SIZLOC_ST(3)
145 INTEGER SIZW_P0
146 INTEGER, DIMENSION(:), ALLOCATABLE :: SIZLOC_WRITE
147 INTEGER OUTPSIZE
148 INTEGER, DIMENSION(:), ALLOCATABLE :: ITABG
149 INTEGER :: LEN_TMP_NAME
150 CHARACTER(len=2148) :: TMP_NAME
151C-----------------------------------------------
152 ALLOCATE(ITABG(LENG))
153 I161=1
154 i16a=i161+lnopt1*nrbody0
155 i16b=i16a+lnopt1*naccelm
156 i16c=i16b+lnopt1*nvolu
157 i16d=i16c+lnopt1*(ninter+nintsub)
158 i16e=i16d+lnopt1*nrwall
159 i16f=i16e !obsolete option removed
160 i16g=i16f+lnopt1*njoint
161 i16h=i16g+lnopt1*nsect
162 i16i=i16h+lnopt1*nlink
163C-----------------------------------------------
164C OPEN FILE
165C-----------------------------------------------
166 IF(ispmd == 0) THEN
167 IF(irootyy == 2)THEN
168 IF(ioutp>=1000)ioutp=1
169 WRITE(choutp,'(I3.3)')ioutp
170 filnam=rootnam(1:rootlen)//'Y'//choutp
171 filen = rootlen + 4
172 ELSE
173 IF(ioutp>=10000)ioutp=1
174 WRITE(choutp,'(I4.4)')ioutp
175 filnam=rootnam(1:rootlen)//'_'//choutp//'.sty'
176 filen = rootlen + 9
177 ENDIF
178 len_tmp_name = outfile_name_len + filen
179 tmp_name=outfile_name(1:outfile_name_len)//filnam(1:filen)
180 OPEN(unit=iugeo,file=tmp_name(1:len_tmp_name),
181 . access='SEQUENTIAL',
182 . form='FORMATTED',status='UNKNOWN')
183C-----------------------------------------------
184C WRITE CONTROL
185C-----------------------------------------------
186 WRITE(iugeo,'(2A)')'#RADIOSS OUTPUT FILE V21 ',
187 . filnam(1:filen)
188C-----------------------------------------------
189C GLOBAL SCALARS
190C-----------------------------------------------
191 WRITE(iugeo,'(A)')'/GLOBAL '
192 WRITE(iugeo,'(A)')' '
193 WRITE(iugeo,'(A)')'#FORMAT: (1P5E16.9)'
194 WRITE(iugeo,'(2A)')'# TIME INTERNAL_ENERGY',
195 . ' KINETIC_ENERGY ROT_KINE_ENERGY EXTE_FORCE_WORK'
196 WRITE(iugeo,'(1P5E16.9)') tt,enint,encin,enrot,output%TH%WFEXT
197 ENDIF
198C-----------------------------------------------
199C MATER SCALARS
200C-----------------------------------------------
201 CALL outp_mt(pm,npart,partsav,ipart,ipm)
202C-----------------------------------------------
203C NODAL VECT
204C-----------------------------------------------
205 IF (nspmd > 1)
206 . CALL spmd_outpitab(itab,weight,nodglob,itabg)
207 CALL outp_n_v('COORDINATE',
208 . ' Coordinates ',
209 . 20,x,numnod,itab,
210 . itabg,leng,nodglob,weight)
211
212 IF(outp_v(1) == 1) CALL outp_n_v('VELOCITY ',
213 . ' Velocity ',
214 . 20,v,numnod,itab,
215 . itabg,leng,nodglob,weight)
216
217 IF(outp_v(2) == 1) CALL outp_n_v('DISPLACEME',
218 . ' Displacement ',
219 . 20,d,numnod,itab,
220 . itabg,leng,nodglob,weight)
221
222 IF(outp_v(3) == 1) CALL outp_n_v('ACCELERATI',
223 . ' Acceleration ',
224 . 20,a,numnod,itab,
225 . itabg,leng,nodglob,weight)
226
227 IF(outp_v(4) == 1) CALL outp_n_vc('CONT_FORCE',
228 . ' Contact Forces ',
229 . 20,cont,numnod,itab,
230 . itabg,leng,nodglob,weight)
231
232 IF(outp_v(5) == 1) CALL outp_n_v('INTE_FORCE',
233 . ' Internal Forces ',
234 . 20,fint,numnod,itab,
235 . itabg,leng,nodglob,weight)
236
237 IF(outp_v(6) == 1) CALL outp_n_v('EXTE_FORCE',
238 . ' External Forces ',
239 . 20,fext,numnod,itab,
240 . itabg,leng,nodglob,weight)
241
242 IF(outp_v(9) == 1 .AND. (idrot==1 .OR. isecut>0 .OR. iisrot>0 .OR. impose_dr>0 ) .AND. iroddl/=0)
243 . CALL outp_n_v('VELO_ROT ',
244 . ' Rotational Velocity ',
245 . 20,vr,numnod,itab,
246 . itabg,leng,nodglob,weight)
247
248 IF(outp_v(12) == 1) THEN
249 CALL outp_n_vc('PCONT_NORM',
250 . ' Contact Pressure from Normal Contact Forces ',
251 . 50,fncont,numnod,itab,
252 . itabg,leng,nodglob,weight)
253 CALL outp_n_vc('PCONT_TANG',
254 . ' Contact Pressure from Tangential Contact Forces ',
255 . 50,ftcont,numnod,itab,
256 . itabg,leng,nodglob,weight)
257 END IF
258C-----------------------------------------------
259C NODAL SCALAR
260C JBM LE 23/9/96
261C
262C-----------------------------------------------
263 innoda=0
264 IF(outp_n(1) == 1) CALL outp_no('NODAL_DT ',
265 . ' Nodal Time Step ',
266 . 20,output%DATA%SCAL_DT,numnod,itab,innoda,
267 . itabg,leng,nodglob,weight)
268
269 IF(outp_n(1)+anim_n(1)>=1)
270 . innoda=innoda+numnod
271
272 IF(outp_n(2) == 1) CALL outp_no('NOD_ADMAS ',
273 . ' Added Mass ',
274 . 20,output%DATA%SCAL_DMAS,numnod,itab,0,
275 . itabg,leng,nodglob,weight)
276
277 IF(outp_n(2)+anim_n(2)>=1)
278 . innoda=innoda+numnod
279
280 IF(outp_n(3) == 1) CALL outp_no('NOD_ADINER',
281 . ' Added Inertia ',
282 . 20,output%DATA%SCAL_DINER,numnod,itab,0,
283 . itabg,leng,nodglob,weight)
284
285C IF(OUTP_N(3)+ANIM_N(12)>=1)
286C . INNODA=INNODA+NUMNOD
287
288C-----------------------------------------------
289C SOLID SCALAR
290C-----------------------------------------------
291 IF (nspmd > 1) THEN
292 CALL outp_arsz_ss(iparg,dd_iad,ipm,ixs,sizp0_ss,sizloc_ss,sizwr_ss)
293 ELSE
294 ALLOCATE(sizloc_write(2*nspgroup+2))
295 CALL count_arsz_ss(iparg,dd_iad,ipm,ixs,sizloc_ss,sizloc_write)
296 sizp0_ss(:)=sizloc_ss(:)
297 sizwr_ss(1:2) = -1
298 DO i=1,nspgroup
299 sizwr_ss(1) = max(sizwr_ss(1),sizloc_write(i))
300 sizwr_ss(2) = max(sizwr_ss(2),sizloc_write(nspgroup+i))
301 ENDDO
302 sizwr_ss(1) = sizwr_ss(1)+6
303 sizwr_ss(2) = sizwr_ss(2)+6
304 DEALLOCATE(sizloc_write)
305 ENDIF
306 IF(outp_ss(1) == 1) CALL outp_s_s( 1,'OFF ',
307 . ' Delete flag (1.0:On 0.0:Off) '
308 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
309
310 IF(outp_ss(2) == 1) CALL outp_s_s(10,'EPSP ',
311 . ' Plastic Strain '
312 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
313
314 IF(outp_ss(3) == 1) CALL outp_s_s( 3,'ENER ',
315 . ' Internal Specific Energy '
316 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
317
318 IF(outp_ss(4) == 1) CALL outp_s_s( 4,'DENS ',
319 . ' Density '
320 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
321
322 IF(outp_ss(5) == 1) CALL outp_s_s(11,'TEMP ',
323 . ' Temperature '
324 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
325
326 IF(outp_ss(6) == 1) CALL outp_s_s( 2,'PRES ',
327 . ' Pressure '
328 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
329
330 IF(outp_ss(7) == 1) CALL outp_s_s(-2,'VONM ',
331 . ' Von Mises Stress '
332 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
333
334 IF(outp_ss(25) == 1) CALL outp_s_s(25,'HOUR ',
335 . ' Hourglass Specific Energy '
336 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
337
338 IF(outp_ss(20) == 1) CALL outp_s_s(20,'USER1 ',
339 . ' User Variable 1 '
340 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
341
342 IF(outp_ss(21) == 1) CALL outp_s_s(21,'USER2 ',
343 . ' User Variable 2 '
344 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
345
346 IF(outp_ss(22) == 1) CALL outp_s_s(22,'USER3 ',
347 . ' User Variable 3 '
348 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
349
350 IF(outp_ss(23) == 1) CALL outp_s_s(23,'USER4 ',
351 . ' User Variable 4 '
352 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
353
354 IF(outp_ss(24) == 1) CALL outp_s_s(24,'USER5 ',
355 . ' User Variable 5 '
356 . ,elbuf_tab,iparg,eani,ixs, ipm,dd_iad,sizloc_ss(1),sizp0_ss(1),sizwr_ss(1))
357
358 IF(outp_ss(26) == 1)CALL outp_s_s(26,'USERS ',
359 . ' all user variable '
360 . ,ELBUF_TAB,IPARG,EANI,IXS, IPM,DD_IAD,SIZLOC_SS(2),SIZP0_SS(2),SIZWR_SS(2))
361C-----------------------------------------------
362C element equivalent stress criteria
363C-----------------------------------------------
364 IF(OUTP_SS(27) == 1)CALL OUTP_S_S(26,'eq stress ',
365 . ' equivalent stress '
366 . ,ELBUF_TAB,IPARG,EANI,IXS, IPM,DD_IAD,SIZLOC_SS(2),SIZP0_SS(2),SIZWR_SS(2))
367C-----------------------------------------------
368C SHELL SCALAR
369C-----------------------------------------------
370
371 IF (NSPMD > 1) THEN
372 CALL OUTP_ARSZ_CS(IPARG,IXC,IXTG,IGEO,IPM,DD_IAD,
373 . SIZP0_CS,SIZLOC_CS,SIZWR_CS)
374 ELSE
375 ALLOCATE(SIZLOC_WRITE(2*NSPGROUP+2))
376 CALL COUNT_ARSZ_CS(IPARG,IXC,IXTG,IGEO,IPM,DD_IAD,
377 . SIZLOC_CS,SIZLOC_WRITE)
378 SIZP0_CS(1:2)=SIZLOC_CS(1:2)+6
379 SIZWR_CS(1:2) = -1
380 DO I=1,NSPGROUP
381 SIZWR_CS(1) = MAX(SIZWR_CS(1),SIZLOC_WRITE(I))
382 SIZWR_CS(2) = MAX(SIZWR_CS(2),SIZLOC_WRITE(NSPGROUP+I))
383 ENDDO
384 SIZWR_CS(1) = SIZWR_CS(1)+6
385 SIZWR_CS(2) = SIZWR_CS(2)+6
386 DEALLOCATE(SIZLOC_WRITE)
387 ENDIF
388
389 IF (OUTP_CS(1) == 1) CALL OUTP_C_S( 6,'off ',
390 . ' delete flag(1.0:on 0.0:off) '
391 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,
392 . IXC ,IXTG ,DD_IAD,SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1))
393
394 IF (OUTP_CS(2) == 1) THEN
395 CALL OUTP_C_S(15,'epsp ',
396 . ' plastic strain '
397 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,SIZLOC_CS(1),
398 . SIZP0_CS(1),THKE,SIZWR_CS(1) )
399
400 CALL OUTP_C_S(-15,'epsp_max ',
401 . ' maximum plastic strain '
402 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,SIZLOC_CS(1),
403 . SIZP0_CS(1),THKE,SIZWR_CS(1) )
404
405 ENDIF
406
407 IF (OUTP_CS(3) == 1) CALL OUTP_C_S( 5,'ener ',
408 . ' internal energy '
409 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,
410 . IXC ,IXTG ,DD_IAD,SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
411
412 IF (OUTP_CS(4) == 1) CALL OUTP_C_S( 3,'thick ',
413 . ' thickness '
414 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,
415 . IXC ,IXTG ,DD_IAD,SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
416
417 IF (OUTP_CS(7) == 1) CALL OUTP_C_S( 1,'vonm ',
418 . ' membrane von mises stress '
419 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,
420 . IXC ,IXTG ,DD_IAD,SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1))
421
422 IF (OUTP_CS(25) == 1) CALL OUTP_C_S(25,'hour ',
423 . ' hourglass energy '
424 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,
425 . IXC ,IXTG ,DD_IAD,SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1))
426
427 IF (OUTP_CS(20) == 1) THEN
428 CALL OUTP_C_S(20,'user1 ',
429 . ' user variable 1 '
430 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
431 . SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
432
433 CALL OUTP_C_S(-20,'user1_max ',
434 . ' Maximum User Variable 1 '
435 . ,elbuf_tab,iparg,eani ,ipm ,igeo,ixc ,ixtg ,dd_iad,
436 . sizloc_cs(1),sizp0_cs(1),thke,sizwr_cs(1))
437
438 ENDIF
439
440 IF (outp_cs(21) == 1) THEN
441 CALL outp_c_s(21,'USER2 ',
442 . ' User Variable 2 '
443 . ,elbuf_tab,iparg,eani ,ipm ,igeo,ixc ,ixtg ,dd_iad,
444 . sizloc_cs(1),sizp0_cs(1),thke,sizwr_cs(1) )
445
446 CALL outp_c_s(-21,'user2_max ',
447 . ' maximum user variable 2 '
448 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
449 . SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
450
451 ENDIF
452
453 IF (OUTP_CS(22) == 1) THEN
454 CALL OUTP_C_S(22,'user3 ',
455 . ' user variable 3 '
456 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
457 . SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1))
458
459 CALL OUTP_C_S(-22,'user3_max ',
460 . ' maximum user variable 3 '
461 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
462 . SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
463
464 ENDIF
465
466 IF (OUTP_CS(23) == 1) THEN
467 CALL OUTP_C_S(23,'user4 ',
468 . ' user variable 4 '
469 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
470 . SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
471
472 CALL OUTP_C_S(-23,'user4_max ',
473 . ' maximum user variable 4 '
474 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
475 . SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
476
477 ENDIF
478
479 IF (OUTP_CS(24) == 1) THEN
480 CALL OUTP_C_S(24,'user5 ',
481 . ' user variable 5 '
482 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
483 . SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
484
485 CALL OUTP_C_S(-24,'user5_max ',
486 . ' maximum user variable 5 '
487 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
488 . SIZLOC_CS(1),SIZP0_CS(1),THKE,SIZWR_CS(1) )
489
490 ENDIF
491
492 IF (OUTP_CS(26) == 1)CALL OUTP_C_S(26,'users ',
493 . ' all user variables '
494 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
495 . SIZLOC_CS(2),SIZP0_CS(2),THKE,SIZWR_CS(2) )
496 DO J = 27, 86
497 JJ = J - 26
498 write(CHAR_JJ,'(i2.2)')JJ
499 VAR_CHAR = 'user'//CHAR_JJ//'_full'
500 TITRE = ' user variable'//CHAR_JJ//' '
501 IF(OUTP_CS(J) == 1)CALL OUTP_C_S(J,VAR_CHAR,
502 . TITRE ,
503 . ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,IXC ,IXTG ,DD_IAD,
504 . SIZLOC_CS(2),SIZP0_CS(2),THKE,SIZWR_CS(2) )
505 ENDDO
506C-----------------------------------------------
507C element equivalent stress criteria
508C-----------------------------------------------
509 IF (OUTP_CS(87) == 1) CALL OUTP_C_S(87,'eq stress ',
510 . ' membrane equivalent stress '
511 . ,ELBUF_TAB,IPARG,EANI ,IPM ,IGEO,
512 . IXC ,IXTG ,DD_IAD,SIZLOC_CS(2),SIZP0_CS(2),THKE,SIZWR_CS(2) )
513C-----------------------------------------------
514C SPH SCALAR
515C-----------------------------------------------
516 IF (NSPMD > 1) THEN
517 CALL OUTP_ARSZ_SPS(IPARG,DD_IAD,SIZLOC,SIZP0,SIZW_P0)
518 ELSE
519 ALLOCATE(SIZLOC_WRITE(NSPGROUP+1))
520 CALL COUNT_ARSZ_SPS(IPARG,DD_IAD,SIZLOC,SIZLOC_WRITE)
521 SIZP0 = SIZLOC
522 SIZW_P0 = -1
523 DO I=1,NSPGROUP
524 SIZW_P0 = MAX(SIZW_P0,SIZLOC_WRITE(I))
525 ENDDO
526 SIZW_P0 = SIZW_P0 + 6
527 DEALLOCATE(SIZLOC_WRITE)
528 ENDIF
529
530 IF(OUTP_SPS(1) == 1) CALL OUTP_SP_S( 1,'off ',
531 . ' delete flag(1.0:on 0.0:off) '
532 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
533
534 IF(OUTP_SPS(2) == 1) CALL OUTP_SP_S(10,'epsp ',
535 . ' plastic strain '
536 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
537
538 IF(OUTP_SPS(3) == 1) CALL OUTP_SP_S( 3,'ener ',
539 . ' internal specific energy '
540 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
541
542 IF(OUTP_SPS(4) == 1) CALL OUTP_SP_S( 4,'dens ',
543 . ' density '
544 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
545
546 IF(OUTP_SPS(5) == 1) CALL OUTP_SP_S(11,'temp ',
547 . ' temperature '
548 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
549
550 IF(OUTP_SPS(6) == 1) CALL OUTP_SP_S( 2,'pres ',
551 . ' pressure '
552 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
553
554 IF(OUTP_SPS(7) == 1) CALL OUTP_SP_S(-2,'vonm ',
555 . ' von mises stress '
556 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
557
558 IF(OUTP_SPS(25) == 1) CALL OUTP_SP_S(25,'slen ',
559 . ' smoothing length '
560 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
561
562 IF(OUTP_SPS(20) == 1) CALL OUTP_SP_S(20,'user1 ',
563 . ' user variable 1 '
564 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
565
566 IF(OUTP_SPS(21) == 1) CALL OUTP_SP_S(21,'user2 ',
567 . ' user variable 2 '
568 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
569
570 IF(OUTP_SPS(22) == 1) CALL OUTP_SP_S(22,'user3 ',
571 . ' user variable 3 '
572 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
573
574 IF(OUTP_SPS(23) == 1) CALL OUTP_SP_S(23,'user4 ',
575 . ' user variable 4 '
576 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
577
578 IF(OUTP_SPS(24) == 1) CALL OUTP_SP_S(24,'user5 ',
579 . ' user variable 5 '
580 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
581C-----------------------------------------------
582C element equivalent stress criteria
583C-----------------------------------------------
584 IF(OUTP_SPS(26) == 1) CALL OUTP_SP_S(24,'eq stress ',
585 . ' equivalent stress '
586 . ,ELBUF_TAB,IPARG,EANI,DD_IAD,KXSP,IPM,SPBUF,SIZLOC,SIZP0,SIZW_P0)
587C-----------------------------------------------
588C SOLID TENSOR(3D)
589C-----------------------------------------------
590 IF (NSPMD > 1) THEN
591 CALL OUTP_ARSZ_ST(IPARG,DD_IAD,SIZLOC_ST,SIZWR_ST,SIZP0_ST)
592 ELSE
593 SIZP0_ST = -1
594 SIZLOC_ST = -1
595 ALLOCATE( SIZLOC_WRITE(3*NSPGROUP+3) )
596 CALL COUNT_ARSZ_ST(IPARG,DD_IAD,SIZLOC_ST,SIZLOC_WRITE)
597 SIZP0_ST(1:3) = SIZLOC_WRITE(3*NSPGROUP+1:3*NSPGROUP+3)+8
598 SIZWR_ST(1:3) = -1
599 DO I=1,NSPGROUP
600 SIZWR_ST(1) = MAX(SIZWR_ST(1),SIZLOC_WRITE(I))
601 SIZWR_ST(2) = MAX(SIZWR_ST(2),SIZLOC_WRITE(NSPGROUP+I))
602 SIZWR_ST(3) = MAX(SIZWR_ST(3),SIZLOC_WRITE(2*NSPGROUP+I))
603 ENDDO
604 SIZWR_ST(1:3)= SIZWR_ST(1:3)+8
605 DEALLOCATE( SIZLOC_WRITE )
606 ENDIF
607
608 IF(OUTP_ST(1) == 1) CALL OUTP_S_T( 2,'stress ',
609 . ' stress ',
610 . ELBUF_TAB,IPARG,DD_IAD,SIZLOC_ST(1),SIZP0_ST(1),SIZWR_ST(1))
611
612 IF(OUTP_ST(2) == 1) CALL OUTP_S_TT( 2,'str_ful ',
613 . ' full stress tensor + plastic strain ',
614 . ELBUF_TAB,IPARG,DD_IAD,IPM, IXS,
615 . SIZLOC_ST(2),SIZP0_ST(2),SIZWR_ST(2))
616
617 IF(OUTP_ST(3) == 1) CALL OUTP_S_TT( 3,'strain_ful ',
618 . ' full strain tensor ',
619 . ELBUF_TAB,IPARG,DD_IAD,IPM, IXS,
620 . SIZLOC_ST(3),SIZP0_ST(3),SIZWR_ST(3))
621C-----------------------------------------------
622C SHELL TENSOR(2D)
623C-----------------------------------------------
624 IF (NSPMD > 1) THEN
625 CALL OUTP_ARSZ_CT(IPARG,DD_IAD,SIZLOC_CT,SIZP0_CT,SIZWR_CT,ELBUF_TAB)
626 ELSE
627 ALLOCATE(SIZLOC_WRITE(3*NSPGROUP+3))
628 CALL COUNT_ARSZ_CT(IPARG,DD_IAD,SIZLOC_CT,SIZLOC_WRITE,ELBUF_TAB)
629 SIZP0_CT(1)=SIZLOC_CT(1)
630 SIZP0_CT(2:3)=SIZLOC_CT(2:3)
631 SIZWR_CT(1:3)=-1
632 DO I=1,NSPGROUP
633 SIZWR_CT(1) = MAX(SIZWR_CT(1),SIZLOC_WRITE(I))
634 SIZWR_CT(2) = MAX(SIZWR_CT(2),SIZLOC_WRITE(NSPGROUP+I))
635 SIZWR_CT(3) = MAX(SIZWR_CT(3),SIZLOC_WRITE(2*NSPGROUP+I))
636 ENDDO
637 SIZWR_CT(1) = SIZWR_CT(1) + 6
638 DEALLOCATE(SIZLOC_WRITE)
639 ENDIF
640
641 IF(OUTP_CT(1) == 1) CALL OUTP_C_T( 1,'stress_mem',
642 . ' membrane stress ',
643 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
644
645 IF(OUTP_CT(2) == 1) CALL OUTP_C_T( 2,'stress_ben',
646 . ' bending stress ',
647 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
648
649 IF(OUTP_CT(3) == 1) CALL OUTP_C_T( 3,'stress_upp',
650 . ' upper surface stress ',
651 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
652
653 IF(OUTP_CT(4) == 1) CALL OUTP_C_T( 4,'stress_low',
654 . ' lower surface stress ',
655 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
656
657 IF(OUTP_CT(5) == 1) CALL OUTP_C_T( 5,'strain_mem',
658 . ' membrane strain ',
659 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
660
661 IF(OUTP_CT(6) == 1) CALL OUTP_C_T( 6,'strain_ben',
662 . ' bending strain ',
663 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
664
665 IF(OUTP_CT(7) == 1) CALL OUTP_C_T( 7,'strain_upp',
666 . ' upper surface strain ',
667 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
668
669 IF(OUTP_CT(8) == 1) CALL OUTP_C_T( 8,'strain_low',
670 . ' lower surface strain ',
671 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
672
673 IF(OUTP_CT(91) == 1) CALL OUTP_C_T(91,'epsdot_mem',
674 . ' membrane strain rate ',
675 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
676
677 IF(OUTP_CT(92) == 1) CALL OUTP_C_T(92,'epsdot_ben',
678 . ' bending strain rate ',
679 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1))
680
681 IF(OUTP_CT(93) == 1) CALL OUTP_C_T(93,'epsdot_upp',
682 . ' upper surface strain rate ',
683 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1))
684
685 IF(OUTP_CT(94) == 1) CALL OUTP_C_T(94,'epsdot_low',
686 . ' lower surface strain rate ',
687 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
688
689 DO I=1,30
690 IF(OUTP_CT(10+I) == 1)THEN
691 WRITE(MES,'(a,i2,a)')
692 . ' stress(layer',I,') '
693 WRITE(T10,'(a,i2)')'stress_l',I
694 CALL OUTP_C_T( 10+I,T10 ,MES,
695 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
696 ENDIF
697
698 IF(OUTP_CT(50+I) == 1)THEN
699 WRITE(MES,'(a,i2,a)')
700 . ' strain(layer',I,') '
701 WRITE(T10,'(a,i2)')'strain_l',I
702 CALL OUTP_C_T( 50+I,T10 ,MES,
703 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
704 ENDIF
705
706 IF(OUTP_CT(100+I) == 1)THEN
707 WRITE(MES,'(a,i2,a)')
708 . ' epsdot(layer',I,') '
709 WRITE(T10,'(a,i2)')'epsdot_l',I
710 CALL OUTP_C_T( 100+I,T10 ,MES,
711 . ELBUF_TAB,IPARG,TANI,DD_IAD,SIZLOC_CT(1),SIZP0_CT(1),SIZWR_CT(1) )
712 ENDIF
713 ENDDO
714 IF(OUTP_CT(95) == 1) CALL OUTP_C_TF(95,'stress_ful',
715 . ' full stress tensor + plastic strain ',
716 . ELBUF_TAB,IPARG ,DD_IAD,SIZLOC_CT(2),SIZP0_CT(2),THKE,GEO,
717 . IGEO ,STACK ,DRAPE_SH4N, DRAPE_SH3N,IXC ,IXTG,SIZWR_CT(2),
718 . DRAPEG )
719
720
721 IF(OUTP_CT(96) == 1) CALL OUTP_C_TF(96,'strain_ful',
722 . ' full strain tensor + plastic strain ',
723 . ELBUF_TAB,IPARG ,DD_IAD,SIZLOC_CT(3),SIZP0_CT(3),THKE,GEO,
724 . IGEO ,STACK ,DRAPE_SH4N, DRAPE_SH3N,IXC ,IXTG,SIZWR_CT(3),
725 . DRAPEG )
726
727C-----------------------------------------------
728C SPH TENSOR(3D)
729C-----------------------------------------------
730 IF (NSPMD > 1) THEN
731 CALL OUTP_ARSZ_SPT(IPARG,DD_IAD,SIZLOC,SIZP0,SIZW_P0)
732 ELSE
733 ALLOCATE(SIZLOC_WRITE(NSPGROUP+1))
734 CALL COUNT_ARSZ_SPT(IPARG,DD_IAD,SIZLOC,SIZLOC_WRITE)
735 SIZP0 = SIZLOC
736 SIZW_P0 = -1
737 DO I=1,NSPGROUP
738 SIZW_P0 = MAX(SIZW_P0,SIZLOC_WRITE(I))
739 ENDDO
740 SIZW_P0 = SIZW_P0 + 6
741 DEALLOCATE(SIZLOC_WRITE)
742 ENDIF
743
744 IF(OUTP_SPT(1) == 1) CALL OUTP_SP_T('stress ',
745 . ' stress ',
746 . ELBUF_TAB,IPARG,DD_IAD,SIZLOC,SIZP0,SIZW_P0)
747
748 IF (NSPMD > 1) THEN
749 CALL OUTP_ARSZ_SPTT(IPARG,DD_IAD,SIZLOC,SIZP0,SIZW_P0)
750 ELSE
751 ALLOCATE(SIZLOC_WRITE(NSPGROUP+1))
752 CALL COUNT_ARSZ_SPTT(IPARG,DD_IAD,SIZLOC,SIZLOC_WRITE)
753 SIZP0 = SIZLOC
754 SIZW_P0 = -1
755 DO I=1,NSPGROUP
756 SIZW_P0 = MAX(SIZW_P0,SIZLOC_WRITE(I))
757 ENDDO
758 DEALLOCATE(SIZLOC_WRITE)
759 ENDIF
760
761 IF(OUTP_SPT(2) == 1) CALL OUTP_SP_TT('str_ful ',
762 . ' full stress tensor + plastic strain ',
763 . ELBUF_TAB,IPARG,DD_IAD,IPM, KXSP,SPBUF,SIZLOC,SIZP0,SIZW_P0)
764C-----------------------------------------------
765C SPRING SCALAR
766C-----------------------------------------------
767 IF (NSPMD > 1) THEN
768 CALL OUTP_ARSZ_RS(IPARG,DD_IAD,SIZLOC,SIZP0,SIZW_P0)
769 ELSE
770 ALLOCATE(SIZLOC_WRITE(NSPGROUP+1))
771 CALL COUNT_ARSZ_RS(IPARG,DD_IAD,SIZLOC,SIZLOC_WRITE)
772 SIZP0 = SIZLOC_WRITE(NSPGROUP+1)
773 SIZW_P0 = -1
774 DO I=1,NSPGROUP
775 SIZW_P0 = MAX(SIZW_P0,SIZLOC_WRITE(I))
776 ENDDO
777 SIZW_P0 = SIZW_P0 + 6
778 DEALLOCATE(SIZLOC_WRITE)
779 ENDIF
780
781 IF(OUTP_RS(1) == 1)CALL OUTP_R_S(1,'off ',
782 . ' delete flag(1.0:on 0.0:off) '
783 . ,ELBUF_TAB,IPARG,DD_IAD,SIZLOC,SIZP0,SIZW_P0)
784
785C-----------------------------------------------
786C springs forces and moments
787C-----------------------------------------------
788 IF (NSPMD > 1) THEN
789 CALL OUTP_ARSZ_RT(IPARG,IGEO,GEO,IXR,DD_IAD,
790 . SIZLOC,SIZP0,SIZW_P0)
791 ELSE
792 ALLOCATE(SIZLOC_WRITE(NSPGROUP+1))
793 CALL COUNT_ARSZ_RT(IPARG,IGEO,GEO,IXR,DD_IAD,
794 . SIZLOC,SIZLOC_WRITE)
795 SIZP0 = SIZLOC
796 SIZW_P0 = -1
797 DO I=1,NSPGROUP
798 SIZW_P0 = MAX(SIZW_P0,SIZLOC_WRITE(I))
799 ENDDO
800 SIZW_P0 = SIZW_P0 + 6
801 DEALLOCATE(SIZLOC_WRITE)
802 ENDIF
803C
804 IF(OUTP_RS(2) == 1)CALL OUTP_R_T('full ',
805 . ' forces and moments ',
806 . ELBUF_TAB,IPARG,DD_IAD,IXR,IGEO,GEO,SIZLOC,SIZP0,SIZW_P0)
807
808C-----------------------------------------------
809C RBODIES+SECTIONS+RWALLS FORCES (FOPT)
810C-----------------------------------------------
811 IF(OUTP_V(8) == 1)
812 . CALL OUTP_N_V2(FOPT ,NPBY ,NOM_OPT(I16G),NOM_OPT(I16D))
813C-----------------------------------------------
814C END
815C-----------------------------------------------
816
817 IF(ISPMD == 0) THEN
818 WRITE(IUGEO,'(a)')'/enddata '
819 CLOSE(UNIT=IUGEO)
820 CALL FILELEN_SYS(FILNAM,FILEN, OUTPSIZE)
821 OUTPFILESIZE=OUTPFILESIZE + OUTPSIZE/1024
822 WRITE (IOUT,1000) FILNAM(1:FILEN)
823 WRITE (ISTDO,1000) FILNAM(1:FILEN)
824 ENDIF
825 1000 FORMAT (4X,' output file:',1X,A,' written')
826C-----------
827 RETURN
828 END
genoutp
subroutine genoutp(x, d, v, a, ixs, ixq, ixc, ixt, ixp, ixr, ixtg, iparg, pm, igeo, ms, cont, itab, partsav, fint, fext, tani, eani, ipart, vr, elbuf_tab, dd_iad, weight, ipm, kxsp, spbuf, nodglob, leng, fopt, nom_opt, npby, fncont, ftcont, geo, thke, stack, drape_sh4n, drape_sh3n, drapeg, output)
Definition genoutp.F:83
max
#define max(a, b)
Definition macros.h:21
drape_mod
Definition drape_mod.F:54
inoutfile_mod
Definition inoutfile_mod.F:96
inoutfile_mod::outfile_name
character(len=outfile_char_len) outfile_name
Definition inoutfile_mod.F:111
inoutfile_mod::outfile_name_len
integer outfile_name_len
Definition inoutfile_mod.F:109
stack_mod
Definition stack_mod.F:75
outp_c_s
subroutine outp_c_s(nbxx, key, text, elbuf_tab, iparg, eani, ipm, igeo, ixc, ixtg, dd_iad, sizloc, sizp0, thke, siz_wr)
Definition outp_c_s.F:36
outp_mt
subroutine outp_mt(pm, npart, partsav, ipart, ipm)
Definition outp_mt.F:32
outp_n_vc
subroutine outp_n_vc(key, text, len, vect, numnod, itab, itabg, leng, nodglob, weight)
Definition outp_n_v.F:244
outp_n_v
subroutine outp_n_v(key, text, len, vect, numnod, itab, itabg, leng, nodglob, weight)
Definition outp_n_v.F:32
outp_no
subroutine outp_no(key, text, len, anin, numnod, itab, innoda, itabg, leng, nodglob, weight)
Definition outp_no.F:32
outp_s_s
subroutine outp_s_s(nbx, key, text, elbuf_tab, iparg, eani, ixs, ipm, dd_iad, sizloc, sizp0, siz_wr)
Definition outp_s_s.F:38
count_arsz_ss
subroutine count_arsz_ss(iparg, dd_iad, ipm, ixs, wasz, siz_write_loc)
Definition outp_s_s.F:401
spmd_outpitab
subroutine spmd_outpitab(v, weight, nodglob, vglob)
Definition spmd_outp.F:1081
outp_arsz_ss
subroutine outp_arsz_ss(iparg, dd_iad, ipm, ixs, p0ars, wasz, wasz_wr)
Definition spmd_outp.F:37
slen
subroutine slen(x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, j, area, aream)
Definition slen.F:31