Télécharger dfnff1.eso

Retour à la liste

Numérotation des lignes :

  1. C DFNFF1 SOURCE GOUNAND 06/08/04 21:15:24 5520
  2. SUBROUTINE DFNFF1(NDDL,IESREF,IESREL,NBPOGO,NBELEV,NBELEF,NBELFV,
  3. $ NPJMI,
  4. $ DFFPG,JMIJAC,SSFACT,
  5. $ JDFFPG,
  6. $ IMPR,IRET)
  7. IMPLICIT REAL*8 (A-H,O-Z)
  8. IMPLICIT INTEGER (I-N)
  9. C***********************************************************************
  10. C NOM : DFNFF1
  11. C DESCRIPTION :
  12. C
  13. C
  14. C
  15. C LANGAGE : Fortran 77 (sauf E/S)
  16. C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  17. C mél : gounand@semt2.smts.cea.fr
  18. C***********************************************************************
  19. C APPELES : -
  20. C APPELE PAR : DFNFR
  21. C***********************************************************************
  22. C ENTREES : * NDDL (type entier) : nombre de ddl (par
  23. C élément) de la transformation géométrique.
  24. C * IESREF (type entier) : dimension de l'espace de
  25. C référence.
  26. C * IESREL (type entier) : dimension de l'espace
  27. C réel (i.e. géométrique).
  28. C * NBPOGO (type entier) : nombre de points
  29. C d'intégration.
  30. C * NBELEM (type entier) : nombre d'éléments du
  31. C maillage élémentaire courant.
  32. C * DFFPG (type réel) : valeurs des dérivées
  33. C premières des fonctions d'interpolation aux
  34. C points de gauss sur l'élément de référence.
  35. C * JMIJAC (type réel) : valeurs de l'inverse de
  36. C la matrice jacobienne aux points de Gauss sur
  37. C le maillage élémentaire.
  38. C ENTREES/SORTIES : * JDFFPG (type réel) : valeurs des dérivées
  39. C premières des fonctions d'interpolation
  40. C aux points de gauss sur l'élément réel.
  41. C SORTIES : -
  42. C***********************************************************************
  43. C VERSION : v2, 03/10/03, refonte complète (modif SMPOUET)
  44. C VERSION : v1, 17/01/03, version initiale
  45. C HISTORIQUE : v1, 17/01/03, création
  46. C HISTORIQUE :
  47. C HISTORIQUE :
  48. C***********************************************************************
  49. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  50. C en cas de modification de ce sous-programme afin de faciliter
  51. C la maintenance !
  52. C***********************************************************************
  53. -INC CCOPTIO
  54. INTEGER NDDL,IESREF,IESREL,NBPOGO,NBELEV,NBELEF,NBELFV,NPJMI
  55. REAL*8 DFFPG (NDDL,IESREF,NBPOGO,NBELFV)
  56. REAL*8 JMIJAC(IESREF,IESREL,NPJMI,NBELEF)
  57. LOGICAL SSFACT(NBELFV,NBELEV)
  58. REAL*8 JDFFPG(NDDL,IESREL,NBPOGO,NBELEF)
  59. *
  60. INTEGER IMPR,IRET
  61. *
  62. INTEGER IDDL,ISREF,ISREL,IBPOGO,IBELEV,IBELEF,IBELFV,IPJMI
  63. *
  64. * Executable statements
  65. *
  66. IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans dfnff1'
  67. IBELEF=0
  68. DO 1 IBELEV=1,NBELEV
  69. DO IBELFV=1,NBELFV
  70. IF (SSFACT(IBELFV,IBELEV)) THEN
  71. IBELEF=IBELEF+1
  72. DO 12 IBPOGO=1,NBPOGO
  73. IF (NPJMI.EQ.1) THEN
  74. IPJMI=1
  75. ELSE
  76. IPJMI=IBPOGO
  77. ENDIF
  78. DO 122 ISREL=1,IESREL
  79. DO 1222 ISREF=1,IESREF
  80. DO 12222 IDDL=1,NDDL
  81. JDFFPG(IDDL,ISREL,IBPOGO,IBELEF)=
  82. $ JDFFPG(IDDL,ISREL,IBPOGO,IBELEF)+
  83. $ (DFFPG(IDDL,ISREF,IBPOGO,IBELFV)
  84. $ *JMIJAC(ISREF,ISREL,IPJMI,IBELEF))
  85. 12222 CONTINUE
  86. 1222 CONTINUE
  87. 122 CONTINUE
  88. 12 CONTINUE
  89. ENDIF
  90. ENDDO
  91. 1 CONTINUE
  92.  
  93. *
  94. * Normal termination
  95. *
  96. IRET=0
  97. RETURN
  98. *
  99. * Format handling
  100. *
  101. *
  102. * Error handling
  103. *
  104. 9999 CONTINUE
  105. IRET=1
  106. WRITE(IOIMP,*) 'An error was detected in subroutine dfnff1'
  107. RETURN
  108. *
  109. * End of subroutine DFNFF1
  110. *
  111. END
  112.  
  113.  
  114.  

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