Télécharger locara.eso

Retour à la liste

Numérotation des lignes :

locara
  1. C LOCARA SOURCE OF166741 26/02/23 21:15:24 12446
  2.  
  3. C-----------------------------------------------------------------------
  4. C DESCRIPTION FONCTIONNELLE :
  5. C -------------------------
  6. C Calcul de la longueur caracteristique d'un element, determinee
  7. C comme la distance maximale entre deux noeuds de l'element
  8. C
  9. C INTERFACE :
  10. C ---------
  11. C ENTREE : IDIM : INTEGER, nombre de dimensions d'espace
  12. C ENTREE : NBNN : INTEGER, nombre de noeuds de l'element
  13. C ENTREE : XEL : tableau de REAL*8 de dimensions (3,NBNN),
  14. C coordonnees des noeuds de l'element
  15. C Ce tableau a ete rempli par la routine DOXE
  16. C appelee au prealable
  17. C SORTIE : LCARAC : REAL*8, longueur caracteristique de l'element
  18. C-----------------------------------------------------------------------
  19. SUBROUTINE LOCARA(IDIM,NBNN,XEL,LCARAC)
  20.  
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23.  
  24. C Arguments de l'interface
  25. C
  26. INTEGER IDIM, NBNN
  27. REAL*8 XEL(3,NBNN), LCARAC
  28. C
  29. C Variables locales
  30. C
  31. INTEGER INO1, INO2, IX
  32. REAL*8 DIST
  33. C-----------------------------------------------------------------------
  34. LCARAC=0.0D0
  35. C
  36. IF (NBNN.GT.1) THEN
  37. C
  38. C Determination de la distance maximale entre deux noeuds de
  39. C l'element
  40. DO 10 INO1=1,NBNN-1
  41. DO 11 INO2=INO1+1,NBNN
  42. DIST=0.0D0
  43. DO 12 IX=1,IDIM
  44. DIST=DIST+(XEL(IX,INO2)-XEL(IX,INO1))**2
  45. 12 CONTINUE
  46. LCARAC = MAX(LCARAC,DIST)
  47. 11 CONTINUE
  48. 10 CONTINUE
  49. C
  50. LCARAC = SQRT(LCARAC)
  51. ENDIF
  52. C
  53. RETURN
  54. END
  55.  
  56.  
  57.  

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