Esope source files

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