Télécharger la plaquette Cast3M


fr en

Télécharger centj1.eso

Retour à la liste

Numérotation des lignes : 

centj1

  1. C CENTJ1    SOURCE    CHAT      05/01/12    21:52:02     5004
  2.       SUBROUTINE CENTJ1(RHO,UX,UY,P,GAM,CNX,CNY,CTX,CTY,
  3.      &                  DFRO,DFRUN,DFRUT,DFRET)
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  CENTJ1
  9. C
  10. C DESCRIPTION       :  Calcul du Jacobien du flux numerique centré
  11. C                      Contribution de l'onde de gauche à droite
  12. C                      Voir CENTJ0
  13. C
  14. C LANGAGE           :  FORTRAN 77
  15. C
  16. C AUTEUR            :  A. BECCANTINI DRN/DMT/SEMT/LTMF
  17. C
  18. C************************************************************************
  19. C
  20. C APPELES           :  MPFUNT
  21. C
  22. C************************************************************************
  23. C
  24. C**** Entrées:
  25. C
  26. C     RHO      =  densité
  27. C
  28. C     UX, UY   =  vitesse
  29. C
  30. C     P        =  pression
  31. C
  32. C     GAM      =  le "gamma"
  33. C
  34. C     CNX, CNY =  normale à la face
  35. C
  36. C     CTX, CTY =  tangente à la face
  37. C
  38. C
  39. C**** Sorties:
  40. C
  41. C     DFRO(4)  = derivés partielles de F_RHO+ par raport a RHO, UX,
  42. C                UY, P
  43. C
  44. C     DFRUN(4) = derivés partielles de F_RHOUN+ par raport a RHO, UX,
  45. C                UY, P
  46. C
  47. C     DFRUT(4) = derivés partielles de F_RHOUT+ par raport a RHO, UX,
  48. C                UY, P
  49. C
  50. C     DFRET(4) = derivés partielles de F_RHOET+ par raport a RHO,
  51. C                UX, UY, P
  52. C
  53. C************************************************************************
  54. C
  55. C HISTORIQUE (Anomalies et modifications éventuelles)
  56. C
  57. C HISTORIQUE :  Créé le 04.09.00
  58. C
  59. C************************************************************************
  60. C
  61. C N.B.: Toutes les variables sont declarées explicitement
  62. C
  63. C
  64.       IMPLICIT INTEGER(I-N)
  65.       REAL*8 RHO,UX,UY,P,GAM,CNX,CNY,CTX,CTY
  66.      &      ,DFRO(4),DFRUN(4),DFRUT(4),DFRET(4)
  67.      &      ,UN,UT,CELL,GSGM1P,ECIN
  68. C
  69.       UN = UX * CNX + UY * CNY
  70.       UT = UX * CTX + UY * CTY
  71.       CELL = RHO * UN
  72. C
  73.       DFRO(1) = 0.5D0 * UN
  74.       DFRO(2) = 0.5D0 * RHO * CNX
  75.       DFRO(3) = 0.5D0 * RHO * CNY
  76.       DFRO(4) = 0.0D0
  77. C
  78.       DFRUN(1) = DFRO(1) * UN
  79.       DFRUN(2) = CELL * CNX
  80.       DFRUN(3) = CELL * CNY
  81.       DFRUN(4) = 0.5D0
  82. C
  83.       DFRUT(1) = DFRO(1) * UT
  84.       DFRUT(2) = (DFRO(2) * UT)  + (0.5D0 * CELL * CTX)
  85.       DFRUT(3) = (DFRO(3) * UT)  + (0.5D0 * CELL * CTY)
  86.       DFRUT(4) = 0.0D0
  87. C
  88.       GSGM1P=(GAM/(GAM-1.0D0))*P
  89.       ECIN=0.5D0*((UX*UX)+(UY*UY))
  90.       DFRET(1) = 0.5D0*ECIN*UN
  91.       DFRET(2) = 0.5D0*((CNX*(GSGM1P+(RHO*ECIN)))+
  92.      &     (CELL*UX))
  93.       DFRET(3) = 0.5D0*((CNY*(GSGM1P+(RHO*ECIN)))+
  94.      &     (CELL*UY))
  95.       DFRET(4) = 0.5D0*UN*(GAM/(GAM-1.D0))
  96. C
  97.       RETURN
  98.       END
  99.  
  100.  
  101.  
  102.  
  103.  

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