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.  
  35. -INC PPARAM
  36. -INC CCOPTIO
  37. CBEGININCLUDE SMCHAEL
  38. SEGMENT MCHAEL
  39. POINTEUR IMACHE(N1).MELEME
  40. POINTEUR ICHEVA(N1).MCHEVA
  41. ENDSEGMENT
  42. SEGMENT MCHEVA
  43. REAL*8 VELCHE(NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM)
  44. ENDSEGMENT
  45. SEGMENT LCHEVA
  46. POINTEUR LISCHE(NBCHE).MCHEVA
  47. ENDSEGMENT
  48. CENDINCLUDE SMCHAEL
  49. INTEGER NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM,N1
  50. POINTEUR FC.MCHEVA
  51. POINTEUR LCOF.LCHEVA
  52. POINTEUR MYCOF.MCHEVA
  53. *
  54. INTEGER IMPR,IRET
  55. *
  56. * Executable statements
  57. *
  58. IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans ccgdiv'
  59. IF (.NOT.(IDIM.GE.1.AND.IDIM.LE.3)) THEN
  60. WRITE(IOIMP,*) 'IDIM=',IDIM,' ?'
  61. GOTO 9999
  62. ENDIF
  63. NLFC=FC.VELCHE(/6)
  64. NPFC=FC.VELCHE(/5)
  65. DO IIDIM=1,IDIM
  66. MYCOF=LCOF.LISCHE(IIDIM)
  67. NLC=MYCOF.VELCHE(/6)
  68. NPC=MYCOF.VELCHE(/5)
  69. DO ILFC=1,NLFC
  70. IF (NLC.EQ.1) THEN
  71. ILC=1
  72. ELSE
  73. ILC=ILFC
  74. ENDIF
  75. DO IPFC=1,NPFC
  76. IF (NPC.EQ.1) THEN
  77. IPC=1
  78. ELSE
  79. IPC=IPFC
  80. ENDIF
  81. FC.VELCHE(1,1,1,1,IPFC,ILFC)=
  82. $ FC.VELCHE(1,1,1,1,IPFC,ILFC)+
  83. $ MYCOF.VELCHE(1,1,1,IIDIM,IPC,ILC)
  84. ENDDO
  85. ENDDO
  86. ENDDO
  87. *
  88. * Normal termination
  89. *
  90. IRET=0
  91. RETURN
  92. *
  93. * Format handling
  94. *
  95. *
  96. * Error handling
  97. *
  98. 9999 CONTINUE
  99. IRET=1
  100. WRITE(IOIMP,*) 'An error was detected in subroutine ccgdiv'
  101. RETURN
  102. *
  103. * End of subroutine CCGDIV
  104. *
  105. END
  106.  
  107.  
  108.  

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