Esope source files

   1 : C PRE4      SOURCE    CHAT      05/01/13    02:25:12     5004           
   2 :       SUBROUTINE PRE4()
   3 : C************************************************************************
   4 : C
   5 : C PROJET            :  CASTEM 2000
   6 : C
   7 : C NOM               :  PRE4
   8 : C
   9 : C DESCRIPTION       :  Transport de scalaires passifs (vitesse aux faces)
  10 : C     
  11 : C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  12 : C      
  13 : C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/TTMF
  14 : C     
  15 : C************************************************************************
  16 : C
  17 : C
  18 : C APPELES (Outils)  : LIRENT, ERREUR 
  19 : C
  20 : C APPELES (Calcul)  :
  21 : C
  22 : C PRE4-------------- PRE41 (1er ordre en espace)
  23 : C
  24 : C************************************************************************
  25 : C
  26 : C PHRASE D'APPEL (GIBIANE) :
  27 : C
  28 : C  a)  1er ordre en espace (1er ordre en temps)
  29 : C
  30 : C       CHAM1S  = 'PRET' 'CLAUDEIS' 'FACE' ENTI1 TAB1 CHPO1  ;
  31 : C
  32 : C
  33 : C ENTREES  : 
  34 : C
  35 : C            ENTI1       : ordre en espace (=1)
  36 : C
  37 : C            MOD1        : OBJET MODELE de TYPE NAVIER_STOKES
  38 : C
  39 : C            CHPO1       : CHPOINT "CENTRE" contenant les scalaires à tranporter
  40 : C
  41 : C SORTIES  : 
  42 : C
  43 : C            CHAM1S      : MCHAML "FACEL" contenant les s scalaires à tranporter
  44 : C                          (a gauche et a droite de chaque face).
  45 : C                          Autant de composantes que CHPO1
  46 : C
  47 : C
  48 : C************************************************************************
  49 : C
  50 : C HISTORIQUE (Anomalies et modifications éventuelles)
  51 : C
  52 : C HISTORIQUE :  Créée le 28.11.01
  53 : C 
  54 : C************************************************************************
  55 : C
  56 : C
  57 : C**** Variables de COOPTIO
  58 : C
  59 : C      INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  60 : C     &        ,IOTER,   IOLEC,  IOIMP,     IOCAR,  IOACQ
  61 : C     &        ,IOPER,   IOSGB,  IOGRA,     IOSAU,  IORES
  62 : C     &        ,IECHO,   IIMPI,  IOSPI
  63 : C     &        ,IDIM
  64 : C     &        ,MCOORD      
  65 : C     &        ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  66 : C     &        ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  67 : C     &        ,NORINC,NORVAL,NORIND,NORVAD
  68 : C     &        ,NUCROU, IPSAUV
  69 : C
  70 : C**** Les variables
  71 : C
  72 :       IMPLICIT INTEGER(I-N)
  73 :       INTEGER   ICOND,  IRETOU
  74 :      &         ,ORDESP,NBOPT, NBOPT0,IP
  75 :       PARAMETER(NBOPT=1)
  76 :       CHARACTER*(8) TOPTIO(NBOPT),MOTOP
  77 : C      
  78 :       DATA TOPTIO/'FACE'/
  79 : C
  80 : C**** Les Includes
  81 : C
  82 : -INC CCOPTIO
  83 : C
  84 : C**** Les mots clés
  85 : C
  86 :       ICOND = 1
  86 :       CALL LIRCHA(MOTOP,ICOND,IRETOU)
  87 :       IF(IERR .NE. 0) GOTO 9999
  88 :       NBOPT0 = NBOPT
  88 :       CALL OPTLI(IP,TOPTIO,MOTOP,NBOPT0)
  89 :       IF(IERR .NE. 0)GOTO 9999
  90 :       IF(IP.NE.1)THEN
  91 : *
  92 : ****     Tentative d'utilisation d'une option non-implementée
  93 : C
  93 :          CALL ERREUR(251)
  94 :          GOTO 9999
  95 :       ENDIF
  96 : C
  97 : C**** Initialisation des parametres d'erreur
  98 : C
  99 :       MOTERR(1:40) = '                                        '
 100 : C
 101 : C**** Lecture ORDESP
 102 : C
 103 :       ICOND = 1
 103 :       CALL LIRENT(ORDESP,ICOND,IRETOU)
 104 :       IF(IERR .NE. 0) GOTO 9999
 105 :       IF((ORDESP .LT. 1) .OR. (ORDESP .GT. 2))THEN
 106 : C
 107 : C******* Message d'erreur standard
 108 : C        -301 0 %m1:40
 109 : C
 110 :          MOTERR(1:40) = 'ENTI1 = ???                             '
 111 :          WRITE(IOIMP,*) MOTERR
 112 : C
 113 : C******* Message d'erreur standard
 114 : C        Entier valant: %i1
 115 : C        -2 0
 116 : C
 117 :          INTERR(1) = ORDESP
 117 :          CALL ERREUR(-2)        
 118 : C
 119 : C******* Message d'erreur standard
 120 : C        21 2
 121 : C        Données incompatibles
 122 : C
 122 :          CALL ERREUR(21)
 123 :          GOTO 9999
 124 :       ENDIF
 125 :       IF(ORDESP .EQ. 1)THEN
 126 : C
 127 : C******* ORDESP = 1
 128 : C
 128 :          CALL PRE41()
 129 :       ELSE
 130 : C
 131 : C******* Message d'erreur standard
 132 : C        5 3
 133 : C        Erreur anormale. Contactez votre support
 134 : C
 134 :          CALL ERREUR(5)
 135 :       ENDIF
 136 : C
 137 :  9999 CONTINUE
 138 : C      
 139 :       RETURN                                                                    
 140 :       END
 141 : 
 142 : 
 143 :  
 144 :  
 145 :  
 146 :  
 147 :  
 148 :