C INITFL SOURCE BP208322 12/01/10 21:15:01 7244 SUBROUTINE INITFL(IPRIGI,IPMASS,IPMV1,IPV1,IFLU) * * ********************************************************************** * * SBR REALISANT L INITAILISATION DES VECTEURS DE PROJECTION * POUR LA RECHERCHE DE MODES PAR LA METHODE SIMULTANEE * * UNE ATTENTION PARTICULIERE POUR LE CAS OU L ON A DU FLUIDE * GENERATON D UN VECTEUR RESPECTANT LAPL.PI = 0 * * * D. BROCHARD (EMSI) 6994 13/07/88 * BP 10/01/2012 : IFLU = nombre d'inconnues PI * * ********************************************************************** * * IMPLICIT INTEGER(I-N) IMPLICIT REAL*8 (A-H,O-Z) -INC CCHAMP -INC CCOPTIO -INC SMCHPOI SEGMENT IDEMEN(0) CHARACTER*4 MOC * IFLU=0 * * write(6,*) ' initfl nomdu ' , nomdu(15) MCHPOI=IPMV1 SEGACT MCHPOI NSOUPO = IPCHP(/1) DO 10 ISOUPO=1,NSOUPO MSOUPO = IPCHP(ISOUPO) SEGACT MSOUPO MPOVAL=IPOVAL segact,MPOVAL NC = NOCOMP(/2) DO 20 IC = 1,NC MOC=NOCOMP(IC) c IF(MOC.EQ.NOMDU(15)) GOTO 110 IF(MOC.EQ.NOMDU(15)) IFLU=IFLU+VPOCHA(/1) 20 CONTINUE segdes,MPOVAL SEGDES MSOUPO 10 CONTINUE SEGDES MCHPOI c GO TO 100 if(IFLU.eq.0) GO TO 100 110 continue C C IL Y A DES TERMES EN PI C * write(6,*) ' on passe dans la mise à zero des termes en PI' IPMV1L=IPMV1 CALL ANCHPO(IPMV1L,NOMDU(15),IPMV1) CALL DTCHPO(IPMV1L) CALL DTCHPO(IPV1) SEGINI IDEMEN IDEMEN(**)=IPMV1 NOID=0 NOEN=0 SEGDES IDEMEN IDEM = IDEMEN CALL RESOU1(IPRIGI,IDEM,NOID,NOEN,1D-18) * recuperer nombre de modes d'ensemble et maillage pts contrainnts CALL LIRENT(NBENT,1,IRETOU) call lirobj('MAILLAGE',ipt8,1,iretou) IDEMEN = IDEM SEGACT IDEMEN IPV1=IDEMEN(1) SEGSUP IDEMEN CALL MUCPRI(IPV1,IPMASS,IPMV1L) CALL ANCHPO(IPMV1L,NOMDU(15),IPMV1) CALL DTCHPO(IPMV1L) c IFLU = 1 * 100 CONTINUE * * RETURN END