Télécharger pre112.eso

Retour à la liste

Numérotation des lignes :

  1. C PRE112 SOURCE PV 09/03/12 21:30:33 6325
  2. C PRE111 SOURCE BECC 98/08/24 21:15:02 3286
  3. SUBROUTINE PRE112(ICEN,IFACE,IFACEL,INORM,
  4. & IROC,IVITC,IPC,IGAMC,
  5. & IROF,IVITF,IPF,IGAMF,
  6. & LOGAN,LOGNEG,LOGBOR,MESERR,VALER,VAL1,VAL2)
  7. C************************************************************************
  8. C
  9. C PROJET : CASTEM 2000
  10. C
  11. C NOM : PRE111
  12. C
  13. C DESCRIPTION : Voir PRE11
  14. C
  15. C Cas Deux Dimensions
  16. C
  17. C Une Espece
  18. C
  19. C 1er ordre en espace, 1re ordre en temps
  20. C
  21. C Creations des objets MCHAML IROF, IVITF, IPF,
  22. C IGAMF
  23. C
  24. C
  25. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  26. C
  27. C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF
  28. C
  29. C************************************************************************
  30. C
  31. C
  32. C APPELES (Outils) : KRIPAD, LICHT
  33. C
  34. C APPELES (Calcul) : AUCUN
  35. C
  36. C
  37. C************************************************************************
  38. C
  39. C ENTREES
  40. C
  41. C 1) Pointeurs de MELEMEs et de CHPOINTs de la table DOMAINE
  42. C
  43. C ICEN : MELEME de 'POI1' SPG des CENTRES
  44. C
  45. C IFACE : MELEME de 'POI1' SPG des FACES
  46. C
  47. C IFACEL : MELEME de 'SEG3' avec
  48. C CENTRE d'Elt "gauche"
  49. C CENTRE de Face
  50. C CENTRE d'Elt "droite"
  51. C
  52. C N.B. = IFACE.NUM(i,1) = IFACEL.NUM(i,2)
  53. C
  54. C INORM : CHPOINT des cosinus directeurs de normales aux faces
  55. C
  56. C 2) Pointeurs des CHPOINTs
  57. C
  58. C IROC : CHPOINT "CENTRE" contenant la masse volumique RHO
  59. C
  60. C IVITC : CHPOINT "CENTRE" contenant la vitesse UX, UY ;
  61. C
  62. C IPC : CHPOINT "CENTRE" contenat la pression P;
  63. C
  64. C IGAMC : CHPOINT "CENTRE" contenat le "Gamma" du gaz
  65. C
  66. C
  67. C SORTIES
  68. C
  69. C
  70. C IROF : MCHAML defini sur le MELEME de pointeur IFACEL,
  71. C contenant la masse volumique RHO
  72. C
  73. C IVITF : MCHAML defini sur le MELEME de pointeur IFACEL,
  74. C contenant la vitesse UN, UT dans le repaire local
  75. C (n,t) et defini sur le MELEME de pointeur IFACE,
  76. C contenant les cosinus directeurs du repere local
  77. C
  78. C IPF : MCHAML defini sur le MELEME de pointeur IFACEL,
  79. C contenant la pression P
  80. C
  81. C IGAMF : MCHAML defini sur le MELEME de pointeur IFACEL,
  82. C contenant le "gamma" du gaz
  83. C
  84. C LOGAN : anomalie detectee (changement de la convention dans
  85. C la table domaine)
  86. C
  87. C LOGNEG : (LOGICAL): si .TRUE. une pression ou une densité
  88. C negative a été detectée -> en interactif le
  89. C programme s'arrete en GIBIANE
  90. C (erreur stocké en MESERR et VALER)
  91. C
  92. C LOGBOR : (LOGICAL): si .TRUE. un gamma a ete detecte
  93. C dehor 1 et 3 (sa valeur stockée en MESERR et VALER;
  94. C en VAL1 et en VAL2 on stocke 1.0 et 3.0)
  95. C
  96. C MESERR
  97. C VALER
  98. C VAL1,
  99. C VAL2 : pour les messages d'erreur
  100. C
  101. C************************************************************************
  102. C
  103. C HISTORIQUE (Anomalies et modifications éventuelles)
  104. C
  105. C HISTORIQUE : Créée le 11.6.98.
  106. C
  107. C************************************************************************
  108. C
  109. C
  110. C ATTENTION: Cet programme marche que si le MAILLAGE est convex;
  111. C si non il faut changer l'argoritme de calcul de
  112. C l'orientation des normales aux faces.
  113. C
  114. C
  115. C************************************************************************
  116. C
  117. C**** Variables de COOPTIO
  118. C
  119. C INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  120. C & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ
  121. C & ,IOPER, IOSGB, IOGRA, IOSAU, IORES
  122. C & ,IECHO, IIMPI, IOSPI
  123. C & ,IDIM
  124. CC & ,MCOORD
  125. C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  126. C & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  127. C & ,NORINC,NORVAL,NORIND,NORVAD
  128. C & ,NUCROU, IPSAUV
  129. C
  130. C**** Les variables
  131. C
  132. IMPLICIT INTEGER(I-N)
  133. INTEGER ICEN, IFACE, IFACEL, IROC, IVITC, IPC ,IGAMC, INORM
  134. & , IROF, IVITF, IPF, IGAMF
  135. & , IGEOM, NFAC
  136. & , N1PTEL, N1EL, N2PTEL, N2EL, N2, N1, N3, L1
  137. & , NLCF, NGCF, NGCEG, NLCEG, NGCED, NLCED, NGCF1
  138. REAL*8 VALER, VAL1, VAL2,XG,YG,ZG,XC,YC,ZC
  139. & ,DXG, DYG, DZG,ORIENT
  140. & , CNX, CNY, CNZ, CTX, CTY, CTZ, CVX, CVY, CVZ
  141. & , ROG, PG, GAMG, UXG, UYG, UZG, UNG, UTG, UVG
  142. & , ROD, PD, GAMD, UXD, UYD, UZD, UND, UTD, UVD
  143. CHARACTER*(40) MESERR
  144. CHARACTER*(8) TYPE
  145. LOGICAL LOGAN,LOGNEG, LOGBOR
  146. C
  147. C**** Les Includes
  148. C
  149. -INC SMCOORD
  150.  
  151. -INC PPARAM
  152. -INC CCOPTIO
  153. -INC SMCHPOI
  154. POINTEUR MPROC.MPOVAL, MPVITC.MPOVAL, MPPC.MPOVAL,
  155. & MPGAMC.MPOVAL, MPNORM.MPOVAL
  156. -INC SMCHAML
  157. POINTEUR MELVNX.MELVAL, MELVNY.MELVAL, MELVNZ.MELVAL,
  158. & MELT1X.MELVAL, MELT1Y.MELVAL, MELT1Z.MELVAL,
  159. & MELT2X.MELVAL, MELT2Y.MELVAL, MELT2Z.MELVAL
  160. POINTEUR MELVUN.MELVAL, MELVUT.MELVAL, MELVUV.MELVAL
  161. POINTEUR MELRO.MELVAL, MELP.MELVAL,
  162. & MELGAM.MELVAL
  163. -INC SMLENTI
  164. -INC SMELEME
  165. C
  166. C**** Initialisation des parametres d'erreur déjà faite, i.e.
  167. C
  168. C LOGNEG = .FALSE.
  169. C LOGBOR = .FALSE.
  170. C MESERR = ' '
  171. C MOTERR(1:40) = MESERR(1:40)
  172. C VALER = 0.0D0
  173. C VAL1 = 0.0D0
  174. C VAL2 = 0.0D0
  175. C
  176. C
  177. C**** KRIPAD pour la correspondance global/local de centre
  178. C
  179. CALL KRIPAD(ICEN,MLENT1)
  180. C
  181. C**** MLENTI1 a MCORD.XCOORD(/1)/(IDIM+1) elements
  182. C
  183. C Si i est le numero global d'un noeud de ICEN,
  184. C MLENT1.LECT(i) contient sa position, i.e.
  185. C
  186. C I = numero global du noeud centre
  187. C MLENT1.LECT(i) = numero local du noeud centre
  188. C
  189. C MLENT1 déjà activé, i.e.
  190. C
  191. C SEGACT MLENT1
  192. C
  193. C**** Activation de CHPOINTs
  194. C
  195. C densité
  196. C vitesse
  197. C pression
  198. C gamma
  199. C cosinus directeurs des normales aux surface
  200. C
  201. CALL LICHT(IROC ,MPROC ,TYPE,IGEOM)
  202. CALL LICHT(IVITC,MPVITC,TYPE,IGEOM)
  203. CALL LICHT(IPC ,MPPC ,TYPE,IGEOM)
  204. CALL LICHT(IGAMC,MPGAMC,TYPE,IGEOM)
  205. CALL LICHT(INORM,MPNORM,TYPE,IGEOM)
  206. C
  207. C**** MPOVA1 - MPOVA5 sont déjà activés i.e.:
  208. C
  209. C SEGACT MPROC
  210. C SEGACT MPVITC
  211. C SEGACT MPPC
  212. C SEGACT MPGAMC
  213. C SEGACT MPNORM
  214. C
  215. C**** Le MELEME FACEL
  216. C
  217. IPT1 = IFACEL
  218. IPT2 = IFACE
  219. SEGACT IPT1
  220. SEGACT IPT2
  221. NFAC = IPT1.NUM(/2)
  222. C
  223. C**** Creation de MCHAMLs contenant les etats gauche et droite,
  224. C
  225. C i.e.:
  226. C
  227. C vitesse + cosinus directors du repere local
  228. C densité
  229. C pression
  230. C gamma
  231. C
  232. C**** Cosinus directors du repere local et vitesse
  233. C
  234. C Les cosinus directeurs
  235. C
  236. N1 = 2
  237. N3 = 6
  238. L1 = 28
  239. SEGINI MCHEL1
  240. IVITF = MCHEL1
  241. MCHEL1.TITCHE = 'U '
  242. MCHEL1.IMACHE(1) = IFACE
  243. MCHEL1.IMACHE(2) = IFACEL
  244. MCHEL1.CONCHE(1) = '(n,t,v)in(x,y,z)'
  245. MCHEL1.CONCHE(2) = ' U in (n,t,v) '
  246. C
  247. C**** Valeurs des cosinus definies par respect au repair global, i.e.
  248. C
  249. MCHEL1.INFCHE(1,1) = 2
  250. MCHEL1.INFCHE(1,3) = NIFOUR
  251. MCHEL1.INFCHE(1,4) = 0
  252. MCHEL1.INFCHE(1,5) = 0
  253. MCHEL1.INFCHE(1,6) = 0
  254. MCHEL1.IFOCHE = IFOUR
  255. C
  256. C**** Valeurs de vitesse definies par respect au repair local, i.e.
  257. C
  258. MCHEL1.INFCHE(2,1) = 1
  259. MCHEL1.INFCHE(2,3) = NIFOUR
  260. MCHEL1.INFCHE(2,4) = 0
  261. MCHEL1.INFCHE(2,5) = 0
  262. MCHEL1.INFCHE(2,6) = 0
  263. C
  264. C**** Le cosinus directeurs
  265. C
  266. N1PTEL = 1
  267. N1EL = NFAC
  268. N2PTEL = 0
  269. N2EL = 0
  270. C
  271. C**** MCHAML a N2 composantes:
  272. C
  273. C cosinus directeurs du repere local (n,t1,t2)=(n,t,v)
  274. C
  275. C IDIM = 3 -> 9 composantes
  276. C
  277. N2 = 9
  278. SEGINI MCHAM1
  279. MCHEL1.ICHAML(1) = MCHAM1
  280. MCHAM1.NOMCHE(1) = 'NX '
  281. MCHAM1.NOMCHE(2) = 'NY '
  282. MCHAM1.NOMCHE(3) = 'NZ '
  283. MCHAM1.NOMCHE(4) = 'TX '
  284. MCHAM1.NOMCHE(5) = 'TY '
  285. MCHAM1.NOMCHE(6) = 'TZ '
  286. MCHAM1.NOMCHE(7) = 'VX '
  287. MCHAM1.NOMCHE(8) = 'VY '
  288. MCHAM1.NOMCHE(9) = 'VZ '
  289. MCHAM1.TYPCHE(1) = 'REAL*8 '
  290. MCHAM1.TYPCHE(2) = 'REAL*8 '
  291. MCHAM1.TYPCHE(3) = 'REAL*8 '
  292. MCHAM1.TYPCHE(4) = 'REAL*8 '
  293. MCHAM1.TYPCHE(5) = 'REAL*8 '
  294. MCHAM1.TYPCHE(6) = 'REAL*8 '
  295. MCHAM1.TYPCHE(7) = 'REAL*8 '
  296. MCHAM1.TYPCHE(8) = 'REAL*8 '
  297. MCHAM1.TYPCHE(9) = 'REAL*8 '
  298. SEGINI MELVNX
  299. SEGINI MELVNY
  300. SEGINI MELVNZ
  301. SEGINI MELT1X
  302. SEGINI MELT1Y
  303. SEGINI MELT1Z
  304. SEGINI MELT2X
  305. SEGINI MELT2Y
  306. SEGINI MELT2Z
  307. MCHAM1.IELVAL(1) = MELVNX
  308. MCHAM1.IELVAL(2) = MELVNY
  309. MCHAM1.IELVAL(3) = MELVNZ
  310. MCHAM1.IELVAL(4) = MELT1X
  311. MCHAM1.IELVAL(5) = MELT1Y
  312. MCHAM1.IELVAL(6) = MELT1Z
  313. MCHAM1.IELVAL(7) = MELT2X
  314. MCHAM1.IELVAL(8) = MELT2Y
  315. MCHAM1.IELVAL(9) = MELT2Z
  316. SEGDES MCHAM1
  317. C
  318. C**** Vitesse
  319. C
  320. N1EL = NFAC
  321. N1PTEL = 3
  322. N2EL = 0
  323. N2PTEL = 0
  324. C
  325. C**** MCHAML a N2 composantes:
  326. C
  327. C IDIM = 3 -> 3 composantes
  328. C
  329. N2 = 3
  330. SEGINI MCHAM1
  331. MCHEL1.ICHAML(2) = MCHAM1
  332. SEGDES MCHEL1
  333. MCHAM1.NOMCHE(1) = 'UN '
  334. MCHAM1.NOMCHE(2) = 'UT '
  335. MCHAM1.NOMCHE(3) = 'UV '
  336. MCHAM1.TYPCHE(1) = 'REAL*8 '
  337. MCHAM1.TYPCHE(2) = 'REAL*8 '
  338. MCHAM1.TYPCHE(3) = 'REAL*8 '
  339. SEGINI MELVUN
  340. SEGINI MELVUT
  341. SEGINI MELVUV
  342. MCHAM1.IELVAL(1) = MELVUN
  343. MCHAM1.IELVAL(2) = MELVUT
  344. MCHAM1.IELVAL(3) = MELVUV
  345. SEGDES MCHAM1
  346. C
  347. C**** Densite
  348. C
  349. N1 = 1
  350. N3 = 6
  351. L1 = 15
  352. SEGINI MCHEL2
  353. IROF = MCHEL2
  354. MCHEL2.IMACHE(1) = IFACEL
  355. MCHEL2.TITCHE = 'RO '
  356. MCHEL2.CONCHE(1) = ' '
  357. C
  358. C**** Valeurs independente du repére, i.e.
  359. C
  360. MCHEL2.INFCHE(1,1) = 0
  361. MCHEL2.INFCHE(1,3) = NIFOUR
  362. MCHEL2.INFCHE(1,4) = 0
  363. MCHEL2.INFCHE(1,5) = 0
  364. MCHEL2.INFCHE(1,6) = 0
  365. MCHEL2.IFOCHE = IFOUR
  366. N2 = 1
  367. SEGINI MCHAM1
  368. MCHEL2.ICHAML(1) = MCHAM1
  369. SEGDES MCHEL2
  370. MCHAM1.NOMCHE(1) = 'SCAL '
  371. MCHAM1.TYPCHE(1) = 'REAL*8 '
  372. SEGINI MELRO
  373. MCHAM1.IELVAL(1) = MELRO
  374. SEGDES MCHAM1
  375. C
  376. C**** Pression
  377. C
  378. MCHEL1 = IROF
  379. SEGINI, MCHEL2 = MCHEL1
  380. IPF = MCHEL2
  381. MCHEL2.TITCHE = 'P '
  382. C
  383. C**** MCHAM1 = MCHAML de la densite
  384. C
  385. SEGINI, MCHAM2 = MCHAM1
  386. MCHEL2.ICHAML(1) = MCHAM2
  387. SEGDES MCHEL2
  388. SEGINI MELP
  389. MCHAM2.IELVAL(1) = MELP
  390. SEGDES MCHAM2
  391. C
  392. C**** Gamma
  393. C
  394. MCHEL1 = IROF
  395. SEGINI, MCHEL2 = MCHEL1
  396. IGAMF = MCHEL2
  397. MCHEL2.TITCHE = 'GAMMA '
  398. C
  399. C**** MCHAM1 = MCHAML de la densite
  400. C
  401. SEGINI, MCHAM2 = MCHAM1
  402. MCHEL2.ICHAML(1) = MCHAM2
  403. SEGDES MCHEL2
  404. SEGINI MELGAM
  405. MCHAM2.IELVAL(1) = MELGAM
  406. SEGDES MCHAM2
  407. C
  408. C**** Recapitulatif: les MELVALs et les MPOVALs actives
  409. C
  410. C MELVNX, MELVNY, MELVNZ
  411. C MELT1X, MELT1Y, MELT1Z
  412. C MELT2X, MELT2Y, MELT2Z
  413. C
  414. C MELVUN, MELVUT , MELVUN-> vitesse
  415. C
  416. C MELRO -> densite
  417. C
  418. C MELP -> pression
  419. C
  420. C MELGAM -> gamma
  421. C
  422. C MPROC -> densite
  423. C
  424. C MPVITC -> vitesse
  425. C
  426. C MPPC -> pression
  427. C
  428. C MPGAMC -> gamma
  429. C
  430. C MPNORM -> normales aux faces
  431. C
  432. C**** Boucle sur le faces
  433. C
  434. IDIMP1=IDIM+1
  435. DO NLCF = 1, NFAC
  436. C
  437. C******* NLCF = numero local du centre de face
  438. C NGCF = numero global du centre de face
  439. C NGCEG = numero global du centre ELT "gauche"
  440. C NLCEG = numero local du centre ELT "gauche"
  441. C NGCED = numero global du centre ELT "droite"
  442. C NLCED = numero local du centre ELT "droite"
  443. C
  444. NGCEG = IPT1.NUM(1,NLCF)
  445. NGCF = IPT1.NUM(2,NLCF)
  446. NGCED = IPT1.NUM(3,NLCF)
  447. NLCEG = MLENT1.LECT(NGCEG)
  448. NLCED = MLENT1.LECT(NGCED)
  449. C
  450. C******* TEST: IPT2.NUM(1,NLCF) = IPT1.NUM(2,NLCF)
  451. C
  452. NGCF1 = IPT2.NUM(1,NLCF)
  453. IF( NGCF1 .NE. NGCF) THEN
  454. LOGAN = .TRUE.
  455. MESERR(1:40) = 'PRET, subroutine pre112.eso '
  456. GOTO 9999
  457. ENDIF
  458. C
  459. C******* Cosinus directeurs des NORMALES aux faces
  460. C
  461. C On impose que les normales sont direct "Gauche" -> "Centre"
  462. C
  463. INDCEL = (NGCEG-1)*IDIMP1
  464. XG = XCOOR(INDCEL+1)
  465. YG = XCOOR(INDCEL+2)
  466. ZG = XCOOR(INDCEL+3)
  467. INDCEL = (NGCF-1)*IDIMP1
  468. XC = XCOOR(INDCEL + 1)
  469. YC = XCOOR(INDCEL + 2)
  470. ZC = XCOOR(INDCEL+3)
  471. DXG = XC - XG
  472. DYG = YC - YG
  473. DZG = ZC - ZG
  474. C
  475. C******* On calcule le sign du pruduit scalare
  476. C (Normales de Castem) * (vecteur "gauche" -> "centre")
  477. C
  478. CNX = MPNORM.VPOCHA(NLCF,7)
  479. CNY = MPNORM.VPOCHA(NLCF,8)
  480. CNZ = MPNORM.VPOCHA(NLCF,9)
  481. ORIENT = CNX * DXG + CNY * DYG + CNZ * DZG
  482. ORIENT = SIGN(1.0D0,ORIENT)
  483. IF(ORIENT .NE. 1.0D0)THEN
  484. LOGAN = .TRUE.
  485. MESERR(1:30)=
  486. & 'PRET , subroutine pre121.eso. '
  487. GOTO 9999
  488. ENDIF
  489. CNX = CNX * ORIENT
  490. CNY = CNY * ORIENT
  491. CNZ = CNZ * ORIENT
  492. C
  493. C********** Cosinus directeurs de tangente 1
  494. C
  495. CTX = MPNORM.VPOCHA(NLCF,1) * ORIENT
  496. CTY = MPNORM.VPOCHA(NLCF,2) * ORIENT
  497. CTZ = MPNORM.VPOCHA(NLCF,3) * ORIENT
  498. C
  499. C********** Cosinus directeurs de tangente 2
  500. C
  501. CVX = MPNORM.VPOCHA(NLCF,4) * ORIENT
  502. CVY = MPNORM.VPOCHA(NLCF,5) * ORIENT
  503. CVZ = MPNORM.VPOCHA(NLCF,6) * ORIENT
  504. C
  505. C******* Les autres MELVALs
  506. C
  507. C
  508. C******* N.B.: On suppose qu'on a déjà controlle RO, P > 0
  509. C GAMMA \in (1,3)
  510. C Si non il faut le faire, en utlisant LOGBOR,
  511. C LOGNEG, VALER, VAL1, VAL2
  512. C
  513. C
  514. C
  515. C******* NGCEG = NGCED -> Mur
  516. C
  517. IF( NGCEG .EQ. NGCED)THEN
  518. ROG = MPROC.VPOCHA(NLCEG , 1)
  519. PG = MPPC.VPOCHA(NLCEG, 1)
  520. GAMG = MPGAMC.VPOCHA(NLCEG, 1)
  521. UXG = MPVITC.VPOCHA(NLCEG , 1)
  522. UYG = MPVITC.VPOCHA(NLCEG , 2)
  523. UZG = MPVITC.VPOCHA(NLCEG , 3)
  524. UNG = UXG * CNX + UYG * CNY + UZG * CNZ
  525. UTG = UXG * CTX + UYG * CTY + UZG * CTZ
  526. UVG = UXG * CVX + UYG * CVY + UZG * CVZ
  527. C
  528. C********** Son etat droite
  529. C
  530. ROD = ROG
  531. PD = PG
  532. GAMD = GAMG
  533. UND = -1.0D0 * UNG
  534. UTD = UTG
  535. UVD = UVG
  536. C
  537. C************* Fin cas mur
  538. C
  539. ELSE
  540. C
  541. C************* Etat gauche
  542. C
  543. ROG = MPROC.VPOCHA(NLCEG, 1)
  544. PG = MPPC.VPOCHA(NLCEG, 1)
  545. GAMG = MPGAMC.VPOCHA(NLCEG, 1)
  546. UXG = MPVITC.VPOCHA(NLCEG , 1)
  547. UYG = MPVITC.VPOCHA(NLCEG , 2)
  548. UZG = MPVITC.VPOCHA(NLCEG , 3)
  549. UNG = UXG * CNX + UYG * CNY + UZG * CNZ
  550. UTG = UXG * CTX + UYG * CTY + UZG * CTZ
  551. UVG = UXG * CVX + UYG * CVY + UZG * CVZ
  552. C
  553. C********** Etat gauche
  554. C
  555. ROD = MPROC.VPOCHA(NLCED,1)
  556. PD = MPPC.VPOCHA(NLCED,1)
  557. GAMD = MPGAMC.VPOCHA(NLCED,1)
  558. C
  559. C************* On suppose qu'on a déjà controlle ROG, PG > 0
  560. C Si non il faut le faire!!!
  561. C
  562. UXD = MPVITC.VPOCHA(NLCED,1)
  563. UYD = MPVITC.VPOCHA(NLCED,2)
  564. UZD = MPVITC.VPOCHA(NLCED,3)
  565. UND = UXD * CNX + UYD * CNY + UZD * CNZ
  566. UTD = UXD * CTX + UYD * CTY + UZD * CTZ
  567. UVD = UXD * CVX + UYD * CVY + UZD * CVZ
  568. ENDIF
  569. C
  570. C************* Les MELVALs
  571. C
  572. MELRO.VELCHE(1,NLCF) = ROG
  573. MELRO.VELCHE(3,NLCF) = ROD
  574. MELP.VELCHE(1,NLCF) = PG
  575. MELP.VELCHE(3,NLCF) = PD
  576. MELGAM.VELCHE(1,NLCF) = GAMG
  577. MELGAM.VELCHE(3,NLCF) = GAMD
  578. MELVUN.VELCHE(1,NLCF) = UNG
  579. MELVUN.VELCHE(3,NLCF) = UND
  580. MELVUT.VELCHE(1,NLCF) = UTG
  581. MELVUT.VELCHE(3,NLCF) = UTD
  582. MELVUV.VELCHE(1,NLCF) = UVG
  583. MELVUV.VELCHE(3,NLCF) = UVD
  584. MELVNX.VELCHE(1,NLCF) = CNX
  585. MELVNY.VELCHE(1,NLCF) = CNY
  586. MELVNZ.VELCHE(1,NLCF) = CNZ
  587. MELT1X.VELCHE(1,NLCF) = CTX
  588. MELT1Y.VELCHE(1,NLCF) = CTY
  589. MELT1Z.VELCHE(1,NLCF) = CTZ
  590. MELT2X.VELCHE(1,NLCF) = CVX
  591. MELT2Y.VELCHE(1,NLCF) = CVY
  592. MELT2Z.VELCHE(1,NLCF) = CVZ
  593. ENDDO
  594. C
  595. C**** Desactivation des SEGMENTs
  596. C
  597. SEGDES IPT1
  598. SEGDES IPT2
  599. C
  600. SEGDES MPROC
  601. SEGDES MPVITC
  602. SEGDES MPPC
  603. SEGDES MPGAMC
  604. SEGDES MPNORM
  605. C
  606. SEGDES MELRO
  607. SEGDES MELP
  608. SEGDES MELGAM
  609. SEGDES MELVUN
  610. SEGDES MELVUT
  611. SEGDES MELVUV
  612. SEGDES MELVNX
  613. SEGDES MELVNY
  614. SEGDES MELVNZ
  615. SEGDES MELT1X
  616. SEGDES MELT1Y
  617. SEGDES MELT1Z
  618. SEGDES MELT2X
  619. SEGDES MELT2Y
  620. SEGDES MELT2Z
  621. C
  622. C**** Destruction du MELNTI correspondance local/global
  623. C
  624. SEGSUP MLENT1
  625. C
  626. 9999 CONTINUE
  627. C
  628. RETURN
  629. END
  630.  
  631.  
  632.  
  633.  
  634.  
  635.  
  636.  
  637.  
  638.  
  639.  

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