Esope source files

   1 : C PRE6      SOURCE    BECC      09/12/03    21:15:08     6574           
   2 :       SUBROUTINE PRE6()
   3 : C************************************************************************
   4 : C
   5 : C PROJET            :  CASTEM 2000
   6 : C
   7 : C NOM               :  PRE6
   8 : C
   9 : C DESCRIPTION       :  Voir PRET 
  10 : C
  11 : C                      Resolution des Equations d'Euler, VF,
  12 : C                     "cell-centred approach". Discrete Equation
  13 : C                      Method for combustion.
  14 : C      
  15 : C                      Calcul des etats "gauche et droit" aux
  16 : C                      interfaces.
  17 : C     
  18 : C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  19 : C      
  20 : C AUTEUR            :  A. BECCANTINI, DEN/DM2S/SFME/LTMF
  21 : C     
  22 : C************************************************************************
  23 : C
  24 : C PHRASE D'APPEL (GIBIANE) :
  25 : C
  26 : C 1) Discrete Equation Method
  27 : C
  28 : C    a)  2nd order in space
  29 : C
  30 : C       CHAMAL1 CHAMAL2 CHAMRN1 CHAMRN2 CHAMV1 CHAMV2 CHAMP1 CHAMP2 = 
  31 : C          'PRET' 'DEM' MOD1   TABPGAS 
  32 : C                               CHPAL1 CHPGRA1 CHPLIA1  
  33 : C                               CHPAL2 CHPGRA1 CHPLIA1
  34 : C                               CHPRN1 CHPGRR1 CHPLIR1
  35 : C                               CHPRN2 CHPGRR2 CHPLIR2
  36 : C                               CHPVN1 CHPGRV1 CHPLIV1
  37 : C                               CHPVN2 CHPGRV2 CHPLIV2
  38 : C                               CHPP1  CHPGRP1 CHPLIP1
  39 : C                               CHPP2  CHPGRP2 CHPLIP2
  40 : C                               EPS ;
  41 : C
  42 : C ENTREES  : 
  43 : C
  44 : C           'PERFTEMP'   : mot clé (gaz "thermally perfect")
  45 : C                    
  46 : C            MOD1        : OBJET MODELE de TYPE NAVIER_STOKES
  47 : C
  48 : C            TABPGAS     : TABLE des proprietes du gaz
  49 : C
  50 : C            CHPAL1 CHPGRA1 CHPLIA1  : CHPOINTs "CENTRE":
  51 : C                          alpha1, its gradient, its limiter
  52 : C
  53 : C            CHPAL2 CHPGRA2 CHPLIA2  : CHPOINTs "CENTRE":
  54 : C                          alpha2, its gradient, its limiter
  55 : C
  56 : C            CHPRN1 CHPGRR1 CHPLIR1  : CHPOINTs "CENTRE":
  57 : C                          rho1, its gradient, its limiter
  58 : C
  59 : C            CHPRN2 CHPGRR2 CHPLIR2  : CHPOINTs "CENTRE":
  60 : C                          rho2, its gradient, its limiter
  61 : C
  62 : C            CHPVN1 CHPGRV1 CHPLIV1  : CHPOINTs "CENTRE":
  63 : C                          v1, its gradient, its limiter
  64 : C
  65 : C            CHPVN2 CHPGRV2 CHPLIV2  : CHPOINTs "CENTRE":
  66 : C                          v2, its gradient, its limiter
  67 : C
  68 : C            CHPVP1 CHPGRP1 CHPLIP1  : CHPOINTs "CENTRE":
  69 : C                          v1, its gradient, its limiter
  70 : C
  71 : C            CHPVP2 CHPGRP2 CHPLIP2  : CHPOINTs "CENTRE":
  72 : C                          v2, its gradient, its limiter
  73 : C
  74 : C            EPS         : FLOTTANT; alpha_i < eps  => no phase_i 
  75 : C     
  76 : C SORTIES  : 
  77 : C
  78 : C            CHAMAL1 CHAMAL2 CHAMRN1 CHAMRN2 CHAMV1 CHAMV2 CHAMP1
  79 : C            CHAMP2      :  MCHAML "FACEL" containing alpha, density
  80 : C                           velocity and pressure on the left and on the
  81 : C                           right
  82 : C
  83 : C************************************************************************
  84 : C
  85 : C HISTORIQUE (Anomalies et modifications éventuelles)
  86 : C
  87 : C HISTORIQUE :  Crée le 03.12.09.
  88 : C 
  89 : C************************************************************************
  90 : C
  91 :       IMPLICIT INTEGER(I-N)
  92 : -INC CCOPTIO
  93 : C
  94 : C******* ORDESP = 2, ORDTEM = 1
  95 : C
  95 :       CALL PRE61()
  96 :       IF(IERR .NE. 0) GOTO 9999
  97 : C
  98 :  9999 CONTINUE
  99 : C
 100 :       RETURN
 101 :       END
 102 :