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

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