Esope source files

   1 : C CENTJ3    SOURCE    CHAT      05/01/12    21:52:09     5004           
   2 :       SUBROUTINE CENTJ3(RHO,UX,UY,CNX,CNY,DFRUN)
   3 : C
   4 : C************************************************************************
   5 : C
   6 : C PROJET            :  CASTEM 2000
   7 : C
   8 : C NOM               :  CENTJ3
   9 : C
  10 : C DESCRIPTION       :  Calcul du Jacobien du flux numerique centré
  11 : C                      Contribution des faces sur les murs
  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           :  MPFUNT
  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     CNX, CNY =  normale à la face
  30 : C
  31 : C
  32 : C**** Sorties:
  33 : C
  34 : C     DFRUN(4)= derivés partielles de F_RHOUN par raport a RHO, UX, 
  35 : C                UY, P
  36 : C
  37 : C************************************************************************
  38 : C
  39 : C HISTORIQUE (Anomalies et modifications éventuelles)
  40 : C
  41 : C HISTORIQUE :  Créé le 04.09.00
  42 : C
  43 : C************************************************************************
  44 : C
  45 : C N.B.: Toutes les variables sont declarées explicitement
  46 : C                          
  47 :       IMPLICIT INTEGER(I-N)
  48 :       REAL*8 RHO,UX,UY,CNX,CNY
  49 :      &     ,DFRUN(4)
  50 :      &     ,UN, CELL
  51 : C
  52 :       UN = UX * CNX + UY * CNY
  53 :       CELL = 2.0D0 * RHO * UN
  54 :       DFRUN(1) = UN * UN
  55 :       DFRUN(2) = CELL * CNX
  56 :       DFRUN(3) = CELL * CNY
  57 :       DFRUN(4) = 1.0D0
  58 : C
  59 :       RETURN
  60 :       END
  61 :       
  62 :  
  63 :  
  64 :