Esope source files

   1 : C VLHJ4     SOURCE    CHAT      05/01/13    04:07:57     5004           
   2 :       SUBROUTINE VLHJ4(RHO,UX,UY,UZ,P,RHOET,GAM,CNX,CNY,CNZ,
   3 :      &     CT1X,CT1Y,CT1Z,CT2X,CT2Y,CT2Z,
   4 :      &     DFROC,DFRUNC,DFRU1C,DFRU2C,DFRETC)
   5 : C************************************************************************
   6 : C
   7 : C PROJET            :  CASTEM 2000
   8 : C
   9 : C NOM               :  VLHJ4
  10 : C
  11 : C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  12 : C                      Contribution de l'onde de gauche à droite
  13 : C                      Cas 3D monoespece
  14 : C
  15 : C LANGAGE           :  FORTRAN 77
  16 : C      
  17 : C AUTEUR            :  A. BECCANTINI SFME/LTMF  
  18 : C
  19 : C************************************************************************
  20 : C
  21 : C APPELES           :  VLHJ5
  22 : C
  23 : C************************************************************************
  24 : C
  25 : C**** Entrées:
  26 : C
  27 : C     RHO             =  densité
  28 : C
  29 : C     UX, UY, UZ      =  vitesse 
  30 : C
  31 : C     P               =  pression
  32 : C
  33 : C     RHOET           =  rho * (ein + ecin)
  34 : C
  35 : C     GAM             =  le "gamma"
  36 : C
  37 : C     CNX,  CNY,  CNZ =  normale à la face
  38 : C
  39 : C     CT1X, CT1Y, CT1Z
  40 : C     CT2X, CT2Y, CT2Z=  tangentes à la face (T1,T2)
  41 : C
  42 : C**** Sorties:
  43 : C
  44 : C     DFROC(5) = derivés partielles de F_RHO+ par raport a RHO, RHOUX, 
  45 : C                RHOUY, RHOUZ, RHOET
  46 : C
  47 : C     DFRUNC(5)= derivés partielles de F_RHOUN+ par raport a RHO, RHOUX, 
  48 : C                RHOUY, RHOUZ, RHOET
  49 : C
  50 : C     DFRU1C(5)= derivés partielles de F_RHOUT1+ par raport a RHO, RHOUX, 
  51 : C                 RHOUY, RHOUZ, RHOET (tangente T1)
  52 : C
  53 : C     DFRU2C(5)= derivés partielles de F_RHOUT2+ par raport a RHO, RHOUX, 
  54 : C                 RHOUY, RHOUZ, RHOET (tangente T2)
  55 : C
  56 : C     DFRETC(5)= derivés partielles de F_RHOET+ par raport a RHO,
  57 : C                RHOUX, RHOUY, RHOUZ, RHOET
  58 : C
  59 : C************************************************************************
  60 : C
  61 : C HISTORIQUE (Anomalies et modifications éventuelles)
  62 : C
  63 : C HISTORIQUE :  Créé le 23.07.01
  64 : C
  65 : C************************************************************************
  66 : C
  67 : C N.B.: Toutes les variables sont declarées explicitement
  68 : C                          
  69 :       IMPLICIT INTEGER(I-N)
  70 :       REAL*8 RHO,UX,UY,UZ,P,RHOET,GAM,CNX,CNY,CNZ,CT1X,CT1Y,CT1Z
  71 :      &     ,CT2X,CT2Y,CT2Z
  72 :      &     ,DFROC(5),DFRUNC(5),DFRU1C(5),DFRU2C(5),DFRETC(5)
  73 :      &     ,DFROP(5),DFRUNP(5),DFRU1P(5),DFRU2P(5),DFRETP(5)
  74 :      &     ,USRO,UXSRO,UYSRO,UZSRO,GM1,GM1EC,GM1UX,GM1UY,GM1UZ
  75 : C
  76 : C**** Jacobian with respect to the primitive variables
  77 : C    
  77 :       CALL VLHJ5(RHO,UX,UY,UZ,P,RHOET,GAM,CNX,CNY,CNZ,
  78 :      &     CT1X,CT1Y,CT1Z,CT2X,CT2Y,CT2Z,
  79 :      &     DFROP,DFRUNP,DFRU1P,DFRU2P,DFRETP)
  80 : C
  81 : C**** Jacobian with respect to the conservative variables
  82 : C    
  83 :       USRO = 1.0D0 / RHO
  84 :       UXSRO = UX * USRO
  85 :       UYSRO = UY * USRO
  86 :       UZSRO = UZ * USRO
  87 :       GM1 = GAM - 1.0D0
  88 :       GM1EC = GM1 * 0.5D0 * (UX * UX + UY * UY + UZ * UZ)
  89 :       GM1UX = GM1 * UX
  90 :       GM1UY = GM1 * UY
  91 :       GM1UZ = GM1 * UZ
  92 : C
  93 :       DFROC(1) = DFROP(1) - UXSRO * DFROP(2) -  UYSRO * DFROP(3)
  94 :      &     - UZSRO  * DFROP(4) + GM1EC * DFROP(5)
  95 :       DFROC(2) =  USRO * DFROP(2) -  GM1UX * DFROP(5)
  96 :       DFROC(3) =  USRO * DFROP(3) -  GM1UY * DFROP(5)
  97 :       DFROC(4) =  USRO * DFROP(4) -  GM1UZ * DFROP(5)
  98 :       DFROC(5) =  GM1 * DFROP(5)
  99 : C
 100 :       DFRUNC(1) = DFRUNP(1) - UXSRO * DFRUNP(2) -  UYSRO * DFRUNP(3)
 101 :      &     - UZSRO  * DFRUNP(4) + GM1EC * DFRUNP(5)
 102 :       DFRUNC(2) =  USRO * DFRUNP(2) -  GM1UX * DFRUNP(5)
 103 :       DFRUNC(3) =  USRO * DFRUNP(3) -  GM1UY * DFRUNP(5)
 104 :       DFRUNC(4) =  USRO * DFRUNP(4) -  GM1UZ * DFRUNP(5)
 105 :       DFRUNC(5) =  GM1 * DFRUNP(5)
 106 : C
 107 :       DFRU1C(1) = DFRU1P(1) - UXSRO * DFRU1P(2) -  UYSRO * DFRU1P(3)
 108 :      &     - UZSRO  * DFRU1P(4) + GM1EC * DFRU1P(5)
 109 :       DFRU1C(2) =  USRO * DFRU1P(2) -  GM1UX * DFRU1P(5)
 110 :       DFRU1C(3) =  USRO * DFRU1P(3) -  GM1UY * DFRU1P(5)
 111 :       DFRU1C(4) =  USRO * DFRU1P(4) -  GM1UZ * DFRU1P(5)
 112 :       DFRU1C(5) =  GM1 * DFRU1P(5)
 113 : C
 114 :       DFRU2C(1) = DFRU2P(1) - UXSRO * DFRU2P(2) -  UYSRO * DFRU2P(3)
 115 :      &     - UZSRO  * DFRU2P(4) + GM1EC * DFRU2P(5)
 116 :       DFRU2C(2) =  USRO * DFRU2P(2) -  GM1UX * DFRU2P(5)
 117 :       DFRU2C(3) =  USRO * DFRU2P(3) -  GM1UY * DFRU2P(5)
 118 :       DFRU2C(4) =  USRO * DFRU2P(4) -  GM1UZ * DFRU2P(5)
 119 :       DFRU2C(5) =  GM1 * DFRU2P(5)
 120 : C
 121 :       DFRETC(1) = DFRETP(1) - UXSRO * DFRETP(2) -  UYSRO * DFRETP(3)
 122 :      &     - UZSRO  * DFRETP(4) + GM1EC * DFRETP(5)
 123 :       DFRETC(2) =  USRO * DFRETP(2) -  GM1UX * DFRETP(5)
 124 :       DFRETC(3) =  USRO * DFRETP(3) -  GM1UY * DFRETP(5)
 125 :       DFRETC(4) =  USRO * DFRETP(4) -  GM1UZ * DFRETP(5)
 126 :       DFRETC(5) =  GM1 * DFRETP(5)
 127 : C
 128 :       RETURN
 129 :       END
 130 :       
 131 :  
 132 :  
 133 :