Esope source files

   1 : C VLHJ1     SOURCE    CHAT      05/01/13    04:07:48     5004           
   2 :       SUBROUTINE VLHJ1(RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY,
   3 :      &                  DFRO,DFRUN,DFRUT,DFRET)
   4 : C************************************************************************
   5 : C
   6 : C PROJET            :  CASTEM 2000
   7 : C
   8 : C NOM               :  VLHJ1
   9 : C
  10 : C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  11 : C                      Contribution de l'onde de gauche à droite
  12 : C                      Voir VLHJ0
  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     RHOET    =  rho * (ein + ecin)
  33 : C
  34 : C     GAM      =  le "gamma"
  35 : C
  36 : C     CNX, CNY =  normale à la face
  37 : C
  38 : C     CTX, CTY =  tangente à la face
  39 : C
  40 : C
  41 : C**** Sorties:
  42 : C
  43 : C     DFRO(4)  = derivés partielles de F_RHO+ par raport a RHO, UX, 
  44 : C                UY, P
  45 : C
  46 : C     DFRUN(4) = derivés partielles de F_RHOUN+ par raport a RHO, UX, 
  47 : C                UY, P
  48 : C
  49 : C     DFRUT(4) = derivés partielles de F_RHOUT+ par raport a RHO, UX, 
  50 : C                UY, P
  51 : C
  52 : C     DFRET(4) = derivés partielles de F_RHOET+ par raport a RHO,
  53 : C                UX, UY, P
  54 : C
  55 : C************************************************************************
  56 : C
  57 : C HISTORIQUE (Anomalies et modifications éventuelles)
  58 : C
  59 : C HISTORIQUE :  Créé le 04.09.00
  60 : C
  61 : C************************************************************************
  62 : C
  63 : C N.B.: Toutes les variables sont declarées explicitement
  64 : C                          
  65 : C
  66 :       IMPLICIT INTEGER(I-N)
  67 :       REAL*8 RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY
  68 :      &      ,DFRO(4),DFRUN(4),DFRUT(4),DFRET(4)
  69 :      &      ,UN,UT,CSON2,CSON,HT
  70 :      &      ,MACH,MNUM,DERMNU,PNUM,DERPNU
  71 :      &      ,CELL, FRHO
  72 : C
  73 :       UN = UX * CNX + UY * CNY
  74 :       UT = UX * CTX + UY * CTY
  75 :       CSON2 = GAM * P / RHO
  76 :       CSON = SQRT(CSON2)
  77 :       MACH = UN / CSON
  78 :       HT = (RHOET + P) / RHO
  78 :       CALL MPFUNT(MACH,MNUM,DERMNU,PNUM,DERPNU)
  79 : C
  80 :       FRHO = MNUM * RHO * CSON
  81 : C
  82 :       DFRO(1) = ((DERMNU * UN) + (MNUM * CSON))  / 2.0D0
  83 :       DFRO(2) = DERMNU * RHO * CNX
  84 :       DFRO(3) = DERMNU * RHO * CNY
  85 :       DFRO(4) = (MNUM - DERMNU * MACH) * RHO * CSON / (2.0D0 * P)
  86 : C
  87 :       DFRUN(1) = (DFRO(1) * UN) + (DERPNU * MACH * P / (2.0D0 * RHO))
  88 :       CELL = (MNUM * RHO * CSON) + (DERPNU  * P / CSON)
  89 :       DFRUN(2) = (DFRO(2) * UN) + (CELL * CNX)
  90 :       DFRUN(3) = (DFRO(3) * UN) + (CELL * CNY)
  91 :       DFRUN(4) = (DFRO(4) * UN) - (0.5D0 * MACH * DERPNU) + PNUM
  92 : C     
  93 :       DFRUT(1) = DFRO(1) * UT
  94 :       DFRUT(2) = (DFRO(2) * UT)  + (FRHO * CTX)
  95 :       DFRUT(3) = (DFRO(3) * UT)  + (FRHO * CTY)
  96 :       DFRUT(4) = DFRO(4) * UT
  97 : C     
  98 :       CELL = CSON2 / (GAM - 1.0D0)
  99 :       DFRET(1) = (DFRO(1) * HT) - (FRHO * CELL / RHO)
 100 :       DFRET(2) = (DFRO(2) * HT)  + (FRHO * UX)
 101 :       DFRET(3) = (DFRO(3) * HT)  + (FRHO * UY)
 102 :       DFRET(4) = (DFRO(4) * HT) + (FRHO * CELL / P)
 103 : C
 104 :       RETURN
 105 :       END
 106 :       
 107 :  
 108 :