Télécharger pre211.eso

Retour à la liste

Numérotation des lignes :

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

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