Télécharger la plaquette Cast3M


fr en

Télécharger bgchel.eso

Retour à la liste

Numérotation des lignes : 

bgchel

  1. C BGCHEL    SOURCE    PV090527  25/01/07    14:42:26     12115          
  2.  
  3. C=====================================================================
  4. C=====================================================================
  5. C BGCHEL    SOURCE   MARINO ARROYO   99/05/21
  6.       SUBROUTINE BGCHEL(IPO1,IEPS,T0,X0,IRET)
  7. C=====================================================================
  8. C
  9. C     EFFECTUE DIVERSES OPERATIONS SUR DES MCHAML
  10. C  ENTREES
  11. C     IPO1=POINTEUR SUR LE CHAMPS PAR ELEMENT
  12. C     IEPS= 1 CONDUCTIVITE
  13. C         = 2 DERIVEE CONDUCTIVITE
  14. C         = 3 CAPACITE
  15. C         = 4 DERIVEE CAPACITE
  16. C     T0= 1 TEMPERATURE DE REFERENCE
  17. C     X0= CARACTERISTIQUE DE REFERENCE
  18. C  SORTIES
  19. C     IRET=POINTEUR SUR LE CHAMP RESULTANT
  20. C         =0  SI OPERATION IMPOSSSIBLE
  21. C
  22. C=====================================================================
  23. C
  24.       IMPLICIT INTEGER(I-N)
  25.       IMPLICIT REAL*8(A-H,O-Z)
  26. C
  27. -INC SMCHAML
  28. -INC SMCOORD
  29.  
  30. -INC PPARAM
  31. -INC CCOPTIO
  32. C
  33.       IRET=0
  34.       MCHEL1 = IPO1
  35.       SEGINI,MCHELM=MCHEL1
  36.       IRET = MCHELM
  37.       NSOUS = IMACHE(/1)
  38. *
  39.       DO 72 IA=1,NSOUS
  40. *
  41.          MCHAM1=ICHAML(IA)
  42.          SEGINI,MCHAML=MCHAM1
  43.          ICHAML(IA)=MCHAML
  44. *
  45.          DO 75 ICOMP=1,IELVAL(/1)
  46. *
  47.             MELVA1 = IELVAL(ICOMP)
  48.             SEGINI,MELVAL=MELVA1
  49.             IELVAL(ICOMP) = MELVAL
  50.             SEGACT MELVA1
  51. *
  52.             IF (TYPCHE(ICOMP).EQ.'REAL*8') GOTO 1666
  53. C
  54.             GOTO 166
  55. C
  56.  1666       N1PTEL=VELCHE(/1)
  57.             N1EL  =VELCHE(/2)
  58.             GOTO( 1, 2, 3, 4, 5, 6),IEPS
  59.  66         IRET=0
  60.             GOTO 666
  61.  1          CONTINUE
  62. C
  63. C         CONDUCTIVITE
  64. C
  65.             DO 1000 IB=1,N1PTEL
  66.                DO 1001 ID=1,N1EL
  67.                   XX=MELVA1.VELCHE(IB,ID)
  68.                   XX=XX*T0
  69.                   VELCHE(IB,ID)=(FBG1(XX))/X0
  70.  1001          CONTINUE
  71.  1000       CONTINUE
  72.             GOTO 166
  73. C
  74. C         DERIVEE CONDUCTIVITE
  75. C
  76.  2          CONTINUE
  77.             DO 20 IB=1,N1PTEL
  78.                DO 21 ID=1,N1EL
  79.                   XX =MELVA1.VELCHE(IB,ID)
  80.                   XX=XX*T0
  81.                   VELCHE(IB,ID)=(FBG2(XX))*T0/X0
  82.  21            CONTINUE
  83.  20         CONTINUE
  84.             GOTO 166
  85. C
  86. C         CAPACITE
  87. C
  88.  3          CONTINUE
  89.             DO 60 IB=1,N1PTEL
  90.                DO 61 ID=1,N1EL
  91.                   XX=MELVA1.VELCHE(IB,ID)
  92.                   XX=XX*T0
  93.                   VELCHE(IB,ID)=(FBG3(XX))/X0
  94.  61            CONTINUE
  95.  60         CONTINUE
  96.             GOTO 166
  97. C
  98. C         DERIVEE CAPACITE
  99. C
  100.  4          CONTINUE
  101.             DO 70 IB=1,N1PTEL
  102.                DO 71 ID=1,N1EL
  103.                   XX=MELVA1.VELCHE(IB,ID)
  104.                   XX=XX*T0
  105.                   VELCHE(IB,ID)=(FBG4(XX))*T0/X0
  106.  71            CONTINUE
  107.  70         CONTINUE
  108.             GOTO 166
  109. C
  110. C         ENDOMMAGEMENT CHIMIQUE
  111. C
  112.  5          CONTINUE
  113.             DO 80 IB=1,N1PTEL
  114.                DO 81 ID=1,N1EL
  115.                   XX=MELVA1.VELCHE(IB,ID)
  116.                   XX=XX*T0
  117.                   VELCHE(IB,ID)=FBG5(XX)
  118.  81            CONTINUE
  119.  80         CONTINUE
  120.             GOTO 166
  121. C
  122. C         DIFFUSIVITE MECANIQUE
  123. C
  124.  6          CONTINUE
  125.             DO 90 IB=1,N1PTEL
  126.                DO 91 ID=1,N1EL
  127.                   XX=MELVA1.VELCHE(IB,ID)
  128.                   XX=XX
  129.                   VELCHE(IB,ID)=(FBG6(XX))/X0
  130.  91            CONTINUE
  131.  90         CONTINUE
  132.             GOTO 166
  133.  
  134.  166        CONTINUE
  135.             SEGDES MELVAL
  136.             SEGDES MELVA1
  137. *
  138.  75      CONTINUE
  139.          SEGDES MCHAML
  140. *
  141.  72   CONTINUE
  142.       SEGDES MCHELM
  143.  666  CONTINUE
  144.       RETURN
  145.       END
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  

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