Esope source files

   1 : C FIMPVF    SOURCE    CHAT      05/01/13    00:00:49     5004           
   2 :       SUBROUTINE FIMPVF
   3 : C************************************************************************
   4 : C
   5 : C PROJET            :  CASTEM 2000
   6 : C
   7 : C NOM               :  FIMPVF
   8 : C
   9 : C DESCRIPTION       :  Subroutine appellée par YFIMP
  10 : C
  11 : C                      Modelisation 2D/3D des equations d'Euler/NS
  12 : C                      Termes source
  13 : C
  14 : C                      Calcul du flux/residu
  15 : C
  16 : C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  17 : C
  18 : C AUTEUR            :  A. BECCANTINI, DEN/DM2S/SEMT/LTMF
  19 : C
  20 : C************************************************************************
  21 : C
  22 : C
  23 : C************************************************************************
  24 : C
  25 : C HISTORIQUE (Anomalies et modifications éventuelles)
  26 : C
  27 : C HISTORIQUE : Décembre 2003. The option 'GRAVKEPS' is added
  28 : C              in order to account for the term source of
  29 : C              the Multi-Component NS equations with k-\eps model.
  30 : C
  31 : C            : February 2004. Mode 2D AXIS: contribution of the pressure
  32 : C              to the residual
  33 : C
  34 : C************************************************************************
  35 :       IMPLICIT INTEGER(I-N)
  36 : -INC CCOPTIO
  37 : C
  38 : C     Variables from CCOPTIO     
  39 : C
  40 : C      INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  41 : C     &        ,IOTER,   IOLEC,  IOIMP,     IOCAR,  IOACQ
  42 : C     &        ,IOPER,   IOSGB,  IOGRA,     IOSAU,  IORES
  43 : C     &        ,IECHO,   IIMPI,  IOSPI
  44 : C     &        ,IDIM,    IFICLE, IPREFI     
  45 : C     &        ,MCOORD
  46 : C     &        ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  47 : C     &        ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  48 : C     &        ,NORINC,NORVAL,NORIND,NORVAD
  49 : C     &        ,NUCROU, IPSAUV
  50 : C
  51 :       INTEGER NBOPT, ICELL
  52 : C
  53 :       PARAMETER (NBOPT=3)
  54 :       CHARACTER*8 LOPT(NBOPT)
  55 : C      
  56 :       DATA LOPT/'GRAVMONO','GRAVKEPS','AXIS'/
  57 : C
  57 :       CALL LIRMOT(LOPT,NBOPT,ICELL,1)
  58 :       IF(IERR .NE. 0) GOTO 9999
  59 :       IF(ICELL .EQ. 1)THEN
  60 : C        'GRAV    '
  61 : C        Gravité
  61 :          CALL FIMVF1()
  62 :          IF(IERR .NE. 0) GOTO 9999
  63 :       ELSEIF(ICELL .EQ. 2) THEN
  63 :          CALL FIMTR1()
  64 :          IF(IERR .NE. 0) GOTO 9999
  65 :       ELSEIF(ICELL .EQ. 3) THEN
  65 :          CALL FIMAXI()
  66 :          IF(IERR .NE. 0) GOTO 9999
  67 :       ENDIF
  68 :  9999 CONTINUE
  69 :       RETURN
  70 :       END
  71 :  
  72 :  
  73 :  
  74 :