Esope source files

   1 : C CENTJ0    SOURCE    CHAT      05/01/12    21:51:59     5004           
   2 :       SUBROUTINE CENTJ0(RHO,UX,UY,P,GAM,CNX,CNY,CTX,CTY,
   3 :      &                  DFROC,DFRUNC,DFRUTC,DFRETC)
   4 : C************************************************************************
   5 : C
   6 : C PROJET            :  CASTEM 2000
   7 : C
   8 : C NOM               :  CENTJ0
   9 : C
  10 : C DESCRIPTION       :  Calcul du Jacobien du flux numerique centré
  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           :  CENTJ1
  20 : C                      Appellée par KONJR1, KONJC1
  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     DFROC(4) = derivés partielles de F_RHO+ par raport a RHO, RHOUX, 
  42 : C                RHOUY, RHOET
  43 : C
  44 : C     DFRUNC(4)= derivés partielles de F_RHOUN+ par raport a RHO, RHOUX, 
  45 : C                RHOUY, RHOET
  46 : C
  47 : C     DFRUNT(4)= derivés partielles de F_RHOUT+ par raport a RHO, RHOUX, 
  48 : C                RHOUY, RHOET
  49 : C
  50 : C     DFRETC(4)= derivés partielles de F_RHOET+ par raport a RHO,
  51 : C                RHOUX, RHOUY, RHOET
  52 : C
  53 : C************************************************************************
  54 : C
  55 : C HISTORIQUE (Anomalies et modifications éventuelles)
  56 : C
  57 : C HISTORIQUE :
  58 : C
  59 : C************************************************************************
  60 : C
  61 : C N.B.: Toutes les variables sont declarées explicitement
  62 : C                          
  63 :       IMPLICIT INTEGER(I-N)
  64 :       REAL*8 RHO,UX,UY,P,GAM,CNX,CNY,CTX,CTY
  65 :      &      ,DFROC(4),DFRUNC(4),DFRUTC(4),DFRETC(4)
  66 :      &      ,DFROP(4),DFRUNP(4),DFRUTP(4),DFRETP(4)
  67 :      &      ,USRO,UXSRO,UYSRO,GM1,GM1EC,GM1UX,GM1UY
  68 : C
  69 : C**** Jacobian with respect to the primitive variables
  70 : C    
  70 :       CALL CENTJ1(RHO,UX,UY,P,GAM,CNX,CNY,CTX,CTY,
  71 :      &                  DFROP,DFRUNP,DFRUTP,DFRETP)
  72 : C
  73 : C**** Jacobian with respect to the conservative variables
  74 : C    
  75 :       USRO = 1.0D0 / RHO
  76 :       UXSRO = UX * USRO
  77 :       UYSRO = UY * USRO
  78 :       GM1 = GAM - 1.0D0
  79 :       GM1EC = GM1 * 0.5D0 * (UX * UX + UY * UY)
  80 :       GM1UX = GM1 * UX
  81 :       GM1UY = GM1 * UY
  82 : C
  83 :       DFROC(1) = DFROP(1) - UXSRO * DFROP(2) -  UYSRO * DFROP(3)
  84 :      &           + GM1EC * DFROP(4)
  85 :       DFROC(2) =  USRO * DFROP(2) -  GM1UX * DFROP(4)
  86 :       DFROC(3) =  USRO * DFROP(3) -  GM1UY * DFROP(4)
  87 :       DFROC(4) =  GM1 * DFROP(4)
  88 : C
  89 :       DFRUNC(1) = DFRUNP(1) - UXSRO * DFRUNP(2) -  UYSRO * DFRUNP(3)
  90 :      &           + GM1EC * DFRUNP(4)
  91 :       DFRUNC(2) =  USRO * DFRUNP(2) -  GM1UX * DFRUNP(4)
  92 :       DFRUNC(3) =  USRO * DFRUNP(3) -  GM1UY * DFRUNP(4)
  93 :       DFRUNC(4) =  GM1 * DFRUNP(4)
  94 : C
  95 :       DFRUTC(1) = DFRUTP(1) - UXSRO * DFRUTP(2) -  UYSRO * DFRUTP(3)
  96 :      &           + GM1EC * DFRUTP(4)
  97 :       DFRUTC(2) =  USRO * DFRUTP(2) -  GM1UX * DFRUTP(4)
  98 :       DFRUTC(3) =  USRO * DFRUTP(3) -  GM1UY * DFRUTP(4)
  99 :       DFRUTC(4) =  GM1 * DFRUTP(4)
 100 : C
 101 :       DFRETC(1) = DFRETP(1) - UXSRO * DFRETP(2) -  UYSRO * DFRETP(3)
 102 :      &           + GM1EC * DFRETP(4)
 103 :       DFRETC(2) =  USRO * DFRETP(2) -  GM1UX * DFRETP(4)
 104 :       DFRETC(3) =  USRO * DFRETP(3) -  GM1UY * DFRETP(4)
 105 :       DFRETC(4) =  GM1 * DFRETP(4)
 106 : C
 107 :       RETURN
 108 :       END
 109 :       
 110 :  
 111 :  
 112 :