1 : C MEULTI SOURCE PV 13/04/15 21:15:05 7764 2 : SUBROUTINE MEULTI(MELEME,MELETI,TYPE) 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8 (A-H,O-Z) 5 : 6 : CHARACTER*8 TYPE 7 : DIMENSION AK(4) 8 : DATA AK/1.D30,-1.D30,2.D24,0.D0/ 9 : -INC CCOPTIO 10 : -INC SMELEME 11 : POINTEUR MELETI.MELEME 12 : -INC SMCOORD 13 : 14 : MELETI=0 14 : CALL ECROBJ('MAILLAGE',MELEME) 14 : CALL ECRCHA('POI1') 14 : CALL PRCHAN 14 : CALL LIROBJ('MAILLAGE',MELETI,1,IRETOU) 15 : IF(IRETOU.EQ.0)RETURN 16 : segini,ipt1=meleti 17 : meleti=ipt1 18 : SEGACT MELETI*MOD 19 : NPT=MELETI.NUM(/2) 20 : 21 : NBPT0=XCOOR(/1)/(IDIM+1) 22 : NBPTS=NBPT0+NPT 23 : SEGADJ MCOORD 24 : 25 : IP=NBPT0 26 : DO 1 I=1,NPT 27 : IP=IP+1 28 : I1=MELETI.NUM(1,I) 29 : DO 2 K=1,(IDIM+1) 30 : XCOOR((IP-1)*(IDIM+1) +K)=XCOOR((I1-1)*(IDIM+1) +K)+AK(K) 31 : 2 CONTINUE 32 : MELETI.NUM(1,I)=IP 33 : 1 CONTINUE 34 : SEGDES MELETI 35 : RETURN 36 : END 37 : 38 : 39 :
© Cast3M 2003 - All rights reserved.
Disclaimer