Sources Esope | Cast3M

pre1
C PRE1      SOURCE    OF166741  24/12/13    21:17:00     12097          
      SUBROUTINE PRE1()
C************************************************************************
C
C PROJET            :  CASTEM 2000
C
C NOM               :  PRE1
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                      mono-espece "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 PRE1-------------- PRE11 (1er ordre en espace)
C   |
C   |
C    --------------- PRE12 (2eme ordre en espace)
C
C************************************************************************
C
C PHRASE D'APPEL (GIBIANE) :
C
C 1) gaz ideal mono-espece
C
C    a)  1er ordre en espace (1er ordre en temps)
C
C       CHAM1S CHAM2S CHAM3S CHAM4S  = 'PRET' 'PERFMONO' ENTI1 ENTI2 MOD1
C                                       CHPO1  CHPO2  CHPO3  CHPO4 ;
C
C
C ENTREES  :
C
C           'PERFMONO'   : 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 "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 "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  = 'PRET' 'PERFMONO' ENTI1 ENTI2 TAB1
C                                       CHPO1  CHPO2  CHPO3
C                                       CHPO4  CHPO5  CHPO6
C                                       CHPO7  CHPO8  CHPO9
C                                       CHPO10 (FLOT1);
C
C
C ENTREES  :
C
C           'PERFMONO'   : 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            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 de
C                          gradient de la pression
C                          (1 composantes, 'P1' ).
C
C            CHPO10      : CHPOINT "CENTRE" contenant le "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
C SORTIES  : voir le cas a)
C
C************************************************************************
C
C HISTORIQUE (Anomalies et modifications éventuelles)
C
C HISTORIQUE :  Créée le 11.6.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**** Les mots clés
C
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 PRE11()
      ELSEIF((ORDESP .EQ. 2))THEN
C
C******* ORDESP = 2, ORDTEM = 1 ou 2
C
         CALL PRE12(ORDTEM)
      ELSE
C
C******* Message d'erreur standard
C        5 3
C        Erreur anormale. Contactez votre support
C
         CALL ERREUR(5)
      ENDIF
C
 9999 CONTINUE
C
      RETURN
      END