demch1
C DEMCH1 SOURCE PV 20/03/24 21:17:06 10554
C UTILITAIRE DE DEMCHA (GESTION NOEUD MILIEU)
C
IMPLICIT INTEGER(I-N)
IMPLICIT REAL*8(A-H,O-Z)
-INC PPARAM
-INC CCOPTIO
-INC SMCOORD
SEGMENT KONPOS(NBPTS)
SEGMENT KONFIN(IKOUR)
SEGMENT KONMIL(IKOUR)
SEGMENT KONSUI(IKOUR)
J1=MIN(I1,I3)
J3=MAX(I1,I3)
ITFA=KONPOS(J1)
IF (ITFA.EQ.0) GOTO 90
85 CONTINUE
ITF=KONSUI(ITFA)
IF (KONFIN(ITFA).EQ.J3) GOTO 80
IF (ITF.EQ.0) GOTO 90
ITFA=ITF
GOTO 85
90 KONCOU=KONCOU+1
IF (KONCOU.GE.KONFIN(/1)) THEN
IKOUR=KONCOU+500
SEGADJ KONFIN,KONMIL,KONSUI
ENDIF
IF (ITFA.EQ.0) THEN
KONPOS(J1)=KONCOU
ELSE
KONSUI(ITFA)=KONCOU
ENDIF
KONFIN(KONCOU)=J3
* CREATION DE NOEUD
segact mcoord*mod
NBPTS=nbpts+1
SEGADJ MCOORD
DO 100 ICOOR=1,IDIM+1
XCOOR((NBPTS-1)*(IDIM+1)+ICOOR)=
$ (XCOOR((J1-1)*(IDIM+1)+ICOOR)
$ +XCOOR((J3-1)*(IDIM+1)+ICOOR))/2.D0
100 CONTINUE
KONMIL(KONCOU)=NBPTS
ITFA=KONCOU
80 CONTINUE
END
					© Cast3M 2003 - Tous droits réservés.
					Mentions légales