Sources Esope | Cast3M

pre3
C PRE3      SOURCE    OF166741  24/12/13    21:17:06     12097          
      SUBROUTINE PRE3()
C************************************************************************
C
C PROJET            :  CASTEM 2000
C
C NOM               :  PRE3
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
C                      mono/multiespeces "thermally 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 PRE3-------------- PRE31 (1er ordre en espace)
C   |
C   |
C    --------------- PRE32 (2eme ordre en espace)
C
C************************************************************************
C
C PHRASE D'APPEL (GIBIANE) :
C
C 1) gaz "thermally perfect"
C
C    a)  1er ordre en espace (1er ordre en temps)
C
C       CHAM1S CHAM2S CHAM3S (CHAM4S) = 'PRET' 'PERFTEMP'
C                     ENTI1  ENTI2 MOD1 TAB2
C                     CHPO1 CHPO2 CHPO3 (CHPO4);
C
C
C ENTREES  :
C
C           'PERFTEMP'   : mot clé (gaz "thermally perfect")
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            TAB2        : TABLE des proprietes du gaz
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 (dans le cas multi-especes,
C                          compatible avec TAB2)
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 (si existe CHPO4)
C
C
C    b) 2eme ordre en espace
C       1er  ordre  en temps
C
C       CHAM1S CHAM2S CHAM3S (CHAM4S) = 'PRET' 'PERFMULT' ENTI1 ENTI2
C                                        TAB1 TAB2
C                                        CHPO1  CHPO2  CHPO3
C                                        CHPO4  CHPO5  CHPO6
C                                        CHPO7  CHPO8  CHPO9
C                                        (CHPO10 CHPO11 CHPO12) ;
C
C
C ENTREES  :
C
C           'PERFTEMP'   : mot clé (gaz "thermally prefect")
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            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    c) 2eme ordre en espace
C       2er  ordre  en temps
C
C       CHAM1S CHAM2S CHAM3S (CHAM4S) = '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 DT ;
C
C
C ENTREES  :
C
C           'PERFTEMP'   : mot clé (gaz "thermally prefect")
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 le "gamma" du gaz
C                          une seule composante ('SCAL')
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 18.12.98.
C
C************************************************************************
C
C**** Les variables
C
      IMPLICIT INTEGER(I-N)
      INTEGER   ICOND,  IRETOU
     &         ,ORDESP, ORDTEM
C
C**** Les Includes
C
 
-INC PPARAM
-INC CCOPTIO
-INC SMCOORD
 
      segact mcoord
 
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
      IF((ORDESP .EQ. 1) .AND. (ORDTEM .EQ. 1))THEN
C
C******* ORDESP = 1, ORDTEM = 1
C
         CALL PRE31()
      ELSEIF((ORDESP .EQ. 2))THEN
C
C******* ORDESP = 2, ORDTEM = 1 ou 2
C
         CALL PRE32(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
C
 9999 CONTINUE
C
      RETURN
      END