Télécharger centj3.eso

Retour à la liste

Numérotation des lignes :

  1. C CENTJ3 SOURCE CHAT 05/01/12 21:52:09 5004
  2. SUBROUTINE CENTJ3(RHO,UX,UY,CNX,CNY,DFRUN)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : CENTJ3
  9. C
  10. C DESCRIPTION : Calcul du Jacobien du flux numerique centré
  11. C Contribution des faces sur les murs
  12. C
  13. C LANGAGE : FORTRAN 77
  14. C
  15. C AUTEUR : A. BECCANTINI DRN/DMT/SEMT/LTMF
  16. C
  17. C************************************************************************
  18. C
  19. C APPELES : MPFUNT
  20. C
  21. C************************************************************************
  22. C
  23. C**** Entrées:
  24. C
  25. C RHO = densité
  26. C
  27. C UX, UY = vitesse
  28. C
  29. C CNX, CNY = normale à la face
  30. C
  31. C
  32. C**** Sorties:
  33. C
  34. C DFRUN(4)= derivés partielles de F_RHOUN par raport a RHO, UX,
  35. C UY, P
  36. C
  37. C************************************************************************
  38. C
  39. C HISTORIQUE (Anomalies et modifications éventuelles)
  40. C
  41. C HISTORIQUE : Créé le 04.09.00
  42. C
  43. C************************************************************************
  44. C
  45. C N.B.: Toutes les variables sont declarées explicitement
  46. C
  47. IMPLICIT INTEGER(I-N)
  48. REAL*8 RHO,UX,UY,CNX,CNY
  49. & ,DFRUN(4)
  50. & ,UN, CELL
  51. C
  52. UN = UX * CNX + UY * CNY
  53. CELL = 2.0D0 * RHO * UN
  54. DFRUN(1) = UN * UN
  55. DFRUN(2) = CELL * CNX
  56. DFRUN(3) = CELL * CNY
  57. DFRUN(4) = 1.0D0
  58. C
  59. RETURN
  60. END
  61.  
  62.  
  63.  
  64.  
  65.  

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