Télécharger etoilm.eso

Retour à la liste

Numérotation des lignes :

etoilm
  1. C ETOILM SOURCE GOUNAND 21/03/31 21:15:06 10931
  2. SUBROUTINE ETOILM(NODE,IPT1,
  3. $ IPT3)
  4. IMPLICIT REAL*8 (A-H,O-Z)
  5. IMPLICIT INTEGER (I-N)
  6. C***********************************************************************
  7. C NOM : ETOILM
  8. C
  9. C DESCRIPTION : Etant donné un maillage complexe IPT1 constitue
  10. C d'elements de type POI1, SEG2, TRI3 ou QUA4 et un noeud NODE, on
  11. C construit IPT3 le maillage obtenu par étoilement de IPT1 avec
  12. C NODE. Interface à ETOIL1.
  13. C L'étoilement est fait avec les éléments de IPT1 qui ne
  14. C contiennent pas NODE.
  15. C IPT3 et ses sous-maillages actifs*MOD.
  16. C
  17. C
  18. C LANGAGE : ESOPE
  19. C AUTEUR : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  20. C mél : gounand@semt2.smts.cea.fr
  21. C***********************************************************************
  22. C VERSION : v1, 27/04/2020, version initiale
  23. C HISTORIQUE : v1, 27/04/2020, création
  24. C HISTORIQUE :
  25. C***********************************************************************
  26. -INC PPARAM
  27. -INC CCOPTIO
  28. -INC SMELEME
  29. *
  30. * Executable statements
  31. *
  32. SEGACT IPT1
  33. NSOUS=IPT1.LISOUS(/1)
  34. IF (NSOUS.NE.0) THEN
  35. NBNN=0
  36. NBELEM=0
  37. NBSOUS=0
  38. NBREF=0
  39. SEGINI IPT3
  40. ENDIF
  41. *
  42. DO ISOUS=1,MAX(1,NSOUS)
  43. IF (NSOUS.EQ.0) THEN
  44. IPT2=IPT1
  45. ELSE
  46. IPT2=IPT1.LISOUS(ISOUS)
  47. SEGACT IPT2
  48. ENDIF
  49. CALL ETOIL1(NODE,IPT2,IPT4)
  50. IF (IERR.NE.0) RETURN
  51. IF (NSOUS.EQ.0) THEN
  52. IPT3=IPT4
  53. ELSE
  54. IPT3.LISOUS(ISOUS)=IPT4
  55. ENDIF
  56. ENDDO
  57. RETURN
  58. *
  59. * End of subroutine ETOILM
  60. *
  61. END
  62.  
  63.  
  64.  

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