C FIMPVF SOURCE CHAT 05/01/13 00:00:49 5004 SUBROUTINE FIMPVF C************************************************************************ C C PROJET : CASTEM 2000 C C NOM : FIMPVF C C DESCRIPTION : Subroutine appellée par YFIMP C C Modelisation 2D/3D des equations d'Euler/NS C Termes source C C Calcul du flux/residu C C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec estensions CISI) C C AUTEUR : A. BECCANTINI, DEN/DM2S/SEMT/LTMF C C************************************************************************ C C C************************************************************************ C C HISTORIQUE (Anomalies et modifications éventuelles) C C HISTORIQUE : Décembre 2003. The option 'GRAVKEPS' is added C in order to account for the term source of C the Multi-Component NS equations with k-\eps model. C C : February 2004. Mode 2D AXIS: contribution of the pressure C to the residual C C************************************************************************ IMPLICIT INTEGER(I-N) -INC CCOPTIO C C Variables from CCOPTIO 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, IFICLE, IPREFI C & ,MCOORD C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE C & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU C & ,NORINC,NORVAL,NORIND,NORVAD C & ,NUCROU, IPSAUV C INTEGER NBOPT, ICELL C PARAMETER (NBOPT=3) CHARACTER*8 LOPT(NBOPT) C DATA LOPT/'GRAVMONO','GRAVKEPS','AXIS'/ C CALL LIRMOT(LOPT,NBOPT,ICELL,1) IF(IERR .NE. 0) GOTO 9999 IF(ICELL .EQ. 1)THEN C 'GRAV ' C Gravité CALL FIMVF1() IF(IERR .NE. 0) GOTO 9999 ELSEIF(ICELL .EQ. 2) THEN CALL FIMTR1() IF(IERR .NE. 0) GOTO 9999 ELSEIF(ICELL .EQ. 3) THEN CALL FIMAXI() IF(IERR .NE. 0) GOTO 9999 ENDIF 9999 CONTINUE RETURN END