Télécharger capap1.eso

Retour à la liste

Numérotation des lignes :

capap1
  1. C CAPAP1 SOURCE CB215821 21/03/03 21:15:07 10910
  2. SUBROUTINE CAPAP1(IPMAIL,IVAMAT,IPMATR,INFOR)
  3.  
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6.  
  7. C=======================================================================
  8. C= C A P A P 1 =
  9. C= ----------- =
  10. C= =
  11. C= Fonction : =
  12. C= ---------- =
  13. C= Calcul de la matrice de CAPACITE CALORIFIQUE d'un element POI1 =
  14. C= =
  15. C= Clement BERTHINIER, le 23 fevrier 2021. =
  16. C=======================================================================
  17.  
  18. -INC PPARAM
  19. -INC CCOPTIO
  20. -INC SMELEME
  21. -INC SMRIGID
  22. -INC SMCHAML
  23.  
  24. SEGMENT MPTVAL
  25. INTEGER IPOS(NS),NSOF(NS),IVAL(NCOSOU)
  26. CHARACTER*16 TYVAL(NCOSOU)
  27. ENDSEGMENT
  28.  
  29. MELEME = IPMAIL
  30. NBELEM = NUM(/2)
  31.  
  32. XMATRI = IPMATR
  33. MPTVAL = IVAMAT
  34.  
  35.  
  36. C BOUCLE SUR LES ELEMENTS DU MAILLAGE ELEMENTAIRE IPMAIL
  37. C ======================================================
  38. DO iElt = 1, NBELEM
  39. IF (INFOR .EQ. 1)THEN
  40. C Cas THERMIQUE : Calcul du terme M * Cp
  41. MELVAL = MPTVAL.IVAL(1)
  42. IEMN = MIN(iElt,MELVAL.VELCHE(/2))
  43. XM = MELVAL.VELCHE(1,IEMN)
  44.  
  45. MELVAL = MPTVAL.IVAL(2)
  46. IEMN = MIN(iElt,MELVAL.VELCHE(/2))
  47. XC = MELVAL.VELCHE(1,IEMN)
  48.  
  49. XRIG = XM * XC
  50.  
  51. ELSEIF(INFOR .EQ. 2)THEN
  52. C Cas DIFFUSION : Directement le CDG (CD généralise)
  53. MELVAL = MPTVAL.IVAL(1)
  54. IEMN = MIN(iElt,MELVAL.VELCHE(/2))
  55. XRIG = MELVAL.VELCHE(1,IEMN)
  56. ENDIF
  57.  
  58. C Remplissage direct de la matrice
  59. RE(1,1,iElt) = XRIG
  60. ENDDO
  61.  
  62. END
  63.  
  64.  

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