Télécharger la plaquette Cast3M


fr en

Télécharger cli122.eso

Retour à la liste

Numérotation des lignes : 

cli122

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

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