Télécharger la plaquette Cast3M


fr en

Télécharger kopidr.eso

Retour à la liste

Numérotation des lignes : 

kopidr

  1. C KOPIDR    SOURCE    FANDEUR   22/01/19    21:15:08     11256          
  2.       SUBROUTINE KOPIDR(IGEOM,LPRIM,MRIGID)
  3.       IMPLICIT REAL*8 (A-H,O-Z)
  4.       IMPLICIT INTEGER (I-N)
  5. C***********************************************************************
  6. C NOM         : KOPIDR
  7. C DESCRIPTION : Transforme un CHPOINt MCHPOI en matrice
  8. C               diagonale MRIGID
  9. C
  10. C
  11. C
  12. C LANGAGE     : ESOPE
  13. C AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  14. C               mél : gounand@semt2.smts.cea.fr
  15. C***********************************************************************
  16. C ENTREES            : IGEOM,LPRIM
  17. C ENTREES/SORTIES    :
  18. C SORTIES            : MRIGID
  19. C***********************************************************************
  20. C VERSION    : v1, 10/05/2011, version initiale
  21. C HISTORIQUE : v1, 10/05/2011, création
  22. C HISTORIQUE :
  23. C HISTORIQUE :
  24. C***********************************************************************
  25.  
  26. -INC PPARAM
  27. -INC CCOPTIO
  28. -INC SMRIGID
  29. -INC SMELEME
  30. -INC SMLMOTS
  31. -INC CCHAMP
  32. *
  33. * Executable statements
  34. *
  35. C
  36. C
  37. C**** On controle que le MELEME soit de type POI1
  38. C     Si non changer
  39. C
  40.       MELEME = IGEOM
  41.       SEGACT MELEME
  42.       NBSOUS = MELEME.LISOUS(/1)
  43.       NTYP = MELEME.ITYPEL
  44. C
  45.       IF ((NBSOUS.NE.0).OR.(NTYP.NE.1)) THEN
  46. * In CHANGE : SEGINI MELEME
  47.          CALL CHANGE(MELEME,1)
  48.          IPT1 = IGEOM
  49.          SEGDES IPT1
  50.          IF (IERR.NE.0) RETURN
  51.       ENDIF
  52. *
  53.       MLMOTS=LPRIM
  54.       SEGACT MLMOTS
  55.       NRIGEL=MOTS(/2)
  56.       SEGINI MRIGID
  57.       MTYMAT='DIAGONAL'
  58.       IFORIG=IFOUR
  59. C
  60.       NBEL=NUM(/2)
  61.       DO I=1,NRIGEL
  62.          COERIG(I)=1.D0
  63.          IRIGEL(1,I)=MELEME
  64.          NLIGRP=1
  65.          NLIGRD=1
  66.          SEGINI DESCR
  67.          LISINC(1)=MOTS(I)
  68.          CALL PLACE(NOMDD,LNOMDD,idx,MOTS(I))
  69.          IF (idx.NE.0) THEN
  70.             LISDUA(1)=NOMDU(idx)
  71.          ELSE
  72.             LISDUA(1)=MOTS(I)
  73.          ENDIF
  74.          NOELEP(1)=1
  75.          NOELED(1)=1
  76.          SEGDES DESCR
  77.          IRIGEL(3,I)=DESCR
  78.          NELRIG=NBEL
  79.          SEGINI XMATRI
  80.          DO IBEL=1,NBEL
  81.             RE(1,1,IBEL)=1.D0
  82.          ENDDO
  83.          SEGDES XMATRI
  84.          IRIGEL(4,I)=XMATRI
  85.       ENDDO
  86.       SEGDES MRIGID
  87.       SEGDES MLMOTS
  88.       SEGDES MELEME
  89. *
  90. * Normal termination
  91. *
  92.       RETURN
  93. *
  94. * End of subroutine KOPIDR
  95. *
  96.       END
  97.  
  98.  
  99.  

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