1 : C DEVERS SOURCE AM 08/02/14 21:20:19 6050 2 : SUBROUTINE DEVERS 3 : IMPLICIT INTEGER(I-N) 4 : IMPLICIT REAL*8(A-H,O-Z) 5 : C ==================================================================== 6 : C = L'OPERATEUR DEVE FABRIQUE UN OBJET ATTACHE CORRESPONDANT A UNE = 7 : C = LIAISON DE TYPE DEVERSEMENT DE FLUIDE = 8 : C = SYNTAXE : = 9 : C = ATTA1=DEVE ELS1 ELS2 GRAV G RHO RHO1 HAUT H RAYO R ZINI Z0 ; = 10 : C = ATTA1 OBJET ATTACHE = 11 : C = ELS1 OBJET ELEMENT-STRUCTURE RASSEMBLANT LES POINTS DU HAUT = 12 : C = DU DEVERSOIR = 13 : C = ELS2 OBJET ELEMENT-STRUCTURE RASSEMBLANT LES POINTS DU BAS = 14 : C = DU DEVERSOIR = 15 : C = ( LES POINTS DU HAUT ET DU BAS DOIVENT SE CORRESPONDRE 2 A 2 ) = 16 : C = MOTS CLEFS SUIVIS D'UN FLOTTANT : = 17 : C = GRAV ACCELERATION DE LA PESANTEUR = 18 : C = RHO MASSE VOLUMIQUE DU FLUIDE = 19 : C = HAUT HAUTEUR DU DEVERSOIR = 20 : C = RAYO RAYON DU DEVERSOIR = 21 : C = DEBC DEBIT CRITIQUE EN M***/S = 22 : C = EPLA EPAISSEUR DE LA LAME FLUIDE EN HAUT = 23 : C = EPLB EPAISSEUR DE LA LAME FLUIDE EN BAS = 24 : C = CCRI COEFFICIENT DU DEBIT CRITIQUE 25 : C = ASPE HAUTEUR DES ASPERITES = 26 : C = CREATION 24/05/85 = 27 : C = PROGRAMMEUR GUILBAUD = 28 : C ==================================================================== 29 : C 30 : -INC CCOPTIO 31 : -INC CCHAMP 32 : -INC SMCOORD 33 : -INC SMATTAC 34 : -INC SMELEME 35 : -INC SMELSTR 36 : -INC SMCHPOI 37 : CHARACTER*4 NOMDEV(9) 38 : CHARACTER*4 MOVALI(6) 39 : DATA NOMDEV/'GRAV','RHO ','HAUT','RAYO','DEBC','EPLA','EPLB', 40 : *'CCRI','ASPE'/ 41 : DATA LOMDEV/9/ 42 : DATA MOVALI/'ZMHA','ZMBA','DEBI','SOHA','SOBA','FOBA'/ 43 : DATA LOVALI/6/ 43 : CALL LIROBJ('ELEMSTRU',IRET1,1,IRETOU) 43 : CALL LIROBJ('ELEMSTRU',IRET2,1,IRETOU) 44 : IF(IERR.NE.0) RETURN 45 : 5 CONTINUE 45 : CALL LIRMOT(NOMDEV,LOMDEV,IMOT,0) 46 : IF(IMOT.EQ.0) GOTO 60 46 : CALL LIRREE(RRET,1,IRETOU) 47 : IF(IERR.NE.0) RETURN 48 : GOTO(10,20,30,40,50,55,57,58,59), IMOT 49 : 10 G=RRET 50 : GOTO 5 51 : 20 RO=RRET 52 : GOTO 5 53 : 30 H=RRET 54 : GOTO 5 55 : 40 R=RRET 56 : GOTO 5 57 : 50 Q=RRET 58 : GOTO 5 59 : 55 EP1=RRET 60 : GOTO 5 61 : 57 EP2=RRET 62 : GOTO 5 63 : 58 CR=RRET 64 : GOTO 5 65 : 59 EPS=RRET 66 : GOTO 5 67 : 60 CONTINUE 68 : MELSTR=IRET1 69 : SEGACT MELSTR 70 : MSOSTU=ISOSTU(1) 71 : IPT1=IMELEM(1) 72 : SEGDES MELSTR 73 : MELSTR=IRET2 74 : SEGACT MELSTR 75 : IPT2=IMELEM(1) 76 : SEGACT IPT1,IPT2 77 : NJON=IPT1.NUM(/2) 78 : N=1 79 : SEGINI MATTAC 80 : M=0 81 : N=NJON 82 : SEGINI MSOUMA 83 : LISATT(1)=MSOUMA 84 : ITYATT='DEVE' 85 : IPHYCH=0 86 : NI=0 87 : NJ=LOMDEV 88 : N1=0 89 : NP=0 90 : NT=0 91 : SEGINI MGEOCH 92 : IGEOCH=MGEOCH 93 : RJEUCH(1)=G 94 : RJEUCH(2)=RO 95 : RJEUCH(3)=H 96 : RJEUCH(4)=R 97 : RJEUCH(5)=Q 98 : RJEUCH(6)=EP1 99 : RJEUCH(7)=EP2 100 : RJEUCH(8)=CR 101 : RJEUCH(9)=EPS 102 : SEGDES MGEOCH 103 : SEGACT MCOORD 104 : C 105 : C CREATION DU CHAMPOINT DES VARIABLES DE LIAISONS 106 : C 107 : NBNO=XCOOR(/1)/(IDIM+1) 108 : XCOOR(**)=0.D0 109 : XCOOR(**)=0.D0 110 : IF(IDIM.EQ.3) XCOOR(**)=0.D0 111 : XCOOR(**)=0.D0 112 : NSOUPO=1 113 : NAT=1 114 : SEGINI MCHPOI 115 : NC=LOVALI 116 : SEGINI MSOUPO 117 : IPCHP(1)=MSOUPO 118 : IFOPOI=2 119 : NBNN=1 120 : NBELEM=1 121 : NBSOUS=0 122 : NBREF=0 123 : SEGINI MELEME 124 : IGEOC=MELEME 125 : ITYPEL=1 126 : NUM(1,1)=NBNO+1 127 : N=1 128 : SEGINI MPOVAL 129 : IPOVAL=MPOVAL 130 : DO 65 NCC=1,NC 131 : NOCOMP(NCC)=MOVALI(NCC) 132 : NOHARM(NCC)=0 133 : VPOCHA(1,NCC)=0.D0 134 : 65 CONTINUE 135 : SEGDES MPOVAL,MELEME,MSOUPO,MCHPOI 135 : CALL ECCHPO(MCHPOI,0) 136 : MJOP=MCHPOI 137 : C 138 : N=2 139 : DO 80 NJ=1,NJON 140 : SEGINI MJONCT 141 : IATREL(NJ)=MJONCT 142 : MJOTYP='DEVE' 143 : MJOPOI=MJOP 144 : MJODDL=' ' 145 : DO 70 IP=1,2 146 : ISTRJO(IP)=MSOSTU 147 : NSOUPO=1 148 : NAT=1 149 : SEGINI MCHPOI 150 : IPCHJO(IP)=MCHPOI 151 : IPOSJO(IP)=0 152 : NC=3 153 : SEGINI MSOUPO 154 : IPCHP(1)=MSOUPO 155 : IFOPOI=2 156 : NBNN=1 157 : NBELEM=1 158 : NBSOUS=0 159 : NBREF=0 160 : SEGINI MELEME 161 : IGEOC=MELEME 162 : ITYPEL=1 163 : IF(IP.EQ.1) NUM(1,1)=IPT1.NUM(1,NJ) 164 : IF(IP.EQ.2) NUM(1,1)=IPT2.NUM(1,NJ) 165 : DO 75 I=1,NC 166 : READ (NOMDD(I),FMT='(A4)') NOCOMP(I) 167 : NOHARM(I)=0 168 : 75 CONTINUE 169 : SEGDES MELEME,MSOUPO,MCHPOI 170 : 70 CONTINUE 171 : SEGDES MJONCT 172 : 80 CONTINUE 173 : SEGDES IPT1,IPT2,MSOUMA,MATTAC 173 : CALL ECROBJ('ATTACHE ',MATTAC) 174 : RETURN 175 : END 176 : 177 : 178 : 179 :
© Cast3M 2003 - All rights reserved.
Disclaimer