Esope source files

   1 : C VLHJ0     SOURCE    CHAT      05/01/13    04:07:44     5004           
   2 :       SUBROUTINE VLHJ0(RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY,
   3 :      &                  DFROC,DFRUNC,DFRUTC,DFRETC)
   4 : C************************************************************************
   5 : C
   6 : C PROJET            :  CASTEM 2000
   7 : C
   8 : C NOM               :  VLHJ0
   9 : C
  10 : C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  11 : C                      Contribution de l'onde de gauche à droite
  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           :  VLHJ1
  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     P        =  pression
  30 : C
  31 : C     RHOET    =  rho * (ein + ecin)
  32 : C
  33 : C     GAM      =  le "gamma"
  34 : C
  35 : C     CNX, CNY =  normale à la face
  36 : C
  37 : C     CTX, CTY =  tangente à la face
  38 : C
  39 : C
  40 : C**** Sorties:
  41 : C
  42 : C     DFROC(4) = derivés partielles de F_RHO+ par raport a RHO, RHOUX, 
  43 : C                RHOUY, RHOET
  44 : C
  45 : C     DFRUNC(4)= derivés partielles de F_RHOUN+ par raport a RHO, RHOUX, 
  46 : C                RHOUY, RHOET
  47 : C
  48 : C     DFRUNT(4)= derivés partielles de F_RHOUT+ par raport a RHO, RHOUX, 
  49 : C                RHOUY, RHOET
  50 : C
  51 : C     DFRETC(4)= derivés partielles de F_RHOET+ par raport a RHO,
  52 : C                RHOUX, RHOUY, RHOET
  53 : C
  54 : C************************************************************************
  55 : C
  56 : C HISTORIQUE (Anomalies et modifications éventuelles)
  57 : C
  58 : C HISTORIQUE :  Créé le 04.09.00
  59 : C
  60 : C************************************************************************
  61 : C
  62 : C N.B.: Toutes les variables sont declarées explicitement
  63 : C                          
  64 :       IMPLICIT INTEGER(I-N)
  65 :       REAL*8 RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY
  66 :      &      ,DFROC(4),DFRUNC(4),DFRUTC(4),DFRETC(4)
  67 :      &      ,DFROP(4),DFRUNP(4),DFRUTP(4),DFRETP(4)
  68 :      &      ,USRO,UXSRO,UYSRO,GM1,GM1EC,GM1UX,GM1UY
  69 : C
  70 : C**** Jacobian with respect to the primitive variables
  71 : C    
  71 :       CALL VLHJ1(RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY,
  72 :      &                  DFROP,DFRUNP,DFRUTP,DFRETP)
  73 : C
  74 : C**** Jacobian with respect to the conservative variables
  75 : C    
  76 :       USRO = 1.0D0 / RHO
  77 :       UXSRO = UX * USRO
  78 :       UYSRO = UY * USRO
  79 :       GM1 = GAM - 1.0D0
  80 :       GM1EC = GM1 * 0.5D0 * (UX * UX + UY * UY)
  81 :       GM1UX = GM1 * UX
  82 :       GM1UY = GM1 * UY
  83 : C
  84 :       DFROC(1) = DFROP(1) - UXSRO * DFROP(2) -  UYSRO * DFROP(3)
  85 :      &           + GM1EC * DFROP(4)
  86 :       DFROC(2) =  USRO * DFROP(2) -  GM1UX * DFROP(4)
  87 :       DFROC(3) =  USRO * DFROP(3) -  GM1UY * DFROP(4)
  88 :       DFROC(4) =  GM1 * DFROP(4)
  89 : C
  90 :       DFRUNC(1) = DFRUNP(1) - UXSRO * DFRUNP(2) -  UYSRO * DFRUNP(3)
  91 :      &           + GM1EC * DFRUNP(4)
  92 :       DFRUNC(2) =  USRO * DFRUNP(2) -  GM1UX * DFRUNP(4)
  93 :       DFRUNC(3) =  USRO * DFRUNP(3) -  GM1UY * DFRUNP(4)
  94 :       DFRUNC(4) =  GM1 * DFRUNP(4)
  95 : C
  96 :       DFRUTC(1) = DFRUTP(1) - UXSRO * DFRUTP(2) -  UYSRO * DFRUTP(3)
  97 :      &           + GM1EC * DFRUTP(4)
  98 :       DFRUTC(2) =  USRO * DFRUTP(2) -  GM1UX * DFRUTP(4)
  99 :       DFRUTC(3) =  USRO * DFRUTP(3) -  GM1UY * DFRUTP(4)
 100 :       DFRUTC(4) =  GM1 * DFRUTP(4)
 101 : C
 102 :       DFRETC(1) = DFRETP(1) - UXSRO * DFRETP(2) -  UYSRO * DFRETP(3)
 103 :      &           + GM1EC * DFRETP(4)
 104 :       DFRETC(2) =  USRO * DFRETP(2) -  GM1UX * DFRETP(4)
 105 :       DFRETC(3) =  USRO * DFRETP(3) -  GM1UY * DFRETP(4)
 106 :       DFRETC(4) =  GM1 * DFRETP(4)
 107 : C
 108 :       RETURN
 109 :       END
 110 :       
 111 :  
 112 :