1 : C DEBIT SOURCE CHAT 06/08/24 21:25:56 5529 2 : SUBROUTINE DEBIT 3 : *_______________________________________________________________________ 4 : * 5 : * OPERATEUR DEBIT ( Inspire de PRESSI ) 6 : * 7 : * FQ = DEBIT MODEL I D MAIL I (MOT1) 8 : * I CHPDEB I 9 : * 10 : * FQ CHPOINT CONTENANT LES FORCES NODALES EQUIVALENTES 11 : * MODEL OBJET MMODEL ,MASSIF , SUR LEQUEL S'APPLIQUE 12 : * LA CONDITION DE DEBIT 13 : * D VALEUR ALGEBRIQUE DU DEBIT 14 : * MAIL POUR LES MASSIFS ,OBJET MAILLAGE REPRESENTANT 15 : * LA FACE SUR LAQUELLE S'APPLIQUE LA CONDITION 16 : * DE DEBIT 17 : * CHPOT CHPOINT CONTENANT LES VALEURS ALGEBRIQUES 18 : * DES DEBITS AUX NOEUDS 19 : * MOT1 NOM DONNE AU RESULAT ( PAR DEFAUT , CELUI DU MODELE) 20 : *______________________________________________________________________ 21 : * 22 : IMPLICIT INTEGER(I-N) 23 : -INC CCOPTIO 24 : -INC SMCHPOI 25 : * 26 : CHARACTER*(4) MOSCAL,MOPOI1,MOT1 27 : REAL*8 XXX 28 : * 29 : DATA MOSCAL/'SCAL'/,MOPOI1/'POI1'/,MOT1/' '/ 30 : * 31 : * ON LIT UN MMODEL 32 : * 32 : CALL LIROBJ('MMODEL ',IPMODL,1,IRETM) 33 : IF (IERR.NE.0) RETURN 34 : 35 : ICONV=0 36 : * 37 : * ON LIT SOIT UN FLOTTANT ET UN MAILLAGE ,SOIT UN CHPOINT 38 : * 38 : CALL LIRREE(XXX,0,IRETOU) 39 : IF (IRETOU.EQ.1) THEN 39 : CALL LIROBJ('MAILLAGE',IPTGEO,1,IRETOU) 40 : IF (IERR.NE.0) GOTO 666 41 : * 41 : CALL ECRREE (XXX) 41 : CALL PROG 41 : CALL ECRCHA(MOSCAL) 41 : CALL ECRENT(1) 41 : CALL ECROBJ('MAILLAGE',IPTGEO) 41 : CALL ECRCHA(MOPOI1) 41 : CALL PRCHAN 42 : IF (IERR.NE.0) GOTO 666 42 : CALL MANUCH 42 : CALL LIROBJ('CHPOINT ',ICHPOP,1,IRETOU) 43 : IF (IERR.NE.0) THEN 43 : CALL DTCHPO(ICHPOP) 44 : GOTO 666 45 : ENDIF 46 : ELSE 46 : CALL LIROBJ('CHPOINT ',ICHPOP,1,IRETOU) 47 : IF (IERR.NE.0) GOTO 666 48 : ENDIF 49 : * 49 : CALL LIRCHA(MOT1,0,IRETOU) 50 : * 50 : CALL FDEBIT(ICHPOP,IPMODL,IPTFP,MOT1,IRET) 51 : IF (IRET.EQ.0) GOTO 666 52 : * 53 : * LE NUMERO DE L HARMONIQUE EST PRIS EGAL A NIFOUR 54 : * POUR TOUTES LES COMPOSANTES DU CHPOINT 55 : * 56 : MCHPOI=IPTFP 57 : SEGACT MCHPOI 58 : * ATTRIBUTION D'UNE NATURE DISCRETE AU CHPO SORTANT 59 : NAT = MAX ( 1 , JATTRI(/1) ) 60 : NSOUPO = IPCHP(/1) 61 : SEGADJ MCHPOI 62 : JATTRI(1) = 2 63 : * 64 : DO 100 IA=1,IPCHP(/1) 65 : MSOUPO=IPCHP(IA) 66 : SEGACT MSOUPO 67 : DO 101 NC=1,NOHARM(/1) 68 : NOHARM(NC)=NIFOUR 69 : 101 CONTINUE 70 : SEGDES MSOUPO 71 : 100 CONTINUE 72 : SEGDES MCHPOI 73 : * 73 : CALL ECROBJ('CHPOINT ',IPTFP) 74 : * 75 : 666 IF (ICONV.EQ.1) THEN 75 : CALL DTMODL(IPMODL) 76 : ENDIF 77 : RETURN 78 : END 79 : 80 : 81 :
© Cast3M 2003 - All rights reserved.
Disclaimer