Esope source files

   1 : C FSTRAR    SOURCE    BECC      11/05/26    21:15:30     6981           
   2 :       SUBROUTINE FSTRAR(
   3 :      &     PC, GAM,
   4 :      &     P, RHO_R, P_R, U_R,
   5 :      &     U, RHO, D,
   6 :      &     LOGDEB, LOGAN)
   7 : *
   8 : *************************************************************************
   9 : *
  10 : * project           :  CAST3M, EUROPLEXUS...
  11 : *
  12 : * name              :  fstrar
  13 : *
  14 : * description       :  euler equations for a mixture of stiffened gases
  15 : *                      flux in the non-reactive case.
  16 : *
  17 : *                      right rarefaction curve as function of p
  18 : *
  19 : * language          :  fortran 77
  20 : *
  21 : * author            :  a. beccantini den/dm2s/sfme/ltmf
  22 : *
  23 : *************************************************************************
  24 : *
  24 : * called by         :
  25 : *
  26 : *
  27 : *************************************************************************
  28 : *
  29 : ***** input
  30 : *
  31 : *     pc, gam         = properties of the gas
  32 : *
  33 : *     p               = pressure in the left side of the right
  34 : *                       rarefaction
  35 : *
  36 : *     rho_r, p_r, u_r
  37 : *                     = density, pressure, velocity on the right
  38 : *
  39 : *     logdeb          = debugging ?
  40 : *
  41 : ***** output
  42 : *
  43 : *     u, rho          = velocity and density in the left side of the
  44 : *                       right rarefaction
  45 : *
  46 : *     d               = rarefaction speed
  47 : *
  48 : *     logan           = anomaly ?
  49 : *
  50 : *************************************************************************
  51 : *
  52 : *     26/11/2009      created
  53 : *     25/05/2011      evolution in CAST3M
  54 : *
  55 : *************************************************************************
  56 : *
  57 : * n.b.: all variables are declared
  58 : *
  59 : C      IMPLICIT NONE
  60 :       IMPLICIT INTEGER(I-N)
  61 :       REAL*8
  62 :      &     PC, GAM,
  63 :      &     P, RHO_R, P_R, U_R,
  64 :      &     U, D
  65 : *
  66 :       REAL*8 A_R, RHO, A, COEF
  67 : *
  68 : *     debugging ?
  69 : *
  70 :       LOGICAL LOGDEB, LOGAN
  71 : *
  71 :       CALL ASTIFF (RHO_R, P_R, GAM, PC, A_R, LOGDEB, LOGAN)
  72 : *      write(*,*) 'ciao', logan, a_r
  73 :       IF (LOGAN) GOTO 9999
  74 : *
  75 : *     density (isentropic transformation)
  76 : *
  77 : *     a negative p here is catched by astiff after...
  78 : *
  79 :       RHO = (P + PC) / (P_R + PC)
  80 :       RHO = RHO ** (1.0D0 / GAM)
  81 :       RHO = RHO * RHO_R
  82 : *
  82 :       CALL ASTIFF (RHO, P, GAM, PC, A, LOGDEB, LOGAN)
  83 :       IF (LOGAN) GOTO 9999
  84 : *
  85 :       COEF =  2.0D0 / (GAM - 1.0D0)
  86 :       U = U_R + (COEF * (A - A_R))
  87 : *
  88 :       D = U + A
  89 : *
  90 :  9999 RETURN
  91 :       END
  92 : 
  93 :