Esope source files

   1 : C VLHJ7     SOURCE    CHAT      05/01/13    04:08:07     5004           
   2 :       SUBROUTINE VLHJ7(RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ,DFRUN)
   3 : C
   4 : C************************************************************************
   5 : C
   6 : C PROJET            :  CASTEM 2000
   7 : C
   8 : C NOM               :  VLHJ7
   9 : C
  10 : C DESCRIPTION       :  Calcul du Jacobien du flux numerique VLH
  11 : C                      Contribution des faces sur les murs
  12 : C                      Voir VLHJ6
  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           :  MPFUNT
  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 =  normale à la face
  35 : C
  36 : C
  37 : C**** Sorties:
  38 : C
  39 : C     DFRUN(5)= derivés partielles de F_RHOUN par raport a RHO, UX, 
  40 : C                UY, UZ, P
  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 :       IMPLICIT INTEGER(I-N)
  53 :       REAL*8 RHO,UX,UY,UZ,P,GAM,CNX,CNY,CNZ
  54 :      &     ,DFRUN(5)
  55 :      &     ,UN,PSRHO,CSON2,CSON
  56 :      &     ,MACH,MNUM,DERMNU,PNUM,DERPNU
  57 :      &     ,CELL
  58 : C
  59 :       UN = UX * CNX + UY * CNY + UZ * CNZ
  60 :       PSRHO = P / RHO
  61 :       CSON2 = GAM * PSRHO
  62 :       CSON = SQRT(CSON2)
  63 :       MACH = UN / CSON
  63 :       CALL MPFUNT(MACH,MNUM,DERMNU,PNUM,DERPNU)
  64 : C
  65 :       CELL = DERMNU * UN + MNUM * CSON
  66 :       DFRUN(1) = CELL * UN + DERPNU * MACH * PSRHO
  67 :       CELL = (CELL *  RHO + P / CSON * DERPNU) * 2.0D0
  68 :       DFRUN(2) = CELL * CNX
  69 :       DFRUN(3) = CELL * CNY
  70 :       DFRUN(4) = CELL * CNZ
  71 :       DFRUN(5) = (CSON * MNUM - UN * DERMNU) * UN / PSRHO +
  72 :      &     2.0D0 * PNUM - MACH * DERPNU
  73 : C
  74 :       RETURN
  75 :       END
  76 :       
  77 :  
  78 :  
  79 :