Télécharger demch1.eso

Retour à la liste

Numérotation des lignes :

demch1
  1. C DEMCH1 SOURCE PV 20/03/24 21:17:06 10554
  2. C UTILITAIRE DE DEMCHA (GESTION NOEUD MILIEU)
  3. C
  4. SUBROUTINE DEMCH1(I1,I3,I2,KONCOU,KONFIN,KONMIL,KONSUI,KONPOS)
  5. IMPLICIT INTEGER(I-N)
  6. IMPLICIT REAL*8(A-H,O-Z)
  7.  
  8. -INC PPARAM
  9. -INC CCOPTIO
  10. -INC SMCOORD
  11. SEGMENT KONPOS(NBPTS)
  12. SEGMENT KONFIN(IKOUR)
  13. SEGMENT KONMIL(IKOUR)
  14. SEGMENT KONSUI(IKOUR)
  15. J1=MIN(I1,I3)
  16. J3=MAX(I1,I3)
  17. ITFA=KONPOS(J1)
  18. IF (ITFA.EQ.0) GOTO 90
  19. 85 CONTINUE
  20. ITF=KONSUI(ITFA)
  21. IF (KONFIN(ITFA).EQ.J3) GOTO 80
  22. IF (ITF.EQ.0) GOTO 90
  23. ITFA=ITF
  24. GOTO 85
  25. 90 KONCOU=KONCOU+1
  26. IF (KONCOU.GE.KONFIN(/1)) THEN
  27. IKOUR=KONCOU+500
  28. SEGADJ KONFIN,KONMIL,KONSUI
  29. ENDIF
  30. IF (ITFA.EQ.0) THEN
  31. KONPOS(J1)=KONCOU
  32. ELSE
  33. KONSUI(ITFA)=KONCOU
  34. ENDIF
  35. KONFIN(KONCOU)=J3
  36. * CREATION DE NOEUD
  37. segact mcoord*mod
  38. NBPTS=nbpts+1
  39. SEGADJ MCOORD
  40. DO 100 ICOOR=1,IDIM+1
  41. XCOOR((NBPTS-1)*(IDIM+1)+ICOOR)=
  42. $ (XCOOR((J1-1)*(IDIM+1)+ICOOR)
  43. $ +XCOOR((J3-1)*(IDIM+1)+ICOOR))/2.D0
  44. 100 CONTINUE
  45. KONMIL(KONCOU)=NBPTS
  46. ITFA=KONCOU
  47. 80 CONTINUE
  48. I2=KONMIL(ITFA)
  49. END
  50.  
  51.  
  52.  
  53.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales