Télécharger dcoga1.eso

Retour à la liste

Numérotation des lignes :

dcoga1
  1. C DCOGA1 SOURCE GOUNAND 05/12/21 21:17:58 5281
  2. SUBROUTINE DCOGA1(NDDL,IESREL,NBPOGO,NLDF,NLCO,NBELEM,
  3. $ JCOEFF,JDFFPG,
  4. $ JDCOFG,
  5. $ IMPR,IRET)
  6. IMPLICIT REAL*8 (A-H,O-Z)
  7. IMPLICIT INTEGER (I-N)
  8. C***********************************************************************
  9. C NOM : DCOGA1
  10. C PROJET :
  11. C DESCRIPTION : Calcul des derivees du coefficient aux points de Gauss
  12. C d'un élément
  13. C de référence pour chaque élément réel.
  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 : DCOGAU
  21. C***********************************************************************
  22. C ENTREES :
  23. C ENTREES/SORTIES :
  24. C SORTIES : -
  25. C CODE RETOUR (IRET) : = 0 si tout s'est bien passé
  26. C***********************************************************************
  27. C VERSION : v3.1, 30/07/04, possiblité de travailler
  28. C dans l'espace de référence
  29. C VERSION : v1, 26/08/03, version initiale
  30. C HISTORIQUE : v1, 26/08/03, création
  31. C HISTORIQUE :
  32. C HISTORIQUE :
  33. C***********************************************************************
  34. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  35. C en cas de modification de ce sous-programme afin de faciliter
  36. C la maintenance !
  37. C***********************************************************************
  38.  
  39. -INC PPARAM
  40. -INC CCOPTIO
  41. INTEGER NDDL,IESREL,NBPOGO,NBELEM
  42. REAL*8 JCOEFF(NDDL,NLCO)
  43. REAL*8 JDFFPG(NDDL,IESREL,NBPOGO,NLDF)
  44. REAL*8 JDCOFG(IESREL,NBPOGO,NBELEM)
  45. *
  46. INTEGER IMPR,IRET
  47. INTEGER IELEM,IPOGO,JESREL,IDDL
  48. *
  49. * Executable statements
  50. *
  51. IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans dcoga1'
  52. DO IELEM=1,NBELEM
  53. IF (NLCO.EQ.1) THEN
  54. ILCO=1
  55. ELSE
  56. ILCO=IELEM
  57. ENDIF
  58. IF (NLDF.EQ.1) THEN
  59. ILDF=1
  60. ELSE
  61. ILDF=IELEM
  62. ENDIF
  63. DO IPOGO=1,NBPOGO
  64. DO JESREL=1,IESREL
  65. DO IDDL=1,NDDL
  66. JDCOFG(JESREL,IPOGO,IELEM)=
  67. $ JDCOFG(JESREL,IPOGO,IELEM)
  68. $ + (JCOEFF(IDDL,ILCO)
  69. $ * JDFFPG(IDDL,JESREL,IPOGO,ILDF))
  70. ENDDO
  71. ENDDO
  72. ENDDO
  73. ENDDO
  74. *
  75. * Normal termination
  76. *
  77. IRET=0
  78. RETURN
  79. *
  80. * Format handling
  81. *
  82. *
  83. * Error handling
  84. *
  85. 9999 CONTINUE
  86. IRET=1
  87. WRITE(IOIMP,*) 'An error was detected in subroutine dcoga1'
  88. RETURN
  89. *
  90. * End of subroutine DCOGA1
  91. *
  92. END
  93.  
  94.  
  95.  

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