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. -INC CCOPTIO
  151. -INC SMCHPOI
  152. POINTEUR MPROC.MPOVAL, MPVITC.MPOVAL, MPPC.MPOVAL,
  153. & MPGAMC.MPOVAL, MPNORM.MPOVAL
  154. -INC SMCHAML
  155. POINTEUR MELVNX.MELVAL, MELVNY.MELVAL, MELVNZ.MELVAL,
  156. & MELT1X.MELVAL, MELT1Y.MELVAL, MELT1Z.MELVAL,
  157. & MELT2X.MELVAL, MELT2Y.MELVAL, MELT2Z.MELVAL
  158. POINTEUR MELVUN.MELVAL, MELVUT.MELVAL, MELVUV.MELVAL
  159. POINTEUR MELRO.MELVAL, MELP.MELVAL,
  160. & MELGAM.MELVAL
  161. -INC SMLENTI
  162. -INC SMELEME
  163. C
  164. C**** Initialisation des parametres d'erreur déjà faite, i.e.
  165. C
  166. C LOGNEG = .FALSE.
  167. C LOGBOR = .FALSE.
  168. C MESERR = ' '
  169. C MOTERR(1:40) = MESERR(1:40)
  170. C VALER = 0.0D0
  171. C VAL1 = 0.0D0
  172. C VAL2 = 0.0D0
  173. C
  174. C
  175. C**** KRIPAD pour la correspondance global/local de centre
  176. C
  177. CALL KRIPAD(ICEN,MLENT1)
  178. C
  179. C**** MLENTI1 a MCORD.XCOORD(/1)/(IDIM+1) elements
  180. C
  181. C Si i est le numero global d'un noeud de ICEN,
  182. C MLENT1.LECT(i) contient sa position, i.e.
  183. C
  184. C I = numero global du noeud centre
  185. C MLENT1.LECT(i) = numero local du noeud centre
  186. C
  187. C MLENT1 déjà activé, i.e.
  188. C
  189. C SEGACT MLENT1
  190. C
  191. C**** Activation de CHPOINTs
  192. C
  193. C densité
  194. C vitesse
  195. C pression
  196. C gamma
  197. C cosinus directeurs des normales aux surface
  198. C
  199. CALL LICHT(IROC ,MPROC ,TYPE,IGEOM)
  200. CALL LICHT(IVITC,MPVITC,TYPE,IGEOM)
  201. CALL LICHT(IPC ,MPPC ,TYPE,IGEOM)
  202. CALL LICHT(IGAMC,MPGAMC,TYPE,IGEOM)
  203. CALL LICHT(INORM,MPNORM,TYPE,IGEOM)
  204. C
  205. C**** MPOVA1 - MPOVA5 sont déjà activés i.e.:
  206. C
  207. C SEGACT MPROC
  208. C SEGACT MPVITC
  209. C SEGACT MPPC
  210. C SEGACT MPGAMC
  211. C SEGACT MPNORM
  212. C
  213. C**** Le MELEME FACEL
  214. C
  215. IPT1 = IFACEL
  216. IPT2 = IFACE
  217. SEGACT IPT1
  218. SEGACT IPT2
  219. NFAC = IPT1.NUM(/2)
  220. C
  221. C**** Creation de MCHAMLs contenant les etats gauche et droite,
  222. C
  223. C i.e.:
  224. C
  225. C vitesse + cosinus directors du repere local
  226. C densité
  227. C pression
  228. C gamma
  229. C
  230. C**** Cosinus directors du repere local et vitesse
  231. C
  232. C Les cosinus directeurs
  233. C
  234. N1 = 2
  235. N3 = 6
  236. L1 = 28
  237. SEGINI MCHEL1
  238. IVITF = MCHEL1
  239. MCHEL1.TITCHE = 'U '
  240. MCHEL1.IMACHE(1) = IFACE
  241. MCHEL1.IMACHE(2) = IFACEL
  242. MCHEL1.CONCHE(1) = '(n,t,v)in(x,y,z)'
  243. MCHEL1.CONCHE(2) = ' U in (n,t,v) '
  244. C
  245. C**** Valeurs des cosinus definies par respect au repair global, i.e.
  246. C
  247. MCHEL1.INFCHE(1,1) = 2
  248. MCHEL1.INFCHE(1,3) = NIFOUR
  249. MCHEL1.INFCHE(1,4) = 0
  250. MCHEL1.INFCHE(1,5) = 0
  251. MCHEL1.INFCHE(1,6) = 0
  252. MCHEL1.IFOCHE = IFOUR
  253. C
  254. C**** Valeurs de vitesse definies par respect au repair local, i.e.
  255. C
  256. MCHEL1.INFCHE(2,1) = 1
  257. MCHEL1.INFCHE(2,3) = NIFOUR
  258. MCHEL1.INFCHE(2,4) = 0
  259. MCHEL1.INFCHE(2,5) = 0
  260. MCHEL1.INFCHE(2,6) = 0
  261. C
  262. C**** Le cosinus directeurs
  263. C
  264. N1PTEL = 1
  265. N1EL = NFAC
  266. N2PTEL = 0
  267. N2EL = 0
  268. C
  269. C**** MCHAML a N2 composantes:
  270. C
  271. C cosinus directeurs du repere local (n,t1,t2)=(n,t,v)
  272. C
  273. C IDIM = 3 -> 9 composantes
  274. C
  275. N2 = 9
  276. SEGINI MCHAM1
  277. MCHEL1.ICHAML(1) = MCHAM1
  278. MCHAM1.NOMCHE(1) = 'NX '
  279. MCHAM1.NOMCHE(2) = 'NY '
  280. MCHAM1.NOMCHE(3) = 'NZ '
  281. MCHAM1.NOMCHE(4) = 'TX '
  282. MCHAM1.NOMCHE(5) = 'TY '
  283. MCHAM1.NOMCHE(6) = 'TZ '
  284. MCHAM1.NOMCHE(7) = 'VX '
  285. MCHAM1.NOMCHE(8) = 'VY '
  286. MCHAM1.NOMCHE(9) = 'VZ '
  287. MCHAM1.TYPCHE(1) = 'REAL*8 '
  288. MCHAM1.TYPCHE(2) = 'REAL*8 '
  289. MCHAM1.TYPCHE(3) = 'REAL*8 '
  290. MCHAM1.TYPCHE(4) = 'REAL*8 '
  291. MCHAM1.TYPCHE(5) = 'REAL*8 '
  292. MCHAM1.TYPCHE(6) = 'REAL*8 '
  293. MCHAM1.TYPCHE(7) = 'REAL*8 '
  294. MCHAM1.TYPCHE(8) = 'REAL*8 '
  295. MCHAM1.TYPCHE(9) = 'REAL*8 '
  296. SEGINI MELVNX
  297. SEGINI MELVNY
  298. SEGINI MELVNZ
  299. SEGINI MELT1X
  300. SEGINI MELT1Y
  301. SEGINI MELT1Z
  302. SEGINI MELT2X
  303. SEGINI MELT2Y
  304. SEGINI MELT2Z
  305. MCHAM1.IELVAL(1) = MELVNX
  306. MCHAM1.IELVAL(2) = MELVNY
  307. MCHAM1.IELVAL(3) = MELVNZ
  308. MCHAM1.IELVAL(4) = MELT1X
  309. MCHAM1.IELVAL(5) = MELT1Y
  310. MCHAM1.IELVAL(6) = MELT1Z
  311. MCHAM1.IELVAL(7) = MELT2X
  312. MCHAM1.IELVAL(8) = MELT2Y
  313. MCHAM1.IELVAL(9) = MELT2Z
  314. SEGDES MCHAM1
  315. C
  316. C**** Vitesse
  317. C
  318. N1EL = NFAC
  319. N1PTEL = 3
  320. N2EL = 0
  321. N2PTEL = 0
  322. C
  323. C**** MCHAML a N2 composantes:
  324. C
  325. C IDIM = 3 -> 3 composantes
  326. C
  327. N2 = 3
  328. SEGINI MCHAM1
  329. MCHEL1.ICHAML(2) = MCHAM1
  330. SEGDES MCHEL1
  331. MCHAM1.NOMCHE(1) = 'UN '
  332. MCHAM1.NOMCHE(2) = 'UT '
  333. MCHAM1.NOMCHE(3) = 'UV '
  334. MCHAM1.TYPCHE(1) = 'REAL*8 '
  335. MCHAM1.TYPCHE(2) = 'REAL*8 '
  336. MCHAM1.TYPCHE(3) = 'REAL*8 '
  337. SEGINI MELVUN
  338. SEGINI MELVUT
  339. SEGINI MELVUV
  340. MCHAM1.IELVAL(1) = MELVUN
  341. MCHAM1.IELVAL(2) = MELVUT
  342. MCHAM1.IELVAL(3) = MELVUV
  343. SEGDES MCHAM1
  344. C
  345. C**** Densite
  346. C
  347. N1 = 1
  348. N3 = 6
  349. L1 = 15
  350. SEGINI MCHEL2
  351. IROF = MCHEL2
  352. MCHEL2.IMACHE(1) = IFACEL
  353. MCHEL2.TITCHE = 'RO '
  354. MCHEL2.CONCHE(1) = ' '
  355. C
  356. C**** Valeurs independente du repére, i.e.
  357. C
  358. MCHEL2.INFCHE(1,1) = 0
  359. MCHEL2.INFCHE(1,3) = NIFOUR
  360. MCHEL2.INFCHE(1,4) = 0
  361. MCHEL2.INFCHE(1,5) = 0
  362. MCHEL2.INFCHE(1,6) = 0
  363. MCHEL2.IFOCHE = IFOUR
  364. N2 = 1
  365. SEGINI MCHAM1
  366. MCHEL2.ICHAML(1) = MCHAM1
  367. SEGDES MCHEL2
  368. MCHAM1.NOMCHE(1) = 'SCAL '
  369. MCHAM1.TYPCHE(1) = 'REAL*8 '
  370. SEGINI MELRO
  371. MCHAM1.IELVAL(1) = MELRO
  372. SEGDES MCHAM1
  373. C
  374. C**** Pression
  375. C
  376. MCHEL1 = IROF
  377. SEGINI, MCHEL2 = MCHEL1
  378. IPF = MCHEL2
  379. MCHEL2.TITCHE = 'P '
  380. C
  381. C**** MCHAM1 = MCHAML de la densite
  382. C
  383. SEGINI, MCHAM2 = MCHAM1
  384. MCHEL2.ICHAML(1) = MCHAM2
  385. SEGDES MCHEL2
  386. SEGINI MELP
  387. MCHAM2.IELVAL(1) = MELP
  388. SEGDES MCHAM2
  389. C
  390. C**** Gamma
  391. C
  392. MCHEL1 = IROF
  393. SEGINI, MCHEL2 = MCHEL1
  394. IGAMF = MCHEL2
  395. MCHEL2.TITCHE = 'GAMMA '
  396. C
  397. C**** MCHAM1 = MCHAML de la densite
  398. C
  399. SEGINI, MCHAM2 = MCHAM1
  400. MCHEL2.ICHAML(1) = MCHAM2
  401. SEGDES MCHEL2
  402. SEGINI MELGAM
  403. MCHAM2.IELVAL(1) = MELGAM
  404. SEGDES MCHAM2
  405. C
  406. C**** Recapitulatif: les MELVALs et les MPOVALs actives
  407. C
  408. C MELVNX, MELVNY, MELVNZ
  409. C MELT1X, MELT1Y, MELT1Z
  410. C MELT2X, MELT2Y, MELT2Z
  411. C
  412. C MELVUN, MELVUT , MELVUN-> vitesse
  413. C
  414. C MELRO -> densite
  415. C
  416. C MELP -> pression
  417. C
  418. C MELGAM -> gamma
  419. C
  420. C MPROC -> densite
  421. C
  422. C MPVITC -> vitesse
  423. C
  424. C MPPC -> pression
  425. C
  426. C MPGAMC -> gamma
  427. C
  428. C MPNORM -> normales aux faces
  429. C
  430. C**** Boucle sur le faces
  431. C
  432. IDIMP1=IDIM+1
  433. DO NLCF = 1, NFAC
  434. C
  435. C******* NLCF = numero local du centre de face
  436. C NGCF = numero global du centre de face
  437. C NGCEG = numero global du centre ELT "gauche"
  438. C NLCEG = numero local du centre ELT "gauche"
  439. C NGCED = numero global du centre ELT "droite"
  440. C NLCED = numero local du centre ELT "droite"
  441. C
  442. NGCEG = IPT1.NUM(1,NLCF)
  443. NGCF = IPT1.NUM(2,NLCF)
  444. NGCED = IPT1.NUM(3,NLCF)
  445. NLCEG = MLENT1.LECT(NGCEG)
  446. NLCED = MLENT1.LECT(NGCED)
  447. C
  448. C******* TEST: IPT2.NUM(1,NLCF) = IPT1.NUM(2,NLCF)
  449. C
  450. NGCF1 = IPT2.NUM(1,NLCF)
  451. IF( NGCF1 .NE. NGCF) THEN
  452. LOGAN = .TRUE.
  453. MESERR(1:40) = 'PRET, subroutine pre112.eso '
  454. GOTO 9999
  455. ENDIF
  456. C
  457. C******* Cosinus directeurs des NORMALES aux faces
  458. C
  459. C On impose que les normales sont direct "Gauche" -> "Centre"
  460. C
  461. INDCEL = (NGCEG-1)*IDIMP1
  462. XG = XCOOR(INDCEL+1)
  463. YG = XCOOR(INDCEL+2)
  464. ZG = XCOOR(INDCEL+3)
  465. INDCEL = (NGCF-1)*IDIMP1
  466. XC = XCOOR(INDCEL + 1)
  467. YC = XCOOR(INDCEL + 2)
  468. ZC = XCOOR(INDCEL+3)
  469. DXG = XC - XG
  470. DYG = YC - YG
  471. DZG = ZC - ZG
  472. C
  473. C******* On calcule le sign du pruduit scalare
  474. C (Normales de Castem) * (vecteur "gauche" -> "centre")
  475. C
  476. CNX = MPNORM.VPOCHA(NLCF,7)
  477. CNY = MPNORM.VPOCHA(NLCF,8)
  478. CNZ = MPNORM.VPOCHA(NLCF,9)
  479. ORIENT = CNX * DXG + CNY * DYG + CNZ * DZG
  480. ORIENT = SIGN(1.0D0,ORIENT)
  481. IF(ORIENT .NE. 1.0D0)THEN
  482. LOGAN = .TRUE.
  483. MESERR(1:30)=
  484. & 'PRET , subroutine pre121.eso. '
  485. GOTO 9999
  486. ENDIF
  487. CNX = CNX * ORIENT
  488. CNY = CNY * ORIENT
  489. CNZ = CNZ * ORIENT
  490. C
  491. C********** Cosinus directeurs de tangente 1
  492. C
  493. CTX = MPNORM.VPOCHA(NLCF,1) * ORIENT
  494. CTY = MPNORM.VPOCHA(NLCF,2) * ORIENT
  495. CTZ = MPNORM.VPOCHA(NLCF,3) * ORIENT
  496. C
  497. C********** Cosinus directeurs de tangente 2
  498. C
  499. CVX = MPNORM.VPOCHA(NLCF,4) * ORIENT
  500. CVY = MPNORM.VPOCHA(NLCF,5) * ORIENT
  501. CVZ = MPNORM.VPOCHA(NLCF,6) * ORIENT
  502. C
  503. C******* Les autres MELVALs
  504. C
  505. C
  506. C******* N.B.: On suppose qu'on a déjà controlle RO, P > 0
  507. C GAMMA \in (1,3)
  508. C Si non il faut le faire, en utlisant LOGBOR,
  509. C LOGNEG, VALER, VAL1, VAL2
  510. C
  511. C
  512. C
  513. C******* NGCEG = NGCED -> Mur
  514. C
  515. IF( NGCEG .EQ. NGCED)THEN
  516. ROG = MPROC.VPOCHA(NLCEG , 1)
  517. PG = MPPC.VPOCHA(NLCEG, 1)
  518. GAMG = MPGAMC.VPOCHA(NLCEG, 1)
  519. UXG = MPVITC.VPOCHA(NLCEG , 1)
  520. UYG = MPVITC.VPOCHA(NLCEG , 2)
  521. UZG = MPVITC.VPOCHA(NLCEG , 3)
  522. UNG = UXG * CNX + UYG * CNY + UZG * CNZ
  523. UTG = UXG * CTX + UYG * CTY + UZG * CTZ
  524. UVG = UXG * CVX + UYG * CVY + UZG * CVZ
  525. C
  526. C********** Son etat droite
  527. C
  528. ROD = ROG
  529. PD = PG
  530. GAMD = GAMG
  531. UND = -1.0D0 * UNG
  532. UTD = UTG
  533. UVD = UVG
  534. C
  535. C************* Fin cas mur
  536. C
  537. ELSE
  538. C
  539. C************* Etat gauche
  540. C
  541. ROG = MPROC.VPOCHA(NLCEG, 1)
  542. PG = MPPC.VPOCHA(NLCEG, 1)
  543. GAMG = MPGAMC.VPOCHA(NLCEG, 1)
  544. UXG = MPVITC.VPOCHA(NLCEG , 1)
  545. UYG = MPVITC.VPOCHA(NLCEG , 2)
  546. UZG = MPVITC.VPOCHA(NLCEG , 3)
  547. UNG = UXG * CNX + UYG * CNY + UZG * CNZ
  548. UTG = UXG * CTX + UYG * CTY + UZG * CTZ
  549. UVG = UXG * CVX + UYG * CVY + UZG * CVZ
  550. C
  551. C********** Etat gauche
  552. C
  553. ROD = MPROC.VPOCHA(NLCED,1)
  554. PD = MPPC.VPOCHA(NLCED,1)
  555. GAMD = MPGAMC.VPOCHA(NLCED,1)
  556. C
  557. C************* On suppose qu'on a déjà controlle ROG, PG > 0
  558. C Si non il faut le faire!!!
  559. C
  560. UXD = MPVITC.VPOCHA(NLCED,1)
  561. UYD = MPVITC.VPOCHA(NLCED,2)
  562. UZD = MPVITC.VPOCHA(NLCED,3)
  563. UND = UXD * CNX + UYD * CNY + UZD * CNZ
  564. UTD = UXD * CTX UYD * CTY + UZD * CTZ
  565. UVD = UXD * CVX + UYD * CVY + UZD * CVZ
  566. ENDIF
  567. C
  568. C************* Les MELVALs
  569. C
  570. MELRO.VELCHE(1,NLCF) = ROG
  571. MELRO.VELCHE(3,NLCF) = ROD
  572. MELP.VELCHE(1,NLCF) = PG
  573. MELP.VELCHE(3,NLCF) = PD
  574. MELGAM.VELCHE(1,NLCF) = GAMG
  575. MELGAM.VELCHE(3,NLCF) = GAMD
  576. MELVUN.VELCHE(1,NLCF) = UNG
  577. MELVUN.VELCHE(3,NLCF) = UND
  578. MELVUT.VELCHE(1,NLCF) = UTG
  579. MELVUT.VELCHE(3,NLCF) = UTD
  580. MELVUV.VELCHE(1,NLCF) = UVG
  581. MELVUV.VELCHE(3,NLCF) = UVD
  582. MELVNX.VELCHE(1,NLCF) = CNX
  583. MELVNY.VELCHE(1,NLCF) = CNY
  584. MELVNZ.VELCHE(1,NLCF) = CNZ
  585. MELT1X.VELCHE(1,NLCF) = CTX
  586. MELT1Y.VELCHE(1,NLCF) = CTY
  587. MELT1Z.VELCHE(1,NLCF) = CTZ
  588. MELT2X.VELCHE(1,NLCF) = CVX
  589. MELT2Y.VELCHE(1,NLCF) = CVY
  590. MELT2Z.VELCHE(1,NLCF) = CVZ
  591. ENDDO
  592. C
  593. C**** Desactivation des SEGMENTs
  594. C
  595. SEGDES IPT1
  596. SEGDES IPT2
  597. C
  598. SEGDES MPROC
  599. SEGDES MPVITC
  600. SEGDES MPPC
  601. SEGDES MPGAMC
  602. SEGDES MPNORM
  603. C
  604. SEGDES MELRO
  605. SEGDES MELP
  606. SEGDES MELGAM
  607. SEGDES MELVUN
  608. SEGDES MELVUT
  609. SEGDES MELVUV
  610. SEGDES MELVNX
  611. SEGDES MELVNY
  612. SEGDES MELVNZ
  613. SEGDES MELT1X
  614. SEGDES MELT1Y
  615. SEGDES MELT1Z
  616. SEGDES MELT2X
  617. SEGDES MELT2Y
  618. SEGDES MELT2Z
  619. C
  620. C**** Destruction du MELNTI correspondance local/global
  621. C
  622. SEGSUP MLENT1
  623. C
  624. 9999 CONTINUE
  625. C
  626. RETURN
  627. END
  628.  
  629.  
  630.  
  631.  
  632.  
  633.  
  634.  
  635.  
  636.  
  637.  

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