Télécharger la plaquette Cast3M


fr en

Télécharger cli163.eso

Retour à la liste

Numérotation des lignes : 

cli163

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

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