OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
rdresa.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!|| rdresa ../engine/source/output/restart/rdresa.F
25!||--- called by ------------------------------------------------------
26!|| radioss2 ../engine/source/engine/radioss2.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| arret ../engine/source/system/arret.F
30!|| cur_fil_c ../common_source/tools/input_output/write_routtines.c
31!|| rdcomi ../engine/source/output/restart/rdcomm.F
32!|| read_i_c ../common_source/tools/input_output/write_routtines.c
33!||--- uses -----------------------------------------------------
34!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
35!|| h3d_mod ../engine/share/modules/h3d_mod.F
36!|| interfaces_mod ../common_source/modules/interfaces/interfaces_mod.F90
37!|| loads_mod ../common_source/modules/loads/loads_mod.F90
38!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
39!|| message_mod ../engine/share/message_module/message_mod.F
40!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
41!|| names_and_titles_mod ../common_source/modules/names_and_titles_mod.F
42!|| pblast_mod ../common_source/modules/loads/pblast_mod.F90
43!|| pinchtype_mod ../common_source/modules/pinchtype_mod.F
44!|| rbe3_mod ../common_source/modules/constraints/rbe3_mod.F90
45!|| sensor_mod ../common_source/modules/sensor_mod.F90
46!|| skew_mod ../common_source/modules/skew_mod.F90
47!|| state_mod ../common_source/modules/state_mod.F
48!||====================================================================
49 SUBROUTINE rdresa(ITSTV1 , LHEADER, LVARINT, LVARREA, PRADIOS,
50 . PRESFIL,MULTI_FVM,H3D_DATA,
51 . PINCH_DATA,DYNAIN_DATA,INTERFACES,SENSORS,LOADS,
52 . MAT_ELEM, NAMES_AND_TITLES,SKEWS ,GLOB_THERM,PBLAST,
53 . RBE3)
54C-----------------------------------------------
55C M o d u l e s
56C-----------------------------------------------
57 USE message_mod
58 USE multi_fvm_mod
59 USE h3d_mod
60 USE pinchtype_mod
61 USE h3d_mod
62 USE state_mod
63 USE interfaces_mod
64 USE sensor_mod
65 USE loads_mod
66 USE mat_elem_mod
67 USE names_and_titles_mod !< names_and_titles host the input deck names and titles for outputs
68 USE skew_mod !< skew type
69 use glob_therm_mod !< global parameters and flags for thermal analysis
70 use pblast_mod !< pblast parameters
71 use rbe3_mod !< rbe3 type
72C-----------------------------------------------
73C I m p l i c i t T y p e s
74C-----------------------------------------------
75#include "implicit_f.inc"
76C-----------------------------------------------
77C C o m m o n B l o c k s
78C-----------------------------------------------
79#include "r4r8_p.inc"
80#include "com01_c.inc"
81#include "com04_c.inc"
82#include "scr03_c.inc"
83#include "scr14_c.inc"
84#include "warn_c.inc"
85#include "units_c.inc"
86C-----------------------------------------------------------------
87C D u m m y A r g u m e n t s
88C-----------------------------------------------
89 INTEGER ITSTV1, LHEADER, LVARINT, LVARREA, PRADIOS,
90 . PRESFIL
91 TYPE(MULTI_FVM_STRUCT) :: MULTI_FVM
92 TYPE(H3D_DATABASE) :: H3D_DATA
93 TYPE(pinch) :: PINCH_DATA
94 TYPE (DYNAIN_DATABASE),INTENT(INOUT) :: DYNAIN_DATA
95 TYPE (INTERFACES_) ,INTENT(INOUT) :: INTERFACES
96 TYPE (SENSORS_) ,INTENT(INOUT) :: SENSORS
97 TYPE (LOADS_) ,INTENT(INOUT) :: LOADS
98 TYPE (MAT_ELEM_) ,INTENT(INOUT) :: MAT_ELEM
99 TYPE (NAMES_AND_TITLES_),INTENT(INOUT) :: NAMES_AND_TITLES
100 TYPE (SKEW_) ,INTENT(INOUT) :: SKEWS
101 type (glob_therm_) ,intent(inout) :: glob_therm
102 TYPE (PBLAST_) ,INTENT(INOUT) :: PBLAST
103 TYPE (RBE3_) ,INTENT(INOUT) :: RBE3
104C-----------------------------------------------
105C L o c a l V a r i a b l e s
106C-----------------------------------------------
107 INTEGER I
108 INTEGER LMXHEAD,LMXVINT
109 PARAMETER (LMXHEAD=100)
110 parameter (lmxvint=4000)
111 INTEGER TABHEAD(LMXHEAD),TABVINT(LMXVINT)
112 my_real
113 . probint
114 INTEGER,DIMENSION(LTITLE) :: ITITLE
115C-----------------------------------------------
116 DO i=1,lmxvint
117 tabvint(i)=0
118 ENDDO
119C-----------------------------------------------
120C Read Restart File
121C-----------------------------------------------
122 CALL cur_fil_c(0)
123C-----------------------------------------------
124 CALL read_i_c(ititle,ltitle) !< input deck title
125 DO i=1,ltitle
126 names_and_titles%TITLE(i:i) = char(ititle(i))
127 ENDDO
128C-----------------------------------------------
129C Read Header
130C-----------------------------------------------
131 CALL read_i_c(lheader,1)
132
133 CALL read_i_c(tabhead,lheader)
134C RADIOSS FILE R/W
135 pradios=tabhead(1)
136 IF(pradios/=iradios)RETURN
137C RESTART FILE
138 presfil=tabhead(2)
139 IF(presfil/=iresfil)RETURN
140C ITESTV
141 itestv =tabhead(3)
142 itestv = iabs(itestv)
143C itestv +/- depending on R4R8 (single/double precision)
144 IF (ir4r8==2) itestv = -itestv
145 IF(itestv/=itstv1.AND.debug(2)<2)RETURN
146C-----
147C PREVIOUS RUN NUMBER
148 pirun =tabhead(4)
149C CODE OF PREVIOUS RUN (-1:UNKNOWN,0:STARTER,1:ENGINE,2:MODIF)
150 pcodrun=tabhead(5)
151C CODVERS FOR PREVIOUS RUN
152 pcodver=tabhead(6)
153C MINOR VERSION CODE OF PREVIOUS RUN
154 pminver=tabhead(7)
155C SOURCE VERSION CODE OF PREVIOUS RUN
156 psrcver=tabhead(8)
157C INPUT VERSION OF PREVIOUS RUN
158 pinvers=tabhead(9)
159C STARTER CODVERS
160 scodver=tabhead(10)
161C STARTER MINOR VERSION CODE.
162 sminver=tabhead(11)
163C STARTER SOURCE VERSION CODE.
164 ssrcver=tabhead(12)
165C STARTER INPUT VERSION
166 invstr =tabhead(13)
167C
168 lvarint= tabhead(14)
169 lvarrea= tabhead(15)
170C ------------------------------------------
171
172 IF(lvarint>lmxvint)THEN
173 CALL ancmsg(msgid=31,anmode=aninfo,i1=lvarint,i2=lmxvint)
174 CALL arret(2)
175 ENDIF
176
177C----------------------------
178C Read integer variables
179C----------------------------
180 CALL read_i_c(tabvint,lvarint)
181 CALL rdcomi(tabvint,multi_fvm,h3d_data,dynain_data,
182 . interfaces%PARAMETERS,sensors,loads ,mat_elem,skews,
183 . glob_therm,pblast,rbe3 )
184
185C
186 IF(anim_ply > 0 .AND.iplyxfem == 0) THEN
187 WRITE(istdo,*) '**WARNING** /ANIM/SHELL/PLY '
188 WRITE(istdo,*)
189 . 'IS AVAILABLE JUST WITH NEW DELAMINATION SHELL FORMULATION'
190 WRITE(iout,*) '**WARNING** /ANIM/SHELL/PLY '
191 WRITE(iout,*)
192 . 'IS AVAILABLE JUST WITH NEW DELAMINATION SHELL FORMULATION'
193 anim_ply = 0
194 ENDIF
195C
196 RETURN
197 END
198
199!||====================================================================
200!|| iniresa ../engine/source/output/restart/rdresa.F
201!||--- called by ------------------------------------------------------
202!|| radioss2 ../engine/source/engine/radioss2.F
203!||--- uses -----------------------------------------------------
204!|| ale_mod ../common_source/modules/ale/ale_mod.F
205!|| aleanim_mod ../engine/share/modules/aleanim_mod.F
206!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
207!|| h3d_mod ../engine/share/modules/h3d_mod.f
208!|| message_mod ../engine/share/message_module/message_mod.F
209!|| pinchtype_mod ../common_source/modules/pinchtype_mod.F
210!||====================================================================
211 SUBROUTINE iniresa(H3D_DATA,PINCH_DATA,GLOB_THERM)
212C-----------------------------------------------
213C M o d u l e s
214C-----------------------------------------------
215 USE message_mod
216 USE h3d_mod
217 USE pinchtype_mod
218 USE ale_mod
219 USE aleanim_mod
220 USE glob_therm_mod
221C-----------------------------------------------
222C I m p l i c i t T y p e s
223C-----------------------------------------------
224#include "implicit_f.inc"
225C-----------------------------------------------
226C G l o b a l P a r a m e t e r s
227C-----------------------------------------------
228#include "mvsiz_p.inc"
229C-----------------------------------------------
230C C o m m o n B l o c k s
231C-----------------------------------------------
232#include "com01_c.inc"
233#include "com04_c.inc"
234#include "com_xfem1.inc"
235#include "param_c.inc"
236#include "units_c.inc"
237#include "scr03_c.inc"
238#include "scr06_c.inc"
239#include "scr07_c.inc"
240#include "scr14_c.inc"
241#include "scr16_c.inc"
242#include "parit_c.inc"
243#include "task_c.inc"
244#include "scrnoi_c.inc"
245#include "scrfs_c.inc"
246#include "scrcut_c.inc"
247#include "scr_fac_c.inc"
248#include "sphcom.inc"
249#include "spmd_c.inc"
250#include "tabsiz_c.inc"
251#include "scr17_c.inc"
252#include "sms_c.inc"
253C-----------------------------------------------
254C D u m m y A r g u m e n t s
255C-----------------------------------------------
256 TYPE(h3d_database) :: H3D_DATA
257 TYPE(PINCH) :: PINCH_DATA
258 type (glob_therm_) ,intent(inout) :: GLOB_THERM
259C-----------------------------------------------
260C L o c a l V a r i a b l e s
261C-----------------------------------------------
262 INTEGER I,NVAR,M16STAK,NPE,NIPMAX,NPSAV_OLD
263 my_real PROBINT
264C-----------------------------------------------
265C
266 pinch_data%SAPINCH = 0
267C ipari0 from starter
268 IF(iparit==-1)iparit=ipari0
269C
270 IF(icrack3d > 0)THEN
271 IF(iparit /= 1) iparit = ipari0
272 ENDIF
273C
274 IF(ipari0==0.AND.iparit==1) THEN
275 IF(ispmd==0)THEN
276 WRITE(istdo,*)' **WARNING PARALLEL ARITHMETIC DISABLED'
277 WRITE(iout,*)' **WARNING PREVIOUS RUN WITHOUT PARITH/ON'
278 WRITE(iout,*)' **PARALLEL ARITHMETIC IS NO MORE ACTIVABLE'
279 ENDIF
280 iparit=0
281 ENDIF
282
283 nisky0 = 8*numels + 4*numelq + 4*numelc + 2*numelt +
284 . 2*numelp + 3*numelr + 3*numeltg +
285 . 6*numels10 + 12*numels20 + 8*numels16 + 3*numeltg6
286 lcnel = nisky0
287C adding surfaces for MV as elems for parith/on SPMD
288 nisky0 = nisky0 + 4*nskymv0
289C pressure loads
290 nisky0 = nisky0 + 4*nconld
291C thermal fluxes
292 nisky0 = nisky0 + 4*glob_therm%NUMCONV + 4*glob_therm%NUMRADIA + 4*glob_therm%NFXFLUX
293 nisky0 = nisky0 + slloadp
294C
295 IF(iparit>=1) iparit=1
296 IF( lwamp /= 0 .AND. lwanmp/=0)THEN
297 lwamp = max(lwamp,2*mvsiz*100)
298 lenwa = max( nthread * lwamp,lwanmp)
299 ENDIF
300 swa = lenwa
301 sfv = nfunct
302 sicodt=numnod+2*numnod*max(iale,ieuler,ialelag)
303 sicodr=numnod*iroddl
304C-------------------------------------------
305C IPARITH
306C-------------------------------------------
307C AR=AM(M73) MIGHT BE USABLE IF IRODDL=0
308C
309C CONTENTRATED LOADS AND FORCES
310C
311 npsav_old = npsav
312 IF(npsav==0)npsav_old = 23
313 npsav = 29
314 IF(iparit==3)THEN
315 sa = 3*numnod+9*numnod*nthread
316 sar = 3*numnod+9*numnod*nthread
317 sstifn = numnod+3*numnod*nthread
318 sviscn = numnod+3*numnod*nthread
319 sstifr = numnod*iroddl+3*numnod*nthread*iroddl
320 sdmsph = 0
321 spartsav = npsav*(npart+nthpart)*nthread
322 sfsky = 0
323 sfskym = 0
324 sfskyd = 0
325 sisky = 0
326 probint=0.5
327 lskyi_sms = max(nint(5*max(4*numnod,imaximp)*probint)+4*(numsph+nskyi18)+lcni2,isumnx)
328 sadsky = 0
329 sprocne = 0
330 saddcni2 = 0
331 sprocni2 = 0
332 siadsdp = 0
333 siadrcp = 0
334 siads = 0
335 siadwal = 0
336 siadrbk = 0
337 siadi2 = 0
338 siadmv2 = 0
339 siadmv3 = 0
340 siadll = 0
341 siadrbm = 0
342 siadi18 = 0
343 siadrbmk = 0
344 saddcncnd = 0
345 sprocncnd = 0
346 siadcnd = 0
347 ELSEIF(iparit/=0)THEN
348 sa = 3*numnod
349 sar = 3*numnod
350 sstifn = numnod
351 sviscn = numnod
352 sstifr = numnod
353 IF(sol2sph_flag/=0) THEN
354 sdmsph = numnod
355 ELSE
356 sdmsph = 0
357 END IF
358 spartsav = npsav*(npart+nthpart)*nthread
359C LCNE0 = NISKY0+voisins (pre calcule dans starter)
360 lsky = lcne0
361 sfsky = 8*lcne0
362C pxfem
363 lskypxfem = lcnepxfem
364C xfem for layered shell (cracks)
365 lskycrkxfem = lcnecrkxfem
366C FSKY MASSE
367 IF(iale+ieuler+glob_therm%ITHERM+ialelag/=0) THEN
368 sfskym = lcne0
369 ELSE
370 sfskym = 0
371 ENDIF
372C FSKY DMASS
373 IF(sol2sph_flag/=0) THEN
374 sfskyd = lcne0
375 ELSE
376 sfskyd = 0
377 ENDIF
378C PROBINT egalement initialisee dans STARTER (lectur)
379 probint=0.5
380
381 lskyi = max(nint(5*max(4*numnod,imaximp)*probint)
382 + +4*(numsph+nskyi18),isumnx)
383 sfskyi = nfskyi*lskyi
384 sisky = lskyi
385
386 lskyi_sms = lskyi+lcni2
387
388 sadsky = numnod + 1
389C INDEX VECTOR
390 sadsky = sadsky + numnod
391 sprocne = lcne0
392 sprocne_pxfe = lcnepxfem
393 sprocne_crkxfe = lcnecrkxfem
394C si I2NSNT>0 (tableau addresse connectivite interface type 2)
395 IF (i2nsnt>0) THEN
396 saddcni2 = numnod+1
397C Rajout numnod pour index machine vectorielle
398 IF(ivector==1) saddcni2 = saddcni2 + numnod
399 ELSE
400 saddcni2 = 0
401 ENDIF
402C LCNI2 tableau de connectivite
403 sprocni2 = lcni2
404 siadsdp = nspmd+1
405 siadrcp = nspmd+1
406C IADS (correspond bien au NISKY0 : connections locales)
407C IADWAL : NSKYRW0
408C IADRBK : NSKYRBK0
409C IADI2 : NSKYI2
410C IADMV2 : NSKYMV0
411C IADMV3 : NSKYIMVC0
412C IADLL : NSKYLL0
413C IADRBM : NSKYRBM0
414C IADI18 : NSKYI18
415 siads = nisky0
416 siadwal = nskyrw0
417 siadrbk = nskyrbk0
418 siadi2 = niskyi2
419 siadmv2 = nskymv0
420 siadmv3 = nskymvc0
421 siadll = nskyll0
422 siadrbm = nskyrbm0
423 siadi18 = nskyi18
424 siadrbmk = nskyrbmk0
425 ELSE
426 sa = 3*numnod*nthread
427 sar = 3*numnod*nthread
428 sstifn = numnod*nthread
429C pinching
430 IF(npinch > 0) THEN
431 pinch_data%SAPINCH = nthread*numnod
432 ELSE
433 pinch_data%SAPINCH = 0
434 ENDIF
435C
436 IF(sol2sph_flag/=0) THEN
437 sdmsph = numnod*nthread
438 ELSE
439 sdmsph = 0
440 END IF
441C j ai un doute sur SVISCN - et du flag KDTINT
442 sviscn = numnod*nthread
443 sstifr = numnod*iroddl*nthread
444 spartsav = npsav*(npart+nthpart)*nthread
445 sfsky = 0
446 sfskym = 0
447 sfskyd = 0
448 sisky = 0
449
450 probint=0.5
451 lskyi_sms = max(nint(5*max(4*numnod,imaximp)*probint)
452 + +4*(numsph+nskyi18)+lcni2,isumnx)
453
454 sadsky = 0
455 sprocne = 0
456 saddcni2 = 0
457 sprocni2 = 0
458 siadsdp = 0
459 siadrcp = 0
460 siads = 0
461 siadwal = 0
462 siadrbk = 0
463 siadi2 = 0
464 siadmv2 = 0
465 siadmv3 = 0
466 siadll = 0
467 siadrbm = 0
468 siadi18 = 0
469 sprocne_pxfe = 0
470 sprocne_crkxfe = 0
471 ENDIF
472C
473 IF (ns10e >0 .AND. ipari0 /= 0) THEN
474 saddcncnd = numnod + 1
475 sprocncnd = lcncnd
476 siadcnd = 2*ns10e
477 ELSE
478 saddcncnd = 0
479 sprocncnd = 0
480 siadcnd = 0
481 END IF
482C
483 IF(ale%GLOBAL%INCOMP==0)THEN
484 sval2=max(iale,glob_therm%ITHERM,ieuler,ialelag)*(numels+numelq+numeltg+nsvois+nqvois+ntgvois)*max(1,nmult)
485 ! working array phi & segsurf
486 sphi = max(iale,glob_therm%ITHERM,ieuler,ialelag)*( numels + numelq + numeltg + nsvois + nqvois + ntgvois + nsegflu )
487 ELSE
488 sval2 = 0
489 sphi = 0
490 ENDIF
491 sr = 0
492 sestif = 0
493
494C-------------------------------------------
495C LINK
496C-------------------------------------------
497 silink=4*nrlink
498 sfr_rl = (nspmd+2)*nrlink
499
500
501 sllink=lllink
502 siadrl = lllink
503
504C---------------------------------------------------------
505C /FLSEC (FLOW SECTION, obsolete option)
506C---------------------------------------------------------
507 scrflsw = 9*nsflsw
508 sflsw = 6*nsflsw
509 sneflsw = 1*nsflsw
510 snnflsw = 8*ntflsw
511C---------------------------------------------------------
512C ANIMATION (vectors for CONTACT,INTERNAL FORCES, EXTERNAL FORCES)
513C CUTS NOISE
514C---------------------------------------------------------
515 sfani = 3*numnod*(
516 . min(1,anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT)
517 . + min(1,anim_v(5)+outp_v(5)+h3d_data%N_VECT_FINT)
518 . + min(1,anim_v(6)+outp_v(6)+h3d_data%N_VECT_FEXT) )
519 . + 6*(nsect+nrbody+nrwall)
520 . + 6*numnod*min(1,anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT)
521 . + 3*numnod*min(1,anim_v(13)+h3d_data%N_VECT_CONT2)
522 . + 6*numnod*min(1,anim_v(27)+h3d_data%N_VECT_PCONT2)
523
524 sxcut = 7*ncuts
525 sanin = numnod*(
526 . min(1,anim_n(1)+outp_n(1)+h3d_data%N_SCAL_DT)
527 . + min(1,anim_n(2)+outp_n(2)+h3d_data%N_SCAL_DMAS)
528 . + min(1,anim_n(12)+outp_n(3)+h3d_data%N_SCAL_DINER))
529 . + numelr*(anim_fe(11)+anim_fe(12)+anim_fe(13)) + 2*numnod*min(1,anim_n(15)*anim_n(16)+h3d_data%N_SCAL_DAMA2)
530 stani = 6*(numelc+numeltg)*iepsdot + 15*(numelt+numelp+numelr)*(max(anim_ft(1),h3d_data%UND_FORC))
531C
532 ssecfcum = 7*numnod*nsect*isecut
533C Work Array private to sph : WAPSH(LWASPH)
534 swasph = lwasph
535C
536 !/H3D/ELEM/VORT,X,Y,Z
537 fani_cell%IS_VORT_X_REQUESTED = .false.
538 fani_cell%IS_VORT_Y_REQUESTED = .false.
539 fani_cell%IS_VORT_Z_REQUESTED = .false.
540 IF(anim_se(10) == 1 .OR. h3d_data%SOL_SCAL_VORTX == 1)fani_cell%IS_VORT_X_REQUESTED=.true.
541 IF(anim_se(4960) == 1 .OR. h3d_data%SOL_SCAL_VORTY == 1)fani_cell%IS_VORT_Y_REQUESTED=.true.
542 IF(anim_se(4961) == 1 .OR. h3d_data%SOL_SCAL_VORTZ == 1)fani_cell%IS_VORT_Z_REQUESTED=.true.
543
544 ! /H3D/ELEM/VECT/CONT
545 ! INTER18 REACTION FORCES ON CELL CENTROIDS
546 fani_cell%IS_F18_FVM_REQUESTED = .false.
547 IF(h3d_data%N_VECT_CONT == 1)fani_cell%IS_F18_FVM_REQUESTED = .true.
548
549 m16stak = 0
550 IF(numels16/=0)THEN
551 npe = 16
552 nipmax = 81
553 m16stak = (11*nipmax+15*npe+3*nipmax*npe)*mvsiz
554 ENDIF
555 IF(numels20/=0)THEN
556 npe = 20
557 nipmax = 81
558 m16stak = max(m16stak,(10*nipmax+15*npe+3*nipmax*npe)*mvsiz)
559 ENDIF
560 npsav = npsav_old
561C NPSAV remis a la bonne valeur dans rdresb
562 IF(numels8a>0)THEN
563 nipmax = 729
564 m16stak = max(m16stak,12*nipmax*mvsiz)
565 ENDIF
566 sw16 = m16stak * nthread
567 sicut= 44*ncuts
568C---------------------------------------------------------
569C tableaux spmd additionnels alloues dans ma en v40
570C---------------------------------------------------------
571C---------------------------------------------------------
572C tableaux spmd additionnels alloues dans am en v40
573C---------------------------------------------------------
574 sdretrio = 5*ninter
575 slbvrs = 21*nibvel
576C---------------------------------------------------------
577C Stockages facultatifs
578C---------------------------------------------------------
579 IF(nnoise==0)THEN
580 nnoise=nnoiser*rnoi
581 nvar=ncnois*rnoi
582 ELSE
583 IF((noisev+noisea+noisep)==0) noisev = 1
584 nvar=3*noisev+3*noisea+noisep
585 ENDIF
586 if01=1
587 mf01=1
588 ifif=if01
589 mfif=mf01
590 sinoise = 0
591 sfnoise = 0
592 IF(nnoise>0)THEN
593 ifif = (2*nnoise+10)
594 mfif = (6*nvar*nnoise+1)
595 sinoise = 2*nnoise+10
596 sfnoise = 6*nvar*nnoise+1
597 ENDIF
598 nnoise_sav = nnoise
599
600C---------------------------------------------------------
601C
602 RETURN
603 END
604C
my_real
#define my_real
Definition cppsort.cpp:32
the
end diagonal values have been computed in the(sparse) matrix id.SOL
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
ale_mod
Definition ale_mod.F:149
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:249
aleanim_mod
Definition aleanim_mod.F:40
aleanim_mod::fani_cell
type(fani_cell_) fani_cell
Definition aleanim_mod.F:55
h3d_mod
Definition h3d_mod.F:308
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::ltitle
integer, parameter ltitle
Definition names_and_titles_mod.F:1007
pinchtype_mod
Definition pinchtype_mod.F:69
pinchtype_mod::npinch
integer npinch
Definition pinchtype_mod.F:73
state_mod
Definition state_mod.F:54
rdcomi
subroutine rdcomi(tabvint, multi_fvm, h3d_data, dynain_data, inter_parameters, sensors, loads, mat_elem, skews, glob_therm, pblast, rbe3)
Definition rdcomm.F:58
rdresa
subroutine rdresa(itstv1, lheader, lvarint, lvarrea, pradios, presfil, multi_fvm, h3d_data, pinch_data, dynain_data, interfaces, sensors, loads, mat_elem, names_and_titles, skews, glob_therm, pblast, rbe3)
Definition rdresa.F:54
iniresa
subroutine iniresa(h3d_data, pinch_data, glob_therm)
Definition rdresa.F:212
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
arret
subroutine arret(nn)
Definition arret.F:87
h3d_mod::h3d_database
Definition h3d_mod.F:435
pinchtype_mod::pinch
Definition pinchtype_mod.F:76
read_i_c
void read_i_c(int *w, int *len)
Definition write_routtines.c:805
cur_fil_c
void cur_fil_c(int *nf)
Definition write_routtines.c:477