Télécharger la plaquette Cast3M


fr en

Télécharger dfnfr1.eso

Retour à la liste

Numérotation des lignes : 

dfnfr1

  1. C DFNFR1    SOURCE    GOUNAND   21/06/02    21:15:48     11022          
  2.       SUBROUTINE DFNFR1(NDDL,IESREF,IESREL,NPFF,NDPOGO,NBELEM,
  3.      $     NPJMI,
  4.      $     DFFPG,JMIJAC,
  5.      $     JDFFPG,
  6.      $     IMPR,IRET)
  7.       IMPLICIT REAL*8 (A-H,O-Z)
  8.       IMPLICIT INTEGER (I-N)
  9. C***********************************************************************
  10. C NOM         : DFNFR1
  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                      * NDPOGO (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, 25/09/03, refonte complète (modif SMTNLIN)
  44. C VERSION    : v1, 19/05/00, version initiale
  45. C HISTORIQUE : v1, 19/05/00, 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.  
  54. -INC PPARAM
  55. -INC CCOPTIO
  56.       INTEGER NDDL,IESREF,IESREL,NDPOGO,NBELEM
  57.       REAL*8 DFFPG (NDDL,IESREF,NPFF)
  58.       REAL*8 JMIJAC(IESREF,IESREL,NPJMI,NBELEM)
  59.       REAL*8 JDFFPG(NDDL,IESREL,NDPOGO,NBELEM)
  60. *
  61.       INTEGER IMPR,IRET
  62. *
  63.       INTEGER IDDL,ISREF,ISREL,IPOGO,IELEM
  64.       REAL*8 INVJAC
  65. *
  66. * Executable statements
  67. *
  68.       IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans dfnfr1'
  69.       DO 1 IELEM=1,NBELEM
  70.          DO 12 IPOGO=1,NDPOGO
  71.             IF (NPFF.EQ.1) THEN
  72.                IPFF=1
  73.             ELSE
  74.                IPFF=IPOGO
  75.             ENDIF
  76.             IF (NPJMI.EQ.1) THEN
  77.                IPJMI=1
  78.             ELSE
  79.                IPJMI=IPOGO
  80.             ENDIF
  81.             DO 122 ISREL=1,IESREL
  82.                DO 1222 ISREF=1,IESREF
  83.                   INVJAC=JMIJAC(ISREF,ISREL,IPJMI,IELEM)
  84.                   DO 12222 IDDL=1,NDDL
  85. *                     JDFFPG(IDDL,ISREL,IPOGO,IELEM)=
  86. *     $                    JDFFPG(IDDL,ISREL,IPOGO,IELEM)+
  87. *     $                    (DFFPG(IDDL,ISREF,IPOGO)
  88. *     $                    *JMIJAC(ISREF,ISREL,IPOGO,IELEM)
  89.                      JDFFPG(IDDL,ISREL,IPOGO,IELEM)=
  90.      $                    JDFFPG(IDDL,ISREL,IPOGO,IELEM)+
  91.      $                    (DFFPG(IDDL,ISREF,IPFF)
  92.      $                    *INVJAC)
  93. 12222             CONTINUE
  94.  1222          CONTINUE
  95.  122        CONTINUE
  96.  12      CONTINUE
  97.  1    CONTINUE
  98.  
  99. *
  100. * Normal termination
  101. *
  102.       IRET=0
  103.       RETURN
  104. *
  105. * Format handling
  106. *
  107. *
  108. * Error handling
  109. *
  110.  9999 CONTINUE
  111.       IRET=1
  112.       WRITE(IOIMP,*) 'An error was detected in subroutine dfnfr1'
  113.       RETURN
  114. *
  115. * End of subroutine DFNFR1
  116. *
  117.       END
  118.  
  119.  
  120.  
  121.  
  122.  

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