Télécharger devers.eso

Retour à la liste

Numérotation des lignes :

  1. C DEVERS SOURCE BP208322 15/06/22 21:17:06 8543
  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/
  44. CALL LIROBJ('ELEMSTRU',IRET1,1,IRETOU)
  45. CALL LIROBJ('ELEMSTRU',IRET2,1,IRETOU)
  46. IF(IERR.NE.0) RETURN
  47. 5 CONTINUE
  48. CALL LIRMOT(NOMDEV,LOMDEV,IMOT,0)
  49. IF(IMOT.EQ.0) GOTO 60
  50. CALL LIRREE(RRET,1,IRETOU)
  51. IF(IERR.NE.0) RETURN
  52. GOTO(10,20,30,40,50,55,57,58,59), IMOT
  53. 10 G=RRET
  54. GOTO 5
  55. 20 RO=RRET
  56. GOTO 5
  57. 30 H=RRET
  58. GOTO 5
  59. 40 R=RRET
  60. GOTO 5
  61. 50 Q=RRET
  62. GOTO 5
  63. 55 EP1=RRET
  64. GOTO 5
  65. 57 EP2=RRET
  66. GOTO 5
  67. 58 CR=RRET
  68. GOTO 5
  69. 59 EPS=RRET
  70. GOTO 5
  71. 60 CONTINUE
  72. MELSTR=IRET1
  73. SEGACT MELSTR
  74. MSOSTU=ISOSTU(1)
  75. IPT1=IMELEM(1)
  76. SEGDES MELSTR
  77. MELSTR=IRET2
  78. SEGACT MELSTR
  79. IPT2=IMELEM(1)
  80. SEGACT IPT1,IPT2
  81. NJON=IPT1.NUM(/2)
  82. N=1
  83. SEGINI MATTAC
  84. M=0
  85. N=NJON
  86. SEGINI MSOUMA
  87. LISATT(1)=MSOUMA
  88. ITYATT='DEVE'
  89. IPHYCH=0
  90. NI=0
  91. NJ=LOMDEV
  92. N1=0
  93. NP=0
  94. NT=0
  95. SEGINI MGEOCH
  96. IGEOCH=MGEOCH
  97. RJEUCH(1)=G
  98. RJEUCH(2)=RO
  99. RJEUCH(3)=H
  100. RJEUCH(4)=R
  101. RJEUCH(5)=Q
  102. RJEUCH(6)=EP1
  103. RJEUCH(7)=EP2
  104. RJEUCH(8)=CR
  105. RJEUCH(9)=EPS
  106. SEGDES MGEOCH
  107. SEGACT MCOORD
  108. C
  109. C CREATION DU CHAMPOINT DES VARIABLES DE LIAISONS
  110. C
  111. NBNO=XCOOR(/1)/(IDIM+1)
  112. XCOOR(**)=0.D0
  113. XCOOR(**)=0.D0
  114. IF(IDIM.EQ.3) XCOOR(**)=0.D0
  115. XCOOR(**)=0.D0
  116. NSOUPO=1
  117. NAT=1
  118. SEGINI MCHPOI
  119. NC=LOVALI
  120. SEGINI MSOUPO
  121. IPCHP(1)=MSOUPO
  122. IFOPOI=2
  123. NBNN=1
  124. NBELEM=1
  125. NBSOUS=0
  126. NBREF=0
  127. SEGINI MELEME
  128. IGEOC=MELEME
  129. ITYPEL=1
  130. NUM(1,1)=NBNO+1
  131. N=1
  132. SEGINI MPOVAL
  133. IPOVAL=MPOVAL
  134. DO 65 NCC=1,NC
  135. NOCOMP(NCC)=MOVALI(NCC)
  136. NOHARM(NCC)=0
  137. VPOCHA(1,NCC)=0.D0
  138. 65 CONTINUE
  139. SEGDES MPOVAL,MELEME,MSOUPO,MCHPOI
  140. CALL ECCHPO(MCHPOI,0)
  141. MJOP=MCHPOI
  142. C
  143. N=2
  144. DO 80 NJ=1,NJON
  145. SEGINI MJONCT
  146. IATREL(NJ)=MJONCT
  147. MJOTYP='DEVE'
  148. MJOPOI=MJOP
  149. MJODDL=' '
  150. DO 70 IP=1,2
  151. ISTRJO(IP)=MSOSTU
  152. NSOUPO=1
  153. NAT=1
  154. SEGINI MCHPOI
  155. IPCHJO(IP)=MCHPOI
  156. IPOSJO(IP)=0
  157. NC=3
  158. SEGINI MSOUPO
  159. IPCHP(1)=MSOUPO
  160. IFOPOI=2
  161. NBNN=1
  162. NBELEM=1
  163. NBSOUS=0
  164. NBREF=0
  165. SEGINI MELEME
  166. IGEOC=MELEME
  167. ITYPEL=1
  168. IF(IP.EQ.1) NUM(1,1)=IPT1.NUM(1,NJ)
  169. IF(IP.EQ.2) NUM(1,1)=IPT2.NUM(1,NJ)
  170. DO 75 I=1,NC
  171. READ (NOMDD(I),FMT='(A4)') NOCOMP(I)
  172. NOHARM(I)=0
  173. 75 CONTINUE
  174. SEGDES MELEME,MSOUPO,MCHPOI
  175. 70 CONTINUE
  176. SEGDES MJONCT
  177. 80 CONTINUE
  178. SEGDES IPT1,IPT2,MSOUMA,MATTAC
  179. CALL ECROBJ('ATTACHE ',MATTAC)
  180. RETURN
  181. END
  182.  
  183.  
  184.  
  185.  
  186.  
  187.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales