Sources Esope | Cast3M

pre2
C PRE2      SOURCE    CB215821  23/01/25    21:15:29     11573          
      SUBROUTINE PRE2()
C************************************************************************
C
C PROJET            :  CASTEM 2000
C
C NOM               :  PRE2
C
C DESCRIPTION       :  Voir PRET
C
C                      Resolution des Equations d'Euler, VF,
C                     "cell-centred approach"
C
C                      Calcul des etats "gauche et droit" aux
C                      interfaces pour le modele de gaz ideal
C                      multiespeces "calorically perfect".
C
C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
C
C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/TTMF
C
C************************************************************************
C
C
C APPELES (Outils)  : LIRENT, ERREUR
C
C APPELES (Calcul)  :
C
C PRE2-------------- PRE21 (1er ordre en espace)
C   |
C   |
C    --------------- PRE22 (2eme ordre en espace)
C
C************************************************************************
C
C PHRASE D'APPEL (GIBIANE) :
C
C 1) gaz ideal multi-espece
C
C    a)  1er ordre en espace (1er ordre en temps)
C
C       CHAM1S CHAM2S CHAM3S CHAM4S CHAM5S = 'PRET' 'PERFMULT'
C                     ENTI1  ENTI2 MOD1 CHPO1 CHPO2 CHPO3  CHPO4 CHPO5;
C
C
C ENTREES  :
C
C           'PERFMULT'   : mot clé (gaz ideal)
C
C            ENTI1       : ordre en espace (=1)
C
C            ENTI2       : ordre en temps (=1 si ENTI1 = 1)
C
C            MOD1        : OBJET MODELE de TYPE NAVIER_STOKES
C
C            CHPO1       : CHPOINT "CENTRE" contenant la masse volumique.
C
C            CHPO2       : CHPOINT "CENTRE" contenant la vitesse
C                         ( NDIM  composantes);
C
C            CHPO3       : CHPOINT "CENTRE" contenant la pression
C
C            CHPO4       : CHPOINT "CENTRE" contenant les fractiones
C                          massiques du gaz
C
C            CHPO5       : CHPOINT "CENTRE" contenant les "gamma" du gaz
C
C SORTIES  :
C
C            CHAM1S      : MCHAML "FACEL" contenant la masse
C                          volumique
C                          (a gauche et a droite de chaque face).
C                          Une seule composante ('SCAL')
C
C            CHAM2S      : MCHAML  "FACEL"  contenant la vitesse et
C                          les cosinus directeurs de (n,t) dans le
C                          repere (x,y);
C                          dans le cas 2D 6 composantes:
C                          'UN' = vitesse normale
C                          'UT' = vitesse tangentiel
C                          'NX' = n.x
C                          'NY' = n.y
C                          'TX' = t.x
C                          'TY' = t.y
C
C            CHAM3S      : MCHAML "FACEL" contenant la pression;
C                          Une seule composante ('SCAL')
C
C            CHAM4S      : MCHAML "FACEL" contenant les fractiones massiques
C                          du gaz
C
C            CHAM5S      : MCHAML "FACEL" contenant les "gamma" du gaz
C                          Une seule composante ('SCAL')
C
C
C    b) 2eme ordre en espace
C       1er  ordre  en temps ou 2eme ordre en temps
C
C       CHAM1S CHAM2S CHAM3S CHAM4S CHAM5S = 'PRET' 'PERFMULT' ENTI1 ENTI2
C                                       TAB1 TAB2
C                                       CHPO1  CHPO2  CHPO3
C                                       CHPO4  CHPO5  CHPO6
C                                       CHPO7  CHPO8  CHPO9
C                                       CHPO10 CHPO11 CHPO12
C                                       CHPO13 (FLOT1) ;
C
C
C ENTREES  :
C
C           'PERFMULT'   : mot clé (gas ideal)
C
C            ENTI1       : ordre en espace (=2)
C
C            ENTI2       : ordre en temps  (=1 ou 2)
C
C            TAB1        : TABLE de sous type "DOMAINE"
C
C            TAB2        : TABLE des proprietes du gaz
C
C            CHPO1       : CHPOINT "CENTRE" contenant la masse
C                          volumique;  une seule composante ('SCAL')
C
C            CHPO2       : CHPOINT "CENTRE" contenant le gradient de
C                          la  masse volumique (2 composantes en 2D,
C                          'P1DX', 'P1DY').
C
C            CHPO3       : CHPOINT "CENTRE"  contenant le limiteur du
C                          gradient de la masse volumique
C                          (une seule composante 'P1  ')
C
C            CHPO4       : CHPOINT "CENTRE" contenant la vitesse
C                          (2 composantes en 2D, 'UX' ey 'UY' );
C
C            CHPO5       : CHPOINT "CENTRE" contenant le gradient de la
C                          vitesse (4 composantes en 2D, 'P1DX','P1DY'
C                          'P2DX','P2DY').
C
C            CHPO6       : CHPOINT "CENTRE"  contenant le limiteur du
C                          gradient de la vitesse (2 composantes en 2D
C                          'P1', 'P2').
C
C            CHPO7       : CHPOINT "CENTRE" contenant la pression;
C                           une seule composante ('SCAL')
C
C            CHPO8       : CHPOINT "CENTRE" contenant le gradient de
C                          la  pression (2 composantes en 2D,
C                          'P1DX', 'P1DY').
C
C            CHPO9       : CHPOINT "CENTRE" contenant le limiteur du
C                          gradient de la pression
C                          (1 composantes, 'P1' ).
C
C            CHPO10      : CHPOINT "CENTRE" contenant les fractiones massiques;
C                          n composantes (n >= 1), compatibles avec
C                          TAB2 . 'ESPEULE'
C
C            CHPO11      : CHPOINT "CENTRE" contenant le gradient des fractiones
C                          massique (2n composantes en 2D)
C
C            CHPO12      : CHPOINT "CENTRE" contenant le limiteur du
C                          gradient des fractiones massiques
C                          (n composantes ).
C
C            CHPO13      : CHPOINT "CENTRE" contenant les "gamma" du gaz
C
C            FLOT1       : FLOTTANT qui contient l'increment du temps
C                          pour la prediction, pour avoir le deuxieme
C                          ordre en temps (ENTI2 = 2)
C                          FLOT1 = (increment du temps  / 2)
C
C SORTIES  : voir le cas a)
C
C************************************************************************
C
C HISTORIQUE (Anomalies et modifications éventuelles)
C
C HISTORIQUE :  Créée le 10.7.98.
C
C************************************************************************
C
C
C**** Variables de COOPTIO
C
C      INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
C     &        ,IOTER,   IOLEC,  IOIMP,     IOCAR,  IOACQ
C     &        ,IOPER,   IOSGB,  IOGRA,     IOSAU,  IORES
C     &        ,IECHO,   IIMPI,  IOSPI
C     &        ,IDIM
C     &        ,MCOORD
C     &        ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
C     &        ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
C     &        ,NORINC,NORVAL,NORIND,NORVAD
C     &        ,NUCROU, IPSAUV
C
C**** Les variables
C
      IMPLICIT INTEGER(I-N)
      INTEGER   ICOND,  IRETOU
     &         ,ORDESP, ORDTEM
C
C**** Les Includes
C
 
-INC SMCOORD
-INC PPARAM
-INC CCOPTIO
C
C**** Initialisation des parametres d'erreur
C
      MOTERR(1:40) = '                                        '
C
C**** Lecture ORDESP
C
      ICOND = 1
      CALL LIRENT(ORDESP,ICOND,IRETOU)
      IF(IERR .NE. 0) GOTO 9999
      IF((ORDESP .LT. 1) .OR. (ORDESP .GT. 2))THEN
C
C******* Message d'erreur standard
C        -301 0 %m1:40
C
         MOTERR(1:40) = 'ENTI1 = ???                             '
         CALL ERREUR(-301)
C
C******* Message d'erreur standard
C        Entier valant: %i1
C        -2 0
C
         INTERR(1) = ORDESP
         CALL ERREUR(-2)
C
C******* Message d'erreur standard
C        21 2
C        Données incompatibles
C
         CALL ERREUR(21)
         GOTO 9999
      ENDIF
C
C**** Lecture ORDTEM
C
      ICOND = 1
      CALL LIRENT(ORDTEM,ICOND,IRETOU)
      IF(IERR .NE. 0) GOTO 9999
      IF((ORDTEM .LT. 1) .OR. (ORDTEM .GT. ORDESP))THEN
C
C******* Message d'erreur standard
C        -301 0 %m1:40
C
         MOTERR(1:40) = 'ENTI2 = ???                             '
         CALL ERREUR(-301)
C
C******* Message d'erreur standard
C        Entier valant: %i1
C        -2 0
C
         INTERR(1) = ORDTEM
         CALL ERREUR(-2)
C
C******* Message d'erreur standard
C        21 2
C        Données incompatibles
C
         CALL ERREUR(21)
         GOTO 9999
      ENDIF
C
      SEGACT,MCOORD
      IF((ORDESP .EQ. 1) .AND. (ORDTEM .EQ. 1))THEN
C
C******* ORDESP = 1, ORDTEM = 1
C
         CALL PRE21
      ELSEIF((ORDESP .EQ. 2))THEN
C
C******* ORDESP = 2, ORDTEM = 1 ou 2
C
         CALL PRE22(ORDTEM)
      ELSE
C
C******* Message d'erreur standard
C        251 2
C        Tentative d'utilisation d'une option non implémentée
C
         CALL ERREUR(251)
      ENDIF
      SEGDES,MCOORD
C
 9999 CONTINUE
C
      RETURN
      END