Télécharger la plaquette Cast3M


fr en

Télécharger ccgadv.eso

Retour à la liste

Numérotation des lignes : 

ccgadv

  1. C CCGADV    SOURCE    GOUNAND   26/01/09    21:15:03     12441          
  2.       SUBROUTINE CCGADV(LCOF,NOMLOI,
  3.      $     FC,
  4.      $     IMPR,IRET)
  5.       IMPLICIT REAL*8 (A-H,O-Z)
  6.       IMPLICIT INTEGER (I-N)
  7. C***********************************************************************
  8. C NOM         : CCGADV
  9. C DESCRIPTION : Calcul de la loi de comportement aux points de Gauss :
  10. C               VOLORI
  11. C
  12. C LANGAGE     : ESOPE
  13. C AUTEUR      : Stéphane GOUNAND (CEA/DRN/DMT/SEMT/LTMF)
  14. C               mél : gounand@semt2.smts.cea.fr
  15. C***********************************************************************
  16. C APPELES          :
  17. C APPELE PAR       :
  18. C***********************************************************************
  19. C ENTREES            :
  20. C ENTREES/SORTIES    :
  21. C SORTIES            : -
  22. C TRAVAIL            :
  23. C***********************************************************************
  24. C VERSION    : v1, 04/08/04, version initiale
  25. C HISTORIQUE : v1, 04/08/04, création
  26. C HISTORIQUE :
  27. C HISTORIQUE :
  28. C***********************************************************************
  29. C Prière de PRENDRE LE TEMPS de compléter les commentaires
  30. C en cas de modification de ce sous-programme afin de faciliter
  31. C la maintenance !
  32. C***********************************************************************
  33. -INC PPARAM
  34. -INC CCOPTIO
  35. -INC TNLIN
  36. *     -INC SMCHAEL
  37.       INTEGER NBLIG,NBCOL,N2LIG,N2COL,NBPOI,NBELM,N1
  38.       POINTEUR FC.MCHEVA
  39.       POINTEUR LCOF.LCHEVA
  40.       POINTEUR JMAJAC.MCHEVA
  41.       POINTEUR JMIJAC.MCHEVA
  42.       POINTEUR JDTJAC.MCHEVA
  43.       CHARACTER*8 NOMLOI
  44.       INTEGER ICOF
  45. *
  46. -INC TMXMAT
  47.       POINTEUR JAC.MXMAT
  48.       POINTEUR JM1.MXMAT
  49.       POINTEUR M1.MXMAT
  50.       POINTEUR M2.MXMAT
  51.       POINTEUR M3.MXMAT
  52. *
  53.       INTEGER IMPR,IRET
  54. *
  55. * Executable statements
  56. *
  57.       IF (IMPR.GT.1) WRITE(IOIMP,*) 'Entrée dans ccgadv'
  58. C      IF (.NOT.(IDIM.EQ.1)) THEN
  59. C         WRITE(IOIMP,*) 'IDIM=',IDIM,' ?'
  60. C         GOTO 9999
  61. C      ENDIF
  62.       NLFC=FC.WELCHE(/6)
  63.       NPFC=FC.WELCHE(/5)
  64.       ICOF=0
  65. *
  66.       ICOF=ICOF+1
  67.       JMAJAC=LCOF.LISCHE(ICOF)
  68.  
  69.       ICOF=ICOF+1
  70.       JMIJAC=LCOF.LISCHE(ICOF)
  71.  
  72.       ICOF=ICOF+1
  73.       JDTJAC=LCOF.LISCHE(ICOF)
  74.       NLJD=JDTJAC.WELCHE(/6)
  75.       NPJD=JDTJAC.WELCHE(/5)
  76.       DO ILFC=1,NLFC
  77.          IF (NLJD.EQ.1) THEN
  78.             ILJD=1
  79.          ELSE
  80.             ILJD=ILFC
  81.          ENDIF
  82.          DO IPFC=1,NPFC
  83.             IF (NPJD.EQ.1) THEN
  84.                IPJD=1
  85.             ELSE
  86.                IPJD=IPFC
  87.             ENDIF
  88.             DET  =JDTJAC.WELCHE(1,1,1,1,IPJD,ILJD)
  89. *            SDET =SIGN(1.D0,DET)
  90.             IF (NOMLOI.EQ.'VOLORI  ') THEN
  91.                CONTRI=DET
  92.             ELSE
  93.                WRITE(IOIMP,*) 'Erreur grave'
  94.                GOTO 9999
  95.             ENDIF
  96.             FC.WELCHE(1,1,1,1,IPFC,ILFC)=
  97.      $           FC.WELCHE(1,1,1,1,IPFC,ILFC)+
  98.      $           CONTRI
  99.          ENDDO
  100.       ENDDO
  101. *
  102. * Normal termination
  103. *
  104.       IRET=0
  105.       RETURN
  106. *
  107. * Format handling
  108. *
  109. *
  110. * Error handling
  111. *
  112.  9999 CONTINUE
  113.       IRET=1
  114.       WRITE(IOIMP,*) 'An error was detected in subroutine ccgadv'
  115.       RETURN
  116. *
  117. * End of subroutine CCGADV
  118. *
  119.       END
  120.  
  121.  

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