Télécharger ccgdiv.eso

Retour à la liste

Numérotation des lignes :

  1. C CCGDIV SOURCE GOUNAND 05/12/21 21:16:26 5281
  2. SUBROUTINE CCGDIV(LCOF,
  3. $ FC,
  4. $ IMPR,IRET)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6. IMPLICIT INTEGER (I-N)
  7. C***********************************************************************
  8. C NOM : CCGDIV
  9. C DESCRIPTION : Calcul de la loi de comportement aux points de Gauss :
  10. C une divergence
  11. C
  12. C
  13. C LANGAGE : ESOPE
  14. C AUTEUR : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  15. C mél : gounand@semt2.smts.cea.fr
  16. C***********************************************************************
  17. C APPELES :
  18. C APPELE PAR :
  19. C***********************************************************************
  20. C ENTREES :
  21. C ENTREES/SORTIES :
  22. C SORTIES : -
  23. C TRAVAIL :
  24. C***********************************************************************
  25. C VERSION : v1, 11/05/04, version initiale
  26. C HISTORIQUE : v1, 11/05/04, création
  27. C HISTORIQUE :
  28. C HISTORIQUE :
  29. C***********************************************************************
  30. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  31. C en cas de modification de ce sous-programme afin de faciliter
  32. C la maintenance !
  33. C***********************************************************************
  34. -INC CCOPTIO
  35. CBEGININCLUDE SMCHAEL
  36. SEGMENT MCHAEL
  37. POINTEUR IMACHE(N1).MELEME
  38. POINTEUR ICHEVA(N1).MCHEVA
  39. ENDSEGMENT
  40. SEGMENT MCHEVA
  41. REAL*8 VELCHE(NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM)
  42. ENDSEGMENT
  43. SEGMENT LCHEVA
  44. POINTEUR LISCHE(NBCHE).MCHEVA
  45. ENDSEGMENT
  46. CENDINCLUDE SMCHAEL
  47. INTEGER NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM,N1
  48. POINTEUR FC.MCHEVA
  49. POINTEUR LCOF.LCHEVA
  50. POINTEUR MYCOF.MCHEVA
  51. *
  52. INTEGER IMPR,IRET
  53. *
  54. * Executable statements
  55. *
  56. IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans ccgdiv'
  57. IF (.NOT.(IDIM.GE.1.AND.IDIM.LE.3)) THEN
  58. WRITE(IOIMP,*) 'IDIM=',IDIM,' ?'
  59. GOTO 9999
  60. ENDIF
  61. NLFC=FC.VELCHE(/6)
  62. NPFC=FC.VELCHE(/5)
  63. DO IIDIM=1,IDIM
  64. MYCOF=LCOF.LISCHE(IIDIM)
  65. NLC=MYCOF.VELCHE(/6)
  66. NPC=MYCOF.VELCHE(/5)
  67. DO ILFC=1,NLFC
  68. IF (NLC.EQ.1) THEN
  69. ILC=1
  70. ELSE
  71. ILC=ILFC
  72. ENDIF
  73. DO IPFC=1,NPFC
  74. IF (NPC.EQ.1) THEN
  75. IPC=1
  76. ELSE
  77. IPC=IPFC
  78. ENDIF
  79. FC.VELCHE(1,1,1,1,IPFC,ILFC)=
  80. $ FC.VELCHE(1,1,1,1,IPFC,ILFC)+
  81. $ MYCOF.VELCHE(1,1,1,IIDIM,IPC,ILC)
  82. ENDDO
  83. ENDDO
  84. ENDDO
  85. *
  86. * Normal termination
  87. *
  88. IRET=0
  89. RETURN
  90. *
  91. * Format handling
  92. *
  93. *
  94. * Error handling
  95. *
  96. 9999 CONTINUE
  97. IRET=1
  98. WRITE(IOIMP,*) 'An error was detected in subroutine ccgdiv'
  99. RETURN
  100. *
  101. * End of subroutine CCGDIV
  102. *
  103. END
  104.  
  105.  
  106.  

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