Esope source files

   1 : C VLHJ6     SOURCE    CHAT      05/01/13    04:08:04     5004           
   2 :       SUBROUTINE VLHJ6(RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ,
   3 :      &                  DFRUNC)
   4 : C************************************************************************
   5 : C
   6 : C PROJET            :  CASTEM 2000
   7 : C
   8 : C NOM               :  VLHJ6
   9 : C
  10 : C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  11 : C                      Contribution des faces sur les murs 
  12 : C                      Cas 3D
  13 : C
  14 : C LANGAGE           :  FORTRAN 77
  15 : C      
  16 : C AUTEUR            :  A. BECCANTINI SFME/LTMF  
  17 : C
  18 : C************************************************************************
  19 : C
  20 : C APPELES           :  VLHJ3
  21 : C
  22 : C************************************************************************
  23 : C
  24 : C**** Entrées:
  25 : C
  26 : C     RHO           =  densité
  27 : C
  28 : C     UX, UY, UZ    =  vitesse 
  29 : C
  30 : C     P             =  pression
  31 : C
  32 : C     GAM           =  le "gamma"
  33 : C
  34 : C     CNX, CNY, CNZ =  normales à la face
  35 : C
  36 : C
  37 : C**** Sorties:
  38 : C
  39 : C     DFRUNC(4)= derivés partielles de F_RHOUN par raport a RHO, RHOUX, 
  40 : C                RHOUY, RHOUZ, RHOET
  41 : C
  42 : C************************************************************************
  43 : C
  44 : C HISTORIQUE (Anomalies et modifications éventuelles)
  45 : C
  46 : C HISTORIQUE :  Créé le 23.07.01
  47 : C
  48 : C************************************************************************
  49 : C
  50 : C N.B.: Toutes les variables sont declarées explicitement
  51 : C                          
  52 : C
  53 :       IMPLICIT INTEGER(I-N)
  54 :       REAL*8 RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ
  55 :      &      ,DFRUNC(5),DFRUNP(5)
  56 :      &      ,USRO,UXSRO,UYSRO,UZSRO,GM1,GM1EC,GM1UX,GM1UY,GM1UZ
  57 : C
  58 : C**** Jacobian with respect to the primitive variables
  59 : C    
  59 :       CALL VLHJ7(RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ,DFRUNP)
  60 : C
  61 : C**** Jacobian with respect to the conservative variables
  62 : C    
  63 :       USRO = 1.0D0 / RHO 
  64 :       UXSRO = UX * USRO
  65 :       UYSRO = UY * USRO
  66 :       UZSRO = UZ * USRO
  67 :       GM1 = GAM - 1.0D0
  68 :       GM1EC = GM1 * 0.5D0 * (UX * UX + UY * UY + UZ * UZ)
  69 :       GM1UX = GM1 * UX
  70 :       GM1UY = GM1 * UY
  71 :       GM1UZ = GM1 * UZ
  72 : C
  73 :       DFRUNC(1) = DFRUNP(1) - UXSRO * DFRUNP(2) -  UYSRO * DFRUNP(3)
  74 :      &           - UZSRO * DFRUNP(4) + GM1EC * DFRUNP(5)
  75 :       DFRUNC(2) =  USRO * DFRUNP(2) -  GM1UX * DFRUNP(5)
  76 :       DFRUNC(3) =  USRO * DFRUNP(3) -  GM1UY * DFRUNP(5)
  77 :       DFRUNC(4) =  USRO * DFRUNP(4) -  GM1UZ * DFRUNP(5)
  78 :       DFRUNC(5) =  GM1 * DFRUNP(5)
  79 : C
  80 :       RETURN
  81 :       END
  82 :       
  83 :  
  84 :  
  85 :