Télécharger la plaquette Cast3M


fr en

Télécharger cli153.eso

Retour à la liste

Numérotation des lignes : 

cli153

  1. C CLI153    SOURCE    CB215821  20/11/25    13:20:11     10792          
  2.       SUBROUTINE CLI153(MELEMF,MELEMC,MELECB,MELEFC,MELRES,INORM,ICHPVO,
  3.      &     ICHPSU,IROC,IVITC,IPC,IGAMC,ICHLIM,ILIINC,ILIINP,IJAC,IJACO)
  4. C************************************************************************
  5. C
  6. C PROJET            :  CASTEM 2000
  7. C
  8. C NOM               :  CLI153
  9. C
  10. C DESCRIPTION       :  Subroutine appellée par CLIM11
  11. C
  12. C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  13. C
  14. C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/LTMF
  15. C
  16. C************************************************************************
  17. C
  18. C APPELES (Calcul)  :
  19. C
  20. C************************************************************************
  21. C
  22. C HISTORIQUE (Anomalies et modifications éventuelles)
  23. C
  24. C HISTORIQUE :
  25. C
  26. C************************************************************************
  27. C
  28.       IMPLICIT INTEGER(I-N)
  29. -INC PPARAM
  30. -INC CCOPTIO
  31. -INC SMLMOTS
  32. -INC SMELEME
  33.       POINTEUR MELEFC.MELEME
  34. -INC SMLENTI
  35.       POINTEUR MLEMC.MLENTI, MLEMCB.MLENTI,MLEMF.MLENTI
  36. -INC SMCHPOI
  37.       POINTEUR MPNORM.MPOVAL, MPVOL.MPOVAL, MPSURF.MPOVAL, MPRC.MPOVAL,
  38.      &     MPVC.MPOVAL, MPPC.MPOVAL, MPGAMC.MPOVAL, MPLIM.MPOVAL
  40.       POINTEUR RR.IZAFM, RUX.IZAFM, RUY.IZAFM, RRET.IZAFM,
  41.      &         UXR.IZAFM, UXUX.IZAFM, UXUY.IZAFM, UXRET.IZAFM,
  42.      &         UYR.IZAFM, UYUX.IZAFM, UYUY.IZAFM, UYRET.IZAFM,
  43.      &         RETR.IZAFM, RETUX.IZAFM, RETUY.IZAFM, RETRET.IZAFM,
  44.      &         RUZ.IZAFM, UXUZ.IZAFM, UYUZ.IZAFM, RETUZ.IZAFM,
  45.      &         UZR.IZAFM, UZUX.IZAFM, UZUY.IZAFM, UZRET.IZAFM,
  46.      &         UZUZ.IZAFM
  47. C
  48. C**** Variables de COOPTIO
  49. C
  50. C      INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  51. C     &        ,IOTER,   IOLEC,  IOIMP,     IOCAR,  IOACQ
  52. C     &        ,IOPER,   IOSGB,  IOGRA,     IOSAU,  IORES
  53. C     &        ,IECHO,   IIMPI,  IOSPI
  54. C     &        ,IDIM,    IFICLE, IPREFI
  55. C     &        ,MCOORD
  56. C     &        ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  57. C     &        ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  58. C     &        ,NORINC,NORVAL,NORIND,NORVAD
  59. C     &        ,NUCROU, IPSAUV
  60. C
  61.       INTEGER MELEMF,MELEMC,MELECB,INORM,ICHPVO,ICHPSU, IROC,IVITC,IPC
  62.      &     ,IGAMC,ICHLIM,ICEL,NFAC,IFAC,MELRES,IJACO
  63.      &     ,NGF,NGC,NLF,NLC,NLCB
  64.      &     ,ILIINC,ILIINP,IJAC
  65.      &     ,MP, NBEL, NBME, NBSOUS, NKID, NKMT, NMATRI, NP, NRIGE
  66.       REAL*8 VOLU,SURF,RC,PC,UXC,UYC,GAMC,CNX,CNY,CTX,CTY,PF
  67.      &     ,CNZ,CTZ,CT2X,CT2Y,CT2Z,UZC
  68.      &     ,VECG(5),VECD(5),DFDUG(5,5),DFDUD(5,5)
  69.      &     ,NVEC(3),TVEC(3),TVEC2(3)
  70.      &     ,DRORO,DROUX,DROUY,DROP
  71.      &     ,DUXRO,DUXUX,DUXUY,DUXP
  72.      &     ,DUYRO,DUYUX,DUYUY,DUYP
  73.      &     ,DPRO,DPUX,DPUY,DPP
  74.      &     ,DROUZ,DUXUZ,DUYUZ,DPUZ
  75.      &     ,DUZRO,DUZUX,DUZUY,DUZUZ,DUZP
  76.      &     ,COEF1,COEF2,COEF3,GM1,COEF,COEF4
  77.       CHARACTER*(8) TYPE
  78. C
  79. C
  80. C**** KRIPAD pour la correspondance global/local
  81. C
  82.       CALL KRIPAD(MELEMC,MLEMC)
  83. C     SEGINI MLEMC
  84.       CALL KRIPAD(MELECB,MLEMCB)
  85. C     SEGINI MLEMCB
  86.       CALL KRIPAD(MELEMF,MLEMF)
  87. C     SEGINI MLEMF
  88. C
  89. C**** CHPOINTs de la table DOMAINE
  90. C
  91.       CALL LICHT(INORM,MPNORM,TYPE,ICEL)
  92.       CALL LICHT(ICHPVO,MPVOL,TYPE,ICEL)
  93.       CALL LICHT(ICHPSU,MPSURF,TYPE,ICEL)
  94. C
  95. C**** LICHT active les MPOVALs en *MOD
  96. C
  97. C     SEGACT MPNORM*MOD
  98. C     SEGACT MPOVSU*MOD
  99. C     SEGACT MPOVOL*MOD
  100. C
  101. C
  102. C**** CHPOINTs des variables
  103. C
  104.       CALL LICHT(IROC,MPRC,TYPE,ICEL)
  105.       CALL LICHT(IVITC,MPVC,TYPE,ICEL)
  106.       CALL LICHT(IPC,MPPC,TYPE,ICEL)
  107.       CALL LICHT(IGAMC,MPGAMC,TYPE,ICEL)
  108.       CALL LICHT(ICHLIM,MPLIM,TYPE,ICEL)
  109. C
  110. C     SEGACT *MOD
  111. C     SEGACT *MOD
  112. C     SEGACT *MOD
  113. C     SEGACT *MOD
  114. C     SEGACT *MOD
  115. C
  116. C
  117. C**** Boucle sur le face pour le calcul des invariants de
  118. C     Riemann et du flux
  119. C
  120.       SEGACT MELEFC
  121.       NFAC=MELEFC.NUM(/2)
  122. C
  123. C**** Objet MATRIK
  124. C
  125.       NRIGE = 7
  126.       NMATRI = 1
  127.       NKID = 9
  128.       NKMT = 7
  129. C
  130.       SEGINI MATRIK
  131.       IJACO = MATRIK
  132.       MATRIK.IRIGEL(1,1) = MELRES
  133.       MATRIK.IRIGEL(2,1) = MELRES
  134. C
  135. C**** Matrice non symetrique
  136. C
  137.       MATRIK.IRIGEL(7,1) = 2
  138. C
  139.       NBME = 25
  140.       NBSOUS = 1
  141.       SEGINI IMATRI
  142.       IF(IJAC.EQ.1)THEN
  143.          MLMOTS=ILIINC
  144.       ELSEIF(IJAC.EQ.2)THEN
  145.          MLMOTS=ILIINP
  146.       ENDIF
  147.       SEGACT MLMOTS
  148.       MATRIK.IRIGEL(4,1) = IMATRI
  149. C
  150.       IMATRI.LISPRI(1) = MLMOTS.MOTS(1)
  151.       IMATRI.LISPRI(2) = MLMOTS.MOTS(2)
  152.       IMATRI.LISPRI(3) = MLMOTS.MOTS(3)
  153.       IMATRI.LISPRI(4) = MLMOTS.MOTS(4)
  154.       IMATRI.LISPRI(5) = MLMOTS.MOTS(5)
  155. C
  156.       IMATRI.LISPRI(6) = MLMOTS.MOTS(1)
  157.       IMATRI.LISPRI(7) = MLMOTS.MOTS(2)
  158.       IMATRI.LISPRI(8) = MLMOTS.MOTS(3)
  159.       IMATRI.LISPRI(9) = MLMOTS.MOTS(4)
  160.       IMATRI.LISPRI(10) = MLMOTS.MOTS(5)
  161. C
  162.       IMATRI.LISPRI(11) = MLMOTS.MOTS(1)
  163.       IMATRI.LISPRI(12) = MLMOTS.MOTS(2)
  164.       IMATRI.LISPRI(13) = MLMOTS.MOTS(3)
  165.       IMATRI.LISPRI(14) = MLMOTS.MOTS(4)
  166.       IMATRI.LISPRI(15) = MLMOTS.MOTS(5)
  167. C
  168.       IMATRI.LISPRI(16) = MLMOTS.MOTS(1)
  169.       IMATRI.LISPRI(17) = MLMOTS.MOTS(2)
  170.       IMATRI.LISPRI(18) = MLMOTS.MOTS(3)
  171.       IMATRI.LISPRI(19) = MLMOTS.MOTS(4)
  172.       IMATRI.LISPRI(20) = MLMOTS.MOTS(5)
  173. C
  174.       IMATRI.LISPRI(21) = MLMOTS.MOTS(1)
  175.       IMATRI.LISPRI(22) = MLMOTS.MOTS(2)
  176.       IMATRI.LISPRI(23) = MLMOTS.MOTS(3)
  177.       IMATRI.LISPRI(24) = MLMOTS.MOTS(4)
  178.       IMATRI.LISPRI(25) = MLMOTS.MOTS(5)
  179. C
  180.       SEGDES MLMOTS
  181.       MLMOTS=ILIINC
  182.       SEGACT MLMOTS
  183. C
  184. C
  185.       IMATRI.LISDUA(1) = MLMOTS.MOTS(1)
  186.       IMATRI.LISDUA(2) = MLMOTS.MOTS(1)
  187.       IMATRI.LISDUA(3) = MLMOTS.MOTS(1)
  188.       IMATRI.LISDUA(4) = MLMOTS.MOTS(1)
  189.       IMATRI.LISDUA(5) = MLMOTS.MOTS(1)
  190. C
  191.       IMATRI.LISDUA(6) = MLMOTS.MOTS(2)
  192.       IMATRI.LISDUA(7) = MLMOTS.MOTS(2)
  193.       IMATRI.LISDUA(8) = MLMOTS.MOTS(2)
  194.       IMATRI.LISDUA(9) = MLMOTS.MOTS(2)
  195.       IMATRI.LISDUA(10) = MLMOTS.MOTS(2)
  196. C
  197.       IMATRI.LISDUA(11) = MLMOTS.MOTS(3)
  198.       IMATRI.LISDUA(12) = MLMOTS.MOTS(3)
  199.       IMATRI.LISDUA(13) = MLMOTS.MOTS(3)
  200.       IMATRI.LISDUA(14) = MLMOTS.MOTS(3)
  201.       IMATRI.LISDUA(15) = MLMOTS.MOTS(3)
  202. C
  203.       IMATRI.LISDUA(16) = MLMOTS.MOTS(4)
  204.       IMATRI.LISDUA(17) = MLMOTS.MOTS(4)
  205.       IMATRI.LISDUA(18) = MLMOTS.MOTS(4)
  206.       IMATRI.LISDUA(19) = MLMOTS.MOTS(4)
  207.       IMATRI.LISDUA(20) = MLMOTS.MOTS(4)
  208. C
  209.       IMATRI.LISDUA(21) = MLMOTS.MOTS(5)
  210.       IMATRI.LISDUA(22) = MLMOTS.MOTS(5)
  211.       IMATRI.LISDUA(23) = MLMOTS.MOTS(5)
  212.       IMATRI.LISDUA(24) = MLMOTS.MOTS(5)
  213.       IMATRI.LISDUA(25) = MLMOTS.MOTS(5)
  214. C
  215.       SEGDES MLMOTS
  216.       NBEL = NFAC
  217.       NBSOUS = 1
  218.       NP = 1
  219.       MP = 1
  220.       SEGINI   RR  , RUX , RUY , RUZ ,RRET ,
  221.      &         UXR , UXUX , UXUY , UXUZ ,UXRET ,
  222.      &         UYR , UYUX , UYUY , UYUZ , UYRET ,
  223.      &         UZR , UZUX , UZUY , UZUZ , UZRET ,
  224.      &         RETR , RETUX , RETUY, RETUZ , RETRET
  225. C
  226. C**** Duale = IMATRI.LISDUA(1) = 'RN'
  227. C     Primale = IMATRI.LISPRI(1) = 'RN'
  228. C     -> IMATRI.LIZAFM(1,1) = RR
  229. C
  230.       IMATRI.LIZAFM(1,1) = RR
  231.       IMATRI.LIZAFM(1,2) = RUX
  232.       IMATRI.LIZAFM(1,3) = RUY
  233.       IMATRI.LIZAFM(1,4) = RUZ
  234.       IMATRI.LIZAFM(1,5) = RRET
  235. C
  236.       IMATRI.LIZAFM(1,6) = UXR
  237.       IMATRI.LIZAFM(1,7) = UXUX
  238.       IMATRI.LIZAFM(1,8) = UXUY
  239.       IMATRI.LIZAFM(1,9) = UXUZ
  240.       IMATRI.LIZAFM(1,10) = UXRET
  241. C
  242.       IMATRI.LIZAFM(1,11) = UYR
  243.       IMATRI.LIZAFM(1,12) = UYUX
  244.       IMATRI.LIZAFM(1,13) = UYUY
  245.       IMATRI.LIZAFM(1,14) = UYUZ
  246.       IMATRI.LIZAFM(1,15) = UYRET
  247. C
  248.       IMATRI.LIZAFM(1,16) = UZR
  249.       IMATRI.LIZAFM(1,17) = UZUX
  250.       IMATRI.LIZAFM(1,18) = UZUY
  251.       IMATRI.LIZAFM(1,19) = UZUZ
  252.       IMATRI.LIZAFM(1,20) = UZRET
  253. C
  254.       IMATRI.LIZAFM(1,21) = RETR
  255.       IMATRI.LIZAFM(1,22) = RETUX
  256.       IMATRI.LIZAFM(1,23) = RETUY
  257.       IMATRI.LIZAFM(1,24) = RETUZ
  258.       IMATRI.LIZAFM(1,25) = RETRET
  259. C
  260.       SEGDES MATRIK
  261.       SEGDES IMATRI
  262. C
  263. C**** Fin definition MATRIK
  264. C
  265.       DO IFAC=1,NFAC,1
  266.          NGF=MELEFC.NUM(1,IFAC)
  267.          NGC=MELEFC.NUM(2,IFAC)
  268.          NLF=MLEMF.LECT(NGF)
  269.          NLC=MLEMC.LECT(NGC)
  270.          NLCB=MLEMCB.LECT(NGF)
  271.          VOLU=MPVOL.VPOCHA(NLC,1)
  272.          SURF=MPSURF.VPOCHA(NLF,1)
  273. C        In CASTEM les normales sont sortantes
  274.          CNX=MPNORM.VPOCHA(NLF,1)
  275.          CNY=MPNORM.VPOCHA(NLF,2)
  276.          CNZ=MPNORM.VPOCHA(NLF,3)
  277.          CTX=MPNORM.VPOCHA(NLF,4)
  278.          CTY=MPNORM.VPOCHA(NLF,5)
  279.          CTZ=MPNORM.VPOCHA(NLF,6)
  280.          CT2X=MPNORM.VPOCHA(NLF,7)
  281.          CT2Y=MPNORM.VPOCHA(NLF,8)
  282.          CT2Z=MPNORM.VPOCHA(NLF,9)
  283. C        Variables au centre
  284.          RC=MPRC.VPOCHA(NLC,1)
  285.          PC=MPPC.VPOCHA(NLC,1)
  286.          UXC=MPVC.VPOCHA(NLC,1)
  287.          UYC=MPVC.VPOCHA(NLC,2)
  288.          UZC=MPVC.VPOCHA(NLC,3)
  289.          GAMC=MPGAMC.VPOCHA(NLC,1)
  290. C        Variables à la face
  291.          PF=MPLIM.VPOCHA(NLCB,1)
  292. C
  293. C******* Jacobian with respect to primitive variables
  294. C
  295.          VECG(1)=RC
  296.          VECG(2)=UXC
  297.          VECG(3)=UYC
  298.          VECG(4)=UZC
  299.          VECG(5)=PC
  300. C
  301.          VECD(1)=RC
  302.          VECD(2)=UXC
  303.          VECD(3)=UYC
  304.          VECD(4)=UZC
  305.          VECD(5)=PF
  306. C
  307.          NVEC(1)=CNX
  308.          NVEC(2)=CNY
  309.          NVEC(3)=CNZ
  310.          TVEC(1)=CTX
  311.          TVEC(2)=CTY
  312.          TVEC(3)=CTZ
  313.          TVEC2(1)=CT2X
  314.          TVEC2(2)=CT2Y
  315.          TVEC2(3)=CT2Z
  316. C
  317.          CALL CONJP4(DFDUG,DFDUD,VECG,VECD,NVEC,TVEC,TVEC2,GAMC)
  318. C
  319.          COEF=-1*SURF/VOLU
  320.          DRORO=(DFDUG(1,1)+DFDUD(1,1))*COEF
  321.          DROUX=(DFDUG(1,2)+DFDUD(1,2))*COEF
  322.          DROUY=(DFDUG(1,3)+DFDUD(1,3))*COEF
  323.          DROUZ=(DFDUG(1,4)+DFDUD(1,4))*COEF
  324.          DROP=DFDUG(1,5)*COEF
  325. C
  326.          DPRO=(DFDUG(5,1)+DFDUD(5,1))*COEF
  327.          DPUX=(DFDUG(5,2)+DFDUD(5,2))*COEF
  328.          DPUY=(DFDUG(5,3)+DFDUD(5,3))*COEF
  329.          DPUZ=(DFDUG(5,4)+DFDUD(5,4))*COEF
  330.          DPP=DFDUG(5,5)*COEF
  331. C
  332.          DUXRO=(DFDUG(2,1)+DFDUD(2,1))*COEF
  333.          DUXUX=(DFDUG(2,2)+DFDUD(2,2))*COEF
  334.          DUXUY=(DFDUG(2,3)+DFDUD(2,3))*COEF
  335.          DUXUZ=(DFDUG(2,4)+DFDUD(2,4))*COEF
  336.          DUXP =DFDUG(2,5)*COEF
  337. C
  338.          DUYRO=(DFDUG(3,1)+DFDUD(3,1))*COEF
  339.          DUYUX=(DFDUG(3,2)+DFDUD(3,2))*COEF
  340.          DUYUY=(DFDUG(3,3)+DFDUD(3,3))*COEF
  341.          DUYUZ=(DFDUG(3,4)+DFDUD(3,4))*COEF
  342.          DUYP =DFDUG(3,5)*COEF
  343. C
  344.          DUZRO=(DFDUG(4,1)+DFDUD(4,1))*COEF
  345.          DUZUX=(DFDUG(4,2)+DFDUD(4,2))*COEF
  346.          DUZUY=(DFDUG(4,3)+DFDUD(4,3))*COEF
  347.          DUZUZ=(DFDUG(4,4)+DFDUD(4,4))*COEF
  348.          DUZP =DFDUG(4,5)*COEF
  349. C
  350. C******* Jacobian with respect to conservative variables
  351. C
  352.          IF(IJAC.EQ.1)THEN
  353.             GM1=GAMC-1.0D0
  354. C
  355.             COEF1=-1.0D0*UXC/RC
  356.             COEF2=-1.0D0*UYC/RC
  357.             COEF3=-1.0D0*UZC/RC
  358.             COEF4=0.5D0*GM1*((UXC*UXC)+(UYC*UYC)+(UZC*UZC))
  359. C
  360.             RR.AM(IFAC,1,1)=DRORO+(DROUX*COEF1)+(DROUY*COEF2)+(DROUZ
  361.      $           *COEF3)+(DROP*COEF4)
  362.             RUX.AM(IFAC,1,1)=(DROUX/RC)-((UXC*GM1)*DROP)
  363.             RUY.AM(IFAC,1,1)=(DROUY/RC)-((UYC*GM1)*DROP)
  364.             RUZ.AM(IFAC,1,1)=(DROUZ/RC)-((UZC*GM1)*DROP)
  365.             RRET.AM(IFAC,1,1)=GM1*DROP
  366. C
  367.             UXR.AM(IFAC,1,1)=DUXRO+(DUXUX*COEF1)+(DUXUY*COEF2)+(DUXUZ
  368.      $           *COEF3)+(DUXP*COEF4)
  369.             UXUX.AM(IFAC,1,1)=(DUXUX/RC)-((UXC*GM1)*DUXP)
  370.             UXUY.AM(IFAC,1,1)=(DUXUY/RC)-((UYC*GM1)*DUXP)
  371.             UXUZ.AM(IFAC,1,1)=(DUXUZ/RC)-((UZC*GM1)*DUXP)
  372.             UXRET.AM(IFAC,1,1)=GM1*DUXP
  373. C
  374.             UYR.AM(IFAC,1,1)=DUYRO+(DUYUX*COEF1)+(DUYUY*COEF2)+(DUYUZ
  375.      $           *COEF3)+(DUYP*COEF4)
  376.             UYUX.AM(IFAC,1,1)=(DUYUX/RC)-((UXC*GM1)*DUYP)
  377.             UYUY.AM(IFAC,1,1)=(DUYUY/RC)-((UYC*GM1)*DUYP)
  378.             UYUZ.AM(IFAC,1,1)=(DUYUZ/RC)-((UZC*GM1)*DUYP)
  379.             UYRET.AM(IFAC,1,1)=GM1*DUYP
  380. C
  381.             UZR.AM(IFAC,1,1)=DUZRO+(DUZUX*COEF1)+(DUZUY*COEF2)+(DUZUZ
  382.      $           *COEF3)+(DUZP*COEF4)
  383.             UZUX.AM(IFAC,1,1)=(DUZUX/RC)-((UXC*GM1)*DUZP)
  384.             UZUY.AM(IFAC,1,1)=(DUZUY/RC)-((UYC*GM1)*DUZP)
  385.             UZUZ.AM(IFAC,1,1)=(DUZUZ/RC)-((UZC*GM1)*DUZP)
  386.             UZRET.AM(IFAC,1,1)=GM1*DUZP
  387. C
  388.             RETR.AM(IFAC,1,1)=DPRO+(DPUX*COEF1)+(DPUY*COEF2)+(DPUZ
  389.      $           *COEF3)+(DPP*COEF4)
  390.             RETUX.AM(IFAC,1,1)=(DPUX/RC)-((UXC*GM1)*DPP)
  391.             RETUY.AM(IFAC,1,1)=(DPUY/RC)-((UYC*GM1)*DPP)
  392.             RETUZ.AM(IFAC,1,1)=(DPUZ/RC)-((UZC*GM1)*DPP)
  393.             RETRET.AM(IFAC,1,1)=GM1*DPP
  394. C
  395.          ELSEIF(IJAC.EQ.2)THEN
  396.             RR.AM(IFAC,1,1)=DRORO
  397.             RUX.AM(IFAC,1,1)=DROUX
  398.             RUY.AM(IFAC,1,1)=DROUY
  399.             RUZ.AM(IFAC,1,1)=DROUZ
  400.             RRET.AM(IFAC,1,1)=DROP
  401. C
  402.             UXR.AM(IFAC,1,1)=DUXRO
  403.             UXUX.AM(IFAC,1,1)=DUXUX
  404.             UXUY.AM(IFAC,1,1)=DUXUY
  405.             UXUZ.AM(IFAC,1,1)=DUXUZ
  406.             UXRET.AM(IFAC,1,1)=DUXP
  407. C
  408.             UYR.AM(IFAC,1,1)=DUYRO
  409.             UYUX.AM(IFAC,1,1)=DUYUX
  410.             UYUY.AM(IFAC,1,1)=DUYUY
  411.             UYUZ.AM(IFAC,1,1)=DUYUZ
  412.             UYRET.AM(IFAC,1,1)=DUYP
  413. C
  414.             UZR.AM(IFAC,1,1)=DUZRO
  415.             UZUX.AM(IFAC,1,1)=DUZUX
  416.             UZUY.AM(IFAC,1,1)=DUZUY
  417.             UZUZ.AM(IFAC,1,1)=DUZUZ
  418.             UZRET.AM(IFAC,1,1)=DUZP
  419. C
  420.             RETR.AM(IFAC,1,1)=DPRO
  421.             RETUX.AM(IFAC,1,1)=DPUX
  422.             RETUY.AM(IFAC,1,1)=DPUY
  423.             RETUZ.AM(IFAC,1,1)=DPUZ
  424.             RETRET.AM(IFAC,1,1)=DPP
  425.          ENDIF
  426.       ENDDO
  427. C
  428.       SEGDES MELEFC
  429. C
  430.       SEGSUP MLEMC
  431.       SEGSUP MLEMCB
  432.       SEGSUP MLEMF
  433. C
  434.       SEGDES MPNORM
  435.       SEGDES MPVOL
  436.       SEGDES MPSURF
  437.       SEGDES MPRC
  438.       SEGDES MPPC
  439.       SEGDES MPVC
  440.       SEGDES MPGAMC
  441.       SEGDES MPLIM
  442. C
  443.       SEGDES   RR  , RUX , RUY , RUZ ,RRET ,
  444.      &         UXR , UXUX , UXUY , UXUZ ,UXRET ,
  445.      &         UYR , UYUX , UYUY , UYUZ , UYRET ,
  446.      &         UZR , UZUX , UZUY , UZUZ , UZRET ,
  447.      &         RETR , RETUX , RETUY, RETUZ , RETRET
  448. C
  449.  9999 CONTINUE
  450.       RETURN
  451.       END
  452.  
  453.  
  454.  
  455.  
  456.  
  457.  
  458.  
  459.  
  460.  
  461.  
  462.  
  463.  
  464.  
  465.  

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