71
72
73
84 USE multi_fvm_mod
98 USE interfaces_mod
99 USE skew_mod
100 USE my_alloc_mod
101 use glob_therm_mod
102 USE connectivity_mod
103 USE nodal_arrays_mod
104 use rbe3_mod
105
106
107
108#include "implicit_f.inc"
109
110
111
112#include "com01_c.inc"
113#include "com04_c.inc"
114#include "com_xfem1.inc"
115#include "param_c.inc"
116#include "remesh_c.inc"
117#include "scr06_c.inc"
118#include "scr07_c.inc"
119#include "scr14_c.inc"
120#include "scr16_c.inc"
121#include "sms_c.inc"
122#include "tabsiz_c.inc"
123#include "task_c.inc"
124#include "couple_c.inc"
125#include "warn_c.inc"
126#include "parit_c.inc"
127#include "intstamp_c.inc"
128#include "spmd_c.inc"
129#include "sphcom.inc"
130#include "fxbcom.inc"
131#include "eigcom.inc"
132#include "inter18.inc"
133#include "stati_c.inc"
134
135 TYPE(connectivity_) :: ELEMENT
136 TYPE(NODAL_ARRAYS_) :: NODES
137 TYPE(MULTI_FVM_STRUCT) :: MULTI_FVM
138 TYPE(H3D_DATABASE) :: H3D_DATA
139 TYPE(PINCH) :: PINCH_DATA
140 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
141 TYPE(t_SEGVAR) :: SEGVAR
142 TYPE(INTERFACES_) :: INTERFACES
143 TYPE(SKEW_) , INTENT(INOUT) :: SKEWS
144 type(glob_therm_) ,intent(inout) :: GLOB_THERM
145 type(rbe3_) ,intent(inout) :: RBE3
146
147
148
149 INTEGER IERR0,IERROR, NBMAT,I,NPT,ISIZ1,ISIZ2
150
151 ierror = 0
152
153 ierr0 = ierr
154
155
156
157
158 CALL allocate_nodal_arrays(nodes, numnod, nthread, iroddl, iparit,
159 . isecut, iisrots, impose_dr, idrot, nrcvvois, sicodt, glob_therm%ITHERM_FE)
160
161
162
163
164
165
166
167
168
169
170 ALLOCATE (
iskwn(siskwn),stat=ierr)
171 IF (ierr/=0) GOTO 1000
173
174 ALLOCATE (
iframe(siframe),stat=ierr)
175 IF (ierr/=0) GOTO 1000
177
178 ALLOCATE (
neth(sneth),stat=ierr)
179 IF (ierr/=0) GOTO 1000
181
182 ALLOCATE (
ibcslag(sibcslag),stat=ierr)
183 IF (ierr/=0) GOTO 1000
185
186 ALLOCATE (
ipart(sipart),stat=ierr)
187 IF (ierr/=0) GOTO 1000
189
190 ALLOCATE (
nom_opt(snom_opt),stat=ierr)
191 IF (ierr/=0) GOTO 1000
193
194 ALLOCATE (
npc(snpc),stat=ierr)
195 IF (ierr/=0) GOTO 1000
197
198 ALLOCATE (
ixtg(sixtg),stat=ierr)
199 IF (ierr/=0) GOTO 1000
201
202 ALLOCATE (
ixtg1(sixtg1),stat=ierr)
203 IF (ierr/=0) GOTO 1000
205
206 ALLOCATE (
ixs(sixs),stat=ierr)
207 IF (ierr/=0) GOTO 1000
209
210 ALLOCATE (
ixq(sixq),stat=ierr)
211 IF (ierr/=0) GOTO 1000
213
214 ALLOCATE (
ixc(sixc),stat=ierr)
215 ALLOCATE(element%shell%ixc(nixc,sixc/nixc))
216 ALLOCATE(element%shell%nodes(4,sixc/nixc))
217 ALLOCATE(element%shell%pid(sixc/nixc))
218 ALLOCATE(element%shell%matid(sixc/nixc))
219 ALLOCATE(element%shell%user_id(sixc/nixc))
221 element%shell%ixc = 0
222
223 ALLOCATE (
ixt(sixt),stat=ierr)
224 IF (ierr/=0) GOTO 1000
226
227 ALLOCATE (
ixp(sixp),stat=ierr)
228 IF (ierr/=0) GOTO 1000
230
231 ALLOCATE (
ixr(sixr),stat=ierr)
232 IF (ierr/=0) GOTO 1000
234
235
236
237
238
239
240
241
242
243 ALLOCATE (
gjbufi(sgjbufi),stat=ierr)
244 IF (ierr/=0) GOTO 1000
246
247 ALLOCATE (ale_connectivity%NALE(
ale%GLOBAL%SNALE),stat=ierr)
248 IF (ierr/=0) GOTO 1000
249 ale_connectivity%NALE(1:
ale%GLOBAL%SNALE) = 0
250
252 IF (ierr/=0) GOTO 1000
254
256 IF (ierr/=0) GOTO 1000
258
259 ALLOCATE (
ifill(sifill),stat=ierr)
260 IF (ierr/=0) GOTO 1000
262
263 ALLOCATE (
ims(sims),stat=ierr)
264 IF (ierr/=0) GOTO 1000
266
267 ALLOCATE (
kxx(skxx),stat=ierr)
268 IF (ierr/=0) GOTO 1000
270
271 ALLOCATE (
ixx(sixx),stat=ierr)
272 IF (ierr/=0) GOTO 1000
274
275 ALLOCATE (
kxsp(skxsp),stat=ierr)
276 IF (ierr/=0) GOTO 1000
278
279 ALLOCATE (
ixsp(sixsp),stat=ierr)
280 IF (ierr/=0) GOTO 1000
282
283 ALLOCATE (
nod2sp(snod2sp),stat=ierr)
284 IF (ierr/=0) GOTO 1000
286
287 ALLOCATE (
ispsym(sispsym),stat=ierr)
288 IF (ierr/=0) GOTO 1000
290
291 ALLOCATE (
ispcond(sispcond),stat=ierr)
292 IF (ierr/=0) GOTO 1000
294
295 ALLOCATE (
isphio(sisphio),stat=ierr)
296 IF (ierr/=0) GOTO 1000
298
299 ALLOCATE (
lprtsph(slprtsph),stat=ierr)
300 IF (ierr/=0) GOTO 1000
302
303 ALLOCATE (
lonfsph(slonfsph),stat=ierr)
304 IF (ierr/=0) GOTO 1000
306
307 ALLOCATE (
ibcl(sibcl),stat=ierr)
308 IF (ierr/=0) GOTO 1000
310
311 ALLOCATE (
ibfv(sibfv),stat=ierr)
312 IF (ierr/=0) GOTO 1000
314
315 ALLOCATE (
iecran(siecran),stat=ierr)
316 IF (ierr/=0) GOTO 1000
318
319 ALLOCATE (
ilas(silas),stat=ierr)
320 IF (ierr/=0) GOTO 1000
322
323 ALLOCATE (
laccelm(slaccelm),stat=ierr)
324 IF (ierr/=0) GOTO 1000
326
327 ALLOCATE (
lgauge(3*nbgauge),stat=ierr)
328 IF (ierr/=0) GOTO 1000
330
331 ALLOCATE (
nnlink(snnlink),stat=ierr)
332 IF (ierr/=0) GOTO 1000
334
335 ALLOCATE (
lnlink(slnlink),stat=ierr)
336 IF (ierr/=0) GOTO 1000
338
339 ALLOCATE (
iparg(siparg),stat=ierr)
340 IF (ierr/=0) GOTO 1000
342
343 ALLOCATE (
dd_iad(sdd_iad),stat=ierr)
344 IF (ierr/=0) GOTO 1000
346
347 ALLOCATE (
igrv(sigrv),stat=ierr)
348 IF (ierr/=0) GOTO 1000
350
351 ALLOCATE (
lgrav(slgrav),stat=ierr)
352 IF (ierr/=0) GOTO 1000
354
355 ALLOCATE (
iexlnk(siexlnk),stat=ierr)
356 IF (ierr/=0) GOTO 1000
358
359 ALLOCATE (
ibvel(sibvel),stat=ierr)
360 IF (ierr/=0) GOTO 1000
362
363 ALLOCATE (
lbvel(slbvel),stat=ierr)
364 IF (ierr/=0) GOTO 1000
366
367 ALLOCATE (
irbe2(sirbe2),stat=ierr)
368 IF (ierr/=0) GOTO 1000
370
371 ALLOCATE (
lrbe2(slrbe2),stat=ierr)
372 IF (ierr/=0) GOTO 1000
374
375 call allocate_rbe3(rbe3, nspmd)
376
377 ALLOCATE (
iactiv(siactiv),stat=ierr)
378 IF (ierr/=0) GOTO 1000
380
381 ALLOCATE (
ibmpc(sibmpc),stat=ierr)
382 IF (ierr/=0) GOTO 1000
384
385 ALLOCATE (
kinet(skinet),stat=ierr)
386 IF (ierr/=0) GOTO 1000
388
389 ALLOCATE (
fasolfr(sfasolfr),stat=ierr)
390 IF (ierr/=0) GOTO 1000
392
393 ALLOCATE (
segquadfr(ssegquadfr),stat=ierr)
394 IF (ierr/=0) GOTO 1000
396
397 ALLOCATE (
ipari(sipari),stat=ierr)
398 IF (ierr/=0) GOTO 1000
400
401
402
403
404
405
406 ALLOCATE (interfaces%INTBUF_TAB(ninter),stat=ierr)
407
408
409
410 ALLOCATE(interfaces%INTBUF_FRIC_TAB(ninterfric),stat=ierr)
411 IF (ierr/=0) GOTO 1000
412
413
414
415
416
417 ALLOCATE (
nprw(snprw),stat=ierr)
418 IF (ierr/=0) GOTO 1000
420
421 ALLOCATE (
lprw(slprw),stat=ierr)
422 IF (ierr/=0) GOTO 1000
424
425 ALLOCATE (
iconx(siconx),stat=ierr)
426 IF (ierr/=0) GOTO 1000
428
429 ALLOCATE (
madprt(smadprt),stat=ierr)
430 IF (ierr/=0) GOTO 1000
432
433 ALLOCATE (
madsh4(smadsh4),stat=ierr)
434 IF (ierr/=0) GOTO 1000
436
437 ALLOCATE (
madsh3(smadsh3),stat=ierr)
438 IF (ierr/=0) GOTO 1000
440
441 ALLOCATE (
madsol(smadsol),stat=ierr)
442 IF (ierr/=0) GOTO 1000
444
445 ALLOCATE (
madnod(smadnod),stat=ierr)
446 IF (ierr/=0) GOTO 1000
448
449 ALLOCATE (
madclnod(madclnods),stat=ierr)
450 IF (ierr/=0) GOTO 1000
452
453 ALLOCATE (
madfail(smadfail),stat=ierr)
454 IF (ierr/=0) GOTO 1000
456 ALLOCATE (
madidx(smadidx),stat=ierr)
457 IF (ierr/=0) GOTO 1000
459
460
461 ALLOCATE (
npby(snpby),stat=ierr)
462 IF (ierr/=0) GOTO 1000
464
465 ALLOCATE (
npbyl(snpbyl),stat=ierr)
466 IF (ierr/=0) GOTO 1000
468
469 ALLOCATE (
lpby(slpby),stat=ierr)
470 IF (ierr/=0) GOTO 1000
472
473 ALLOCATE (
lpbyl(slpbyl),stat=ierr)
474 IF (ierr/=0) GOTO 1000
476
477 ALLOCATE (
lrivet(slrivet),stat=ierr)
478 IF (ierr/=0) GOTO 1000
480
481 ALLOCATE (
nstrf(snstrf),stat=ierr)
482 IF (ierr/=0) GOTO 1000
484
485 ALLOCATE (
nom_sect(snom_sect),stat=ierr)
486 IF (ierr/=0) GOTO 1000
488
489 ALLOCATE (
ljoint(sljoint),stat=ierr)
490 IF (ierr/=0) GOTO 1000
492
493 ALLOCATE (
nodpor(snodpor),stat=ierr)
494 IF (ierr/=0) GOTO 1000
496
497 ALLOCATE (
monvol(smonvol),stat=ierr)
498 IF (ierr/=0) GOTO 1000
500
501 ALLOCATE (
lagbuf(slagbuf),stat=ierr)
502 IF (ierr/=0) GOTO 1000
504
505 ALLOCATE (rthbuf(srthbuf),stat=ierr)
506 IF (ierr/=0) GOTO 1000
507 rthbuf = zero
508
509 ALLOCATE (nodes%BOUNDARY_ADD(2,siad_elem/2),stat=ierr)
510 IF (ierr/=0) GOTO 1000
511 nodes%BOUNDARY_ADD = 0
512
513 ALLOCATE (
iad_rby(siad_rby),stat=ierr)
514 IF (ierr/=0) GOTO 1000
516
517 ALLOCATE (
iad_rby2(siad_rby2),stat=ierr)
518 IF (ierr/=0) GOTO 1000
520
521 ALLOCATE (
iad_i2m(siad_i2m),stat=ierr)
522 IF (ierr/=0) GOTO 100
524
525 ALLOCATE (
iadcj(siadcj),stat=ierr)
526 IF (ierr/=0) GOTO 1000
528
529 ALLOCATE (
iad_rbm(siad_rbm),stat=ierr)
530 IF (ierr/=0) GOTO 1000
532
534 IF (ierr/=0) GOTO 1000
536
537 ALLOCATE (
iad_rbe2(siad_rbe2),stat=ierr)
538 IF (ierr/=0) GOTO 1000
540
541 ALLOCATE (
iad_sec(siad_sec),stat=ierr
542 IF (ierr/=0) GOTO 1000
544
545 ALLOCATE (
iad_cut(siad_cut),stat=ierr)
546 IF (ierr/=0) GOTO 1000
548
549 siad_rbym = siad_rby
550 siad_rbym2 = siad_rby2
551 ALLOCATE (
iad_rbym(siad_rbym),stat=ierr)
552 IF (ierr/=0) GOTO 1000
554
555 ALLOCATE (
iad_rbym2(siad_rbym2),stat=ierr)
556 IF (ierr/=0) GOTO 1000
558
560 IF (ierr/=0) GOTO
562
563 ALLOCATE (
iad_frnor((nspmd+1)*ninter25),stat=ierr)
564 IF (ierr/=0) GOTO 1000
566
567 ALLOCATE (interfaces%SPMD_ARRAYS%FR_EDG(2*nbddedgt),stat=ierr)
568 IF (ierr/=0) GOTO 1000
569 interfaces%SPMD_ARRAYS%FR_EDG = 0
570
571 ALLOCATE (interfaces%SPMD_ARRAYS%IAD_FREDG((nspmd+1)*ninter25),stat=ierr)
572 IF (ierr/=0) GOTO 1000
573 interfaces%SPMD_ARRAYS%IAD_FREDG = 0
574
575 ALLOCATE (
fr_nor(nbddnort),stat=ierr)
576 IF (ierr/=0) GOTO 1000
578
579 ALLOCATE (
procnor(nbccnor),stat=ierr)
580 IF (ierr/=0) GOTO 1000
582 nodes%BOUNDARY_SIZE = sfr_elem
583 ALLOCATE (nodes%BOUNDARY(sfr_elem),stat=ierr)
584 IF (ierr/=0) GOTO 1000
585 nodes%BOUNDARY = 0
586
587 ALLOCATE (
fr_rby(sfr_rby),stat=ierr)
588 IF (ierr/=0) GOTO 1000
590
591 ALLOCATE (
fr_wall(sfr_wall),stat=ierr)
592 IF (ierr/=0) GOTO 1000
594
595 ALLOCATE (
fr_rby2(sfr_rby2),stat=ierr)
596 IF (ierr/=0) GOTO 1000
598
599 ALLOCATE (
fr_i2m(sfr_i2m),stat=ierr)
600 IF (ierr/=0) GOTO 1000
602
603 ALLOCATE (
fr_mv(sfr_mv),stat=ierr)
604 IF (ierr/=0) GOTO 1000
606
607 ALLOCATE (
fr_ll(sfr_ll),stat=ierr)
608 IF (ierr/=0) GOTO 1000
610
611 ALLOCATE (
fr_cj(sfr_cj),stat=ierr)
612 IF (ierr/=0) GOTO 1000
614
615 ALLOCATE (
fr_rbm(sfr_rbm),stat=ierr)
616 IF (ierr/=0) GOTO 1000
618
619 ALLOCATE (
fr_rbm2(sfr_rbm2),stat=ierr)
620 IF (ierr/=0) GOTO 1000
622
623 ALLOCATE (
fr_rbe2(sfr_rbe2),stat=ierr)
624 IF (ierr/=0) GOTO 1000
626
627 ALLOCATE (
fr_sec(sfr_sec),stat=ierr)
628 IF (ierr/=0) GOTO 1000
630
631 ALLOCATE (
fr_cut(sfr_cut),stat=ierr)
632 IF (ierr/=0) GOTO 1000
634
635 ALLOCATE (
rg_cut(srg_cut),stat=ierr)
636 IF (ierr/=0) GOTO 1000
638
639 ALLOCATE (
fr_mad(sfr_mad),stat=ierr)
640 IF (ierr/=0) GOTO 1000
642
643 ALLOCATE (
fr_i18(sfr_i18),stat=ierr)
644 IF (ierr/=0) GOTO 1000
646
647 ALLOCATE (
dd_r2r(sdd_r2r),stat=ierr)
648 IF (ierr/=0) GOTO 1000
650
652 IF (ierr/=0) GOTO 1000
654
655 ALLOCATE (
fr_rbym(sfr_rbym),stat=ierr)
656 IF (ierr/=0) GOTO 1000
658
659 ALLOCATE (
fr_rbym2(sfr_rbym2),stat=ierr)
660 IF (ierr/=0) GOTO 1000
662 ALLOCATE (
iskwp(siskwp),stat=ierr)
663 IF (ierr/=0) GOTO 1000
665
666 ALLOCATE (
nskwp(snskwp),stat=ierr)
667 IF (ierr/=0) GOTO 1000
669
670 ALLOCATE (
iskwp_l(siskwp),stat=ierr)
671 IF (ierr/=0) GOTO 1000
673
674 ALLOCATE (
isensp(sisensp),stat=ierr)
675 IF (ierr/=0) GOTO 1000
677
678 ALLOCATE (
nsensp(snsensp),stat=ierr)
679 IF (ierr/=0) GOTO 1000
681
682 ALLOCATE (
iaccp(siaccp),stat=ierr)
683 IF (ierr/=0) GOTO 1000
685
686 ALLOCATE (
naccp(snaccp),stat=ierr)
687 IF (ierr/=0) GOTO 1000
689
690 ALLOCATE (
igaup(sigaup),stat=ierr)
691 IF (ierr/=0) GOTO 1000
693
694 ALLOCATE (
ngaup(sngaup),stat=ierr)
695 IF (ierr/=0) GOTO 1000
697
698 ALLOCATE (
fr_lagf(sfr_lagf),stat=ierr)
699 IF (ierr/=0) GOTO 1000
701
702
703
704
705
706
707
708
709
710
711
712
713
714 ALLOCATE (
newfront(snewfront),stat=ierr)
715 IF (ierr/=0) GOTO 1000
717
718
719
720
721
722 ALLOCATE (
nbrcvois(snbrcvois),stat=ierr)
723 IF (ierr/=0) GOTO 1000
725
726 ALLOCATE (
lnrcvois(slnrcvois),stat=ierr)
727 IF (ierr/=0) GOTO 1000
729
730 ALLOCATE (
nbsdvois(snbsdvois),stat=ierr)
731 IF (ierr/=0) GOTO 1000
733
734 ALLOCATE (
lnsdvois(slnsdvois),stat=ierr)
735 IF (ierr/=0) GOTO 1000
737
738 ALLOCATE (
nercvois(snercvois),stat=ierr)
739 IF (ierr/=0) GOTO 1000
741
742 ALLOCATE (
lercvois(slercvois),stat=ierr)
743 IF (ierr/=0) GOTO 1000
745
746 ALLOCATE (
nesdvois(snesdvois),stat=ierr)
747 IF (ierr/=0) GOTO 1000
749
750 ALLOCATE (
lesdvois(slesdvois),stat=ierr)
751 IF (ierr/=0) GOTO 1000
753
754 IF(ale_connectivity%has_ne_connect)THEN
756 IF (ierr/=0) GOTO 1000
758
760 IF (ierr/=0) GOTO 1000
762 ELSE
765 ENDIF
766
767 ALLOCATE (
npsegcom(snpsegcom),stat=ierr)
768 IF (ierr/=0) GOTO 1000
770
771 ALLOCATE (
lsegcom(slsegcom),stat=ierr)
772 IF (ierr/=0) GOTO 1000
774
775 ALLOCATE (
nporgeo(snporgeo),stat=ierr)
776 IF (ierr/=0) GOTO 1000
778
779 ALLOCATE (
lnodpor(slnodpor),stat=ierr)
780 IF (ierr/=0) GOTO 1000
782
783 ALLOCATE (
llagf(sllagf),stat=ierr)
784 IF (ierr/=0) GOTO 1000
786
787
788
789
790
791
792
793
794
795
796 ALLOCATE (interfaces%PON%ISKY(sisky),stat=ierr)
797 IF (ierr/=0) GOTO 1000
798 interfaces%PON%ISKY = 0
799
800 ALLOCATE (element%PON%ADSKY(sadsky),stat=ierr
801 IF (ierr/=0) GOTO 1000
802 element%PON%ADSKY = 0
803 element%PON%SADSKY = sadsky
804
805 ALLOCATE (element%PON%PROCNE(sprocne),stat=ierr)
806 IF (ierr/=0) GOTO 1000
807 element%PON%PROCNE = 0
808
809 ALLOCATE (
addcni2(saddcni2),stat=ierr)
810 IF (ierr/=0) GOTO 1000
812
813 ALLOCATE (
procni2(sprocni2),stat=ierr)
814 IF (ierr/=0) GOTO 1000
816
817 ALLOCATE (element%PON%IADSDP(siadsdp),stat=ierr)
818 IF (ierr/=0) GOTO 1000
819 element%PON%IADSDP = 0
820
821 ALLOCATE (element%PON%IADRCP(siadrcp),stat=ierr)
822 IF (ierr/=0) GOTO 1000
823 element%PON%IADRCP = 0
824
825 CALL my_alloc(element%PON%IADS,8,numels,msg="IADS")
826 CALL my_alloc(element%PON%IADS10,6,numels10,msg="iads10")
827 CALL MY_ALLOC(ELEMENT%PON%IADS20,12,NUMELS20,MSG="iads20")
828 CALL my_alloc(element%PON%IADS16,8,numels16,msg="IADS16")
829 CALL my_alloc(element%PON%IADQ,4,numelq,msg="IADQ")
830 CALL my_alloc(element%PON%IADC,4,numelc,msg="IADC")
831 CALL my_alloc(element%PON%IAD_TRUSS,2,numelt,msg="IAD_TRUSS")
832 CALL my_alloc(element%PON%IAD_BEAM,2,numelp,msg="IAD_BEAM")
833 CALL my_alloc(element%PON%IAD_SPRING,3,numelr,msg="IAD_SPRING")
834 CALL my_alloc(element%PON%IAD_TG,3,numeltg,msg="IAD_TG")
835 CALL my_alloc(element%PON%IAD_TG6,3,numeltg6,msg="IAD_TG6")
836 CALL my_alloc(element%PON%IAD_MV,4,nskymv0,msg="IAD_MV")
837 CALL my_alloc(element%PON%IAD_CONLD,4,nconld,msg="IAD_CONLD")
838 CALL my_alloc(element%PON%IAD_CONV,4,glob_therm%NUMCONV,msg="IAD_CONV")
839 CALL my_alloc(element%PON%IAD_RADIA,4,glob_therm%NUMRADIA,msg="IAD_RADIA")
840 CALL my_alloc(element%PON%IAD_LOADP,slloadp,msg="IAD_LOADP")
841 CALL my_alloc(element%PON%IAD_FXFLUX,4,glob_therm%NFXFLUX,msg="IAD_FXFLUX")
842
843
844
845
846
847
848 ALLOCATE (
iadwal(siadwal),stat=ierr)
849 IF (ierr/=0) GOTO 1000
851
852 ALLOCATE (
iadrbk(siadrbk),stat=ierr)
853 IF (ierr/=0) GOTO 1000
855
856 ALLOCATE (
iadi2(siadi2),stat=ierr)
857 IF (ierr/=0) GOTO 1000
859
860 ALLOCATE (
iadmv2(siadmv2),stat=ierr)
861 IF (ierr/=0) GOTO 1000
863
864 ALLOCATE (
iadmv3(siadmv3),stat=ierr)
865 IF (ierr/=0) GOTO 1000
867
868 ALLOCATE (
iadll(siadll),stat=ierr)
869 IF (ierr/=0) GOTO 1000
871
872 ALLOCATE (
iadrbm(siadrbm),stat=ierr)
873 IF (ierr/=0) GOTO 1000
875
876 ALLOCATE (
iadi18(siadi18),stat=ierr)
877 IF (ierr/=0) GOTO 1000
879 ALLOCATE (
iadrbmk(siadrbmk),stat=ierr)
880 IF (ierr/=0) GOTO 1000
882 ALLOCATE (
ilink(silink),stat=ierr)
883 IF (ierr/=0) GOTO 1000
885
886 ALLOCATE (
fr_rl(sfr_rl),stat=ierr)
887 IF (ierr/=0) GOTO 1000
889
890 ALLOCATE (
llink(sllink),stat=ierr)
891 IF (ierr/=0) GOTO 1000
893
894 ALLOCATE (
iadrl(siadrl),stat=ierr)
895 IF (ierr/=0) GOTO 1000
897
898 ALLOCATE (
linale(slinale+llinal),stat=ierr)
899 IF (ierr/=0) GOTO 1000
901
902 ALLOCATE (
neflsw(sneflsw),stat=ierr)
903 IF (ierr/=0) GOTO 1000
905
906 ALLOCATE (
nnflsw(snnflsw),stat=ierr)
907 IF (ierr/=0) GOTO 1000
909
910 ALLOCATE (
icut(sicut),stat=ierr)
911 IF (ierr/=0) GOTO 1000
913
914 ALLOCATE (
igeo(sigeo),stat=ierr)
915 IF (ierr/=0) GOTO 1000
917
918 ALLOCATE (
ipm(sipm),stat=ierr)
919 IF (ierr/=0) GOTO 1000
921
923 IF (ierr/=0) GOTO 1000
925
926 ALLOCATE (
icontact(sicontact),stat=ierr)
927 IF (ierr/=0) GOTO 1000
929
930 IF(nadmesh /= 0)THEN
931 ALLOCATE (
sh4tree(ksh4tree,numelc),stat=ierr)
932 ierror = ierror+ierr
933 ALLOCATE (
sh3tree(ksh3tree,numeltg),stat=ierr)
934 ierror = ierror+ierr
935 ALLOCATE (
ipadmesh(kipadmesh,npart),stat=ierr)
936 ierror = ierror+ierr
937 ALLOCATE (
sh4trim(abs(lsh4trim)),stat=ierr)
938 ierror = ierror+ierr
939 ALLOCATE (
sh3trim(abs(lsh3trim)),stat=ierr)
940 ierror = ierror+ierr
941 ELSE
947 END IF
948
949 ALLOCATE (
ithvar(sithvar),stat=ierr)
950 IF (ierr/=0) GOTO 1000
952
953 ALLOCATE (
icfield(sizfield*nloadc),stat=ierr)
954 IF (ierr/=0) GOTO 1000
956
957 ALLOCATE (
lcfield(slcfield),stat=ierr)
958 IF (ierr/=0) GOTO 1000
960
961 ALLOCATE (
iloadp(sizloadp*nloadp),stat=ierr)
962 IF (ierr/=0) GOTO 1000
964
965 ALLOCATE (
lloadp(slloadp),stat=ierr)
966 IF (ierr/=0) GOTO 1000
968
970 IF (ierr/=0) GOTO 1000
972
973 IF(nsphsol/=0)THEN
974
975 ALLOCATE(
sph2sol(numsph),stat=ierr)
976 IF (ierr/=0) GOTO 1000
978
979 ALLOCATE(
sol2sph(2*numels8),stat=ierr)
980 IF (ierr/=0) GOTO 1000
982
983 ALLOCATE(
irst(3*nsphsol),stat=ierr)
984 IF (ierr/=0) GOTO 1000
986
988 IF (ierr/=0) GOTO 1000
990
991 ELSE
992
994 IF (ierr/=0) GOTO 1000
996 IF (ierr/=0) GOTO 1000
997 ALLOCATE(
irst(0),stat=ierr)
998 IF (ierr/=0) GOTO 1000
1000 IF (ierr/=0) GOTO 1000
1001
1002 END IF
1003
1004 ALLOCATE (element%PON%FSKYD(sfskyd),stat=ierr)
1005 IF (ierr/=0) GOTO 1000
1006 element%PON%FSKYD = 0
1007
1008 ALLOCATE (
ibcscyc(sibcscyc),stat=ierr)
1009 IF (ierr/=0) GOTO 1000
1011
1012 ALLOCATE (
lbcscyc(slbcscyc),stat=ierr)
1013 IF (ierr/=0) GOTO 1000
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040 ALLOCATE (thke(sthke),stat=ierr)
1041 IF (ierr/=0) GOTO 1000
1042 thke = 0
1043
1044 ALLOCATE (dampr(sdampr),stat=ierr)
1045 IF (ierr/=0) GOTO 1000
1046 dampr = 0
1047
1048 ALLOCATE (damp(sdamp),stat=ierr)
1049 IF (ierr/=0) GOTO 1000
1050 damp = 0
1051
1052
1053
1054
1055
1056
1057
1058
1059
1060 ALLOCATE (tf(stf),stat=ierr)
1061 IF (ierr/=0) GOTO 1000
1062 tf = 0
1063
1064 ALLOCATE (pm(spm),stat=ierr)
1065 IF (ierr/=0) GOTO 1000
1066 pm = 0
1067
1068 ALLOCATE (xframe(sxframe),stat=ierr)
1069 IF (ierr/=0) GOTO 1000
1070 xframe = 0
1071
1072 ALLOCATE (geo(sgeo),stat=ierr)
1073 IF (ierr/=0) GOTO 1000
1074 geo = 0
1075
1076 ALLOCATE (eani(seani),stat=ierr)
1077 IF (ierr/=0) GOTO 1000
1078 eani = 0
1079
1080 ALLOCATE (bufmat(sbufmat),stat=ierr)
1081 IF (ierr/=0) GOTO 1000
1082 bufmat = 0
1083
1084 ALLOCATE (
bufgeo(sbufgeo),stat=ierr)
1085 IF (ierr/=0) GOTO 1000
1087
1088 ALLOCATE (bufsf(sbufsf),stat=ierr)
1089 IF (ierr/=0) GOTO 1000
1090 bufsf = 0
1091
1092 ALLOCATE (rbmpc(srbmpc),stat=ierr)
1093 IF (ierr/=0) GOTO 1000
1094 rbmpc = 0
1095
1096 ALLOCATE (gjbufr(sgjbufr),stat=ierr)
1097 IF (ierr/=0) GOTO 1000
1098 gjbufr = 0
1099
1100 ALLOCATE (w(sw),stat=ierr)
1101 IF (ierr/=0) GOTO 1000
1102 w = 0
1103
1104 ALLOCATE (veul(sveul),stat=ierr)
1105 IF (ierr/=0) GOTO 1000
1106 veul = 0
1107
1108 ALLOCATE (fill(sfill),stat=ierr)
1109 IF (ierr/=0) GOTO 1000
1110 fill = 0
1111
1112 ALLOCATE (dfill(sdfill),stat=ierr)
1113 IF (ierr/=0) GOTO 1000
1114 dfill = 0
1115
1116 ALLOCATE (alph(salph),stat=ierr)
1117 IF (ierr/=0) GOTO 1000
1118 alph = 0
1119
1120 ALLOCATE (wb(swb),stat=ierr)
1121 IF (ierr/=0) GOTO 1000
1122 wb = 0
1123
1124 ALLOCATE (dsave(sdsave),stat=ierr)
1125 IF (ierr/=0) GOTO 1000
1126 dsave = 0
1127
1128 ALLOCATE (asave(sasave),stat=ierr)
1129 IF (ierr/=0) GOTO 1000
1130 asave = 0
1131
1132 ALLOCATE (spbuf(sspbuf),stat=ierr)
1133 IF (ierr/=0) GOTO 1000
1134 spbuf = 0
1135
1136 ALLOCATE (vsphio(svsphio),stat=ierr)
1137 IF (ierr/=0) GOTO 1000
1138 vsphio = 0
1139
1140 ALLOCATE (sphveln(ssphveln),stat=ierr)
1141 IF (ierr/=0) GOTO 1000
1142 sphveln = 0
1143
1144 ALLOCATE (msnf(smsnf),stat=ierr)
1145 IF (ierr/=0) GOTO 1000
1146 msnf = 0
1147
1148
1149 ALLOCATE (forc(sforc),stat=ierr)
1150 IF (ierr/=0) GOTO 1000
1151 forc = 0
1152
1153 ALLOCATE (vel(svel),stat=ierr)
1154 IF (ierr/=0) GOTO 1000
1155 vel = 0
1156
1157 ALLOCATE (fsav(sfsav),stat=ierr)
1158 IF (ierr/=0) GOTO 1000
1159 fsav = 0
1160
1161 ALLOCATE (
tab_ump(7,taille),stat=ierr)
1162 IF (ierr/=0) GOTO 1000
1164
1165 ALLOCATE (
poin_ump(nummat),stat=ierr)
1166 IF (ierr/=0) GOTO 1000
1168
1169 IF(ispmd==0) THEN
1170 ALLOCATE (cputime_mp_glob(nbr_gpmp,nspmd+1),stat=ierr)
1171 IF (ierr/=0) GOTO 1000
1172 cputime_mp_glob = 0
1173 ENDIF
1174
1175 ALLOCATE (cputime_mp(taille),stat=ierr)
1176 IF (ierr/=0) GOTO 1000
1177 cputime_mp = 0
1178
1179 ALLOCATE (tab_mat(stab_mat),stat=ierr)
1180 IF (ierr/=0) GOTO 1000
1181 tab_mat = zero
1182
1183 ALLOCATE (fzero(sfzero),stat=ierr)
1184 IF (ierr/=0) GOTO 1000
1185 fzero = 0
1186
1187 ALLOCATE (xlas(sxlas
1188 IF (ierr/=0) GOTO 1000
1189 xlas = 0
1190
1191 ALLOCATE (accelm(saccelm),stat=ierr)
1192 IF (ierr/=0) GOTO 1000
1193 accelm = 0
1194
1195 ALLOCATE (gauge(llgauge*nbgauge),stat=ierr)
1196 IF (ierr/=0) GOTO 1000
1197 gauge = 0
1198
1199 ALLOCATE (fbvel(sfbvel),stat=ierr)
1200 IF (ierr/=0) GOTO 1000
1201 fbvel = 0
1202
1203 ALLOCATE (factiv(lractiv*nactiv),stat=ierr)
1204 IF (ierr/=0) GOTO 1000
1205 factiv = 0
1206
1207 ALLOCATE (grav(sgrav),stat=ierr)
1208 IF (ierr/=0) GOTO 1000
1209 grav = 0
1210
1211 ALLOCATE (fr_wave(sfr_wave),stat=ierr)
1212 IF (ierr/=0) GOTO 1000
1213 fr_wave = 0
1214
1215 ALLOCATE (parts0(sparts0),stat=ierr)
1216 IF (ierr/=0) GOTO 1000
1217 parts0 = 0
1218
1219 ALLOCATE (elbuf(selbuf),stat=ierr)
1220 IF (ierr/=0) GOTO 1000
1221 elbuf = 0
1222
1223 ALLOCATE (rwbuf(srwbuf),stat=ierr)
1224 IF (ierr/=0) GOTO 1000
1225 rwbuf = 0
1226
1227 ALLOCATE (rwsav(srwsav),stat=ierr)
1228 IF (ierr/=0) GOTO 1000
1229 rwsav = 0
1230
1231 ALLOCATE (rby(srby),stat=ierr)
1232 IF (ierr/=0) GOTO 1000
1233 rby = 0
1234
1235 ALLOCATE (rbyl(srbyl),stat=ierr)
1236 IF (ierr/=0) GOTO 1000
1237 rbyl = 0
1238
1239 ALLOCATE (rivet(srivet),stat=ierr)
1240 IF (ierr/=0) GOTO 1000
1241 rivet = 0
1242
1243 ALLOCATE (secbuf(ssecbuf),stat
1244 IF (ierr/=0) GOTO 1000
1245 secbuf = 0
1246
1247 ALLOCATE (volmon(svolmon),stat=ierr)
1248 IF (ierr/=0) GOTO 1000
1249 volmon = 0
1250
1251 ALLOCATE (lambda(slambda),stat=ierr)
1252 IF (ierr/=0) GOTO 1000
1253 lambda = 0
1254
1255 ALLOCATE (rconx(srconx),stat=ierr)
1256 IF (ierr/=0) GOTO 1000
1257 rconx = 0
1258
1259 ALLOCATE (rcontact(srcontact),stat=ierr)
1260 IF (ierr/=0) GOTO 1000
1261 rcontact = ep30
1262
1263 ALLOCATE (acontact(srcontact),stat=ierr)
1264 IF (ierr/=0) GOTO 1000
1265 acontact = ep30
1266
1267 ALLOCATE (pcontact(srcontact),stat=ierr)
1268 IF (ierr/=0) GOTO 1000
1269 pcontact = zero
1270
1271 ALLOCATE (wa(swa),stat=ierr)
1272 IF (ierr/=0) GOTO 1000
1273 wa = 0
1274
1275 ALLOCATE (fv(sfv),stat=ierr)
1276 IF (ierr/=0) GOTO 1000
1277 fv = 0
1278
1279
1280
1281
1282
1283
1284! IF (ierr/=0) GOTO 1000
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
1295
1296
1297
1298
1299 ALLOCATE (partsav(spartsav),stat=ierr)
1300 IF (ierr/=0) GOTO 1000
1301 partsav = 0
1302
1303 ALLOCATE (element%PON%FSKY(8,sfsky/8),stat=ierr)
1304 IF (ierr/=0) GOTO 1000
1305 element%PON%FSKY = 0
1306 element%PON%SFSKY = sfsky
1307 element%PON%MAX_SFSKY = sfsky
1308
1309
1310 ALLOCATE (element%PON%FSKYM(sfskym),stat=ierr)
1311 IF (ierr/=0) GOTO 1000
1312 element%PON%FSKYM = 0
1313
1314 ALLOCATE (interfaces%PON%FSKYI(lskyi,nfskyi),stat=ierr)
1315 IF (ierr/=0) GOTO 1000
1316 interfaces%PON%FSKYI = 0
1317
1318 ALLOCATE (uwa(suwa),stat=ierr)
1319 IF (ierr/=0) GOTO 1000
1320 uwa = 0
1321
1322 ALLOCATE (val2(sval2),stat=ierr)
1323 IF (ierr/=0) GOTO 1000
1324 val2 = 0
1325
1326 ALLOCATE (phi(sphi),stat=ierr)
1327 IF (ierr/=0) GOTO 1000
1328 phi = 0
1329
1330 ALLOCATE (phie(sphi),stat=ierr)
1331 IF (ierr/=0) GOTO 1000
1332 phie = zero
1333 ALLOCATE (por(sphi),stat=ierr)
1334 IF (ierr/=0) GOTO 1000
1335 por = zero
1336
1337 ALLOCATE (r(sr),stat=ierr)
1338 IF (ierr/=0) GOTO 1000
1339 r = 0
1340
1341 ALLOCATE (estif(sestif),stat=ierr)
1342 IF (ierr/=0) GOTO 1000
1343 estif = 0
1344
1345 ALLOCATE (crflsw(scrflsw),stat=ierr)
1346 IF (ierr/=0) GOTO 1000
1347 crflsw = 0
1348
1349 ALLOCATE (flsw(sflsw),stat=ierr)
1350 IF (ierr/=0) GOTO 1000
1351 flsw = 0
1352
1353 ALLOCATE (fani(sfani),stat=ierr)
1354 IF (ierr/=0) GOTO 1000
1355 fani = 0
1356
1357 isiz1=2
1358 isiz2=numelq+numeltg
1359 IF(n2d == 0)THEN
1360 isiz1=3
1361 isiz2=numels
1362 ENDIF
1363
1365 ALLOCATE (
fani_cell%VORT_X(isiz2),stat=ierr)
1366 IF (ierr/=0) GOTO 1000
1368 ENDIF
1370 ALLOCATE (
fani_cell%VORT_Y(isiz2),stat=ierr)
1371 IF (ierr/=0) GOTO 1000
1373 ENDIF
1375 ALLOCATE (
fani_cell%VORT_Z(isiz2),stat=ierr)
1376 IF (ierr/=0) GOTO 1000
1378 ENDIF
1379
1381 ALLOCATE (
fani_cell%F18(isiz1,isiz2),stat=ierr)
1382 IF (ierr/=0) GOTO 1000
1384 ENDIF
1385
1386 ALLOCATE (xcut(sxcut),stat=ierr)
1387 IF (ierr/=0) GOTO 1000
1388 xcut = 0
1389
1390 ALLOCATE (anin(sanin),stat=ierr)
1391 IF (ierr/=0) GOTO 1000
1392 anin = 0
1393
1394 ALLOCATE (tani(stani),stat=ierr)
1395 IF (ierr/=0) GOTO 1000
1396 tani = 0
1397
1398 ALLOCATE (secfcum(ssecfcum),stat=ierr)
1399 IF (ierr/=0) GOTO 1000
1400 secfcum = 0
1401
1402 ALLOCATE (wasph(swasph),stat=ierr)
1403 IF (ierr/=0) GOTO
1404 wasph
1405
1406ALLOCATE (w16(sw16),stat=ierr)
1407 IF (ierr/=0) GOTO
1408 w16 = 0
1409
1410 ALLOCATE (dretrio(sdretrio),stat=ierr)
1411 IF (ierr/=0) GOTO 1000
1412
1413
1414 ALLOCATE (lbvrs(slbvrs),stat=ierr)
1415 IF (ierr/=0) GOTO 1000
1416 lbvrs = 0
1417
1418 ALLOCATE (pv(spv),stat=ierr)
1419 IF (ierr/=0) GOTO 1000
1420 pv = 0
1421
1422 IF(nadmesh /= 0 .OR. irest_mselt/=0)THEN
1423
1424 ALLOCATE(msc(numelc),stat=ierr)
1425 ierror = ierror+ierr
1426
1427 ALLOCATE(inc(numelc),stat=ierr)
1428 ierror = ierror+ierr
1429
1430 ALLOCATE(mstg(numeltg),stat=ierr)
1431 ierror = ierror+ierr
1432
1433 ALLOCATE(intg(numeltg),stat=ierr)
1434 ierror = ierror+ierr
1435
1436 ALLOCATE(ptg(3,numeltg),stat=ierr)
1437 ierror = ierror+ierr
1438
1439 ELSE
1440
1441 ALLOCATE(msc(0))
1442 ALLOCATE(inc(0))
1443 ALLOCATE(mstg(0))
1444 ALLOCATE(intg(0))
1445 ALLOCATE(ptg(0,0))
1446
1447 END IF
1448
1449 IF(irest_mselt/=0)THEN
1450
1451 ALLOCATE(mssa(numels),stat=ierr)
1452 ierror = ierror+ierr
1453
1454 ALLOCATE(mstr(numelt
1455 ierror = ierror+ierr
1456
1457 ALLOCATE(msp(numelp),stat=ierr)
1458 ierror = ierror
1459
1460 ALLOCATE(msrt(numelr),stat=ierr)
1461 ierror = ierror+ierr
1462
1463 ELSE
1464
1465 ALLOCATE(mssa(0))
1466 ALLOCATE(mstr(0))
1467 ALLOCATE(msp(0))
1468 ALLOCATE(msrt(0))
1469
1470 END IF
1471
1472 IF(nadmesh /= 0)THEN
1473
1474 ALLOCATE(padmesh(kpadmesh,npart),stat=ierr)
1475 ierror = ierror+ierr
1476
1477 IF(glob_therm%ITHERM_FE > 0)THEN
1478
1479 ALLOCATE(mcpc(numelc),stat=ierr)
1480 ierror = ierror+ierr
1481
1482 ALLOCATE(mcptg(numeltg),stat=ierr)
1483 ierror = ierror+ierr
1484
1485 ELSE
1486
1487 ALLOCATE(mcpc(0),mcptg(0))
1488
1489 END IF
1490
1491 ELSE
1492 ALLOCATE(padmesh(0,0))
1493 ALLOCATE(mcpc(0),mcptg(0))
1494 END IF
1495
1496 IF(istatcnd/=0)THEN
1497 ALLOCATE(mscnd(numnod),incnd(numnod),stat=ierr)
1498 ierror = ierror+ierr
1499 ELSE
1500 ALLOCATE(mscnd(0),incnd(0))
1501 END IF
1502
1503
1504
1505 IF(irigid_mat > 0) THEN
1506 ALLOCATE(rbym(nfrbym*nrbym),stat=ierr)
1507 IF (ierr/=0) GOTO 1000
1508 rbym = zero
1509 ALLOCATE(
irbym(nirbym*nrbym),stat=ierr)
1510 IF (ierr/=0) GOTO 1000
1512 ALLOCATE(
lnrbym(ngslnrbym),stat=ierr)
1513 IF (ierr/=0) GOTO 1000
1516 IF (ierr/=0) GOTO 1000
1518 ALLOCATE(
icodrbym(nrbym*2),stat=ierr)
1519 IF (ierr/=0) GOTO 1000
1521 ELSE
1524 ENDIF
1525
1526
1527 IF(glob_therm%NUMCONV > 0) THEN
1528 ALLOCATE(
ibcv(glob_therm%NICONV*glob_therm%NUMCONV),stat=ierr)
1529 IF (ierr/=0) GOTO 1000
1531 ALLOCATE(fconv(glob_therm%LFACTHER*glob_therm%NUMCONV),stat=ierr)
1532 IF (ierr/=0) GOTO 1000
1533 fconv = zero
1534 ELSE
1535 ALLOCATE(
ibcv(0),fconv(0))
1536 ENDIF
1537
1538 IF(glob_therm%NUMRADIA > 0) THEN
1539 ALLOCATE(
ibcr(glob_therm%NIRADIA*glob_therm%NUMRADIA),stat=ierr)
1540 IF (ierr/=0) GOTO 1000
1542 ALLOCATE(fradia(glob_therm%LFACTHER*glob_therm%NUMRADIA),stat=ierr)
1543 IF (ierr/=0) GOTO 1000
1544 fradia = zero
1545 ELSE
1546 ALLOCATE(
ibcr(0),fradia(0))
1547 ENDIF
1548
1549 IF(glob_therm%NFXFLUX > 0) THEN
1550 ALLOCATE(
ibfflux(glob_therm%NITFLUX*glob_therm%NFXFLUX),stat=ierr)
1551 IF (ierr/=0) GOTO 1000
1553 ALLOCATE(fbfflux(glob_therm%LFACTHER*glob_therm%NFXFLUX),stat=ierr)
1554 IF (ierr/=0) GOTO 1000
1555 fbfflux = zero
1556 ELSE
1557 ALLOCATE(
ibfflux(0),fbfflux(0))
1558 ENDIF
1559
1560 IF(glob_therm%NFXTEMP > 0) THEN
1561 ALLOCATE(
ibftemp(glob_therm%NIFT
1562 IF (ierr/=0) GOTO 1000
1564 ALLOCATE(fbftemp(glob_therm%LFACTHER*glob_therm%NFXTEMP),stat=ierr)
1565 IF(ierr/=0) GOTO 1000
1566 fbftemp = zero
1567 ELSE
1568 ALLOCATE(
ibftemp(0),fbftemp(0))
1569 ENDIF
1570
1571 ALLOCATE (msf(sms),stat=ierr)
1572 IF (ierr/=0) GOTO 1000
1573 IF(ialelag > 0) THEN
1574 msf = zero
1575 ENDIF
1576
1577
1578 IF(iplyxfem > 0 ) THEN
1579 ALLOCATE(ms_ply(nplyxfe*nplymax),stat=ierr)
1580 IF (ierr/=0) GOTO 1000
1581 ms_ply = zero
1582 ALLOCATE(zi_ply(nplyxfe*nplymax),stat=ierr)
1583 IF (ierr/=0) GOTO 1000
1584 zi_ply = zero
1585 ALLOCATE(msz2(nplyxfe),stat=ierr)
1586 IF (ierr/=0) GOTO 1000
1587 msz2 = zero
1589 IF (ierr/=0) GOTO 1000
1592 IF (ierr/=0) GOTO 1000
1595 IF (ierr/=0) GOTO 1000
1599
1601 IF (ierr/=0) GOTO 1000
1603 ALLOCATE
1604 IF (ierr/=0) GOTO 1000
1606
1608 IF (ierr/=0) GOTO 1000
1610 ELSE
1612 ALLOCATE(msz2(0))
1615 ENDIF
1616 IF(iplybcs > 0) THEN
1618 IF (ierr/=0) GOTO 1000
1621 IF (ierr/=0) GOTO 1000
1624 IF (ierr/=0) GOTO 1000
1626 ALLOCATE(
ibc_ply(numnod),stat=ierr)
1627 IF (ierr/=0) GOTO 1000
1629 ELSE
1631 ENDIF
1632 IF(intplyxfem > 0 ) THEN
1633 ALLOCATE(wagap
1634 IF (ierr/=0) GOTO 1000
1635 wagap = zero
1636 ELSE
1637 ALLOCATE(wagap(0),stat=ierr)
1638 ENDIF
1639
1640
1641
1642 IF(idtmins_old==1) THEN
1643 ALLOCATE(admsms(numnod),stat=ierr)
1644 IF (ierr/=0) GOTO 1000
1645 END IF
1646 IF(idtmins_old==2) THEN
1647 ALLOCATE(dmelc(numelc) ,stat=ierr)
1648 IF (ierr/=0) GOTO 1000
1649 ALLOCATE(dmeltg(numeltg),stat=ierr)
1650 IF (ierr/=0) GOTO 1000
1651 ALLOCATE(dmels(numels) ,stat=ierr)
1652 IF (ierr/=0) GOTO 1000
1653 ALLOCATE(dmeltr(numelt) ,stat=ierr)
1654 IF (ierr/=0) GOTO 1000
1655 ALLOCATE(dmelp(numelp) ,stat=ierr)
1656 IF (ierr/=0) GOTO 1000
1657 ALLOCATE(dmelrt(numelr) ,stat=ierr)
1658 IF (ierr/=0) GOTO 1000
1659 END IF
1660 IF(idtmins_old/=0.OR.idtmins_int_old/=0) THEN
1661 ALLOCATE(res_sms(3*numnod),stat=ierr)
1662 IF (ierr/=0) GOTO 1000
1663 END IF
1664 IF(idtmins_old==2.OR.idtmins_int_old/=0) THEN
1665 ALLOCATE(diag_sms(numnod),stat=ierr)
1666 IF (ierr/=0) GOTO 1000
1667 ALLOCATE(dmint2(4*i2nsn25),stat=ierr)
1668 IF (ierr/=0) GOTO 1000
1669 END IF
1670 IF (isms_selec /= 0) THEN
1672 IF (ierr/=0) GOTO 1000
1674 ELSE
1676 ENDIF
1677
1678
1679 IF(nintstamp/=0.AND.anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0)THEN
1680 ALLOCATE(fcontg(3,numnodg),stat=ierr)
1681 IF (ierr/=0) GOTO 1000
1682 fcontg=zero
1683 ELSE
1684 ALLOCATE(fcontg(0,0),stat=ierr)
1685 END IF
1686
1688 IF(anim_v(26)+h3d_data%N_VECT_CONT_MAX >0.OR.
ifcontmax > 0)
THEN
1689 IF(nspmd > 1 .AND. nintstamp/=0) THEN
1690 IF(ispmd == 0) THEN
1691 ALLOCATE(fcont_max(3*numnodg),stat=ierr)
1692 IF (ierr/=0) GOTO 1000
1693 fcont_max=zero
1695 ELSE
1696 ALLOCATE(fcont_max(0),stat=ierr)
1697 ENDIF
1698 ELSE
1699 ALLOCATE(fcont_max(3*numnod),stat=ierr)
1700 IF (ierr/=0) GOTO 1000
1701 fcont_max=zero
1703 ENDIF
1704 ELSE
1705 ALLOCATE(fcont_max(0),stat=ierr)
1706 END IF
1707
1708
1709 IF(nintstamp/=0.AND.anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT>0)THEN
1710 ALLOCATE(fncontg(3,numnodg),ftcontg(3,numnodg),stat=ierr)
1711 IF (ierr/=0) GOTO 1000
1712 fncontg=zero
1713 ftcontg=zero
1714 ELSE
1715 ALLOCATE(fncontg(0,0),ftcontg(0,0),stat=ierr)
1716 END IF
1717
1719 IF(h3d_data%N_VECT_PCONT_MAX >0.OR.
ifcontpmax > 0)
THEN
1720 IF(nspmd > 1 .AND. nintstamp/=0) THEN
1721 IF(ispmd == 0) THEN
1722 ALLOCATE(ftcont_max(3*numnodg),fncont_max(3*numnodg),stat=ierr)
1723 IF (ierr/=0) GOTO 1000
1724 ftcont_max=zero
1725 fncont_max=zero
1727 ELSE
1728 ALLOCATE(fncont_max(0),ftcont_max(0),stat=ierr)
1729 ENDIF
1730 ELSE
1731 ALLOCATE(ftcont_max(3*numnod),fncont_max(3*numnod),stat=ierr)
1732 IF (ierr/=0) GOTO 1000
1733 ftcont_max=zero
1734 fncont_max=zero
1736 ENDIF
1737 ELSE
1738 ALLOCATE(fncont_max(0),stat=ierr)
1739 ALLOCATE(ftcont_max(0),stat=ierr)
1740 END IF
1741
1742 IF(h3d_data%N_VECT_CONT2_MAX >0.OR.
ifcont2max > 0)
THEN
1743 ALLOCATE(fcont2_max(3*numnod),stat=ierr)
1744 IF (ierr/=0) GOTO 1000
1745 fcont2_max=zero
1746 ELSE
1747 ALLOCATE(fcont2_max(0),stat=ierr)
1748 END IF
1749
1750 IF(h3d_data%N_VECT_PCONT2_MAX >0.OR.
ifcontp2max > 0)
THEN
1751 ALLOCATE(ftcont2_max(3*numnod),fncont2_max(3*numnod),stat=ierr)
1752 ALLOCATE(npcont2_max(3*numnod),stat=ierr)
1753 IF (ierr/=0) GOTO 1000
1754 fncont2_max=zero
1755 ftcont2_max=zero
1756 npcont2_max=zero
1757 ELSE
1758 ALLOCATE(fncont2_max(0),stat=ierr)
1759 ALLOCATE(ftcont2_max(0),stat=ierr)
1760 ALLOCATE(npcont2_max(0),stat=ierr)
1761 END IF
1762
1763
1764 IF(h3d_data%N_VECT_CONT2_MIN >0.OR.
ifcont2min > 0)
THEN
1765 ALLOCATE(fcont2_min(3*numnod),stat=ierr)
1766 IF (ierr/=0) GOTO 1000
1767 fcont2_min=zero
1768 ELSE
1769 ALLOCATE(fcont2_min(0),stat=ierr)
1770 END IF
1771
1772 IF(h3d_data%N_VECT_PCONT2_MIN >0.OR.
ifcontp2min > 0)
THEN
1773 ALLOCATE(ftcont2_min(3*numnod),fncont2_min(3*numnod),stat=ierr)
1774 ALLOCATE(npcont2_min(3*numnod),stat=ierr)
1775 IF (ierr/=0) GOTO 1000
1776 fncont2_min=zero
1777 ftcont2_min=zero
1778 npcont2_min=zero
1779 ELSE
1780 ALLOCATE(fncont2_min(0),stat=ierr)
1781 ALLOCATE(ftcont2_min(0),stat=ierr)
1782 ALLOCATE(npcont2_min(0),stat=ierr)
1783 END IF
1784
1785 IF(h3d_data%N_VECT_CONT2M > 0)THEN
1786 ALLOCATE(mcont2(3,numnod),stat=ierr)
1787 IF (ierr/=0) GOTO 1000
1788 mcont2=zero
1789 mcont2=zero
1790 ELSE
1791 ALLOCATE(mcont2(3,0),stat=ierr)
1792 ALLOCATE(mcont2(3,0),stat=ierr)
1793 END IF
1794
1795
1796
1800 ELSE
1801 ALLOCATE(efric(0,0))
1802 ENDIF
1805 ELSE
1806 ALLOCATE(efric_stamp(0,0))
1807 ENDIF
1808 END IF
1809
1811 ALLOCATE(efricg(numnod))
1812 IF(nintstamp/=0) THEN
1813 ALLOCATE(efricg_stamp(numnodg))
1814 ELSE
1815 ALLOCATE(efricg_stamp(0))
1816 ENDIF
1817
1818
1819
1820 ENDIF
1821
1822 IF(nintstamp/=0.AND.(anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT>0.OR.
1823 . anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0))THEN
1825 IF (ierr/=0) GOTO 1000
1828 ELSE
1830 ENDIF
1831
1832 IF(nintskidold >0)THEN
1833 IF(nintstamp/=0) THEN
1834 ALLOCATE(pskids(nintskidold,numnodg),stat=ierr)
1835 IF (ierr/=0) GOTO 1000
1836 ELSE
1837 ALLOCATE(pskids(nintskidold,numnod),stat=ierr)
1838 IF (ierr/=0) GOTO 1000
1839 ENDIF
1840 END IF
1841
1842
1843
1844 IF(nintstamp/=0)THEN
1845 ALLOCATE(
intstamp(nintstamp),stat=ierr)
1846 IF (ierr/=0)THEN
1847 CALL ancmsg(msgid=20,anmode=aninfo)
1849 ENDIF
1851 ELSE
1853 END IF
1854
1855 ALLOCATE (cfield(lfacload*nloadc),stat=ierr)
1856 IF (ierr/=0) GOTO 1000
1857 cfield = 0
1858
1859 ALLOCATE (loadp(lfacload
1860 IF (ierr/=0) GOTO 100
1861 loadp = 0
1862
1863 IF (nintloadp>0) THEN
1864
1866 IF (ierr/=0) GOTO 1000
1868
1869 ALLOCATE (
loadpinter(ninter*nloadp_hyd),stat=ierr)
1870 IF (ierr/=0) GOTO 1000
1872
1873 ALLOCATE (dgaploadint(ninter*nloadp_hyd),stat=ierr)
1874 IF (ierr/=0) GOTO 1000
1875 dgaploadint = zero
1876 ELSE
1879 ALLOCATE (dgaploadint(0),stat=ierr
1880 ENDIF
1881
1882
1883 IF (ncluster > 0) THEN
1884 ALLOCATE(
cluster(ncluster),stat=ierr)
1885 IF (ierr /= 0) THEN
1886 CALL ancmsg(msgid=20,anmode=aninfo)
1888 ENDIF
1889 ELSE
1891 ENDIF
1892
1893
1894 IF(icrack3d > 0) THEN
1896 IF (ierr/=0) GOTO 1000
1899 IF (ierr/=0) GOTO 1000
1901 ALLOCATE(
iadc_crkxfem(4*ecrkxfec+3*ecrkxfetg),stat=ierr)
1902 IF (ierr/=0) GOTO 1000
1907 IF (ierr/=0) GOTO 1000
1910 IF (ierr/=0) GOTO 1000
1913 IF (ierr/=0) GOTO 1000
1916 IF (ierr/=0) GOTO 1000
1918 ALLOCATE (
crknodiad(sprocne_crkxfe),stat=ierr)
1919 IF (ierr/=0) GOTO 1000
1921 ELSE
1927 ENDIF
1928
1929 IF(icrack3d > 0)THEN
1930 ALLOCATE(
ibordnode(sibordnode),stat=ierr)
1931 IF (ierr/=0) GOTO 1000
1933
1934 ALLOCATE(
iedgesh(siedgesh),stat=ierr)
1935 IF (ierr/=0) GOTO 1000
1937
1938 ALLOCATE (
nodedge(2*numedges),stat=ierr)
1939 IF (ierr/GOTO 1000
1941
1942 ALLOCATE (
iedge(numedges),stat=ierr)
1943 IF (ierrGOTO 1000
1945
1946 ALLOCATE (
iad_edge(siad_edge),stat=ierr)
1947 IF (ierr/=0) GOTO 1000
1949
1950 ALLOCATE (
fr_edge(sfr_edge),stat=ierr)
1951 IF (ierr/=0) GOTO 1000
1953
1954 ALLOCATE (
fr_nbedge(sfr_nbedge),stat=ierr)
1955 IF (ierr/=0) GOTO 1000
1957
1958 ALLOCATE (
nodlevxf(ncrkxfe),stat=ierr)
1959 IF (ierr/=0) GOTO 1000
1961
1963 IF (ierr/=0) GOTO 1000
1965 ELSE
1975 END IF
1976
1977 ALLOCATE (
kxig3d(skxig3d),stat=ierr)
1978 IF (ierr/=0) GOTO 1000
1980
1981 ALLOCATE (
ixig3d(sixig3d),stat=ierr)
1982 IF (ierr/=0) GOTO 1000
1984
1985 ALLOCATE (knot(sknot),stat=ierr)
1986 IF (ierr/=0) GOTO 1000
1987 knot = zero
1988
1989
1990 ALLOCATE( wige(1), stat=ierr)
1991 IF (ierr/=0) GOTO 1000
1992 wige = zero
1993
1994 ALLOCATE (knotlocpc(sknotlocpc),stat=ierr)
1995 IF (ierr/=0) GOTO 1000
1996 knotlocpc = zero
1997
1998 ALLOCATE (knotlocel(sknotlocel),stat=ierr)
1999 IF (ierr/=0) GOTO 1000
2000 knotlocel = zero
2001
2002 IF(icrack3d > 0)THEN
2003 ALLOCATE (
elcutc(2*(numelc+numeltg)),stat=ierr)
2004 IF (ierr/=0) GOTO 1000
2006
2007 ALLOCATE (
nodenr(ncrkxfe),stat=ierr)
2008 IF (ierr/=0) GOTO 1000
2010
2012 IF (ierr/=0) GOTO 1000
2014
2015 ALLOCATE (
enrtag(numnod*ienrnod*nlevmax),stat=ierr)
2016 IF (ierr/=0) GOTO 1000
2018
2019 ELSE
2024 END IF
2025
2027
2029 IF (ierr /= 0) GOTO 1000
2031
2033 IF (ierr /= 0) GOTO 1000
2035
2037 IF (ierr /= 0) GOTO 1000
2039
2041 IF (ierr /= 0) GOTO 1000
2043
2045 IF (ierr /= 0) GOTO 1000
2047
2049 IF (ierr /= 0) GOTO 1000
2051
2052
2053
2054
2055
2056
2057
2058
2060
2061 ENDIF
2062
2063 IF (
alemuscl_param%IALEMUSCL > 0 .AND. .NOT. multi_fvm%IS_USED)
THEN
2065 nbmat = trimat
2067 nbmat = 1
2068 ELSE
2069 nbmat = 0
2070 ENDIF
2071 IF(nbmat>0)THEN
2072 ALLOCATE(
alemuscl_buffer%VOLUME_FRACTION(numels + nsvois + numelq + nqvois,nbmat), stat = ierr)
2073 IF (ierr /= 0) GOTO 1000
2076 IF (ierr /= 0) GOTO 1000
2079 IF (ierr /= 0) GOTO 1000
2081 ALLOCATE(
alemuscl_buffer%GRAD(numels + nsvois + numelq + nqvois,3,nbmat), stat = ierr)
2082 IF (ierr /= 0) GOTO 1000
2084 ALLOCATE(
alemuscl_buffer%ELCENTER(numels + nsvois + numelq + nqvois,3), stat = ierr)
2085 IF (ierr /= 0) GOTO 1000
2087 ENDIF
2088 ENDIF
2089
2090
2091 IF(multi_fvm%IS_USED)THEN
2092 nbmat = multi_fvm%NBMAT
2094 nbmat = trimat
2096 nbmat = 1
2097 ELSE
2098 nbmat = 0
2099 ENDIF
2100 CALL segvar%create(nbmat)
2101
2102
2103 IF (multi_fvm%IS_USED) THEN
2106 IF (ierr /= 0) GOTO 1000
2108 ALLOCATE(multi_fvm%FVM_CONNECTIVITY%KVOIS(
ale%GLOBAL%SIELVS), stat = ierr)
2109 IF (ierr /= 0) GOTO 1000
2110 multi_fvm%FVM_CONNECTIVITY%KVOIS = 0
2111 ENDIF
2112
2113
2114
2115 IF (numply > 0) THEN
2116 ALLOCATE (
ply_info(3,numply),stat=ierr)
2117 IF (ierr /= 0) GOTO 1000
2119 ELSE
2121 END IF
2122
2123 IF (ns10e > 0) THEN
2124 ALLOCATE (
itagnd(numnod),
icnds10(3*ns10e),vnd(3*ns10e),stifnd(ns10e),stat=ierr)
2125 IF (ierr /= 0) GOTO 1000
2128 vnd(1:3*ns10e) = zero
2129 stifnd(1:ns10e) = zero
2130 IF (idamp/=0.OR.ndamp>0.OR.istat/=0) THEN
2131 ALLOCATE (vmd(3*numnod),stat=ierr)
2132 IF (ierr /= 0) GOTO 1000
2133 vmd = zero
2134 END IF
2135 IF (iparit/=0.AND.iparit/=3) THEN
2137 IF (ierr /= 0) GOTO 1000
2138 fskycnd(1:4*lcncnd)=zero
2139 ELSE
2141 END IF
2142 ELSE
2144 END IF
2145
2147 IF (ierr/=0) GOTO 1000
2150 IF (ierr/=0) GOTO 1000
2152
2153 ALLOCATE (
addcncnd(saddcncnd),stat=ierr)
2154 IF (ierr/=0) GOTO 1000
2156 ALLOCATE (
procncnd(sprocncnd),stat=ierr)
2157 IF (ierr/=0) GOTO 1000
2159 ALLOCATE (
iadcnd(siadcnd),stat=ierr)
2160 IF (ierr/=0) GOTO 1000
2162
2163 IF(n2d/=0) THEN
2164 ALLOCATE(ms_2d(numnod),stat=ierr)
2165 IF (ierr/=0) GOTO 1000
2166 ELSE
2167 ALLOCATE(ms_2d(0))
2168 ENDIF
2169
2170
2171 ALLOCATE(pinch_data%XPINCH(3,
max(1,
npinch)))
2172 pinch_data%XPINCH = zero
2173 ALLOCATE(pinch_data%DPINCH(3,
max(1,
npinch)))
2174 pinch_data%DPINCH = zero
2175 ALLOCATE(pinch_data%VPINCH(3,
max(1,
npinch)))
2176 pinch_data%VPINCH = zero
2177 ALLOCATE(pinch_data%APINCH(3,
max(1,pinch_data%SAPINCH)))
2178 pinch_data%APINCH = zero
2179 ALLOCATE(pinch_data%STIFPINCH(
max(1,pinch_data%SAPINCH)))
2180 pinch_data%STIFPINCH = zero
2181 ALLOCATE(pinch_data%MSPINCH(
max(1,pinch_data%SAPINCH)))
2182 pinch_data%MSPINCH = 0
2183
2184
2185
2187 sxcell = 0
2188
2189 IF ((iale > 0 .AND.
ale%GRID%NWALE_RST == 6) )
THEN
2190 sxcell = numels + numelq + ne_nsvois + ne_nqvois
2191 ALLOCATE(xcell(3, sxcell),stat=ierr)
2192 IF (ierr/=0) GOTO 1000
2193 ALLOCATE(xface(3, 6, numels + numelq + ne_nsvois + ne_nqvois),stat=ierr)
2194 IF (ierr/=0) GOTO 1000
2195 ENDIF
2196
2197 IF (inter18_is_variable_gap_defined) THEN
2198 sxcell = numels + numelq + ne_nsvois + ne_nqvois
2199 IF(.NOT.ALLOCATED(xcell))ALLOCATE(xcell(3, sxcell),stat=ierr)
2200 IF (ierr/=0) GOTO 1000
2201 ENDIF
2202 IF(.NOT.ALLOCATED(xcell)) ALLOCATE(xcell(3,sxcell))
2203 IF(.NOT.ALLOCATED(xface)) ALLOCATE(xface(3,6,0))
2204 xcell(1:3,1:sxcell) = zero
2205
2206 ALLOCATE (
fxbipm(nbipm,nfxbody),stat=ierr)
2207 IF (ierr/=0) GOTO 1000
2209 ALLOCATE (
fxbnod(lennod),stat=ierr)
2210 IF (ierr/=0) GOTO 1000
2212 ALLOCATE (fxbmod(6,lenmod),stat=ierr)
2213 IF (ierr/=0) GOTO 1000
2214 fxbmod = 0
2215 ALLOCATE (fxbglm(lenglm),stat=ierr)
2216 IF (ierr/=0) GOTO 1000
2217 fxbglm = 0
2218 ALLOCATE (fxbcpm(lencp),stat=ierr)
2219 IF (ierr/=0) GOTO 1000
2220 fxbcpm = 0
2221 ALLOCATE (fxbcps(lencp),stat=ierr)
2222 IF (ierr/=0) GOTO 1000
2223 fxbcps = 0
2224 ALLOCATE (fxblm(lenlm),stat=ierr)
2225 IF (ierr/=0) GOTO 1000
2226 fxblm = 0
2227 ALLOCATE (fxbfls(lenfls),stat=ierr)
2228 IF (ierr/=0) GOTO 1000
2229 fxbfls = 0
2230 ALLOCATE (fxbdls(lendls),stat
2231 IF (ierr/=0) GOTO 1000
2232 fxbdls = 0
2233 ALLOCATE (fxbdep(lenvar),stat=ierr)
2234 IF (ierr/=0) GOTO 1000
2235 fxbdep = 0
2236 ALLOCATE (fxbvit(lenvar),stat=ierr)
2237 IF (ierr/=0) GOTO 1000
2238 fxbvit = 0
2239 ALLOCATE (fxbacc(lenvar),stat=ierr)
2240 IF (ierr/=0) GOTO 1000
2241 fxbacc = 0
2242 ALLOCATE (fxbrpm(lenrpm),stat=ierr)
2243 IF (ierr/=0) GOTO 1000
2244 fxbrpm = 0
2245 ALLOCATE (
fxbelm(lenelm),stat=ierr)
2246 IF (ierr/=0) GOTO 1000
2248 ALLOCATE (fxbsig(lensig),stat=ierr)
2249 IF (ierr/=0) GOTO 1000
2250 fxbsig = 0
2251 ALLOCATE (
fxbgrvi(lengrvi),stat=ierr)
2252 IF (ierr/=0) GOTO 1000
2254 ALLOCATE (fxbgrvr(lengrvr),stat=ierr)
2255 IF (ierr/=0) GOTO 1000
2256 fxbgrvr = 0
2257
2258 ALLOCATE (
eigipm(neipm,neig),stat=ierr)
2259 IF (ierrGOTO 1000
2261 ALLOCATE (eigrpm(nerpm,neig),stat=ierr)
2262 IF (ierr/=0) GOTO 1000
2263 eigrpm = 0
2264 ALLOCATE (
eigibuf(leibuf),stat=ierr)
2265 IF (ierr/=0) GOTO 1000
2267
2268 ALLOCATE (
slipring(nslipring),stat=ierr
2269 IF (ierr/=0) GOTO 1000
2270
2271 ALLOCATE (
retractor(nretractor),stat=ierr)
2272 IF (ierr/=0) GOTO 1000
2273
2274 IF (n_anchor_remote > 0) THEN
2276 IF (ierr/=0) GOTO 1000
2278 IF (ierr/=0) GOTO 1000
2279 ENDIF
2280 IF (n_anchor_remote_send > 0) THEN
2282 IF (ierr/=0) GOTO 1000
2284 IF (ierr/=0) GOTO 1000
2285 ENDIF
2286
2287 ALLOCATE (tm_dis(
lmax_dis*numnod),stat=ierr)
2288 IF (ierr/=0) GOTO 1000
2289ALLOCATE
2290 IF (ierr/=0) GOTO 1000
2291 ALLOCATE (tm_nsig1(
lmax_nsig*numnod),stat=ierr)
2292 IF (ierr/=0) GOTO 1000
2293 ALLOCATE (tm_nstra1(
lmax_nstra*numnod),stat=ierr)
2294 IF (ierr/=0) GOTO 1000
2295 ALLOCATE (tm_nsig3(
lmax_nsig*numnod),stat=ierr)
2296 IF (ierr/=0) GOTO 1000
2297 ALLOCATE (tm_nstra3(
lmax_nstra*numnod),stat=ierr)
2298 IF (ierr/=0) GOTO 1000
2300 ALLOCATE (tm_pdis(numnod),stat=ierr)
2301 IF (ierr/=0) GOTO 1000
2302 ENDIF
2304 ALLOCATE (tm_pvel(numnod),stat=ierr)
2305 IF (ierr/=0) GOTO 1000
2306 ENDIF
2308 ALLOCATE (tm_pnsig(2*numnod),stat=ierr)
2309 IF (ierr/=0) GOTO 1000
2310 ENDIF
2312 ALLOCATE (tm_pnstra(2*numnod),stat=ierr)
2313 IF (ierr/=0) GOTO 1000
2314 ENDIF
2316 ALLOCATE(gpstmp(6*numnod),p2tmp(2*numnod),stat=ierr)
2317 IF (ierr/=0) GOTO 1000
2318 END IF
2319 IF (nconld > 0) THEN
2320 ALLOCATE (dpl0cld(6*nconld),stat=ierr)
2321 IF (ierr/=0) GOTO 1000
2322 dpl0cld = zero
2323 ALLOCATE (vel0cld(6*nconld),stat=ierr)
2324 IF (ierr/=0) GOTO 1000
2325 vel0cld = zero
2326 ELSE
2327 ALLOCATE(dpl0cld(0),vel0cld(0))
2328 END IF
2329
2331 IF (ierr/=0) GOTO 1000
2333 IF (ierr/=0) GOTO 1000
2334
2335
2336 CALL my_alloc(skews%SKEW,skews%N_SKEW_VAR,skews%TOTAL_SKEW_NUMBER)
2337
2338
2339
2340 IF(debug(6)==1)THEN
2342 END IF
2343
2344 RETURN
2345 1000 CONTINUE
2346 ierror = ierror + ierr
2347 ierr = ierr0
2348 IF (ierror/=0) THEN
2349 CALL ancmsg(msgid=20,anmode=aninfo)
2351 END IF
2352
subroutine multi_allocate(multi_fvm, h3d_data)
type(fani_cell_) fani_cell
type(alefvm_buffer_), target alefvm_buffer
type(alefvm_param_), target alefvm_param
type(alemuscl_param_) alemuscl_param
type(alemuscl_buffer_) alemuscl_buffer
integer, dimension(:), allocatable id_damp_vrel
integer, dimension(:), allocatable fr_damp_vrel
integer, dimension(:), pointer iad_cndm1
integer, dimension(:), pointer fr_nbcccnd1
integer, dimension(:), pointer iadcnd
integer, dimension(:), pointer fr_cndm
integer, dimension(:), pointer fr_cndm1
integer, dimension(:), pointer itagnd
integer, dimension(:), pointer procncnd
integer, dimension(:), pointer icnds10
integer, dimension(:), pointer addcncnd
integer, dimension(:), pointer iad_cndm
integer, dimension(:), pointer fr_nbcccnd
integer, dimension(:), allocatable eigibuf
integer, dimension(:,:), allocatable eigipm
integer, dimension(:,:), allocatable fxbipm
integer, dimension(:), allocatable fxbnod
integer, dimension(:), allocatable fxbelm
integer, dimension(:), allocatable fxbgrvi
type(intstamp_data), dimension(:), allocatable intstamp
integer, dimension(:), allocatable indexcont
integer, dimension(:), allocatable tagcont
integer, dimension(:), allocatable poin_ump
integer, dimension(:), allocatable iconx
integer, dimension(:), allocatable, target igrv
integer, dimension(:), allocatable fr_sec
integer, dimension(:), allocatable iad_rby
integer, dimension(:), allocatable id_global_vois
integer, dimension(:), allocatable fr_nbedge
integer, dimension(:), allocatable iadrbmk
integer, dimension(:), allocatable fr_mad
integer, dimension(:), allocatable ibcv
integer, dimension(:), allocatable lagbuf
integer, dimension(:), allocatable ixx
integer, dimension(:), allocatable, target lpby
integer, dimension(:), allocatable fr_rl
integer, dimension(:), allocatable, target ixs
integer, dimension(:), allocatable iad_rbym2
integer, dimension(:), allocatable iad_i2m
integer, dimension(:), allocatable iad_cut
integer, dimension(:), allocatable lgrav
integer, dimension(:), allocatable, target npby
integer, dimension(:), allocatable kxig3d
integer, dimension(:), allocatable iadmv3
integer, dimension(:), allocatable face_vois
integer, dimension(:), pointer iframe
integer, dimension(:), pointer madfail
integer, dimension(:), allocatable lesdvois
integer, dimension(:), allocatable lnrcvois
integer, dimension(:), allocatable nodenr
integer, dimension(:), allocatable nativ0_sms
integer, dimension(:), allocatable newfront
integer, dimension(:), allocatable iadc_crkxfem
integer, dimension(:), allocatable nodpor
integer, dimension(:), allocatable ilink
integer, dimension(:), allocatable iadrcp_pxfem
integer, dimension(:), allocatable ibc_ply
integer, dimension(:), allocatable llink
integer, dimension(:), allocatable madclnod
integer, dimension(:,:), allocatable ipadmesh
integer, dimension(:), allocatable lbvel
integer, dimension(:), allocatable lprtsph
integer, dimension(:), allocatable nbsdvois
integer, dimension(:), allocatable lnodpor
integer, dimension(:), allocatable ibcr
integer, dimension(:), allocatable ne_nercvois
integer, dimension(:), allocatable iadmv2
integer, dimension(:), allocatable neflsw
integer, dimension(:), allocatable ixig3d
integer, dimension(:), allocatable linale
type(cluster_), dimension(:), allocatable cluster
integer, dimension(:), allocatable iactiv
integer, dimension(:), allocatable crknodiad
integer, dimension(:), allocatable ibcslag
integer, dimension(:), allocatable icodt_ply
integer, dimension(:), allocatable ibufssg_io
integer, dimension(:,:), allocatable sh4tree
integer, dimension(:), allocatable fr_lagf
integer, dimension(:), allocatable ispsym
integer, dimension(:), allocatable sh4trim
integer, dimension(:), allocatable addcsrect
integer, dimension(:), allocatable ipm
integer, dimension(:), allocatable, target ipart
integer, dimension(:), allocatable iadwal
integer, dimension(:), allocatable fr_nor
integer, dimension(:), allocatable isphio
integer, dimension(:), allocatable fr_i18
integer, dimension(:), allocatable, target ipari
integer, dimension(:), allocatable fr_rbym2
integer, dimension(:), allocatable igaup
integer, dimension(:), allocatable icodrbym
integer, dimension(:), allocatable iad_rbm
integer, dimension(:), allocatable iskew_ply
integer, dimension(:), allocatable nnflsw
integer, dimension(:), allocatable nercvois
integer, dimension(:), allocatable ispcond
integer, dimension(:), allocatable ibordnode
integer, dimension(:), allocatable sh3trim
integer, dimension(:), allocatable iadrbk
integer, dimension(:), allocatable iecran
integer, dimension(:), allocatable, target iedgesh
integer, dimension(:), allocatable fr_i2m
integer, dimension(:), allocatable ixt
integer, dimension(:), allocatable lnlink
integer, dimension(:), allocatable ibftemp
integer, dimension(:), allocatable ibfv
integer, dimension(:), allocatable iaccp
integer, dimension(:), allocatable, target iel_crkxfem
integer, dimension(:), allocatable iadi18
integer, dimension(:), allocatable lsegcom
integer, dimension(:), allocatable iskwp_l
integer, dimension(:), allocatable inod_pxfem
integer, dimension(:), allocatable dd_r2r_elem
integer, dimension(:), allocatable kloadpinter
integer, dimension(:), allocatable ixr
integer, dimension(:,:), allocatable sh3tree
integer, dimension(:), allocatable lonfsph
integer, dimension(:), pointer madnod
integer, dimension(:), allocatable madidx
integer, dimension(:), allocatable iexlnk
integer, dimension(:), allocatable, target ixtg
integer, dimension(:), pointer lpbyl
integer, dimension(:), allocatable nnlink
integer, dimension(:), allocatable, target ibcl
integer, dimension(:), allocatable adsky_crkxfem
integer, dimension(:), pointer madprt
integer, dimension(:), allocatable fr_rbym
integer, dimension(:), allocatable monvol
integer, dimension(:), allocatable ifill
integer, dimension(:), allocatable kxfenod2elc
integer, dimension(:), allocatable iskwp
integer, dimension(:), allocatable isensp
integer, dimension(:), allocatable fr_rbe2
integer, dimension(:), allocatable irbe2
integer, dimension(:), allocatable inod_crkxfem
integer, dimension(:), allocatable nporgeo
integer, dimension(:), allocatable procne_crkxfem
integer, dimension(:), allocatable fr_rbm2
integer, dimension(:), allocatable iadsdp_pxfem
integer, dimension(:), allocatable kxsp
integer, dimension(:), allocatable neth
integer, dimension(:), allocatable enrtag
integer, dimension(:), allocatable nodlevxf
integer, dimension(:), allocatable fr_wall
integer, dimension(:), allocatable loadpinter
integer, dimension(:), allocatable elcutc
integer, dimension(:), allocatable fr_ll
integer, dimension(:), allocatable iad_sec
integer, dimension(:), allocatable nsensp
integer, dimension(:), allocatable dd_iad
integer, dimension(:), allocatable gjbufi
integer, dimension(:), pointer madsh3
integer, dimension(:), allocatable icut
integer, dimension(:), allocatable fr_cj
integer, dimension(:), allocatable, target iskwn
integer, dimension(:), allocatable nesdvois
integer, dimension(:), allocatable cne_crkxfem
integer, dimension(:), allocatable, target iloadp
integer, dimension(:), allocatable iadcj
integer, dimension(:), allocatable nprw
integer, dimension(:), allocatable iadsdp_crkxfem
integer, dimension(:), allocatable ngaup
integer, dimension(:), allocatable lnrbym
integer, dimension(:), allocatable nod2sp
integer, dimension(:), allocatable adsky_pxfem
integer, dimension(:), allocatable nodglobxfe
integer, dimension(:), allocatable weight_rm
integer, dimension(:), allocatable ixp
integer, dimension(:), allocatable laccelm
integer, dimension(:), allocatable, target nom_opt
integer, dimension(:), allocatable iad_rbe2
double precision, dimension(:), allocatable bufgeo
integer, dimension(:), pointer madsol
integer, dimension(:), allocatable fasolfr
integer, dimension(:), allocatable iadi2
integer, dimension(:), allocatable, target npc
integer, dimension(:), allocatable igeo
integer, dimension(:), allocatable, target ibmpc
integer, dimension(:), allocatable ixtg1
integer, dimension(:), allocatable fr_mv
integer, dimension(:), allocatable ims
integer, dimension(:), allocatable fr_edge
integer, dimension(:), allocatable lercvois
integer, dimension(:), allocatable iadrl
integer, dimension(:), allocatable addcni2
integer, dimension(:), allocatable lbcscyc
integer, dimension(:), allocatable fr_rby2
integer, dimension(:), allocatable iad_frnor
integer, dimension(:), allocatable iadc_pxfem
integer, dimension(:), allocatable ibvel
integer, dimension(:), allocatable iadrbm
integer, dimension(:), allocatable lrivet
integer, dimension(:), pointer npbyl
integer, dimension(:), allocatable iad_edge
integer, dimension(:), allocatable, target icfield
integer, dimension(:), allocatable kinet
integer, dimension(:), allocatable lgauge
integer, dimension(:), allocatable nstrf
integer, dimension(:), allocatable ibcscyc
integer, dimension(:), allocatable procnor
integer, dimension(:), allocatable tag_skins6
integer, dimension(:), allocatable fr_rby
integer, dimension(:), allocatable iad_rbm2
integer, dimension(:), allocatable irbym
integer, dimension(:,:), allocatable ixsp
integer, dimension(:), allocatable iadrcp_crkxfem
integer, dimension(:), pointer madsh4
integer, dimension(:), allocatable iparg
integer, dimension(:), allocatable ixq
integer, dimension(:), allocatable iedge
integer, dimension(:), allocatable ibfflux
integer, dimension(:), allocatable nodedge
integer, dimension(:), allocatable ilas
integer, dimension(:), allocatable iad_rby2
integer, dimension(:), allocatable icode_ply
integer, dimension(:), allocatable lloadp
integer, dimension(:), allocatable fr_cut
integer, dimension(:), allocatable segquadfr
integer, dimension(:), allocatable procni2
integer, dimension(:), allocatable lcfield
integer, dimension(:), allocatable kxx
integer, dimension(:), allocatable nskwp
integer, dimension(:), allocatable ne_nesdvois
integer, dimension(:), allocatable iad_rbym
integer, dimension(:), allocatable nom_sect
integer, dimension(:), allocatable lprw
integer, dimension(:), allocatable npsegcom
integer, dimension(:), allocatable rg_cut
integer, dimension(:), allocatable iel_pxfem
integer, dimension(:), allocatable lrbe2
integer, dimension(:), allocatable fr_rbm
integer, dimension(:), allocatable llagf
integer, dimension(:), allocatable ixc
integer, dimension(:), allocatable nbrcvois
integer, dimension(:), allocatable ljoint
integer, dimension(:), allocatable lnsdvois
integer, dimension(:,:), allocatable tab_ump
integer, dimension(:), allocatable naccp
integer, dimension(:), allocatable procne_pxfem
integer, dimension(:), allocatable ithvar
integer, dimension(:), allocatable dd_r2r
integer, dimension(:), allocatable iadll
integer, dimension(:), allocatable icontact
integer, dimension(:), allocatable ipart_state
type(retractor_struct), dimension(:), allocatable retractor
type(seatbelt_remote_nodes_struct) anchor_remote_send
type(seatbelt_remote_nodes_struct) anchor_remote
type(slipring_struct), dimension(:), allocatable slipring
integer, dimension(:), allocatable sph2sol
integer, dimension(:), allocatable irst
integer, dimension(:), allocatable sol2sph_typ
integer, dimension(:), allocatable sol2sph
integer, dimension(:,:), allocatable ply_info
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)