74
75
76
87 USE multi_fvm_mod
98 USE output_mod
101 USE interfaces_mod
102 USE skew_mod
103 USE my_alloc_mod
104 use glob_therm_mod
105 USE connectivity_mod
106 USE nodal_arrays_mod
107 use rbe3_mod
108 use rwall_mod
109 use element_mod , only : nixs,nixq,nixc,nixt,nixr,nixp,nixtg
110
111
112
113#include "implicit_f.inc"
114
115
116
117#include "com01_c.inc"
118#include "com04_c.inc"
119#include "com_xfem1.inc"
120#include "param_c.inc"
121#include "remesh_c.inc"
122#include "scr06_c.inc"
123#include "scr07_c.inc"
124#include "scr14_c.inc"
125#include "scr16_c.inc"
126#include "sms_c.inc"
127#include "tabsiz_c.inc"
128#include "task_c.inc"
129#include "couple_c.inc"
130#include "warn_c.inc"
131#include "parit_c.inc"
132#include "intstamp_c.inc"
133#include "spmd_c.inc"
134#include "sphcom.inc"
135#include "fxbcom.inc"
136#include "eigcom.inc"
137#include "inter18.inc"
138#include "stati_c.inc"
139
140 type(output_), INTENT(INOUT) :: OUTPUT
141 TYPE(connectivity_) :: ELEMENT
142 TYPE(NODAL_ARRAYS_) :: NODES
143 TYPE(MULTI_FVM_STRUCT) :: MULTI_FVM
144 TYPE(H3D_DATABASE) :: H3D_DATA
145 TYPE(PINCH) :: PINCH_DATA
146 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
147 TYPE(t_SEGVAR) :: SEGVAR
148 TYPE(INTERFACES_) :: INTERFACES
149 TYPE(SKEW_) , INTENT(INOUT) :: SKEWS
150 type(glob_therm_) ,intent(inout) :: GLOB_THERM
151 type(rbe3_) ,intent(inout) :: RBE3
152 type(rwall_),intent(inout) :: RWALL
153
154
155
156 INTEGER IERR0,IERROR, NBMAT,ISIZ1,ISIZ2
157
158 ierror = 0
159
160 ierr0 = ierr
161
162
163
164
165 CALL allocate_nodal_arrays(nodes, numnod, nthread, iroddl, iparit,
166 . isecut, iisrots, impose_dr, idrot, nrcvvois, sicodt, glob_therm%ITHERM_FE)
167
168
169
170
171
172
173
174
175
176
177 ALLOCATE (
iskwn(siskwn),stat=ierr)
178 IF (ierr/=0) GOTO 1000
180
181 ALLOCATE (
iframe(siframe),stat=ierr)
182 IF (ierr/=0) GOTO 1000
184
185 ALLOCATE (
neth(sneth),stat=ierr)
186 IF (ierr/=0) GOTO 1000
188
189 ALLOCATE (
ibcslag(sibcslag),stat=ierr)
190 IF (ierr/=0) GOTO 1000
192
193 ALLOCATE (
ipart(sipart),stat=ierr)
194 IF (ierr/=0) GOTO 1000
196
197 ALLOCATE (
nom_opt(snom_opt),stat=ierr)
198 IF (ierr/=0) GOTO 1000
200
201 ALLOCATE (
npc(snpc),stat=ierr)
202 IF (ierr/=0) GOTO 1000
204
205 ALLOCATE (
ixtg(sixtg),stat=ierr)
206 IF (ierr/=0) GOTO 1000
208
209 ALLOCATE (
ixtg1(sixtg1),stat=ierr)
210 IF (ierr/=0) GOTO 1000
212
213 ALLOCATE (
ixs(sixs),stat=ierr)
214 IF (ierr/=0) GOTO 1000
216
217 ALLOCATE (
ixq(sixq),stat=ierr)
218 IF (ierr/=0) GOTO 1000
220
221 ALLOCATE (
ixc(sixc),stat=ierr)
222 ALLOCATE(element%shell%ixc(nixc,sixc/nixc))
223 ALLOCATE(element%shell%nodes(4,sixc/nixc))
224 ALLOCATE(element%shell%pid(sixc/nixc))
225 ALLOCATE(element%shell%matid(sixc/nixc))
226 ALLOCATE(element%shell%user_id(sixc/nixc))
228 element%shell%ixc = 0
229
230 ALLOCATE (
ixt(sixt),stat=ierr)
231 IF (ierr/=0) GOTO 1000
233
234 ALLOCATE (
ixp(sixp),stat=ierr)
235 IF (ierr/=0) GOTO 1000
237
238 ALLOCATE (
ixr(sixr),stat=ierr)
239 IF (ierr/=0) GOTO 1000
241
242
243
244
245
246
247
248
249
250 ALLOCATE (
gjbufi(sgjbufi),stat=ierr)
251 IF (ierr/=0) GOTO 1000
253
254 ALLOCATE (ale_connectivity%NALE(
ale%GLOBAL%SNALE),stat=ierr)
255 IF (ierr/=0) GOTO 1000
256 ale_connectivity%NALE(1:
ale%GLOBAL%SNALE) = 0
257
259 IF (ierr/=0) GOTO 1000
261
263 IF (ierr/=0) GOTO 1000
265
266 ALLOCATE (
ifill(sifill),stat=ierr)
267 IF (ierr/=0) GOTO 1000
269
270 ALLOCATE (
ims(sims),stat=ierr)
271 IF (ierr/=0) GOTO 1000
273
274 ALLOCATE (
kxx(skxx),stat=ierr)
275 IF (ierr/=0) GOTO 1000
277
278 ALLOCATE (
ixx(sixx),stat=ierr)
279 IF (ierr/=0) GOTO 1000
281
282 ALLOCATE (
kxsp(skxsp),stat=ierr)
283 IF (ierr/=0) GOTO 1000
285
286 ALLOCATE (
ixsp(sixsp),stat=ierr)
287 IF (ierr/=0) GOTO 1000
289
290 ALLOCATE (
nod2sp(snod2sp),stat=ierr)
291 IF (ierr/=0) GOTO 1000
293
294 ALLOCATE (
ispsym(sispsym),stat=ierr)
295 IF (ierr/=0) GOTO 1000
297
298 ALLOCATE (
ispcond(sispcond),stat=ierr)
299 IF (ierr/=0) GOTO 1000
301
302 ALLOCATE (
isphio(sisphio),stat=ierr)
303 IF (ierr/=0) GOTO 1000
305
306 ALLOCATE (
lprtsph(slprtsph),stat=ierr)
307 IF (ierr/=0) GOTO 1000
309
310 ALLOCATE (
lonfsph(slonfsph),stat=ierr)
311 IF (ierr/=0) GOTO 1000
313
314 ALLOCATE (
ibcl(sibcl),stat=ierr)
315 IF (ierr/=0) GOTO 1000
317
318 ALLOCATE (
ibfv(sibfv),stat=ierr)
319 IF (ierr/=0) GOTO 1000
321
322 ALLOCATE (
iecran(siecran),stat=ierr)
323 IF (ierr/=0) GOTO 1000
325
326 ALLOCATE (
ilas(silas),stat=ierr)
327 IF (ierr/=0) GOTO 1000
329
330 ALLOCATE (
laccelm(slaccelm),stat=ierr)
331 IF (ierr/=0) GOTO 1000
333
334 ALLOCATE (
lgauge(3*nbgauge),stat=ierr)
335 IF (ierr/=0) GOTO 1000
337
338 ALLOCATE (
nnlink(snnlink),stat=ierr)
339 IF (ierr/=0) GOTO 1000
341
342 ALLOCATE (
lnlink(slnlink),stat=ierr)
343 IF (ierr/=0) GOTO 1000
345
346 ALLOCATE (
iparg(siparg),stat=ierr)
347 IF (ierr/=0) GOTO 1000
349
350 ALLOCATE (
dd_iad(sdd_iad),stat=ierr)
351 IF (ierr/=0) GOTO 1000
353
354 ALLOCATE (
igrv(sigrv),stat=ierr)
355 IF (ierr/=0) GOTO 1000
357
358 ALLOCATE (
lgrav(slgrav),stat=ierr)
359 IF (ierr/=0) GOTO 1000
361
362 ALLOCATE (
iexlnk(siexlnk),stat=ierr)
363 IF (ierr/=0) GOTO 1
365
366 ALLOCATE (
ibvel(sibvel),stat=ierr)
367 IF (ierr/=0) GOTO 1000
369
370 ALLOCATE (
lbvel(slbvel),stat=ierr)
371 IF (ierr/=0) GOTO 1000
373
374 ALLOCATE (
irbe2(sirbe2),stat=ierr)
375 IF (ierr/=0) GOTO 1000
377
378 ALLOCATE (
lrbe2(slrbe2),stat=ierr)
379 IF (ierr/=0) GOTO 1000
381
382 call allocate_rbe3(rbe3, nspmd)
383
384 ALLOCATE (
iactiv(siactiv),stat=ierr)
385 IF (ierr/=0) GOTO 1000
387
388 ALLOCATE (
ibmpc(sibmpc),stat=ierr)
389 IF (ierr/=0) GOTO 1000
391
392 ALLOCATE (
fasolfr(sfasolfr),stat=ierr)
393 IF (ierr/=0) GOTO 1000
395
396 ALLOCATE (
segquadfr(ssegquadfr),stat=ierr)
397 IF (ierr/=0) GOTO 1000
399
400 ALLOCATE (
ipari(sipari),stat=ierr
401 IF (ierr/=0) GOTO 1000
403
404
405
406
407
408
409 ALLOCATE (interfaces%INTBUF_TAB(ninter),stat=ierr)
410
411
412
413 ALLOCATE(interfaces%INTBUF_FRIC_TAB(ninterfric),stat=ierr)
414 IF (ierr/=0) GOTO 1000
415
416
417
418
419 call allocate_rwall(rwall, nspmd)
420
421
422 ALLOCATE (
iconx(siconx),stat=ierr)
423 IF (ierr/=0) GOTO 1000
425
426 ALLOCATE (
madprt(smadprt),stat=ierr)
427 IF (ierr/=0) GOTO 1000
429
430 ALLOCATE (
madsh4(smadsh4),stat=ierr)
431 IF (ierr/=0) GOTO 1000
433
434 ALLOCATE (
madsh3(smadsh3),stat=ierr)
435 IF (ierr/=0) GOTO 1000
437
438 ALLOCATE (
madsol(smadsol),stat=ierr)
439 IF (ierr/=0) GOTO 1000
441
442 ALLOCATE (
madnod(smadnod),stat=ierr)
443 IF (ierr/=0) GOTO 1000
445
446 ALLOCATE (
madclnod(madclnods),stat=ierr)
447 IF (ierr/=0) GOTO 1000
449
450 ALLOCATE (
madfail(smadfail),stat=ierr)
451 IF (ierr/=0) GOTO 1000
453 ALLOCATE (
madidx(smadidx),stat=ierr)
454 IF (ierr/=0) GOTO 1000
456
457
458 ALLOCATE (
npby(snpby),stat=ierr)
459 IF (ierr/=0) GOTO 1000
461
462 ALLOCATE (
npbyl(snpbyl),stat=ierr)
463 IF (ierr/=0) GOTO 1000
465
466 ALLOCATE (
lpby(slpby),stat=ierr)
467 IF (ierr/=0) GOTO 1000
469
470 ALLOCATE (
lpbyl(slpbyl),stat=ierr)
471 IF (ierr/=0) GOTO 1000
473
474 ALLOCATE (
lrivet(slrivet),stat=ierr)
475 IF (ierr/=0) GOTO 1000
477
478 ALLOCATE (
nstrf(snstrf),stat=ierr)
479 IF (ierr/=0) GOTO 1000
481
483 IF (ierr/=0) GOTO 1000
485
486 ALLOCATE (
ljoint(sljoint),stat=ierr)
487 IF (ierr/=0) GOTO 1000
489
490 ALLOCATE (
nodpor(snodpor),stat=ierr)
491 IF (ierr/=0) GOTO 1000
493
494 ALLOCATE (
monvol(smonvol),stat=ierr)
495 IF (ierr/=0) GOTO 1000
497
498 ALLOCATE (
lagbuf(slagbuf),stat=ierr)
499 IF (ierr/=0) GOTO 1000
501
502 ALLOCATE (rthbuf(srthbuf),stat=ierr)
503 IF (ierr/=0) GOTO 1000
504 rthbuf = zero
505
506 ALLOCATE (nodes%BOUNDARY_ADD(2,siad_elem/2),stat=ierr)
507 IF (ierr/=0) GOTO 1000
508 nodes%BOUNDARY_ADD = 0
509
510 ALLOCATE (
iad_rby(siad_rby),stat=ierr)
511 IF (ierr/=0) GOTO 1000
513
514 ALLOCATE (
iad_rby2(siad_rby2),stat=ierr)
515 IF (ierr/=0) GOTO 1000
517
518 ALLOCATE (
iad_i2m(siad_i2m),stat=ierr)
519 IF (ierr/=0) GOTO 1000
521
522 ALLOCATE (
iadcj(siadcj),stat=ierr)
523 IF (ierr/=0) GOTO 1000
525
526 ALLOCATE (
iad_rbm(siad_rbm),stat=ierr)
527 IF (ierr/=0) GOTO 1000
529
530 ALLOCATE (
iad_rbm2(siad_rbm2),stat=ierr)
531 IF (ierr/=0) GOTO 1000
533
534 ALLOCATE (
iad_rbe2(siad_rbe2),stat=ierr)
535 IF (ierr/=0) GOTO 1000
537
538 ALLOCATE (
iad_sec(siad_sec),stat=ierr)
539 IF (ierr/=0) GOTO 1000
541
542 ALLOCATE (
iad_cut(siad_cut),stat=ierr)
543 IF (ierr/=0) GOTO 1000
545
546 siad_rbym = siad_rby
547 siad_rbym2 = siad_rby2
548 ALLOCATE (
iad_rbym(siad_rbym),stat=ierr)
549 IF (ierr/=0) GOTO 1000
551
552 ALLOCATE (
iad_rbym2(siad_rbym2),stat=ierr)
553 IF (ierr/=0) GOTO 1000
555
557 IF (ierr/=0) GOTO 1000
559
560 ALLOCATE (
iad_frnor((nspmd+1)*ninter25),stat=ierr)
561 IF (ierr/=0) GOTO 1000
563
564 ALLOCATE (interfaces%SPMD_ARRAYS%FR_EDG(2*nbddedgt),stat=ierr)
565 IF (ierr/=0) GOTO 1000
566 interfaces%SPMD_ARRAYS%FR_EDG = 0
567
568 ALLOCATE (interfaces%SPMD_ARRAYS%IAD_FREDG((nspmd+1)*ninter25),stat=ierr)
569 IF (ierr/=0) GOTO 1000
570 interfaces%SPMD_ARRAYS%IAD_FREDG = 0
571
572 ALLOCATE (
fr_nor(nbddnort),stat=ierr)
573 IF (ierr/=0) GOTO 1000
575
576 ALLOCATE (
procnor(nbccnor),stat=ierr)
577 IF (ierr/=0) GOTO 1000
579 nodes%BOUNDARY_SIZE = sfr_elem
580 ALLOCATE (nodes%BOUNDARY(sfr_elem),stat=ierr)
581 IF (ierr/=0) GOTO 1000
582 nodes%BOUNDARY = 0
583
584 ALLOCATE (
fr_rby(sfr_rby),stat=ierr)
585 IF (ierr/=0) GOTO 1000
587
588 ALLOCATE (
fr_rby2(sfr_rby2),stat=ierr)
589 IF (ierr/=0) GOTO 1000
591
592 ALLOCATE (
fr_i2m(sfr_i2m),stat=ierr)
593 IF (ierr/=0) GOTO 1000
595
596 ALLOCATE (
fr_mv(sfr_mv),stat=ierr)
597 IF (ierr/=0) GOTO 1000
599
600 ALLOCATE (
fr_ll(sfr_ll),stat=ierr)
601 IF (ierr/=0) GOTO 1000
603
604 ALLOCATE (
fr_cj(sfr_cj),stat=ierr)
605 IF (ierr/=0) GOTO 1000
607
608 ALLOCATE (
fr_rbm(sfr_rbm),stat=ierr)
609 IF (ierr/=0) GOTO 1000
611
612 ALLOCATE (
fr_rbm2(sfr_rbm2),stat=ierr)
613 IF (ierr/=0) GOTO 1000
615
616 ALLOCATE (
fr_rbe2(sfr_rbe2),stat=ierr)
617 IF (ierr/=0) GOTO 1000
619
620 ALLOCATE (
fr_sec(sfr_sec),stat=ierr)
621 IF (ierr/=0) GOTO 1000
623
624 ALLOCATE (
fr_cut(sfr_cut),stat=ierr)
625 IF (ierr/=0) GOTO 1000
627
628 ALLOCATE (
rg_cut(srg_cut),stat=ierr)
629 IF (ierr/=0) GOTO 1000
631
632 ALLOCATE (
fr_mad(sfr_mad),stat=ierr)
633 IF (ierr/=0) GOTO 1000
635
636 ALLOCATE (
fr_i18(sfr_i18),stat=ierr)
637 IF (ierr/=0) GOTO 1000
639
640 ALLOCATE (
dd_r2r(sdd_r2r),stat=ierr)
641 IF (ierr/=0) GOTO 1000
643
645 IF (ierr/=0) GOTO 1000
647
648 ALLOCATE (
fr_rbym(sfr_rbym),stat=ierr)
649 IF (ierr/=0) GOTO 1000
651
652 ALLOCATE (
fr_rbym2(sfr_rbym2),stat=ierr)
653 IF (ierr/=0) GOTO 1000
655 ALLOCATE (
iskwp(siskwp),stat=ierr)
656 IF (ierr/=0) GOTO 1000
658
659 ALLOCATE (
nskwp(snskwp),stat=ierr)
660 IF (ierr/=0) GOTO 1000
662
663 ALLOCATE (
iskwp_l(siskwp),stat=ierr)
664 IF (ierr/=0) GOTO 1000
666
667 ALLOCATE (
isensp(sisensp),stat=ierr)
668 IF (ierr/=0) GOTO 1000
670
671 ALLOCATE (
nsensp(snsensp),stat=ierr)
672 IF (ierr/=0) GOTO 1000
674
675 ALLOCATE (
iaccp(siaccp),stat=ierr)
676 IF (ierr/=0) GOTO 1000
678
679 ALLOCATE (
naccp(snaccp),stat=ierr)
680 IF (ierr/=0) GOTO 1000
682
683 ALLOCATE (
igaup(sigaup),stat=ierr)
684 IF (ierr/=0) GOTO 1000
686
687 ALLOCATE (
ngaup(sngaup),stat=ierr)
688 IF (ierr/=0) GOTO 1000
690
691 ALLOCATE (
fr_lagf(sfr_lagf),stat=ierr)
692 IF (ierr/=0) GOTO 1000
694
695
696
697
698
699
700
701
702
703
704
705
706
707 ALLOCATE (
newfront(snewfront),stat=ierr)
708 IF (ierr/=0) GOTO 1000
710
711
712
713
714
715 ALLOCATE (
nbrcvois(snbrcvois),stat=ierr)
716 IF (ierr/=0) GOTO 1000
718
719 ALLOCATE (
lnrcvois(slnrcvois),stat=ierr)
720 IF (ierr/=0) GOTO 1000
722
723 ALLOCATE (
nbsdvois(snbsdvois),stat=ierr)
724 IF (ierr/=0) GOTO 1000
726
727 ALLOCATE (
lnsdvois(slnsdvois),stat=ierr)
728 IF (ierr/=0) GOTO 1000
730
731 ALLOCATE (
nercvois(snercvois),stat=ierr)
732 IF (ierr/=0) GOTO 1000
734
735 ALLOCATE (
lercvois(slercvois),stat=ierr)
736 IF (ierr/=0) GOTO 1000
738
739 ALLOCATE (
nesdvois(snesdvois),stat=ierr)
740 IF (ierr/=0) GOTO 1000
742
743 ALLOCATE (
lesdvois(slesdvois),stat=ierr)
744 IF (ierr/=0) GOTO 1000
746
747 IF(ale_connectivity%has_ne_connect)THEN
749 IF (ierr/=0) GOTO 1000
751
753 IF (ierr/=0) GOTO 1000
755 ELSE
758 ENDIF
759
760 ALLOCATE (
npsegcom(snpsegcom),stat=ierr)
761 IF (ierr/=0) GOTO 1000
763
764 ALLOCATE (
lsegcom(slsegcom),stat=ierr)
765 IF (ierr/=0) GOTO 1000
767
768 ALLOCATE (
nporgeo(snporgeo),stat=ierr)
769 IF (ierr/=0) GOTO 1000
771
772 ALLOCATE (
lnodpor(slnodpor),stat=ierr)
773 IF (ierr/=0) GOTO 1000
775
776 ALLOCATE (
llagf(sllagf),stat=ierr)
777 IF (ierr/=0) GOTO 1000
779
780
781
782
783
784
785
786
787
788
789 ALLOCATE (interfaces%PON%ISKY(sisky),stat=ierr)
790 IF (ierr/=0) GOTO 1000
791 interfaces%PON%ISKY = 0
792
793 ALLOCATE (element%PON%ADSKY(sadsky),stat=ierr)
794 IF (ierr/=0) GOTO 1000
795 element%PON%ADSKY = 0
796 element%PON%SADSKY = sadsky
797
798 ALLOCATE (element%PON%PROCNE(sprocne),stat=ierr)
799 IF (ierr/=0) GOTO 1000
800 element%PON%PROCNE = 0
801
802 ALLOCATE (
addcni2(saddcni2),stat=ierr)
803 IF (ierr/=0) GOTO 1000
805
806 ALLOCATE (
procni2(sprocni2),stat=ierr)
807 IF (ierr/=0) GOTO 1000
809
810 ALLOCATE (element%PON%IADSDP(siadsdp),stat=ierr)
811 IF (ierr/=0) GOTO 1000
812 element%PON%IADSDP = 0
813
814 ALLOCATE (element%PON%IADRCP(siadrcp),stat=ierr)
815 IF (ierr/=0) GOTO 1000
816 element%PON%IADRCP = 0
817
818 CALL my_alloc(element%PON%IADS,8,numels)
819 CALL my_alloc(element%PON%IADS10,6,numels10)
820 CALL my_alloc(element%PON%IADS20,12,numels20)
821 CALL my_alloc(element%PON%IADS16,8,numels16)
822 CALL my_alloc(element%PON%IADQ,4,numelq)
823 CALL my_alloc(element%PON%IADC,4,numelc)
824 CALL my_alloc(element%PON%IAD_TRUSS,2,numelt)
825 CALL my_alloc(element%PON%IAD_BEAM,2,numelp)
826 CALL my_alloc(element%PON%IAD_SPRING,3,numelr)
827 CALL my_alloc(element%PON%IAD_TG,3,numeltg)
828 CALL my_alloc(element%PON%IAD_TG6,3,numeltg6)
829 CALL my_alloc(element%PON%IAD_MV,4,nskymv0)
830 CALL my_alloc(element%PON%IAD_CONLD,4,nconld)
831 CALL my_alloc(element%PON%IAD_CONV,4,glob_therm%NUMCONV)
832 CALL my_alloc(element%PON%IAD_RADIA,4,glob_therm%NUMRADIA)
833 CALL my_alloc(element%PON%IAD_LOADP,slloadp)
834 CALL my_alloc(element%PON%IAD_FXFLUX,4,glob_therm%NFXFLUX)
835
836
837
838
839
840
841 ALLOCATE (
iadwal(siadwal),stat=ierr)
842 IF (ierr/=0) GOTO 1000
844
845 ALLOCATE (
iadrbk(siadrbk),stat=ierr)
846 IF (ierr/=0) GOTO 1000
848
849 ALLOCATE (
iadi2(siadi2),stat=ierr)
850 IF (ierr/=0) GOTO 1000
852
853 ALLOCATE (
iadmv2(siadmv2),stat=ierr)
854 IF (ierr/=0) GOTO 1000
856
857 ALLOCATE (
iadmv3(siadmv3),stat=ierr)
858 IF (ierr/=0) GOTO 1000
860
861 ALLOCATE (
iadll(siadll),stat=ierr)
862 IF (ierr/=0) GOTO 1000
864
865 ALLOCATE (
iadrbm(siadrbm),stat=ierr)
866 IF (ierr/=0) GOTO 1000
868
869 ALLOCATE (
iadi18(siadi18),stat=ierr)
870 IF (ierr/=0) GOTO 1000
872 ALLOCATE (
iadrbmk(siadrbmk),stat=ierr)
873 IF (ierr/=0) GOTO 1000
875 ALLOCATE (
ilink(silink),stat=ierr)
876 IF (ierr/=0) GOTO 1000
878
879 ALLOCATE (
fr_rl(sfr_rl),stat=ierr)
880 IF (ierr/=0) GOTO 1000
882
883 ALLOCATE (
llink(sllink),stat=ierr)
884 IF (ierr/=0) GOTO 1000
886
887 ALLOCATE (
iadrl(siadrl),stat=ierr)
888 IF (ierr/=0) GOTO 1000
890
891 ALLOCATE (
linale(slinale+llinal),stat=ierr)
892 IF (ierr/=0) GOTO 1000
894
895 ALLOCATE (
neflsw(sneflsw),stat=ierr)
896 IF (ierr/=0) GOTO 1000
898
899 ALLOCATE (
nnflsw(snnflsw),stat=ierr)
900 IF (ierr/=0) GOTO 1000
902
903 ALLOCATE (
icut(sicut),stat=ierr)
904 IF (ierr/=0) GOTO 1000
906
907 ALLOCATE (
igeo(sigeo),stat=ierr)
908 IF (ierr/=0) GOTO 1000
910
911 ALLOCATE (
ipm(sipm),stat=ierr)
912 IF (ierr/=0) GOTO 1000
914
916 IF (ierr/=0) GOTO 1000
918
919 ALLOCATE (
icontact(sicontact),stat=ierr)
920 IF (ierr/=0) GOTO 1000
922
923 IF(nadmesh /= 0)THEN
924 ALLOCATE (
sh4tree(ksh4tree,numelc),stat=ierr)
925 ierror = ierror+ierr
926 ALLOCATE (
sh3tree(ksh3tree,numeltg),stat=ierr)
927 ierror = ierror+ierr
928 ALLOCATE (
ipadmesh(kipadmesh,npart),stat=ierr)
929 ierror = ierror+ierr
930 ALLOCATE (
sh4trim(abs(lsh4trim)),stat=ierr)
931 ierror = ierror+ierr
932 ALLOCATE (
sh3trim(abs(lsh3trim)),stat=ierr)
933 ierror = ierror+ierr
934 ELSE
940 END IF
941
942 ALLOCATE (
ithvar(sithvar),stat=ierr)
943 IF (ierr/=0) GOTO 1000
945
946 ALLOCATE (
icfield(sizfield*nloadc),stat=ierr)
947 IF (ierr/=0) GOTO 1000
949
950 ALLOCATE (
lcfield(slcfield),stat=ierr)
951 IF (ierr/=0) GOTO 1000
953
954 ALLOCATE (
iloadp(sizloadp*nloadp),stat=ierr)
955 IF (ierr/=0) GOTO 1000
957
958 ALLOCATE (
lloadp(slloadp),stat=ierr)
959 IF (ierr/=0) GOTO 1000
961
963 IF (ierr/=0) GOTO 1000
965
966 IF(nsphsol/=0)THEN
967
968 ALLOCATE(
sph2sol(numsph),stat=ierr)
969 IF (ierr/=0) GOTO 1000
971
972 ALLOCATE(
sol2sph(2*numels8),stat=ierr)
973 IF (ierr/=0) GOTO 1000
975
976 ALLOCATE(
irst(3*nsphsol),stat=ierr)
977 IF (ierr/=0) GOTO 1000
979
981 IF (ierr/=0) GOTO 1000
983
984 ELSE
985
987 IF (ierr/=0) GOTO 1000
989 IF (ierr/=0) GOTO 1000
990 ALLOCATE(
irst(0),stat=ierr)
991 IF (ierr/=0) GOTO 1000
993 IF (ierr/=0) GOTO 1000
994
995 END IF
996
997 ALLOCATE (element%PON%FSKYD(sfskyd),stat=ierr)
998 IF (ierr/=0) GOTO 1000
999 element%PON%FSKYD = 0
1000
1001 ALLOCATE (
ibcscyc(sibcscyc),stat=ierr)
1002 IF (ierr/=0) GOTO 1000
1004
1005 ALLOCATE (
lbcscyc(slbcscyc),stat=ierr)
1006 IF (ierr/=0) GOTO 1000
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033 ALLOCATE (thke(sthke),stat=ierr)
1034 IF (ierr/=0) GOTO 1000
1035 thke = 0
1036
1037 ALLOCATE (dampr(sdampr),stat=ierr)
1038 IF (ierr/=0) GOTO 1000
1039 dampr = 0
1040
1041 ALLOCATE (damp(sdamp),stat=ierr)
1042 IF (ierr/=0) GOTO 1000
1043 damp = 0
1044
1045
1046
1047
1048
1049
1050
1051
1052
1053 ALLOCATE (tf(stf),stat=ierr)
1054 IF (ierr/=0) GOTO 1000
1055 tf = 0
1056
1057 ALLOCATE (pm(spm),stat=ierr)
1058 IF (ierr/=0) GOTO 1000
1059 pm = 0
1060
1061 ALLOCATE (xframe(sxframe),stat=ierr)
1062 IF (ierr/=0) GOTO 1000
1063 xframe = 0
1064
1065 ALLOCATE (geo(sgeo),stat=ierr)
1066 IF (ierr/=0) GOTO 1000
1067 geo = 0
1068
1069 ALLOCATE (eani(seani),stat=ierr)
1070 IF (ierr/=0) GOTO 1000
1071 eani = 0
1072
1073 ALLOCATE (bufmat(sbufmat),stat=ierr)
1074 IF (ierr/=0) GOTO 1000
1075 bufmat = 0
1076
1077 ALLOCATE (
bufgeo(sbufgeo),stat=ierr)
1078 IF (ierr/=0) GOTO 1000
1080
1081 ALLOCATE (bufsf(sbufsf),stat=ierr)
1082 IF (ierr/=0) GOTO 1000
1083 bufsf = 0
1084
1085 ALLOCATE (rbmpc(srbmpc),stat=ierr)
1086 IF (ierr/=0) GOTO 1000
1087 rbmpc = 0
1088
1089 ALLOCATE (gjbufr(sgjbufr),stat=ierr)
1090 IF (ierr/=0) GOTO 1000
1091 gjbufr = 0
1092
1093 ALLOCATE (w(sw),stat=ierr)
1094 IF (ierr/=0) GOTO 1000
1095 w = 0
1096
1097 ALLOCATE (veul(sveul),stat=ierr)
1098 IF (ierr/=0) GOTO 1000
1099 veul = 0
1100
1101 ALLOCATE (fill(sfill),stat=ierr)
1102 IF (ierr/=0) GOTO 1000
1103 fill = 0
1104
1105 ALLOCATE (dfill(sdfill),stat=ierr)
1106 IF (ierr/=0) GOTO 1000
1107 dfill = 0
1108
1109 ALLOCATE (alph(salph),stat=ierr)
1110 IF (ierr/=0) GOTO 1000
1111 alph = 0
1112
1113 ALLOCATE (wb(swb),stat=ierr)
1114 IF (ierr/=0) GOTO 1000
1115 wb = 0
1116
1117 ALLOCATE (dsave(sdsave),stat=ierr)
1118 IF (ierr/=0) GOTO 1000
1119 dsave = 0
1120
1121 ALLOCATE (asave(sasave),stat=ierr)
1122 IF (ierr/=0) GOTO 1000
1123 asave = 0
1124
1125 ALLOCATE (spbuf(sspbuf),stat=ierr)
1126 IF (ierr/=0) GOTO 1000
1127 spbuf = 0
1128
1129 ALLOCATE (vsphio(svsphio),stat=ierr)
1130 IF (ierr/=0) GOTO 1000
1131 vsphio = 0
1132
1133 ALLOCATE (sphveln(ssphveln),stat=ierr)
1134 IF (ierr/=0) GOTO 1000
1135 sphveln = 0
1136
1137 ALLOCATE (msnf(smsnf),stat=ierr)
1138 IF (ierr/=0) GOTO 1000
1139 msnf = 0
1140
1141
1142 ALLOCATE (forc(sforc),stat=ierr)
1143 IF (ierr/=0) GOTO 1000
1144 forc = 0
1145
1146 ALLOCATE (vel(svel),stat=ierr)
1147 IF (ierr/=0) GOTO 1000
1148 vel = 0
1149
1150 ALLOCATE (fsav(sfsav),stat=ierr)
1151 IF (ierr/=0) GOTO 1000
1152 fsav = 0
1153
1154 ALLOCATE (
tab_ump(7,taille),stat=ierr)
1155 IF (ierr/=0) GOTO 1000
1157
1158 ALLOCATE (
poin_ump(nummat),stat=ierr)
1159 IF (ierr/=0) GOTO 1000
1161
1162 IF(ispmd==0) THEN
1163 ALLOCATE (cputime_mp_glob(nbr_gpmp,nspmd+1),stat=ierr)
1164 IF (ierr/=0) GOTO 1000
1165 cputime_mp_glob = 0
1166 ENDIF
1167
1168 ALLOCATE (cputime_mp(taille),stat=ierr)
1169 IF (ierr/=0) GOTO 1000
1170 cputime_mp = 0
1171
1172 ALLOCATE (tab_mat(stab_mat),stat=ierr)
1173 IF (ierr/=0) GOTO 1000
1174 tab_mat = zero
1175
1176 ALLOCATE (fzero(sfzero),stat=ierr)
1177 IF (ierr/=0) GOTO 1000
1178 fzero = 0
1179
1180 ALLOCATE (xlas(sxlas),stat=ierr)
1181 IF (ierr/=0) GOTO 1000
1182 xlas = 0
1183
1184 ALLOCATE (accelm(saccelm),stat=ierr)
1185 IF (ierr/=0) GOTO 1000
1186 accelm = 0
1187
1188 ALLOCATE (gauge(llgauge*nbgauge),stat=ierr)
1189 IF (ierr/=0) GOTO 1000
1190 gauge = 0
1191
1192 ALLOCATE (fbvel(sfbvel),stat=ierr)
1193 IF (ierr/=0) GOTO 1000
1194 fbvel = 0
1195
1196 ALLOCATE (factiv(lractiv*nactiv),stat=ierr)
1197 IF (ierr/=0) GOTO 1000
1198 factiv = 0
1199
1200 ALLOCATE (grav(sgrav),stat=ierr)
1201 IF (ierr/=0) GOTO 1000
1202 grav = 0
1203
1204 ALLOCATE (fr_wave(sfr_wave),stat=ierr)
1205 IF (ierr/=0) GOTO 1000
1206 fr_wave = 0
1207
1208 ALLOCATE (parts0(sparts0),stat=ierr)
1209 IF (ierr/=0) GOTO 1000
1210 parts0 = 0
1211
1212 ALLOCATE (elbuf(selbuf),stat=ierr)
1213 IF (ierr/=0) GOTO 1000
1214 elbuf = 0
1215
1216
1217 ALLOCATE (rby(srby),stat=ierr)
1218 IF (ierr/=0) GOTO 1000
1219 rby = 0
1220
1221 ALLOCATE (rbyl(srbyl),stat=ierr)
1222 IF (ierr/=0) GOTO 1000
1223 rbyl = 0
1224
1225 ALLOCATE (rivet(srivet),stat=ierr)
1226 IF (ierr/=0) GOTO 1000
1227 rivet = 0
1228
1229 ALLOCATE (secbuf(ssecbuf),stat=ierr)
1230 IF (ierr/=0) GOTO 1000
1231 secbuf = 0
1232
1233 ALLOCATE (volmon(svolmon),stat=ierr)
1234 IF (ierr/=0) GOTO 1000
1235 volmon = 0
1236
1237 ALLOCATE (lambda(slambda),stat=ierr)
1238 IF (ierr/=0) GOTO 1000
1239 lambda = 0
1240
1241 ALLOCATE (rconx(srconx),stat=ierr)
1242 IF (ierr/=0) GOTO 1000
1243 rconx = 0
1244
1245 ALLOCATE (rcontact(srcontact),stat=ierr)
1246 IF (ierr/=0) GOTO 1000
1247 rcontact = ep30
1248
1249 ALLOCATE (acontact(srcontact),stat=ierr)
1250 IF (ierr/=0) GOTO 1000
1251 acontact = ep30
1252
1253 ALLOCATE (pcontact(srcontact),stat=ierr)
1254 IF (ierr/=0) GOTO 1000
1255 pcontact = zero
1256
1257 ALLOCATE (wa(swa),stat=ierr)
1258 IF (ierr/=0) GOTO 1000
1259 wa = 0
1260
1261 ALLOCATE (fv(sfv),stat=ierr)
1262 IF (ierr/=0) GOTO 1000
1263 fv = 0
1264
1265
1266
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1283
1284
1285 ALLOCATE (partsav(spartsav),stat=ierr)
1286 IF (ierr/=0) GOTO 1000
1287 partsav = 0
1288
1289 ALLOCATE (element%PON%FSKY(8,sfsky/8),stat=ierr)
1290 IF (ierr/=0) GOTO 1000
1291 element%PON%FSKY = 0
1292 element%PON%SFSKY = sfsky
1293 element%PON%MAX_SFSKY = sfsky
1294
1295
1296 ALLOCATE (element%PON%FSKYM(sfskym),stat=ierr)
1297 IF (ierr/=0) GOTO 1000
1298 element%PON%FSKYM = 0
1299
1300 ALLOCATE (interfaces%PON%FSKYI(lskyi,nfskyi),stat=ierr)
1301 IF (ierr/=0) GOTO 1000
1302 interfaces%PON%FSKYI = 0
1303
1304 ALLOCATE (uwa(suwa),stat=ierr)
1305 IF (ierr/=0) GOTO 1000
1306 uwa = 0
1307
1308 ALLOCATE (val2(sval2),stat=ierr)
1309 IF (ierr/=0) GOTO 1000
1310 val2 = 0
1311
1312 ALLOCATE (phi(sphi),stat=ierr)
1313 IF (ierr/=0) GOTO 1000
1314 phi = 0
1315
1316 ALLOCATE (phie(sphi),stat=ierr)
1317 IF (ierr/=0) GOTO 1000
1318 phie = zero
1319 ALLOCATE (por(sphi),stat=ierr)
1320 IF (ierr/=0) GOTO 1000
1321 por = zero
1322
1323 ALLOCATE (r(sr),stat=ierr)
1324 IF (ierr/=0) GOTO 1000
1325 r = 0
1326
1327 ALLOCATE (estif(sestif),stat=ierr)
1328 IF (ierr/=0) GOTO 1000
1329 estif = 0
1330
1331 ALLOCATE (crflsw(scrflsw),stat=ierr)
1332 IF (ierr/=0) GOTO 1000
1333 crflsw = 0
1334
1335 ALLOCATE (flsw(sflsw),stat=ierr)
1336 IF (ierr/=0) GOTO 1000
1337 flsw = 0
1338
1339
1340
1341
1342
1343 isiz1=2
1344 isiz2=numelq+numeltg
1345 IF(n2d == 0)THEN
1346 isiz1=3
1347 isiz2=numels
1348 ENDIF
1349
1350 IF(output%DATA%FANI_CELL%IS_VORT_X_REQUESTED)THEN
1351 ALLOCATE (output%DATA%FANI_CELL%VORT_X(isiz2),stat=ierr)
1352 IF (ierr/=0) GOTO 1000
1353 output%DATA%FANI_CELL%VORT_X(:) = zero
1354 ENDIF
1355 IF(output%DATA%FANI_CELL%IS_VORT_Y_REQUESTED)THEN
1356 ALLOCATE (output%DATA%FANI_CELL%VORT_Y(isiz2),stat=ierr)
1357 IF (ierr/=0) GOTO 1000
1358 output%DATA%FANI_CELL%VORT_Y(:) = zero
1359 ENDIF
1360 IF(output%DATA%FANI_CELL%IS_VORT_Z_REQUESTED)THEN
1361 ALLOCATE (output%DATA%FANI_CELL%VORT_Z(isiz2),stat=ierr)
1362 IF (ierr/=0) GOTO 1000
1363 output%DATA%FANI_CELL%VORT_Z(:) = zero
1364 ENDIF
1365
1366 IF(output%DATA%FANI_CELL%IS_F18_FVM_REQUESTED)THEN
1367 ALLOCATE (output%DATA%FANI_CELL%F18(isiz1,isiz2),stat=ierr)
1368 IF (ierr/=0) GOTO 1000
1369 output%DATA%FANI_CELL%F18(:,:) = zero
1370 ENDIF
1371
1372 ALLOCATE (xcut(sxcut),stat=ierr)
1373 IF (ierr/=0) GOTO 1000
1374 xcut = 0
1375
1376
1377
1378
1379
1380 ALLOCATE (tani(stani),stat=ierr)
1381 IF (ierr/=0) GOTO 1000
1382 tani = 0
1383
1384 ALLOCATE (secfcum(ssecfcum),stat=ierr)
1385 IF (ierr/=0) GOTO 1000
1386 secfcum = 0
1387
1388 ALLOCATE (wasph(swasph),stat=ierr)
1389 IF (ierr/=0) GOTO 1000
1390 wasph = 0
1391
1392 ALLOCATE (w16(sw16),stat=ierr)
1393 IF (ierr/=0) GOTO 1000
1394 w16 = 0
1395
1396 ALLOCATE (dretrio(sdretrio),stat=ierr)
1397 IF (ierr/=0) GOTO 1000
1398 dretrio = 0
1399
1400 ALLOCATE (lbvrs(slbvrs),stat=ierr)
1401 IF (ierr/=0) GOTO 1000
1402 lbvrs = 0
1403
1404 ALLOCATE (pv(spv),stat=ierr)
1405 IF (ierr/=0) GOTO 1000
1406 pv = 0
1407
1408 IF(nadmesh /= 0 .OR. irest_mselt/=0)THEN
1409
1410 ALLOCATE(msc(numelc),stat=ierr)
1411 ierror = ierror+ierr
1412
1413 ALLOCATE(inc(numelc),stat=ierr)
1414 ierror = ierror+ierr
1415
1416 ALLOCATE(mstg(numeltg),stat=ierr)
1417 ierror = ierror+ierr
1418
1419 ALLOCATE(intg(numeltg),stat=ierr)
1420 ierror = ierror+ierr
1421
1422 ALLOCATE(ptg(3,numeltg),stat=ierr)
1423 ierror = ierror+ierr
1424
1425 ELSE
1426
1427 ALLOCATE(msc(0))
1428 ALLOCATE(inc(0))
1429 ALLOCATE(mstg(0))
1430 ALLOCATE(intg(0))
1431 ALLOCATE(ptg(0,0))
1432
1433 END IF
1434
1435 IF(irest_mselt/=0)THEN
1436
1437 ALLOCATE(mssa(numels),stat=ierr)
1438 ierror = ierror+ierr
1439
1440 ALLOCATE(mstr(numelt),stat=ierr)
1441 ierror = ierror+ierr
1442
1443 ALLOCATE(msp(numelp),stat=ierr)
1444 ierror = ierror+ierr
1445
1446 ALLOCATE(msrt(numelr),stat=ierr)
1447 ierror = ierror+ierr
1448
1449 ELSE
1450
1451 ALLOCATE(mssa(0))
1452 ALLOCATE(mstr(0))
1453 ALLOCATE(msp(0))
1454 ALLOCATE(msrt(0))
1455
1456 END IF
1457
1458 IF(nadmesh /= 0)THEN
1459
1460 ALLOCATE(padmesh(kpadmesh,npart),stat=ierr)
1461 ierror = ierror+ierr
1462
1463 IF(glob_therm%ITHERM_FE > 0)THEN
1464
1465 ALLOCATE(mcpc(numelc),stat=ierr)
1466 ierror = ierror+ierr
1467
1468 ALLOCATE(mcptg(numeltg),stat=ierr)
1469 ierror = ierror+ierr
1470
1471 ELSE
1472
1473 ALLOCATE(mcpc(0),mcptg(0))
1474
1475 END IF
1476
1477 ELSE
1478 ALLOCATE(padmesh(0,0))
1479 ALLOCATE(mcpc(0),mcptg(0))
1480 END IF
1481
1482 IF(istatcnd/=0)THEN
1483 ALLOCATE(mscnd(numnod),incnd(numnod),stat=ierr)
1484 ierror = ierror+ierr
1485 ELSE
1486 ALLOCATE(mscnd(0),incnd(0))
1487 END IF
1488
1489
1490
1491 IF(irigid_mat > 0) THEN
1492 ALLOCATE(rbym(nfrbym*nrbym),stat=ierr)
1493 IF (ierr/=0) GOTO 1000
1494 rbym = zero
1495 ALLOCATE(
irbym(nirbym*nrbym),stat=ierr)
1496 IF (ierr/=0) GOTO 1000
1498 ALLOCATE(
lnrbym(ngslnrbym),stat=ierr)
1499 IF (ierr/=0) GOTO 1000
1502 IF (ierr/=0) GOTO 1000
1504 ALLOCATE(
icodrbym(nrbym*2),stat=ierr)
1505 IF (ierr/=0) GOTO 1000
1507 ELSE
1510 ENDIF
1511
1512
1513 IF(glob_therm%NUMCONV > 0) THEN
1514 ALLOCATE(
ibcv(glob_therm%NICONV*glob_therm%NUMCONV),stat=ierr)
1515 IF (ierr/=0) GOTO 1000
1517 ALLOCATE(fconv(glob_therm%LFACTHER*glob_therm%NUMCONV),stat=ierr)
1518 IF (ierr/=0) GOTO 1000
1519 fconv = zero
1520 ELSE
1521 ALLOCATE(
ibcv(0),fconv(0))
1522 ENDIF
1523
1524 IF(glob_therm%NUMRADIA > 0) THEN
1525 ALLOCATE(
ibcr(glob_therm%NIRADIA*glob_therm%NUMRADIA),stat=ierr)
1526 IF (ierr/=0) GOTO 1000
1528 ALLOCATE(fradia(glob_therm%LFACTHER*glob_therm%NUMRADIA),stat=ierr)
1529 IF (ierr/=0) GOTO 1000
1530 fradia = zero
1531 ELSE
1532 ALLOCATE(
ibcr(0),fradia(0))
1533 ENDIF
1534
1535 IF(glob_therm%NFXFLUX > 0) THEN
1536 ALLOCATE(
ibfflux(glob_therm%NITFLUX*glob_therm%NFXFLUX),stat=ierr)
1537 IF (ierr/=0) GOTO 1000
1539 ALLOCATE(fbfflux(glob_therm%LFACTHER*glob_therm%NFXFLUX),stat=ierr)
1540 IF (ierr/=0) GOTO 1000
1541 fbfflux = zero
1542 ELSE
1543 ALLOCATE(
ibfflux(0),fbfflux(0))
1544 ENDIF
1545
1546 IF(glob_therm%NFXTEMP > 0) THEN
1547 ALLOCATE(
ibftemp(glob_therm%NIFT*glob_therm%NFXTEMP),stat=ierr)
1548 IF (ierr/=0) GOTO 1000
1550 ALLOCATE(fbftemp(glob_therm%LFACTHER*glob_therm%NFXTEMP),stat=ierr)
1551 IF(ierr/=0) GOTO 1000
1552 fbftemp = zero
1553 ELSE
1554 ALLOCATE(
ibftemp(0),fbftemp(0))
1555 ENDIF
1556
1557 ALLOCATE (msf(sms),stat=ierr)
1558 IF (ierr/=0) GOTO 1000
1559 IF(ialelag > 0) THEN
1560 msf = zero
1561 ENDIF
1562
1563
1564 IF(iplyxfem > 0 ) THEN
1565 ALLOCATE(ms_ply(nplyxfe*nplymax),stat=ierr)
1566 IF (ierr/=0) GOTO 1000
1567 ms_ply = zero
1568 ALLOCATE(zi_ply(nplyxfe*nplymax),stat=ierr)
1569 IF (ierr/=0) GOTO 1000
1570 zi_ply = zero
1571 ALLOCATE(msz2(nplyxfe),stat=ierr)
1572 IF (ierr/=0) GOTO 1000
1573 msz2 = zero
1575 IF (ierr/=0) GOTO 1000
1578 IF (ierr/=0) GOTO 1000
1581 IF (ierr/=0) GOTO 1000
1585
1587 IF (ierr/=0) GOTO 1000
1590 IF (ierr/=0) GOTO 1000
1592
1594 IF (ierr/=0) GOTO 1000
1596 ELSE
1598 ALLOCATE(msz2(0))
1601 ENDIF
1602 IF(iplybcs > 0) THEN
1604 IF (ierr/=0) GOTO 1000
1607 IF (ierr/=0) GOTO 1000
1610 IF (ierr/=0) GOTO 1000
1612 ALLOCATE(
ibc_ply(numnod),stat=ierr)
1613 IF (ierr/=0) GOTO 1000
1615 ELSE
1617 ENDIF
1618 IF(intplyxfem > 0 ) THEN
1619 ALLOCATE(wagap(2*numnod),stat=ierr)
1620 IF (ierr/=0) GOTO 1000
1621 wagap = zero
1622 ELSE
1623 ALLOCATE(wagap(0),stat=ierr)
1624 ENDIF
1625
1626
1627
1628 IF(idtmins_old==1) THEN
1629 ALLOCATE(admsms(numnod),stat=ierr)
1630 IF (ierr/=0) GOTO 1000
1631 END IF
1632 IF(idtmins_old==2) THEN
1633 ALLOCATE(dmelc(numelc) ,stat=ierr)
1634 IF (ierr/=0) GOTO 1000
1635 ALLOCATE(dmeltg(numeltg),stat=ierr)
1636 IF (ierr/=0) GOTO 1000
1637 ALLOCATE(dmels(numels) ,stat=ierr)
1638 IF (ierr/=0) GOTO 1000
1639 ALLOCATE(dmeltr(numelt) ,stat=ierr)
1640 IF (ierr/=0) GOTO 1000
1641 ALLOCATE(dmelp(numelp) ,stat=ierr)
1642 IF (ierr/=0) GOTO 1000
1643 ALLOCATE(dmelrt(numelr) ,stat=ierr)
1644 IF (ierr/=0) GOTO 1000
1645 END IF
1646 IF(idtmins_old/=0.OR.idtmins_int_old/=0) THEN
1647 ALLOCATE(res_sms(3*numnod),stat=ierr)
1648 IF (ierr/=0) GOTO 1000
1649 END IF
1650 IF(idtmins_old==2.OR.idtmins_int_old/=0) THEN
1651 ALLOCATE(diag_sms(numnod),stat=ierr)
1652 IF (ierr/=0) GOTO 1000
1653 ALLOCATE(dmint2(4*i2nsn25),stat=ierr)
1654 IF (ierr/=0) GOTO 1000
1655 END IF
1656 IF (isms_selec /= 0) THEN
1658 IF (ierr/=0) GOTO 1000
1660 ELSE
1662 ENDIF
1663
1664
1665 IF(nintstamp/=0.AND.anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0)THEN
1666 ALLOCATE(fcontg(3,numnodg),stat=ierr)
1667 IF (ierr/=0) GOTO 1000
1668 fcontg=zero
1669 ELSE
1670 ALLOCATE(fcontg(0,0),stat=ierr)
1671 END IF
1672
1674 IF(anim_v(26)+h3d_data%N_VECT_CONT_MAX >0.OR.
ifcontmax > 0)
THEN
1675 IF(nspmd > 1 .AND. nintstamp/=0) THEN
1676 IF(ispmd == 0) THEN
1677 ALLOCATE(fcont_max(3*numnodg),stat=ierr)
1678 IF (ierr/=0) GOTO 1000
1679 fcont_max=zero
1681 ELSE
1682 ALLOCATE(fcont_max(0),stat=ierr)
1683 ENDIF
1684 ELSE
1685 ALLOCATE(fcont_max(3*numnod),stat=ierr)
1686 IF (ierr/=0) GOTO 1000
1687 fcont_max=zero
1689 ENDIF
1690 ELSE
1691 ALLOCATE(fcont_max(0),stat=ierr)
1692 END IF
1693
1694
1695 IF(nintstamp/=0.AND.anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT>0)THEN
1696 ALLOCATE(fncontg(3,numnodg),ftcontg(3,numnodg),stat=ierr)
1697 IF (ierr/=0) GOTO 1000
1698 fncontg=zero
1699 ftcontg=zero
1700 ELSE
1701 ALLOCATE(fncontg(0,0),ftcontg(0,0),stat=ierr)
1702 END IF
1703
1705 IF(h3d_data%N_VECT_PCONT_MAX >0.OR.
ifcontpmax > 0)
THEN
1706 IF(nspmd > 1 .AND. nintstamp/=0) THEN
1707 IF(ispmd == 0) THEN
1708 ALLOCATE(ftcont_max(3*numnodg),fncont_max(3*numnodg),stat=ierr)
1709 IF (ierr/=0) GOTO 1000
1710 ftcont_max=zero
1711 fncont_max=zero
1713 ELSE
1714 ALLOCATE(fncont_max(0),ftcont_max(0),stat=ierr)
1715 ENDIF
1716 ELSE
1717 ALLOCATE(ftcont_max(3*numnod),fncont_max(3*numnod),stat=ierr)
1718 IF (ierr/=0) GOTO 1000
1719 ftcont_max=zero
1720 fncont_max=zero
1722 ENDIF
1723 ELSE
1724 ALLOCATE(fncont_max(0),stat=ierr)
1725 ALLOCATE(ftcont_max(0),stat=ierr)
1726 END IF
1727
1728 IF(h3d_data%N_VECT_CONT2_MAX >0.OR.
ifcont2max > 0)
THEN
1729 ALLOCATE(fcont2_max(3*numnod),stat=ierr)
1730 IF (ierr/=0) GOTO 1000
1731 fcont2_max=zero
1732 ELSE
1733 ALLOCATE(fcont2_max(0),stat=ierr)
1734 END IF
1735
1736 IF(h3d_data%N_VECT_PCONT2_MAX >0.OR.
ifcontp2max > 0)
THEN
1737 ALLOCATE(ftcont2_max(3*numnod),fncont2_max(3*numnod),stat=ierr)
1738 ALLOCATE(npcont2_max(3*numnod),stat
1739 IF (ierr/=0) GOTO 1000
1740 fncont2_max=zero
1741 ftcont2_max=zero
1742 npcont2_max=zero
1743 ELSE
1744 ALLOCATE(fncont2_max(0),stat=ierr)
1745 ALLOCATE(ftcont2_max(0),stat=ierr)
1746 ALLOCATE(npcont2_max(0),stat=ierr)
1747 END IF
1748
1749
1750 IF(h3d_data%N_VECT_CONT2_MIN >0.OR.
ifcont2min > 0)
THEN
1751 ALLOCATE(fcont2_min(3*numnod),stat=ierr)
1752 IF (ierr/=0) GOTO 1000
1753 fcont2_min=zero
1754 ELSE
1755 ALLOCATE(fcont2_min(0),stat=ierr)
1756 END IF
1757
1758 IF(h3d_data%N_VECT_PCONT2_MIN >0.OR.
ifcontp2min > 0)
THEN
1759 ALLOCATE(ftcont2_min(3*numnod),fncont2_min(3*numnod),stat=ierr)
1760 ALLOCATE(npcont2_min(3*numnod),stat=ierr)
1761 IF (ierr/=0) GOTO 1000
1762 fncont2_min=zero
1763 ftcont2_min=zero
1764 npcont2_min=zero
1765 ELSE
1766 ALLOCATE(fncont2_min(0),stat=ierr)
1767 ALLOCATE(ftcont2_min(0),stat=ierr)
1768 ALLOCATE(npcont2_min(0),stat=ierr)
1769 END IF
1770
1771 IF(h3d_data%N_VECT_CONT2M > 0)THEN
1772 ALLOCATE(mcont2(3,numnod),stat=ierr)
1773 IF (ierr/=0) GOTO 1000
1774 mcont2=zero
1775 mcont2=zero
1776 ELSE
1777 ALLOCATE(mcont2(3,0),stat=ierr)
1778 ALLOCATE(mcont2(3,0),stat=ierr)
1779 END IF
1780
1781
1782
1783 IF(output%DATA%S_EFRICINT > 0)THEN
1784 IF(output%DATA%NINEFRIC > 0) THEN
1785 ALLOCATE(output%DATA%EFRIC(output%DATA%NINEFRIC,numnod))
1786 ELSE
1787 ALLOCATE(output%DATA%EFRIC(0,0))
1788 ENDIF
1789 IF(output%DATA%NINEFRIC_STAMP > 0) THEN
1790 ALLOCATE(output%DATA%EFRIC_STAMP(output%DATA%NINEFRIC_STAMP,numnodg
1791 ELSE
1792 ALLOCATE(output%DATA%EFRIC_STAMP(0,0))
1793 ENDIF
1794 END IF
1795
1796 IF(output%DATA%S_EFRIC > 0)THEN
1797 ALLOCATE(output%DATA%EFRICG(numnod))
1798 IF(nintstamp/=0) THEN
1799 ALLOCATE(output%DATA%EFRICG_STAMP(numnodg))
1800 ELSE
1801 ALLOCATE(output%DATA%EFRICG_STAMP(0))
1802 ENDIF
1803
1804
1805
1806 ENDIF
1807
1808 IF(nintstamp/=0.AND.(anim_v(12)+outp_v(12)+h3d_data%N_VECT_PCONT>0.OR.
1809 . anim_v(4)+outp_v(4)+h3d_data%N_VECT_CONT >0))THEN
1811 IF (ierr/=0) GOTO 1000
1814 ELSE
1816 ENDIF
1817
1818 IF(nintskidold >0)THEN
1819 IF(nintstamp/=0) THEN
1820 ALLOCATE(pskids(nintskidold,numnodg),stat=ierr)
1821 IF (ierr/=0) GOTO 1000
1822 ELSE
1823 ALLOCATE(pskids(nintskidold,numnod),stat=ierr)
1824 IF (ierr/=0) GOTO 1000
1825 ENDIF
1826 END IF
1827
1828
1829
1830 IF(nintstamp/=0)THEN
1831 ALLOCATE(
intstamp(nintstamp),stat=ierr)
1832 IF (ierr/=0)THEN
1833 CALL ancmsg(msgid=20,anmode=aninfo)
1835 ENDIF
1837 ELSE
1839 END IF
1840
1841 ALLOCATE (cfield(lfacload*nloadc),stat=ierr)
1842 IF (ierr/=0) GOTO 1000
1843 cfield = 0
1844
1845 ALLOCATE (loadp(lfacload*nloadp),stat=ierr)
1846 IF (ierr/=0) GOTO 1000
1847 loadp = 0
1848
1849 IF (nintloadp>0) THEN
1850
1852 IF (ierr/=0) GOTO 1000
1854
1855 ALLOCATE (
loadpinter(ninter*nloadp_hyd),stat=ierr)
1856 IF (ierr/=0) GOTO 1000
1858
1859 ALLOCATE (dgaploadint(ninter*nloadp_hyd),stat=ierr)
1860 IF (ierr/=0) GOTO 1000
1861 dgaploadint = zero
1862 ELSE
1865 ALLOCATE (dgaploadint(0),stat=ierr)
1866 ENDIF
1867
1868
1869 IF (ncluster > 0) THEN
1870 ALLOCATE(
cluster(ncluster),stat=ierr)
1871 IF (ierr /= 0) THEN
1872 CALL ancmsg(msgid=20,anmode=aninfo)
1874 ENDIF
1875 ELSE
1877 ENDIF
1878
1879
1880 IF(icrack3d > 0) THEN
1882 IF (ierr/=0) GOTO 1000
1885 IF (ierr/=0) GOTO 1000
1887 ALLOCATE(
iadc_crkxfem(4*ecrkxfec+3*ecrkxfetg),stat=ierr)
1888 IF (ierr/=0) GOTO 1000
1893 IF (ierr/=0) GOTO 1000
1896 IF (ierr/=0) GOTO 1000
1899 IF (ierr/=0) GOTO 1000
1902 IF (ierr/=0) GOTO 1000
1904 ALLOCATE (
crknodiad(sprocne_crkxfe),stat=ierr)
1905 IF (ierr/=0) GOTO 1000
1907 ELSE
1913 ENDIF
1914
1915 IF(icrack3d > 0)THEN
1916 ALLOCATE(
ibordnode(sibordnode),stat=ierr)
1917 IF (ierr/=0) GOTO 1000
1919
1920 ALLOCATE(
iedgesh(siedgesh),stat=ierr)
1921 IF (ierr/=0) GOTO 1000
1923
1924 ALLOCATE (
nodedge(2*numedges),stat=ierr)
1925 IF (ierr/=0) GOTO 1000
1927
1928 ALLOCATE (
iedge(numedges),stat=ierr)
1929 IF (ierr/=0) GOTO 1000
1931
1932 ALLOCATE (
iad_edge(siad_edge),stat=ierr)
1933 IF (ierr/=0) GOTO 1000
1935
1936 ALLOCATE (
fr_edge(sfr_edge),stat=ierr)
1937 IF (ierr/=0) GOTO 1000
1939
1940 ALLOCATE (
fr_nbedge(sfr_nbedge),stat=ierr)
1941 IF (ierr/=0) GOTO 1000
1943
1944 ALLOCATE (
nodlevxf(ncrkxfe),stat=ierr)
1945 IF (ierr/=0) GOTO 1000
1947
1949 IF (ierr/=0) GOTO 1000
1951 ELSE
1961 END IF
1962
1963 ALLOCATE (
kxig3d(skxig3d),stat=ierr)
1964 IF (ierr/=0) GOTO 1000
1966
1967 ALLOCATE (
ixig3d(sixig3d),stat=ierr)
1968 IF (ierr/=0) GOTO 1000
1970
1971 ALLOCATE (knot(sknot),stat=ierr)
1972 IF (ierr/=0) GOTO 1000
1973 knot = zero
1974
1975
1976 ALLOCATE( wige(1), stat=ierr)
1977 IF (ierr/=0) GOTO 1000
1978 wige = zero
1979
1980 ALLOCATE (knotlocpc(sknotlocpc),stat=ierr)
1981 IF (ierr/=0) GOTO 1000
1982 knotlocpc = zero
1983
1984 ALLOCATE (knotlocel(sknotlocel),stat=ierr)
1985 IF (ierr/=0) GOTO 1000
1986 knotlocel = zero
1987
1988 IF(icrack3d > 0)THEN
1989 ALLOCATE (
elcutc(2*(numelc+numeltg)),stat=ierr)
1990 IF (ierr/=0) GOTO 1000
1992
1993 ALLOCATE (
nodenr(ncrkxfe),stat=ierr)
1994 IF (ierr/=0) GOTO 1000
1996
1998 IF (ierr/=0) GOTO 1000
2000
2001 ALLOCATE (
enrtag(numnod*ienrnod*nlevmax),stat=ierr)
2002 IF (ierr/=0) GOTO 1000
2004
2005 ELSE
2010 END IF
2011
2013
2015 IF (ierr /= 0) GOTO 1000
2017
2019 IF (ierr /= 0) GOTO 1000
2021
2023 IF (ierr /= 0) GOTO 1000
2025
2027 IF (ierr /= 0) GOTO 1000
2029
2031 IF (ierr /= 0) GOTO 1000
2033
2035 IF (ierr /= 0) GOTO 1000
2037
2038
2039
2040
2041
2042
2043
2044
2046
2047 ENDIF
2048
2049 IF (
alemuscl_param%IALEMUSCL > 0 .AND. .NOT. multi_fvm%IS_USED)
THEN
2051 nbmat = trimat
2053 nbmat = 1
2054 ELSE
2055 nbmat = 0
2056 ENDIF
2057 IF(nbmat>0)THEN
2058 ALLOCATE(
alemuscl_buffer%VOLUME_FRACTION(numels + nsvois + numelq + nqvois,nbmat), stat = ierr)
2059 IF (ierr /= 0) GOTO 1000
2062 IF (ierr /= 0) GOTO 1000
2065 IF (ierr /= 0) GOTO 1000
2067 ALLOCATE(
alemuscl_buffer%GRAD(numels + nsvois + numelq + nqvois,3,nbmat), stat = ierr)
2068 IF (ierr /= 0) GOTO
2070 ALLOCATE(
alemuscl_buffer%ELCENTER(numels + nsvois + numelq + nqvois,3), stat = ierr)
2071 IF (ierr /= 0) GOTO 1000
2073 ENDIF
2074 ENDIF
2075
2076
2077 IF(multi_fvm%IS_USED)THEN
2078 nbmat = multi_fvm%NBMAT
2080
2082 nbmat = 1
2083 ELSE
2084 nbmat = 0
2085 ENDIF
2086 CALL segvar%create(nbmat)
2087
2088
2089 IF (multi_fvm%IS_USED) THEN
2092 IF (ierr /= 0) GOTO 1000
2094 ALLOCATE(multi_fvm%FVM_CONNECTIVITY%KVOIS(
ale%GLOBAL%SIELVS), stat
2095 IF (ierr /= 0) GOTO 1000
2096 multi_fvm%FVM_CONNECTIVITY%KVOIS = 0
2097 ENDIF
2098
2099
2100
2101 IF (numply > 0) THEN
2102 ALLOCATE (
ply_info(3,numply),stat=ierr)
2103 IF (ierr /= 0) GOTO 1000
2105 ELSE
2107 END IF
2108
2109 IF (ns10e > 0) THEN
2110 ALLOCATE (
itagnd(numnod),
icnds10(3*ns10e),vnd(3*ns10e),stifnd(ns10e),stat=ierr)
2111 IF (ierr /= 0) GOTO 1000
2114 vnd(1:3*ns10e) = zero
2115 stifnd(1:ns10e) = zero
2116 IF (idamp/=0.OR.ndamp>0.OR.istat/=0) THEN
2117 ALLOCATE (vmd(3*numnod),stat=ierr)
2118 IF (ierr /= 0) GOTO 1000
2119 vmd = zero
2120 END IF
2121 IF (iparit/=0.AND.iparit/=3) THEN
2123 IF (ierr /= 0) GOTO 1000
2124 fskycnd(1:4*lcncnd)=zero
2125 ELSE
2127 END IF
2128 ELSE
2130 END IF
2131
2133 IF (ierr/=0) GOTO 1000
2136 IF (ierr/=0) GOTO 1000
2138
2139 ALLOCATE (
addcncnd(saddcncnd),stat=ierr)
2140 IF (ierr/=0) GOTO 1000
2142 ALLOCATE (
procncnd(sprocncnd),stat=ierr)
2143 IF (ierr/=0) GOTO 1000
2145 ALLOCATE (
iadcnd(siadcnd),stat=ierr)
2146 IF (ierr/=0) GOTO 1000
2148
2149 IF(n2d/=0) THEN
2150 ALLOCATE(ms_2d(numnod),stat=ierr)
2151 IF (ierr/=0) GOTO 1000
2152 ELSE
2153 ALLOCATE(ms_2d(0))
2154 ENDIF
2155
2156
2157 ALLOCATE(pinch_data%XPINCH(
2158 pinch_data%XPINCH = zero
2159 ALLOCATE(pinch_data%DPINCH(3,
max(1,
npinch)))
2160 pinch_data%DPINCH = zero
2161 ALLOCATE(pinch_data%VPINCH(3,
max(1,
npinch)))
2162 pinch_data%VPINCH = zero
2163 ALLOCATE(pinch_data%APINCH(3,
max(1,pinch_data%SAPINCH)))
2164 pinch_data%APINCH = zero
2165 ALLOCATE(pinch_data%STIFPINCH(
max(1,pinch_data%SAPINCH)))
2166 pinch_data%STIFPINCH = zero
2167 ALLOCATE(pinch_data%MSPINCH(
max(1,pinch_data%SAPINCH)))
2168 pinch_data%MSPINCH = 0
2169
2170
2171
2173 sxcell = 0
2174
2175 IF ((iale > 0 .AND.
ale%GRID%NWALE_RST == 6) )
THEN
2176 sxcell = numels + numelq + ne_nsvois
2177 ALLOCATE(xcell(3, sxcell),stat=ierr)
2178 IF (ierr/=0) GOTO 1000
2179 ALLOCATE(xface(3,
2180 IF (ierr/=0) GOTO 1000
2181 ENDIF
2182
2183 IF (inter18_is_variable_gap_defined) THEN
2184 sxcell = numels + numelq + ne_nsvois + ne_nqvois
2185 IF(.NOT.ALLOCATED(xcell))ALLOCATE(xcell(3, sxcell),stat=ierr)
2186 IF (ierr/=0) GOTO 1000
2187 ENDIF
2188 IF(.NOT.ALLOCATED(xcellALLOCATE
2189 IF(.NOT.ALLOCATED(xface)) ALLOCATE(xface(3,6,0))
2190 xcell(1:3,1:sxcell) = zero
2191
2193 IF (ierr/=0) GOTO 1000
2194
2195 ALLOCATE (
fxbnod(lennod),stat
2196 IF (ierr/=0) GOTO 1000
2198 ALLOCATE (fxbmod(6,lenmod),stat=ierr)
2199 IF (ierr/=0) GOTO 1000
2200
2201 ALLOCATE (fxbglm(lenglm),stat=ierr)
2202 IF (ierrGOTO 1000
2203 fxbglm = 0
2204 ALLOCATE
2205 IF (ierr/=0) GOTO 1000
2206 fxbcpm = 0
2207 ALLOCATE (fxbcps(lencp),stat=ierr)
2208 IF (ierr/=0) GOTO 1000
2209 fxbcps = 0
2210 ALLOCATE (fxblm
2211 IF (ierr/=0) GOTO 1000
2212 fxblm = 0
2213 ALLOCATE (fxbfls(lenfls
2214 IF (ierr/=0) GOTO 1000
2215 fxbfls = 0
2216 ALLOCATE (fxbdls(lendls),stat=ierr)
2217 IF (ierr/=0) GOTO 1000
2218 fxbdls = 0
2219 ALLOCATE (fxbdep(lenvar),stat=ierr)
2220 IF (ierr/=0) GOTO 1000
2221 fxbdep = 0
2222 ALLOCATE (fxbvit(lenvar),stat=ierr)
2223 IF (ierr/=0) GOTO 1000
2224 fxbvit = 0
2225 ALLOCATE (fxbacc(lenvar),stat=ierr)
2226 IF (ierr/=0) GOTO 1000
2227 fxbacc = 0
2228 ALLOCATE (fxbrpm(lenrpm),stat=ierr)
2229 IF (ierr/=0) GOTO 1000
2230 fxbrpm = 0
2231 ALLOCATE (
fxbelm(lenelm),stat=ierr)
2232 IF (ierr/=0) GOTO 1000
2234 ALLOCATE (fxbsig(lensig),stat=ierr)
2235 IF (ierr/=0) GOTO 1000
2236 fxbsig = 0
2237 ALLOCATE (
fxbgrvi(lengrvi),stat=ierr)
2238 IF (ierr/=0) GOTO 1000
2240 ALLOCATE (fxbgrvr(lengrvr),stat=ierr)
2241 IF (ierr/=0) GOTO 1000
2242 fxbgrvr = 0
2243
2244 ALLOCATE (
eigipm(neipm,neig),stat=ierr)
2245 IF (ierr/=0) GOTO 1000
2247 ALLOCATE (eigrpm(nerpm,neig),stat
2248 IFGOTO 1000
2249 eigrpm = 0
2251 IF (ierr/=0) GOTO 1000
2253
2255 IF (ierr/=0) GOTO 1000
2256
2257 ALLOCATE (
retractor(nretractor),stat=ierr)
2258 IF (ierr/=0) GOTO 1000
2259
2260 IF (n_anchor_remote > 0) THEN
2262 IF (ierr/=0) GOTO 1000
2264 IF (ierr/=0) GOTO 1000
2265 ENDIF
2266 IF (n_anchor_remote_send > 0) THEN
2268 IF (ierr/=0) GOTO 1000
2270 IF (ierr/=0) GOTO 1000
2271 ENDIF
2272
2273 ALLOCATE (tm_dis(
lmax_dis*numnod),stat=ierr)
2274 IF (ierr/=0) GOTO 1000
2275 ALLOCATE (tm_vel(
lmax_vel*numnod),stat=ierr)
2276 IF (ierr/=0) GOTO 1000
2277 ALLOCATE (tm_nsig1(
lmax_nsig*numnod),stat=ierr)
2278 IF (ierrGOTO
2279 ALLOCATE (tm_nstra1(
lmax_nstra*numnod),stat=ierr)
2280 IF (ierr/=0) GOTO 1000
2282 IF (ierr/=0) GOTO 1000
2283 ALLOCATE (tm_nstra3(
lmax_nstra*numnod),stat=ierr)
2284 IF (ierr/=0) GOTO 1000
2286 ALLOCATE (tm_pdis(numnod),stat=ierr)
2287 IF (ierr/=0) GOTO 1000
2288 ENDIF
2290 ALLOCATE (tm_pvel
2291 IF (ierr/=0) GOTO 1000
2292 ENDIF
2294 ALLOCATE (tm_pnsig(2*numnod),stat=ierr)
2295 IF (ierr/=0) GOTO 1000
2296 ENDIF
2298 ALLOCATE (tm_pnstra(2*numnod
2299 IF (ierr/=0) GOTO 1000
2300 ENDIF
2302 ALLOCATE(gpstmp(6*numnod),p2tmp(2*numnod),stat=ierr
2303 IF (ierr/=0) GOTO 1000
2304 END IF
2305 IF (nconld > 0) THEN
2306 ALLOCATE (dpl0cld(6*nconld),stat=ierr)
2307 IF (ierr/=0) GOTO 1000
2308 dpl0cld = zero
2309 ALLOCATE (vel0cld(6*nconld),stat=ierr)
2310 IF (ierr/=0) GOTO 1000
2311 vel0cld = zero
2312 ELSE
2313 ALLOCATE(dpl0cld(0),vel0cld(0))
2314 END IF
2315
2317 IF (ierr/=0) GOTO 1000
2319 IF (ierr/=0) GOTO 1000
2320
2321
2322 CALL my_alloc(skews%SKEW,skews%N_SKEW_VAR,skews%TOTAL_SKEW_NUMBER)
2323
2324
2325
2326 IF(debug(6)==1)THEN
2328 END IF
2329
2330 RETURN
2331 1000 CONTINUE
2332 ierror = ierror + ierr
2333 ierr = ierr0
2334 IF (ierror/=0) THEN
2335 CALL ancmsg(msgid=20,anmode=aninfo)
2337 END IF
2338
subroutine multi_allocate(multi_fvm, h3d_data)
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 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 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 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 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)