Télécharger la plaquette Cast3M


fr en

Télécharger dfnff1.eso

Retour à la liste

Numérotation des lignes : 

dfnff1

  1. C DFNFF1    SOURCE    GOUNAND   21/06/02    21:15:46     11022          
  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 SMTNLIN)
  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.  
  54. -INC PPARAM
  55. -INC CCOPTIO
  56.       INTEGER NDDL,IESREF,IESREL,NBPOGO,NBELEV,NBELEF,NBELFV,NPJMI
  57.       REAL*8 DFFPG (NDDL,IESREF,NBPOGO,NBELFV)
  58.       REAL*8 JMIJAC(IESREF,IESREL,NPJMI,NBELEF)
  59.       LOGICAL SSFACT(NBELFV,NBELEV)
  60.       REAL*8 JDFFPG(NDDL,IESREL,NBPOGO,NBELEF)
  61. *
  62.       INTEGER IMPR,IRET
  63. *
  64.       INTEGER IDDL,ISREF,ISREL,IBPOGO,IBELEV,IBELEF,IBELFV,IPJMI
  65. *
  66. * Executable statements
  67. *
  68.       IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans dfnff1'
  69.       IBELEF=0
  70.       DO 1 IBELEV=1,NBELEV
  71.          DO IBELFV=1,NBELFV
  72.             IF (SSFACT(IBELFV,IBELEV)) THEN
  73.                IBELEF=IBELEF+1
  74.                DO 12 IBPOGO=1,NBPOGO
  75.                   IF (NPJMI.EQ.1) THEN
  76.                      IPJMI=1
  77.                   ELSE
  78.                      IPJMI=IBPOGO
  79.                   ENDIF
  80.                   DO 122 ISREL=1,IESREL
  81.                      DO 1222 ISREF=1,IESREF
  82.                         DO 12222 IDDL=1,NDDL
  83.                            JDFFPG(IDDL,ISREL,IBPOGO,IBELEF)=
  84.      $                          JDFFPG(IDDL,ISREL,IBPOGO,IBELEF)+
  85.      $                          (DFFPG(IDDL,ISREF,IBPOGO,IBELFV)
  86.      $                          *JMIJAC(ISREF,ISREL,IPJMI,IBELEF))
  87. 12222                   CONTINUE
  88.  1222                CONTINUE
  89.  122              CONTINUE
  90.  12            CONTINUE
  91.             ENDIF
  92.          ENDDO
  93.  1    CONTINUE
  94.  
  95. *
  96. * Normal termination
  97. *
  98.       IRET=0
  99.       RETURN
  100. *
  101. * Format handling
  102. *
  103. *
  104. * Error handling
  105. *
  106.  9999 CONTINUE
  107.       IRET=1
  108.       WRITE(IOIMP,*) 'An error was detected in subroutine dfnff1'
  109.       RETURN
  110. *
  111. * End of subroutine DFNFF1
  112. *
  113.       END
  114.  
  115.  
  116.  
  117.  

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