Télécharger la plaquette Cast3M


fr en

Télécharger konjp3.eso

Retour à la liste

Numérotation des lignes : 

konjp3

  1. C KONJP3    SOURCE    CB215821  20/11/25    13:32:26     10792          
  2.       SUBROUTINE KONJP3(ILINC,ILINP,IRN,IUN,IPN,IGAMN,INORM,ICHPVO
  3.      $     ,ICHPSU,MELEMC,MELEFE,MELLIM,IMAT)
  4. C
  5. C************************************************************************
  6. C
  7. C PROJET            :  CASTEM 2000
  8. C
  9. C NOM               :  KONJA3
  10. C
  11. C DESCRIPTION       :  Voir KON14
  12. C                      Calcul du jacobien du résidu pour la méthode de
  13. C                      VLH (variable primales = variables primitives,
  14. C                      variable duales = variables conservatives)
  15. C
  16. C                      Cas 3D, gaz "calorically perfect"
  17. C
  18. C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  19. C
  20. C AUTEUR            :  A. BECCANTINI, SFME/LTMF
  21. C
  22. C************************************************************************
  23. C
  24. C
  25. C APPELES (Outils
  26. C          CASTEM)  :  KRIPAD, LICHT, ERREUR
  27. C
  28. C APPELES (Calcul)  :  VLHJ5, VLHJ7
  29. C
  30. C************************************************************************
  31. C
  32. C ENTREES
  33. C
  34. C     ILINC  : liste des inconnues duales (pointeur d'un LISTMOTS)
  35. C
  36. C     ILINP  : liste des inconnues primales (pointeur d'un LISTMOTS)
  37. C
  38. C  1) Pointeurs des CHPOINT
  39. C
  40. C     IRN     : CHPOINT CENTRE  contenant la masse volumique ;
  41. C
  42. C     IUN     : CHPOINT CENTRE  contenant la vitesse ;
  43. C
  44. C     IPN     : CHPOINT CENTRE  contenant la pression ;
  45. C
  46. C     IGAMN   : CHPOINT CENTRE  contenant le gamma ;
  47. C
  48. C     INORM   : CHPOINT FACE    contenant les normales aux faces ;
  49. C
  50. C     ICHPOVO : CHPOINT VOLUME  contenant le volume
  51. C
  52. C     ICHPOSU : CHPOINT FACE    contenant la surface des faces
  53. C
  54. C
  55. C  2) Pointeurs de MELEME de la table DOMAINE
  56. C
  57. C     MELEMC  : MELEME 'CENTRE' du SPG des CENTRES
  58. C
  59. C     MELEFE  : MELEME 'FACEL' du connectivité Faces -> Elts
  60. C
  61. C     MELLIM  : MELEME SPG des conditions aux bords
  62. C
  63. C SORTIES
  64. C
  65. C     IMAT    : pointeur de la MATRIK du jacobien du residu
  66. C
  67. C************************************************************************
  68. C
  69. C HISTORIQUE (Anomalies et modifications éventuelles)
  70. C
  71. C HISTORIQUE :
  72. C
  73. C************************************************************************
  74. C
  75. C
  76. C N.B.: On suppose qu'on a déjà controllé RO, P > 0
  77. C                                         GAMMA \in (1,3)
  78. C       Si non il faut le faire!!!
  79. C
  80. C************************************************************************
  81. C
  82. C
  83. C**** Variables de COOPTIO
  84. C
  85. C      INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  86. C     &        ,IOTER,   IOLEC,  IOIMP,     IOCAR,  IOACQ
  87. C     &        ,IOPER,   IOSGB,  IOGRA,     IOSAU,  IORES
  88. C     &        ,IECHO,   IIMPI,  IOSPI
  89. C     &        ,IDIM,    IFICLE, IPREFI
  90. C     &        ,MCOORD
  91. C     &        ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  92. C     &        ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  93. C     &        ,NORINC,NORVAL,NORIND,NORVAD
  94. C     &        ,NUCROU, IPSAUV
  95. C
  96.       IMPLICIT INTEGER(I-N)
  97.       INTEGER ILINC, ILINP, IRN,IUN,IPN,IGAMN,INORM,ICHPVO,ICHPSU
  98.      &        , IMAT, IGEOMC, IGEOMF
  99.      &        , NFAC, NBSOUS, NBREF, NBELEM, NBNN, NRIGE, NMATRI, NKID
  100.      &        , NKMT, NBME, NBEL,  MP, NP
  101.      &        , IFAC, NGCF, NLCF, NGCG, NGCD, NLCG, NLCD, NLFL
  102.       REAL*8 ROG, PG, UXG, UYG, UZG, RETG, GAMG, VOLG
  103.      &       , ROD, PD, UXD, UYD, UZD, RETD, GAMD, VOLD
  104.      &       , SURF, CNX, CNY, CNZ, CT1X, CT1Y, CT1Z, CT2X, CT2Y, CT2Z
  105.      &       , FUNCEL
  106.      &       , DFRO(5), DFRET(5), DFRUN(5), DFRUT1(5), DFRUT2(5)
  107.       CHARACTER*8 TYPE
  108. C
  109. C**** LES INCLUDES
  110. C
  111.  
  112. -INC PPARAM
  113. -INC CCOPTIO
  114. -INC SMCHPOI
  116. -INC SMELEME
  117. -INC SMLMOTS
  118. -INC SMLENTI
  119.       POINTEUR MPRN.MPOVAL, MPUN.MPOVAL, MPPN.MPOVAL, MPGAMN.MPOVAL,
  120.      &         MPNORM.MPOVAL, MPVOLU.MPOVAL, MPOVSU.MPOVAL
  121.       POINTEUR MELEMC.MELEME, MELEMF.MELEME, MELEFE.MELEME,
  122.      &         MELEDU.MELEME, MELLIM.MELEME
  123.       POINTEUR MLENTC.MLENTI, MLENTF.MLENTI, MLELIM.MLENTI
  124.       POINTEUR RR.IZAFM, RUX.IZAFM, RUY.IZAFM, RUZ.IZAFM, RP.IZAFM,
  125.      &     UXR.IZAFM, UXUX.IZAFM, UXUY.IZAFM, UXUZ.IZAFM, UXP.IZAFM,
  126.      &     UYR.IZAFM, UYUX.IZAFM, UYUY.IZAFM, UYUZ.IZAFM, UYP.IZAFM,
  127.      &     UZR.IZAFM, UZUX.IZAFM, UZUY.IZAFM, UZUZ.IZAFM, UZP.IZAFM,
  128.      &     RETR.IZAFM, RETUX.IZAFM, RETUY.IZAFM, RETUZ.IZAFM,
  129.      &     RETP.IZAFM
  130.       POINTEUR MLMINC.MLMOTS
  131. C
  132. C**** KRIPAD pour la correspondance global/local des conditions limits
  133. C
  134.       CALL KRIPAD(MELLIM,MLELIM)
  135. c     SEGACT MELLIM
  136. C
  137. C**** KRIPAD pour la correspondance global/local des centres
  138. C
  139.       CALL KRIPAD(MELEMC,MLENTC)
  140. C
  141. C     SEGACT MLENTC
  142.       SEGACT MELEMC
  143. C
  144.       SEGACT MELEFE
  145. C
  146.       CALL LICHT(ICHPSU,MPOVSU,TYPE,IGEOMF)
  147.       CALL LICHT(INORM,MPNORM,TYPE,IGEOMF)
  148.       CALL LICHT(ICHPVO,MPVOLU,TYPE,IGEOMC)
  149. C
  150. C**** LICHT active les MPOVALs en *MOD
  151. C
  152. C     i.e.
  153. C
  154. C     SEGACT MPOVSU*MOD
  155. C     SEGACT MPOVNO*MOD
  156. C     SEGACT MPVOLU*MOD
  157. C
  158.       MELEMF = IGEOMF
  159.       CALL KRIPAD(MELEMF,MLENTF)
  160. C
  161. C     SEGACT MLENTF
  162.       SEGACT MELEMF
  163. C
  164.       CALL LICHT(IRN,MPRN,TYPE,IGEOMC)
  165.       CALL LICHT(IPN,MPPN,TYPE,IGEOMC)
  166.       CALL LICHT(IUN,MPUN,TYPE,IGEOMC)
  167.       CALL LICHT(IGAMN,MPGAMN,TYPE,IGEOMC)
  168. C
  169. C     SEGACT MPRN*MOD
  170. C     SEGACT MPPN*MOD
  171. C     SEGACT MPUN*MOD
  172. C     SEGACT MPGAMN*MOD
  173. C
  174.       NFAC = MELEFE.NUM(/2)
  175. C
  176. C**** Maillage des inconnues primales
  177. C
  178.       NBSOUS = 0
  179.       NBREF = 0
  180.       NBELEM =  NFAC
  181.       NBNN = 2
  182. C
  183.       SEGINI MELEDU
  184. C     MELEPR = MELEDU
  185. C
  186. C**** MELEDU = 'SEG2'
  187. C
  188.       MELEDU.ITYPEL =  2
  189. C
  190.       NRIGE = 7
  191.       NMATRI = 1
  192.       NKID = 9
  193.       NKMT = 7
  194. C
  195.       SEGINI MATRIK
  196.       IMAT = MATRIK
  197.       MATRIK.IRIGEL(1,1) = MELEDU
  198.       MATRIK.IRIGEL(2,1) = MELEDU
  199. C
  200. C**** Matrice non symetrique
  201. C
  202.       MATRIK.IRIGEL(7,1) = 2
  203. C
  204.       NBME = 25
  205.       NBSOUS = 1
  206.       SEGINI IMATRI
  207.       MATRIK.IRIGEL(4,1) = IMATRI
  208. C
  209. C**** Variables primales (primitives)
  210. C
  211.       MLMINC = ILINP
  212.       SEGACT MLMINC
  213.       IMATRI.LISPRI(1) = MLMINC.MOTS(1)
  214.       IMATRI.LISPRI(2) = MLMINC.MOTS(2)
  215.       IMATRI.LISPRI(3) = MLMINC.MOTS(3)
  216.       IMATRI.LISPRI(4) = MLMINC.MOTS(4)
  217.       IMATRI.LISPRI(5) = MLMINC.MOTS(5)
  218. C
  219.       IMATRI.LISPRI(6)  = MLMINC.MOTS(1)
  220.       IMATRI.LISPRI(7)  = MLMINC.MOTS(2)
  221.       IMATRI.LISPRI(8)  = MLMINC.MOTS(3)
  222.       IMATRI.LISPRI(9)  = MLMINC.MOTS(4)
  223.       IMATRI.LISPRI(10) = MLMINC.MOTS(5)
  224. C
  225.       IMATRI.LISPRI(11) = MLMINC.MOTS(1)
  226.       IMATRI.LISPRI(12) = MLMINC.MOTS(2)
  227.       IMATRI.LISPRI(13) = MLMINC.MOTS(3)
  228.       IMATRI.LISPRI(14) = MLMINC.MOTS(4)
  229.       IMATRI.LISPRI(15) = MLMINC.MOTS(5)
  230. C
  231.       IMATRI.LISPRI(16) = MLMINC.MOTS(1)
  232.       IMATRI.LISPRI(17) = MLMINC.MOTS(2)
  233.       IMATRI.LISPRI(18) = MLMINC.MOTS(3)
  234.       IMATRI.LISPRI(19) = MLMINC.MOTS(4)
  235.       IMATRI.LISPRI(20) = MLMINC.MOTS(5)
  236. C
  237.       IMATRI.LISPRI(21) = MLMINC.MOTS(1)
  238.       IMATRI.LISPRI(22) = MLMINC.MOTS(2)
  239.       IMATRI.LISPRI(23) = MLMINC.MOTS(3)
  240.       IMATRI.LISPRI(24) = MLMINC.MOTS(4)
  241.       IMATRI.LISPRI(25) = MLMINC.MOTS(5)
  242.       SEGDES MLMINC
  243. C
  244. C**** Variables duales (conservatives)
  245. C
  246.       MLMINC = ILINC
  247.       SEGACT MLMINC
  248.       IMATRI.LISDUA(1) = MLMINC.MOTS(1)
  249.       IMATRI.LISDUA(2) = MLMINC.MOTS(1)
  250.       IMATRI.LISDUA(3) = MLMINC.MOTS(1)
  251.       IMATRI.LISDUA(4) = MLMINC.MOTS(1)
  252.       IMATRI.LISDUA(5) = MLMINC.MOTS(1)
  253. C
  254.       IMATRI.LISDUA(6)  = MLMINC.MOTS(2)
  255.       IMATRI.LISDUA(7)  = MLMINC.MOTS(2)
  256.       IMATRI.LISDUA(8)  = MLMINC.MOTS(2)
  257.       IMATRI.LISDUA(9)  = MLMINC.MOTS(2)
  258.       IMATRI.LISDUA(10) = MLMINC.MOTS(2)
  259. C
  260.       IMATRI.LISDUA(11) = MLMINC.MOTS(3)
  261.       IMATRI.LISDUA(12) = MLMINC.MOTS(3)
  262.       IMATRI.LISDUA(13) = MLMINC.MOTS(3)
  263.       IMATRI.LISDUA(14) = MLMINC.MOTS(3)
  264.       IMATRI.LISDUA(15) = MLMINC.MOTS(3)
  265. C
  266.       IMATRI.LISDUA(16) = MLMINC.MOTS(4)
  267.       IMATRI.LISDUA(17) = MLMINC.MOTS(4)
  268.       IMATRI.LISDUA(18) = MLMINC.MOTS(4)
  269.       IMATRI.LISDUA(19) = MLMINC.MOTS(4)
  270.       IMATRI.LISDUA(20) = MLMINC.MOTS(4)
  271. C
  272.       IMATRI.LISDUA(21) = MLMINC.MOTS(5)
  273.       IMATRI.LISDUA(22) = MLMINC.MOTS(5)
  274.       IMATRI.LISDUA(23) = MLMINC.MOTS(5)
  275.       IMATRI.LISDUA(24) = MLMINC.MOTS(5)
  276.       IMATRI.LISDUA(25) = MLMINC.MOTS(5)
  277.       SEGDES MLMINC
  278. C
  279.       NBEL = NBELEM
  280.       NBSOUS = 1
  281.       NP = 2
  282.       MP = 2
  283.       SEGINI   RR  , RUX , RUY , RUZ, RP ,
  284.      &         UXR , UXUX , UXUY , UXUZ, UXP ,
  285.      &         UYR , UYUX , UYUY , UYUZ, UYP ,
  286.      &         UZR , UZUX , UZUY , UZUZ, UZP ,
  287.      &         RETR , RETUX , RETUY , RETUZ, RETP
  288. C
  289. C**** Duale = IMATRI.LISDUA(1) = 'RN'
  290. C     Primale = IMATRI.LISPRI(1) = 'RN'
  291. C     -> IMATRI.LIZAFM(1,1) = RR
  292. C
  293. C     Duale = IMATRI.LISDUA(2) = 'RN'
  294. C     Primale = IMATRI.LISPRI(1) = 'RUXN'
  295. C     -> IMATRI.LIZAFM(1,2) = RUX
  296. C     ...
  297. C
  298.       IMATRI.LIZAFM(1,1) = RR
  299.       IMATRI.LIZAFM(1,2) = RUX
  300.       IMATRI.LIZAFM(1,3) = RUY
  301.       IMATRI.LIZAFM(1,4) = RUZ
  302.       IMATRI.LIZAFM(1,5) = RP
  303. C
  304.       IMATRI.LIZAFM(1,6) = UXR
  305.       IMATRI.LIZAFM(1,7) = UXUX
  306.       IMATRI.LIZAFM(1,8) = UXUY
  307.       IMATRI.LIZAFM(1,9) = UXUZ
  308.       IMATRI.LIZAFM(1,10) = UXP
  309. C
  310.       IMATRI.LIZAFM(1,11) = UYR
  311.       IMATRI.LIZAFM(1,12) = UYUX
  312.       IMATRI.LIZAFM(1,13) = UYUY
  313.       IMATRI.LIZAFM(1,14) = UYUZ
  314.       IMATRI.LIZAFM(1,15) = UYP
  315. C
  316.       IMATRI.LIZAFM(1,16) = UZR
  317.       IMATRI.LIZAFM(1,17) = UZUX
  318.       IMATRI.LIZAFM(1,18) = UZUY
  319.       IMATRI.LIZAFM(1,19) = UZUZ
  320.       IMATRI.LIZAFM(1,20) = UZP
  321. C
  322.       IMATRI.LIZAFM(1,21) = RETR
  323.       IMATRI.LIZAFM(1,22) = RETUX
  324.       IMATRI.LIZAFM(1,23) = RETUY
  325.       IMATRI.LIZAFM(1,24) = RETUZ
  326.       IMATRI.LIZAFM(1,25) = RETP
  327. C
  328.       DO IFAC = 1, NFAC, 1
  329.          NGCF = MELEFE.NUM(2,IFAC)
  330.          NLCF = MLENTF.LECT(NGCF)
  331.          IF(NLCF .NE. IFAC)THEN
  332.             WRITE(IOIMP,*) 'Il ne faut pas jouer avec la table domaine'
  333.             CALL ERREUR(5)
  334.             GOTO 9999
  335.          ENDIF
  336.          NGCG = MELEFE.NUM(1,IFAC)
  337.          NGCD = MELEFE.NUM(3,IFAC)
  338.          NLFL = MLELIM.LECT(NGCF)
  339.          IF(NLFL .NE. 0)THEN
  340. C
  341. C********** The point belongs on BC -> No contribution to jacobian!
  342. C
  343.             MELEDU.NUM(1,IFAC) = NGCG
  344.             MELEDU.NUM(2,IFAC) = NGCD
  345.          ELSEIF(NGCG .NE. NGCD)THEN
  346. C
  347. C********** Les MELEMEs
  348. C
  349.             MELEDU.NUM(1,IFAC) = NGCG
  350.             MELEDU.NUM(2,IFAC) = NGCD
  351. C
  352. C********** Les etats G et D
  353. C
  354.             NLCG = MLENTC.LECT(NGCG)
  355.             NLCD = MLENTC.LECT(NGCD)
  356. C
  357.             ROG = MPRN.VPOCHA(NLCG,1)
  358.             PG = MPPN.VPOCHA(NLCG,1)
  359.             UXG = MPUN.VPOCHA(NLCG,1)
  360.             UYG = MPUN.VPOCHA(NLCG,2)
  361.             UZG = MPUN.VPOCHA(NLCG,3)
  362.             GAMG = MPGAMN.VPOCHA(NLCG,1)
  363.             RETG= PG / (GAMG - 1.0D0) + 0.5D0 * ROG *
  364.      &           (UXG * UXG + UYG * UYG + UZG * UZG)
  365.             VOLG = MPVOLU.VPOCHA(NLCG,1)
  366. C
  367.             ROD = MPRN.VPOCHA(NLCD,1)
  368.             PD = MPPN.VPOCHA(NLCD,1)
  369.             UXD = MPUN.VPOCHA(NLCD,1)
  370.             UYD = MPUN.VPOCHA(NLCD,2)
  371.             UZD = MPUN.VPOCHA(NLCD,3)
  372.             GAMD = MPGAMN.VPOCHA(NLCD,1)
  373.             RETD= PD / (GAMD - 1.0D0) + 0.5D0 * ROD *
  374.      &           (UXD * UXD + UYD * UYD + UZD * UZD)
  375.             VOLD = MPVOLU.VPOCHA(NLCD,1)
  376. C
  377. C********** La normale G->D
  378. C           La tangente
  379. C
  380.             SURF = MPOVSU.VPOCHA(NLCF,1)
  381.             CNX = MPNORM.VPOCHA(NLCF,7)
  382.             CNY = MPNORM.VPOCHA(NLCF,8)
  383.             CNZ = MPNORM.VPOCHA(NLCF,9)
  384. C
  385. C********** Cosinus directeurs de tangente 1
  386. C
  387.             CT1X = MPNORM.VPOCHA(NLCF,1)
  388.             CT1Y = MPNORM.VPOCHA(NLCF,2)
  389.             CT1Z = MPNORM.VPOCHA(NLCF,3)
  390. C
  391. C********** Cosinus directeurs de tangente 2
  392. C
  393.             CT2X = MPNORM.VPOCHA(NLCF,4)
  394.             CT2Y = MPNORM.VPOCHA(NLCF,5)
  395.             CT2Z = MPNORM.VPOCHA(NLCF,6)
  396. C
  397. C********** La contribution de Gauche
  398. C
  399.             CALL VLHJ5(ROG,UXG,UYG,UZG,PG,RETG,GAMG,CNX,CNY,CNZ,
  400.      &           CT1X,CT1Y,CT1Z,CT2X,CT2Y,CT2Z,
  401.      &           DFRO,DFRUN,DFRUT1,DFRUT2,DFRET)
  402. C
  403. C
  404. C********** AB.AM(IFAC,IPRIM,IDUAL)
  405. C           A = nom de l'inconnu duale (Ro,rUX,rUY,RET)
  406. C           B = nom de l'inconnu primale (Ro,rUX,rUY,RET)
  407. C           IPRIM = 1, 2  -> G, D
  408. C           IDUAL = 1, 2  -> G, D
  409. C           i.e.
  410. C           A_IDUAL = AB.AM(IFAC,IPRIM,IDUAL) * B_IPRIM + ...
  411. C
  412. C
  413. C********** Dual RN
  414. C
  415.             FUNCEL = SURF * DFRO(1)
  416.             RR.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  417.             RR.AM(IFAC,1,2) = FUNCEL / VOLD
  418. C
  419.             FUNCEL = SURF * DFRO(2)
  420.             RUX.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  421.             RUX.AM(IFAC,1,2) = FUNCEL / VOLD
  422. C
  423.             FUNCEL = SURF * DFRO(3)
  424.             RUY.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  425.             RUY.AM(IFAC,1,2) = FUNCEL / VOLD
  426. C
  427.             FUNCEL = SURF * DFRO(4)
  428.             RUZ.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  429.             RUZ.AM(IFAC,1,2) = FUNCEL / VOLD
  430. C
  431.             FUNCEL = SURF * DFRO(5)
  432.             RP.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  433.             RP.AM(IFAC,1,2) = FUNCEL / VOLD
  434. C
  435. C********** Dual RUXN
  436. C
  437.             FUNCEL = SURF * (DFRUN(1) * CNX + DFRUT1(1) * CT1X
  438.      &           + DFRUT2(1) * CT2X)
  439.             UXR.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  440.             UXR.AM(IFAC,1,2) = FUNCEL / VOLD
  441. C
  442.             FUNCEL = SURF * (DFRUN(2) * CNX + DFRUT1(2) * CT1X
  443.      &           + DFRUT2(2) * CT2X)
  444.             UXUX.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  445.             UXUX.AM(IFAC,1,2) = FUNCEL / VOLD
  446. C
  447.             FUNCEL = SURF * (DFRUN(3) * CNX + DFRUT1(3) * CT1X
  448.      &           + DFRUT2(3) * CT2X)
  449.             UXUY.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  450.             UXUY.AM(IFAC,1,2) = FUNCEL / VOLD
  451. C
  452.             FUNCEL = SURF * (DFRUN(4) * CNX + DFRUT1(4) * CT1X
  453.      &           + DFRUT2(4) * CT2X)
  454.             UXUZ.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  455.             UXUZ.AM(IFAC,1,2) = FUNCEL / VOLD
  456. C
  457.             FUNCEL = SURF * (DFRUN(5) * CNX + DFRUT1(5) * CT1X
  458.      &           + DFRUT2(5) * CT2X)
  459.             UXP.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  460.             UXP.AM(IFAC,1,2) = FUNCEL / VOLD
  461. C
  462. C********** Dual RUYN
  463. C
  464.             FUNCEL = SURF * (DFRUN(1) * CNY + DFRUT1(1) * CT1Y
  465.      &           + DFRUT2(1) * CT2Y )
  466.             UYR.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  467.             UYR.AM(IFAC,1,2) = FUNCEL / VOLD
  468. C
  469.             FUNCEL = SURF * (DFRUN(2) * CNY + DFRUT1(2) * CT1Y
  470.      &           + DFRUT2(2) * CT2Y)
  471.             UYUX.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  472.             UYUX.AM(IFAC,1,2) = FUNCEL / VOLD
  473. C
  474.             FUNCEL = SURF * (DFRUN(3) * CNY + DFRUT1(3) * CT1Y
  475.      &           + DFRUT2(3) * CT2Y)
  476.             UYUY.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  477.             UYUY.AM(IFAC,1,2) = FUNCEL / VOLD
  478. C
  479.             FUNCEL = SURF * (DFRUN(4) * CNY + DFRUT1(4) * CT1Y
  480.      &           + DFRUT2(4) * CT2Y)
  481.             UYUZ.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  482.             UYUZ.AM(IFAC,1,2) = FUNCEL / VOLD
  483. C
  484.             FUNCEL = SURF * (DFRUN(5) * CNY + DFRUT1(5) * CT1Y
  485.      &           + DFRUT2(5) * CT2Y)
  486.             UYP.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  487.             UYP.AM(IFAC,1,2) = FUNCEL / VOLD
  488. C
  489. C********** Dual RUZN
  490. C
  491.             FUNCEL = SURF * (DFRUN(1) * CNZ + DFRUT1(1) * CT1Z
  492.      &           + DFRUT2(1) * CT2Z )
  493.             UZR.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  494.             UZR.AM(IFAC,1,2) = FUNCEL / VOLD
  495. C
  496.             FUNCEL = SURF * (DFRUN(2) * CNZ + DFRUT1(2) * CT1Z
  497.      &           + DFRUT2(2) * CT2Z)
  498.             UZUX.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  499.             UZUX.AM(IFAC,1,2) = FUNCEL / VOLD
  500. C
  501.             FUNCEL = SURF * (DFRUN(3) * CNZ + DFRUT1(3) * CT1Z
  502.      &           + DFRUT2(3) * CT2Z)
  503.             UZUY.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  504.             UZUY.AM(IFAC,1,2) = FUNCEL / VOLD
  505. C
  506.             FUNCEL = SURF * (DFRUN(4) * CNZ + DFRUT1(4) * CT1Z
  507.      &           + DFRUT2(4) * CT2Z)
  508.             UZUZ.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  509.             UZUZ.AM(IFAC,1,2) = FUNCEL / VOLD
  510. C
  511.             FUNCEL = SURF * (DFRUN(5) * CNZ + DFRUT1(5) * CT1Z
  512.      &           + DFRUT2(5) * CT2Z)
  513.             UZP.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  514.             UZP.AM(IFAC,1,2) = FUNCEL / VOLD
  515. C
  516. C********** Dual RETN
  517. C
  518.             FUNCEL = SURF * DFRET(1)
  519.             RETR.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  520.             RETR.AM(IFAC,1,2) = FUNCEL / VOLD
  521. C
  522.             FUNCEL = SURF * DFRET(2)
  523.             RETUX.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  524.             RETUX.AM(IFAC,1,2) = FUNCEL / VOLD
  525. C
  526.             FUNCEL = SURF * DFRET(3)
  527.             RETUY.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  528.             RETUY.AM(IFAC,1,2) = FUNCEL / VOLD
  529. C
  530.             FUNCEL = SURF * DFRET(4)
  531.             RETUZ.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  532.             RETUZ.AM(IFAC,1,2) = FUNCEL / VOLD
  533. C
  534.             FUNCEL = SURF * DFRET(5)
  535.             RETP.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  536.             RETP.AM(IFAC,1,2) = FUNCEL / VOLD
  537. C
  538. C
  539. C********** La contribution de D
  540. C
  541.             CNX = -1.0D0 * CNX
  542.             CNY = -1.0D0 * CNY
  543.             CNZ = -1.0D0 * CNZ
  544.             CT1X = -1.0D0 * CT1X
  545.             CT1Y = -1.0D0 * CT1Y
  546.             CT1Z = -1.0D0 * CT1Z
  547.             CT2X = -1.0D0 * CT2X
  548.             CT2Y = -1.0D0 * CT2Y
  549.             CT2Z = -1.0D0 * CT2Z
  550. C
  551.             CALL VLHJ5(ROD,UXD,UYD,UZD,PD,RETD,GAMD,CNX,CNY,CNZ,
  552.      &           CT1X,CT1Y,CT1Z,CT2X,CT2Y,CT2Z,
  553.      &           DFRO,DFRUN,DFRUT1,DFRUT2,DFRET)
  554. C
  555. C
  556. C********** Dual RN
  557. C
  558.             FUNCEL = SURF * DFRO(1)
  559.             RR.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  560.             RR.AM(IFAC,2,1) = FUNCEL / VOLG
  561. C
  562.             FUNCEL = SURF * DFRO(2)
  563.             RUX.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  564.             RUX.AM(IFAC,2,1) = FUNCEL / VOLG
  565. C
  566.             FUNCEL = SURF * DFRO(3)
  567.             RUY.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  568.             RUY.AM(IFAC,2,1) = FUNCEL / VOLG
  569. C
  570.             FUNCEL = SURF * DFRO(4)
  571.             RUZ.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  572.             RUZ.AM(IFAC,2,1) = FUNCEL / VOLG
  573. C
  574.             FUNCEL = SURF * DFRO(5)
  575.             RP.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  576.             RP.AM(IFAC,2,1) = FUNCEL / VOLG
  577. C
  578. C********** Dual RUXN
  579. C
  580.             FUNCEL = SURF * (DFRUN(1) * CNX + DFRUT1(1) * CT1X
  581.      &           + DFRUT2(1) * CT2X)
  582.             UXR.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  583.             UXR.AM(IFAC,2,1) = FUNCEL / VOLG
  584. C
  585.             FUNCEL = SURF * (DFRUN(2) * CNX + DFRUT1(2) * CT1X
  586.      &           + DFRUT2(2) * CT2X)
  587.             UXUX.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  588.             UXUX.AM(IFAC,2,1) = FUNCEL / VOLG
  589. C
  590.             FUNCEL = SURF * (DFRUN(3) * CNX + DFRUT1(3) * CT1X
  591.      &           + DFRUT2(3) * CT2X)
  592.             UXUY.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  593.             UXUY.AM(IFAC,2,1) = FUNCEL / VOLG
  594. C
  595.             FUNCEL = SURF * (DFRUN(4) * CNX + DFRUT1(4) * CT1X
  596.      &           + DFRUT2(4) * CT2X)
  597.             UXUZ.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  598.             UXUZ.AM(IFAC,2,1) = FUNCEL / VOLG
  599. C
  600.             FUNCEL = SURF * (DFRUN(5) * CNX + DFRUT1(5) * CT1X
  601.      &           + DFRUT2(5) * CT2X)
  602.             UXP.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  603.             UXP.AM(IFAC,2,1) = FUNCEL / VOLG
  604. C
  605. C********** Dual RUYN
  606. C
  607.             FUNCEL = SURF * (DFRUN(1) * CNY + DFRUT1(1) * CT1Y
  608.      &           + DFRUT2(1) * CT2Y )
  609.             UYR.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  610.             UYR.AM(IFAC,2,1) = FUNCEL / VOLG
  611. C
  612.             FUNCEL = SURF * (DFRUN(2) * CNY + DFRUT1(2) * CT1Y
  613.      &           + DFRUT2(2) * CT2Y)
  614.             UYUX.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  615.             UYUX.AM(IFAC,2,1) = FUNCEL / VOLG
  616. C
  617.             FUNCEL = SURF * (DFRUN(3) * CNY + DFRUT1(3) * CT1Y
  618.      &           + DFRUT2(3) * CT2Y)
  619.             UYUY.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  620.             UYUY.AM(IFAC,2,1) = FUNCEL / VOLG
  621. C
  622.             FUNCEL = SURF * (DFRUN(4) * CNY + DFRUT1(4) * CT1Y
  623.      &           + DFRUT2(4) * CT2Y)
  624.             UYUZ.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  625.             UYUZ.AM(IFAC,2,1) = FUNCEL / VOLG
  626. C
  627.             FUNCEL = SURF * (DFRUN(5) * CNY + DFRUT1(5) * CT1Y
  628.      &           + DFRUT2(5) * CT2Y)
  629.             UYP.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  630.             UYP.AM(IFAC,2,1) = FUNCEL / VOLG
  631. C
  632. C********** Dual RUZN
  633. C
  634.             FUNCEL = SURF * (DFRUN(1) * CNZ + DFRUT1(1) * CT1Z
  635.      &           + DFRUT2(1) * CT2Z )
  636.             UZR.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  637.             UZR.AM(IFAC,2,1) = FUNCEL / VOLG
  638. C
  639.             FUNCEL = SURF * (DFRUN(2) * CNZ + DFRUT1(2) * CT1Z
  640.      &           + DFRUT2(2) * CT2Z)
  641.             UZUX.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  642.             UZUX.AM(IFAC,2,1) = FUNCEL / VOLG
  643. C
  644.             FUNCEL = SURF * (DFRUN(3) * CNZ + DFRUT1(3) * CT1Z
  645.      &           + DFRUT2(3) * CT2Z)
  646.             UZUY.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  647.             UZUY.AM(IFAC,2,1) = FUNCEL / VOLG
  648. C
  649.             FUNCEL = SURF * (DFRUN(4) * CNZ + DFRUT1(4) * CT1Z
  650.      &           + DFRUT2(4) * CT2Z)
  651.             UZUZ.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  652.             UZUZ.AM(IFAC,2,1) = FUNCEL / VOLG
  653. C
  654.             FUNCEL = SURF * (DFRUN(5) * CNZ + DFRUT1(5) * CT1Z
  655.      &           + DFRUT2(5) * CT2Z)
  656.             UZP.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  657.             UZP.AM(IFAC,2,1) = FUNCEL / VOLG
  658. C
  659. C********** Dual RETN
  660. C
  661.             FUNCEL = SURF * DFRET(1)
  662.             RETR.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  663.             RETR.AM(IFAC,2,1) = FUNCEL / VOLG
  664. C
  665.             FUNCEL = SURF * DFRET(2)
  666.             RETUX.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  667.             RETUX.AM(IFAC,2,1) = FUNCEL / VOLG
  668. C
  669.             FUNCEL = SURF * DFRET(3)
  670.             RETUY.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  671.             RETUY.AM(IFAC,2,1) = FUNCEL / VOLG
  672. C
  673.             FUNCEL = SURF * DFRET(4)
  674.             RETUZ.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  675.             RETUZ.AM(IFAC,2,1) = FUNCEL / VOLG
  676. C
  677.             FUNCEL = SURF * DFRET(5)
  678.             RETP.AM(IFAC,2,2) = -1.0D0 * FUNCEL / VOLD
  679.             RETP.AM(IFAC,2,1) = FUNCEL / VOLG
  680. C
  681.          ELSE
  682. C
  683. C********** Murs (NGCG = NGCD)
  684. C
  685. C
  686. C********** Les MELEMEs
  687. C
  688.             MELEDU.NUM(1,IFAC) = NGCG
  689.             MELEDU.NUM(2,IFAC) = NGCD
  690.             NLCG = MLENTC.LECT(NGCG)
  691. C
  692.             ROG = MPRN.VPOCHA(NLCG,1)
  693.             PG = MPPN.VPOCHA(NLCG,1)
  694.             UXG = MPUN.VPOCHA(NLCG,1)
  695.             UYG = MPUN.VPOCHA(NLCG,2)
  696.             UZG = MPUN.VPOCHA(NLCG,3)
  697.             GAMG = MPGAMN.VPOCHA(NLCG,1)
  698.             VOLG = MPVOLU.VPOCHA(NLCG,1)
  699. C
  700. C********** La normale sortante
  701. C
  702.             SURF = MPOVSU.VPOCHA(NLCF,1)
  703.             CNX = MPNORM.VPOCHA(NLCF,7)
  704.             CNY = MPNORM.VPOCHA(NLCF,8)
  705.             CNZ = MPNORM.VPOCHA(NLCF,9)
  706. C
  707.             CALL VLHJ7(ROG,UXG,UYG,UZG,PG,GAMG,CNX,CNY,CNZ,
  708.      &           DFRUN)
  709. C
  710. C********** Dual RN
  711. C
  712.             RR.AM(IFAC,1,1) = 0.0D0
  713.             RR.AM(IFAC,1,2) = 0.0D0
  714. C
  715.             RUX.AM(IFAC,1,1) = 0.0D0
  716.             RUX.AM(IFAC,1,2) = 0.0D0
  717. C
  718.             RUY.AM(IFAC,1,1) = 0.0D0
  719.             RUY.AM(IFAC,1,2) = 0.0D0
  720. C
  721.             RUZ.AM(IFAC,1,1) = 0.0D0
  722.             RUZ.AM(IFAC,1,2) = 0.0D0
  723. C
  724.             RP.AM(IFAC,1,1) = 0.0D0
  725.             RP.AM(IFAC,1,2) = 0.0D0
  726. C
  727. C********** Dual RUXN
  728. C
  729.             FUNCEL = SURF * DFRUN(1) * CNX
  730.             UXR.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  731.             UXR.AM(IFAC,1,2) = 0.0D0
  732. C
  733.             FUNCEL = SURF * DFRUN(2) * CNX
  734.             UXUX.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  735.             UXUX.AM(IFAC,1,2) = 0.0D0
  736. C
  737.             FUNCEL = SURF * DFRUN(3) * CNX
  738.             UXUY.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  739.             UXUY.AM(IFAC,1,2) = 0.0D0
  740. C
  741.             FUNCEL = SURF * DFRUN(4) * CNX
  742.             UXUZ.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  743.             UXUZ.AM(IFAC,1,2) = 0.0D0
  744. C
  745.             FUNCEL = SURF * DFRUN(5) * CNX
  746.             UXP.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  747.             UXP.AM(IFAC,1,2) = 0.0D0
  748. C
  749. C********** Dual RUYN
  750. C
  751.             FUNCEL = SURF * DFRUN(1) * CNY
  752.             UYR.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  753.             UYR.AM(IFAC,1,2) = 0.0D0
  754. C
  755.             FUNCEL = SURF * DFRUN(2) * CNY
  756.             UYUX.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  757.             UYUX.AM(IFAC,1,2) = 0.0D0
  758. C
  759.             FUNCEL = SURF * DFRUN(3) * CNY
  760.             UYUY.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  761.             UYUY.AM(IFAC,1,2) = 0.0D0
  762. C
  763.             FUNCEL = SURF * DFRUN(4) * CNY
  764.             UYUZ.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  765.             UYUZ.AM(IFAC,1,2) = 0.0D0
  766. C
  767.             FUNCEL = SURF * DFRUN(5) * CNY
  768.             UYP.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  769.             UYP.AM(IFAC,1,2) = 0.0D0
  770. C
  771. C********** Dual RUZN
  772. C
  773.             FUNCEL = SURF * DFRUN(1) * CNZ
  774.             UZR.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  775.             UZR.AM(IFAC,1,2) = 0.0D0
  776. C
  777.             FUNCEL = SURF * DFRUN(2) * CNZ
  778.             UZUX.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  779.             UZUX.AM(IFAC,1,2) = 0.0D0
  780. C
  781.             FUNCEL = SURF * DFRUN(3) * CNZ
  782.             UZUY.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  783.             UZUY.AM(IFAC,1,2) = 0.0D0
  784. C
  785.             FUNCEL = SURF * DFRUN(4) * CNZ
  786.             UZUZ.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  787.             UZUZ.AM(IFAC,1,2) = 0.0D0
  788. C
  789.             FUNCEL = SURF * DFRUN(5) * CNZ
  790.             UZP.AM(IFAC,1,1) = -1.0D0 * FUNCEL / VOLG
  791.             UZP.AM(IFAC,1,2) = 0.0D0
  792. C
  793. C********** Dual RETN
  794. C
  795.             RETR.AM(IFAC,1,1) = 0.0D0
  796.             RETR.AM(IFAC,1,2) = 0.0D0
  797. C
  798.             RETUX.AM(IFAC,1,1) = 0.0D0
  799.             RETUX.AM(IFAC,1,2) = 0.0D0
  800. C
  801.             RETUY.AM(IFAC,1,1) = 0.0D0
  802.             RETUY.AM(IFAC,1,2) = 0.0D0
  803. C
  804.             RETUZ.AM(IFAC,1,1) = 0.0D0
  805.             RETUZ.AM(IFAC,1,2) = 0.0D0
  806. C
  807.             RETP.AM(IFAC,1,1) = 0.0D0
  808.             RETP.AM(IFAC,1,2) = 0.0D0
  809. C
  810. C********** Dual RN
  811. C
  812.             RR.AM(IFAC,2,2) = 0.0D0
  813.             RR.AM(IFAC,2,1) = 0.0D0
  814. C
  815.             RUX.AM(IFAC,2,2) = 0.0D0
  816.             RUX.AM(IFAC,2,1) = 0.0D0
  817. C
  818.             RUY.AM(IFAC,2,2) = 0.0D0
  819.             RUY.AM(IFAC,2,1) = 0.0D0
  820. C
  821.             RUZ.AM(IFAC,2,2) = 0.0D0
  822.             RUZ.AM(IFAC,2,1) = 0.0D0
  823. C
  824.             RP.AM(IFAC,2,2) = 0.0D0
  825.             RP.AM(IFAC,2,1) = 0.0D0
  826. C
  827. C********** Dual RUXN
  828. C
  829.             UXR.AM(IFAC,2,2) = 0.0D0
  830.             UXR.AM(IFAC,2,1) = 0.0D0
  831. C
  832.             UXUX.AM(IFAC,2,2) = 0.0D0
  833.             UXUX.AM(IFAC,2,1) = 0.0D0
  834. C
  835.             UXUY.AM(IFAC,2,2) = 0.0D0
  836.             UXUY.AM(IFAC,2,1) = 0.0D0
  837. C
  838.             UXUZ.AM(IFAC,2,2) = 0.0D0
  839.             UXUZ.AM(IFAC,2,1) = 0.0D0
  840. C
  841.             UXP.AM(IFAC,2,2) = 0.0D0
  842.             UXP.AM(IFAC,2,1) = 0.0D0
  843. C
  844. C********** Dual RUYN
  845. C
  846.             UYR.AM(IFAC,2,2) = 0.0D0
  847.             UYR.AM(IFAC,2,1) = 0.0D0
  848. C
  849.             UYUX.AM(IFAC,2,2) = 0.0D0
  850.             UYUX.AM(IFAC,2,1) = 0.0D0
  851. C
  852.             UYUY.AM(IFAC,2,2) = 0.0D0
  853.             UYUY.AM(IFAC,2,1) = 0.0D0
  854. C
  855.             UYUZ.AM(IFAC,2,2) = 0.0D0
  856.             UYUZ.AM(IFAC,2,1) = 0.0D0
  857. C
  858.             UYP.AM(IFAC,2,2) = 0.0D0
  859.             UYP.AM(IFAC,2,1) = 0.0D0
  860. C
  861. C********** Dual RUZN
  862. C
  863.             UZR.AM(IFAC,2,2) = 0.0D0
  864.             UZR.AM(IFAC,2,1) = 0.0D0
  865. C
  866.             UZUX.AM(IFAC,2,2) = 0.0D0
  867.             UZUX.AM(IFAC,2,1) = 0.0D0
  868. C
  869.             UZUY.AM(IFAC,2,2) = 0.0D0
  870.             UZUY.AM(IFAC,2,1) = 0.0D0
  871. C
  872.             UZUZ.AM(IFAC,2,2) = 0.0D0
  873.             UZUZ.AM(IFAC,2,1) = 0.0D0
  874. C
  875.             UZP.AM(IFAC,2,2) = 0.0D0
  876.             UZP.AM(IFAC,2,1) = 0.0D0
  877. C
  878. C********** Dual RETN
  879. C
  880.             RETR.AM(IFAC,2,2) = 0.0D0
  881.             RETR.AM(IFAC,2,1) = 0.0D0
  882. C
  883.             RETUX.AM(IFAC,2,2) = 0.0D0
  884.             RETUX.AM(IFAC,2,1) = 0.0D0
  885. C
  886.             RETUY.AM(IFAC,2,2) = 0.0D0
  887.             RETUY.AM(IFAC,2,1) = 0.0D0
  888. C
  889.             RETUZ.AM(IFAC,2,2) = 0.0D0
  890.             RETUZ.AM(IFAC,2,1) = 0.0D0
  891. C
  892.             RETP.AM(IFAC,2,2) = 0.0D0
  893.             RETP.AM(IFAC,2,1) = 0.0D0
  894. C
  895.          ENDIF
  896.       ENDDO
  897. C
  898.       SEGDES MELEMC
  899.       SEGDES MELEFE
  900.       SEGDES MELEMF
  901. C
  902.       SEGDES MPOVSU
  903.       SEGDES MPVOLU
  904.       SEGDES MPNORM
  905. C
  906.       SEGDES MPRN
  907.       SEGDES MPPN
  908.       SEGDES MPUN
  909.       SEGDES MPGAMN
  910. C
  911.       SEGDES MELEDU
  912.       SEGDES MATRIK
  913.       SEGDES IMATRI
  914. C
  915.       SEGDES   RR  , RUX , RUY , RUZ, RP ,
  916.      &         UXR , UXUX , UXUY , UXUZ, UXP ,
  917.      &         UYR , UYUX , UYUY , UYUZ, UYP ,
  918.      &         UZR , UZUX , UZUY , UZUZ, UYP ,
  919.      &         RETR , RETUX , RETUY , RETUZ, RETP
  920.  
  921.       SEGSUP MLENTC
  922.       SEGSUP MLENTF
  923.       SEGDES MLMINC
  924.       SEGSUP MLELIM
  925.  
  926.  9999 CONTINUE
  927.       RETURN
  928.       END
  929.  
  930.  
  931.  
  932.  
  933.  
  934.  
  935.  
  936.  
  937.  
  938.  
  939.  
  940.  

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