Télécharger la plaquette Cast3M


fr en

Télécharger cli113.eso

Retour à la liste

Numérotation des lignes : 

cli113

  1. C CLI113    SOURCE    CB215821  20/11/25    13:19:58     10792          
  2.       SUBROUTINE CLI113(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               :  CLI113
  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.      &     ,CNZ,CTZ,CT1X,CT1Y,CT1Z,UZC,UZF,UT1F,UT1,UZ
  69.      &     ,DUZDG1,DFMZG1,DG1DUZ,DROUZ,DUXUZ,DUYUZ,DPUZ
  70.      &     ,DUZRO,DUZUX,DUZUY,DUZUZ,DUZP
  71.      &     ,RF,PF,UXF,UYF
  72.      &     ,UNC,UNF,UTF,SF,ASONC,ASONF,ASON
  73. *     &     ,UTC,SC
  74.      &     ,GM1,USGM1,DSGM1,G1,G3,ASON2,S,UT,UN,RHO,P,UX,UY
  75.      &     ,DUNDG1,DASDG1,DRHDG1,DPDG1,DUXDG1,DUYDG1
  76. *     &     ,CELL,EPS,CACCA
  77.      &     ,DFRDG1,DFMXG1,DFMYG1,DFEDG1
  78.      &     ,DG1DR,DG1DP,DG1DUX,DG1DUY,COEF
  79.      &     ,DRORO,DROUX,DROUY,DROP
  80.      &     ,DUXRO,DUXUX,DUXUY,DUXP
  81.      &     ,DUYRO,DUYUX,DUYUY,DUYP
  82.      &     ,DPRO,DPUX,DPUY,DPP
  83.      &     ,COEF1,COEF2,COEF3,COEF4
  84.       CHARACTER*(8) TYPE
  85. C
  86. C
  87. C**** KRIPAD pour la correspondance global/local
  88. C
  89.       CALL KRIPAD(MELEMC,MLEMC)
  90. C     SEGINI MLEMC
  91.       CALL KRIPAD(MELECB,MLEMCB)
  92. C     SEGINI MLEMCB
  93.       CALL KRIPAD(MELEMF,MLEMF)
  94. C     SEGINI MLEMF
  95. C
  96. C**** CHPOINTs de la table DOMAINE
  97. C
  98.       CALL LICHT(INORM,MPNORM,TYPE,ICEL)
  99.       CALL LICHT(ICHPVO,MPVOL,TYPE,ICEL)
  100.       CALL LICHT(ICHPSU,MPSURF,TYPE,ICEL)
  101. C
  102. C**** LICHT active les MPOVALs en *MOD
  103. C
  104. C     SEGACT MPNORM*MOD
  105. C     SEGACT MPOVSU*MOD
  106. C     SEGACT MPOVOL*MOD
  107. C
  108. C
  109. C**** CHPOINTs des variables
  110. C
  111.       CALL LICHT(IROC,MPRC,TYPE,ICEL)
  112.       CALL LICHT(IVITC,MPVC,TYPE,ICEL)
  113.       CALL LICHT(IPC,MPPC,TYPE,ICEL)
  114.       CALL LICHT(IGAMC,MPGAMC,TYPE,ICEL)
  115.       CALL LICHT(ICHLIM,MPLIM,TYPE,ICEL)
  116. C
  117. C     SEGACT *MOD
  118. C     SEGACT *MOD
  119. C     SEGACT *MOD
  120. C     SEGACT *MOD
  121. C     SEGACT *MOD
  122. C
  123. C
  124. C**** Boucle sur le face pour le calcul des invariants de
  125. C     Riemann et du flux
  126. C
  127.       SEGACT MELEFC
  128.       NFAC=MELEFC.NUM(/2)
  129. C
  130. C**** Objet MATRIK
  131. C
  132.       NRIGE = 7
  133.       NMATRI = 1
  134.       NKID = 9
  135.       NKMT = 7
  136. C
  137.       SEGINI MATRIK
  138.       IJACO = MATRIK
  139.       MATRIK.IRIGEL(1,1) = MELRES
  140.       MATRIK.IRIGEL(2,1) = MELRES
  141. C
  142. C**** Matrice non symetrique
  143. C
  144.       MATRIK.IRIGEL(7,1) = 2
  145. C
  146.       NBME = 25
  147.       NBSOUS = 1
  148.       SEGINI IMATRI
  149.       IF(IJAC.EQ.1)THEN
  150.          MLMOTS=ILIINC
  151.       ELSEIF(IJAC.EQ.2)THEN
  152.          MLMOTS=ILIINP
  153.       ENDIF
  154.       SEGACT MLMOTS
  155.       MATRIK.IRIGEL(4,1) = IMATRI
  156. C
  157.       IMATRI.LISPRI(1) = MLMOTS.MOTS(1)
  158.       IMATRI.LISPRI(2) = MLMOTS.MOTS(2)
  159.       IMATRI.LISPRI(3) = MLMOTS.MOTS(3)
  160.       IMATRI.LISPRI(4) = MLMOTS.MOTS(4)
  161.       IMATRI.LISPRI(5) = MLMOTS.MOTS(5)
  162. C
  163.       IMATRI.LISPRI(6) = MLMOTS.MOTS(1)
  164.       IMATRI.LISPRI(7) = MLMOTS.MOTS(2)
  165.       IMATRI.LISPRI(8) = MLMOTS.MOTS(3)
  166.       IMATRI.LISPRI(9) = MLMOTS.MOTS(4)
  167.       IMATRI.LISPRI(10) = MLMOTS.MOTS(5)
  168. C
  169.       IMATRI.LISPRI(11) = MLMOTS.MOTS(1)
  170.       IMATRI.LISPRI(12) = MLMOTS.MOTS(2)
  171.       IMATRI.LISPRI(13) = MLMOTS.MOTS(3)
  172.       IMATRI.LISPRI(14) = MLMOTS.MOTS(4)
  173.       IMATRI.LISPRI(15) = MLMOTS.MOTS(5)
  174. C
  175.       IMATRI.LISPRI(16) = MLMOTS.MOTS(1)
  176.       IMATRI.LISPRI(17) = MLMOTS.MOTS(2)
  177.       IMATRI.LISPRI(18) = MLMOTS.MOTS(3)
  178.       IMATRI.LISPRI(19) = MLMOTS.MOTS(4)
  179.       IMATRI.LISPRI(20) = MLMOTS.MOTS(5)
  180. C
  181.       IMATRI.LISPRI(21) = MLMOTS.MOTS(1)
  182.       IMATRI.LISPRI(22) = MLMOTS.MOTS(2)
  183.       IMATRI.LISPRI(23) = MLMOTS.MOTS(3)
  184.       IMATRI.LISPRI(24) = MLMOTS.MOTS(4)
  185.       IMATRI.LISPRI(25) = MLMOTS.MOTS(5)
  186. C
  187.       SEGDES MLMOTS
  188.       MLMOTS=ILIINC
  189.       SEGACT MLMOTS
  190. C
  191.       IMATRI.LISDUA(1) = MLMOTS.MOTS(1)
  192.       IMATRI.LISDUA(2) = MLMOTS.MOTS(1)
  193.       IMATRI.LISDUA(3) = MLMOTS.MOTS(1)
  194.       IMATRI.LISDUA(4) = MLMOTS.MOTS(1)
  195.       IMATRI.LISDUA(5) = MLMOTS.MOTS(1)
  196. C
  197.       IMATRI.LISDUA(6) = MLMOTS.MOTS(2)
  198.       IMATRI.LISDUA(7) = MLMOTS.MOTS(2)
  199.       IMATRI.LISDUA(8) = MLMOTS.MOTS(2)
  200.       IMATRI.LISDUA(9) = MLMOTS.MOTS(2)
  201.       IMATRI.LISDUA(10) = MLMOTS.MOTS(2)
  202. C
  203.       IMATRI.LISDUA(11) = MLMOTS.MOTS(3)
  204.       IMATRI.LISDUA(12) = MLMOTS.MOTS(3)
  205.       IMATRI.LISDUA(13) = MLMOTS.MOTS(3)
  206.       IMATRI.LISDUA(14) = MLMOTS.MOTS(3)
  207.       IMATRI.LISDUA(15) = MLMOTS.MOTS(3)
  208. C
  209.       IMATRI.LISDUA(16) = MLMOTS.MOTS(4)
  210.       IMATRI.LISDUA(17) = MLMOTS.MOTS(4)
  211.       IMATRI.LISDUA(18) = MLMOTS.MOTS(4)
  212.       IMATRI.LISDUA(19) = MLMOTS.MOTS(4)
  213.       IMATRI.LISDUA(20) = MLMOTS.MOTS(4)
  214. C
  215.       IMATRI.LISDUA(21) = MLMOTS.MOTS(5)
  216.       IMATRI.LISDUA(22) = MLMOTS.MOTS(5)
  217.       IMATRI.LISDUA(23) = MLMOTS.MOTS(5)
  218.       IMATRI.LISDUA(24) = MLMOTS.MOTS(5)
  219.       IMATRI.LISDUA(25) = MLMOTS.MOTS(5)
  220. C
  221.       SEGDES MLMOTS
  222.       NBEL = NFAC
  223.       NBSOUS = 1
  224.       NP = 1
  225.       MP = 1
  226.       SEGINI   RR  , RUX , RUY , RUZ ,RRET ,
  227.      &         UXR , UXUX , UXUY , UXUZ ,UXRET ,
  228.      &         UYR , UYUX , UYUY , UYUZ , UYRET ,
  229.      &         UZR , UZUX , UZUY , UZUZ , UZRET ,
  230.      &         RETR , RETUX , RETUY, RETUZ , RETRET
  231. C
  232. C**** Duale = IMATRI.LISDUA(1) = 'RN'
  233. C     Primale = IMATRI.LISPRI(1) = 'RN'
  234. C     -> IMATRI.LIZAFM(1,1) = RR
  235. C
  236.       IMATRI.LIZAFM(1,1) = RR
  237.       IMATRI.LIZAFM(1,2) = RUX
  238.       IMATRI.LIZAFM(1,3) = RUY
  239.       IMATRI.LIZAFM(1,4) = RUZ
  240.       IMATRI.LIZAFM(1,5) = RRET
  241. C
  242.       IMATRI.LIZAFM(1,6) = UXR
  243.       IMATRI.LIZAFM(1,7) = UXUX
  244.       IMATRI.LIZAFM(1,8) = UXUY
  245.       IMATRI.LIZAFM(1,9) = UXUZ
  246.       IMATRI.LIZAFM(1,10) = UXRET
  247. C
  248.       IMATRI.LIZAFM(1,11) = UYR
  249.       IMATRI.LIZAFM(1,12) = UYUX
  250.       IMATRI.LIZAFM(1,13) = UYUY
  251.       IMATRI.LIZAFM(1,14) = UYUZ
  252.       IMATRI.LIZAFM(1,15) = UYRET
  253. C
  254.       IMATRI.LIZAFM(1,16) = UZR
  255.       IMATRI.LIZAFM(1,17) = UZUX
  256.       IMATRI.LIZAFM(1,18) = UZUY
  257.       IMATRI.LIZAFM(1,19) = UZUZ
  258.       IMATRI.LIZAFM(1,20) = UZRET
  259. C
  260.       IMATRI.LIZAFM(1,21) = RETR
  261.       IMATRI.LIZAFM(1,22) = RETUX
  262.       IMATRI.LIZAFM(1,23) = RETUY
  263.       IMATRI.LIZAFM(1,24) = RETUZ
  264.       IMATRI.LIZAFM(1,25) = RETRET
  265. C
  266.       SEGDES MATRIK
  267.       SEGDES IMATRI
  268. C
  269. C**** Fin definition MATRIK
  270. C
  271.       DO IFAC=1,NFAC,1
  272.          NGF=MELEFC.NUM(1,IFAC)
  273.          NGC=MELEFC.NUM(2,IFAC)
  274.          NLF=MLEMF.LECT(NGF)
  275.          NLC=MLEMC.LECT(NGC)
  276.          NLCB=MLEMCB.LECT(NGF)
  277.          VOLU=MPVOL.VPOCHA(NLC,1)
  278.          SURF=MPSURF.VPOCHA(NLF,1)
  279. C        In CASTEM les normales sont sortantes
  280.          CNX=-1*MPNORM.VPOCHA(NLF,1)
  281.          CNY=-1*MPNORM.VPOCHA(NLF,2)
  282.          CNZ=-1*MPNORM.VPOCHA(NLF,3)
  283.          CTX=-1*MPNORM.VPOCHA(NLF,4)
  284.          CTY=-1*MPNORM.VPOCHA(NLF,5)
  285.          CTZ=-1*MPNORM.VPOCHA(NLF,6)
  286.          CT1X=-1*MPNORM.VPOCHA(NLF,7)
  287.          CT1Y=-1*MPNORM.VPOCHA(NLF,8)
  288.          CT1Z=-1*MPNORM.VPOCHA(NLF,9)
  289. C        Variables au centre
  290.          RC=MPRC.VPOCHA(NLC,1)
  291.          PC=MPPC.VPOCHA(NLC,1)
  292.          UXC=MPVC.VPOCHA(NLC,1)
  293.          UYC=MPVC.VPOCHA(NLC,2)
  294.          UZC=MPVC.VPOCHA(NLC,3)
  295.          GAMC=MPGAMC.VPOCHA(NLC,1)
  296. C        Variables à la face
  297.          RF=MPLIM.VPOCHA(NLCB,1)
  298.          UXF=MPLIM.VPOCHA(NLCB,2)
  299.          UYF=MPLIM.VPOCHA(NLCB,3)
  300.          UZF=MPLIM.VPOCHA(NLCB,4)
  301.          PF=MPLIM.VPOCHA(NLCB,IDIM+2)
  302. C
  303. C******* On calcule UN, UT ASON, S
  304. C
  305.          UNC=(UXC*CNX)+(UYC*CNY)+(UZC*CNZ)
  306. *         UTC=(UXC*CTX)+(UYC*CTY)
  307.          UNF=(UXF*CNX)+(UYF*CNY)+(UZF*CNZ)
  308.          UTF=(UXF*CTX)+(UYF*CTY)+(UZF*CTZ)
  309.          UT1F=(UXF*CT1X)+(UYF*CT1Y)+(UZF*CT1Z)
  310. C
  311.          ASONC=(GAMC*PC/RC)**0.5D0
  312.          ASONF=(GAMC*PF/RF)**0.5D0
  313. C
  314. *         SC=PC/(RC**GAMC)
  315.          SF=PF/(RF**GAMC)
  316. C
  317. C******* Densite, vitesse, pression sur le bord
  318. C
  319.          GM1=(GAMC-1.0D0)
  320.          USGM1=1.0D0/GM1
  321.          DSGM1=2.0D0*USGM1
  322.          G1=UNC-(DSGM1*ASONC)
  323.          G3=UNF+(DSGM1*ASONF)
  324.          UN=0.5D0*(G1+G3)
  325.          ASON=(0.5D0*(G3-G1))
  326.          ASON=ASON/DSGM1
  327.          ASON2=ASON*ASON
  328.          S=SF
  329.          UT=UTF
  330.          UT1=UT1F
  331.          RHO=ASON2/(GAMC*S)
  332.          RHO=RHO**USGM1
  333.          P=RHO*ASON2/GAMC
  334.          UX=(UN*CNX)+(UT*CTX)+(UT1*CT1X)
  335.          UY=(UN*CNY)+(UT*CTY)+(UT1*CT1Y)
  336.          UZ=(UN*CNZ)+(UT*CTZ)+(UT1*CT1Z)
  337. C
  338. C******* Derivatives
  339. C
  340.          DUNDG1=0.5D0
  341.          DASDG1=-0.5D0/DSGM1
  342.          DRHDG1=GAMC*S
  343.          DRHDG1=1.0D0/DRHDG1
  344.          DRHDG1=DRHDG1**USGM1
  345.          DRHDG1=DRHDG1*DSGM1
  346.          DRHDG1=DRHDG1*(ASON**((3.0D0-GAMC)/(GAMC-1.0D0)))
  347.          DRHDG1=DRHDG1*DASDG1
  348.          DPDG1=((ASON2/GAMC)*DRHDG1)+(((2*ASON*RHO)/GAMC)*DASDG1)
  349. CC
  350. CC******* Test
  351. CC
  352. C         DFEDG1=(DUNDG1*GAMC*USGM1*P) + (UN*GAMC*USGM1*DPDG1) +
  353. C     &        (0.5D0*DRHDG1*UN*((UN*UN)+(UT*UT))) +
  354. C     &        (0.5D0*RHO*DUNDG1*((UN*UN)+(UT*UT))) +
  355. C     &        (RHO*UN*UN*DUNDG1)
  356. C         CACCA=UN*GAMC*USGM1*P + 0.5D0*RHO*UN*(UN*UN+UT*UT)
  357. C         EPS=1.0D-6
  358. C         G1=G1*(1+EPS)
  359. C         CELL=UN
  360. C         UN=0.5D0*(G1+G3)
  361. C         write(*,*) DUNDG1, ((UN - CELL)/(EPS*G1))
  362. C         CELL=ASON
  363. C         ASON=(0.5D0*(G3-G1))
  364. C         ASON=ASON/DSGM1
  365. C         write(*,*) DASDG1, ((ASON - CELL)/(EPS*G1))
  366. C         ASON2=ASON*ASON
  367. C         S=SF
  368. C         UT=UTF
  369. C         CELL=RHO
  370. C         RHO=ASON2/(GAMC*S)
  371. C         RHO=RHO**USGM1
  372. C         write(*,*) DRHDG1, ((RHO - CELL)/(EPS*G1))
  373. C         CELL=P
  374. C         P=RHO*ASON2/GAMC
  375. C         write(*,*) DPDG1, ((P - CELL)/(EPS*G1))
  376. C         CELL=CACCA
  377. C         CACCA=UN*GAMC*USGM1*P + 0.5D0*RHO*UN*(UN*UN+UT*UT)
  378. C         write(*,*) DFEDG1, ((CACCA - CELL)/(EPS*G1))
  379. CC
  380. CC************************************************************
  381. CC
  382.          DUXDG1=DUNDG1*CNX
  383.          DUYDG1=DUNDG1*CNY
  384.          DUZDG1=DUNDG1*CNZ
  385. C
  386.          DFRDG1=(DRHDG1*UN)+(RHO*DUNDG1)
  387.          DFMXG1=(DRHDG1*UN*UX)+(RHO*DUNDG1*UX)+
  388.      &        (RHO*UN*DUXDG1)+(DPDG1*CNX)
  389.          DFMYG1=(DRHDG1*UN*UY)+(RHO*DUNDG1*UY)+
  390.      &        (RHO*UN*DUYDG1)+(DPDG1*CNY)
  391.          DFMZG1=(DRHDG1*UN*UZ)+(RHO*DUNDG1*UZ)+
  392.      &        (RHO*UN*DUZDG1)+(DPDG1*CNZ)
  393.          DFEDG1=(DUNDG1*GAMC*USGM1*P) + (UN*GAMC*USGM1*DPDG1) +
  394.      &        (0.5D0*DRHDG1*UN*((UN*UN)+(UT*UT)+(UT1*UT1))) +
  395.      &        (0.5D0*RHO*DUNDG1*((UN*UN)+(UT*UT)+(UT1*UT1))) +
  396.      &        (RHO*UN*UN*DUNDG1)
  397. C
  398. C******* Jacobian with respect to primitive variables
  399. C
  400.          DG1DR=USGM1*ASONC/RC
  401.          DG1DP=-1.0D0*USGM1*ASONC/PC
  402.          DG1DUX=CNX
  403.          DG1DUY=CNY
  404.          DG1DUZ=CNZ
  405. C
  406.          COEF=SURF/VOLU
  407. C
  408.          DRORO=DFRDG1*DG1DR*COEF
  409.          DROUX=DFRDG1*DG1DUX*COEF
  410.          DROUY=DFRDG1*DG1DUY*COEF
  411.          DROUZ=DFRDG1*DG1DUZ*COEF
  412.          DROP=DFRDG1*DG1DP*COEF
  413. C
  414.          DUXRO=DFMXG1*DG1DR*COEF
  415.          DUXUX=DFMXG1*DG1DUX*COEF
  416.          DUXUY=DFMXG1*DG1DUY*COEF
  417.          DUXUZ=DFMXG1*DG1DUZ*COEF
  418.          DUXP=DFMXG1*DG1DP*COEF
  419. C
  420.          DUYRO=DFMYG1*DG1DR*COEF
  421.          DUYUX=DFMYG1*DG1DUX*COEF
  422.          DUYUY=DFMYG1*DG1DUY*COEF
  423.          DUYUZ=DFMYG1*DG1DUZ*COEF
  424.          DUYP=DFMYG1*DG1DP*COEF
  425. C
  426.          DUZRO=DFMZG1*DG1DR*COEF
  427.          DUZUX=DFMZG1*DG1DUX*COEF
  428.          DUZUY=DFMZG1*DG1DUY*COEF
  429.          DUZUZ=DFMZG1*DG1DUZ*COEF
  430.          DUZP=DFMZG1*DG1DP*COEF
  431. C
  432.          DPRO=DFEDG1*DG1DR*COEF
  433.          DPUX=DFEDG1*DG1DUX*COEF
  434.          DPUY=DFEDG1*DG1DUY*COEF
  435.          DPUZ=DFEDG1*DG1DUZ*COEF
  436.          DPP=DFEDG1*DG1DP*COEF
  437. C
  438. C******* Jacobian with respect to conservative variables
  439. C
  440.          IF(IJAC.EQ.1)THEN
  441. C
  442.             COEF1=-1.0D0*UXC/RC
  443.             COEF2=-1.0D0*UYC/RC
  444.             COEF3=-1.0D0*UZC/RC
  445.             COEF4=0.5D0*GM1*((UXC*UXC)+(UYC*UYC)+(UZC*UZC))
  446. C
  447.             RR.AM(IFAC,1,1)=DRORO+(DROUX*COEF1)+(DROUY*COEF2)+(DROUZ
  448.      $           *COEF3)+(DROP*COEF4)
  449.             RUX.AM(IFAC,1,1)=(DROUX/RC)-((UXC*GM1)*DROP)
  450.             RUY.AM(IFAC,1,1)=(DROUY/RC)-((UYC*GM1)*DROP)
  451.             RUZ.AM(IFAC,1,1)=(DROUZ/RC)-((UZC*GM1)*DROP)
  452.             RRET.AM(IFAC,1,1)=GM1*DROP
  453. C
  454.             UXR.AM(IFAC,1,1)=DUXRO+(DUXUX*COEF1)+(DUXUY*COEF2)+(DUXUZ
  455.      $           *COEF3)+(DUXP*COEF4)
  456.             UXUX.AM(IFAC,1,1)=(DUXUX/RC)-((UXC*GM1)*DUXP)
  457.             UXUY.AM(IFAC,1,1)=(DUXUY/RC)-((UYC*GM1)*DUXP)
  458.             UXUZ.AM(IFAC,1,1)=(DUXUZ/RC)-((UZC*GM1)*DUXP)
  459.             UXRET.AM(IFAC,1,1)=GM1*DUXP
  460. C
  461.             UYR.AM(IFAC,1,1)=DUYRO+(DUYUX*COEF1)+(DUYUY*COEF2)+(DUYUZ
  462.      $           *COEF3)+(DUYP*COEF4)
  463.             UYUX.AM(IFAC,1,1)=(DUYUX/RC)-((UXC*GM1)*DUYP)
  464.             UYUY.AM(IFAC,1,1)=(DUYUY/RC)-((UYC*GM1)*DUYP)
  465.             UYUZ.AM(IFAC,1,1)=(DUYUZ/RC)-((UZC*GM1)*DUYP)
  466.             UYRET.AM(IFAC,1,1)=GM1*DUYP
  467. C
  468.             UZR.AM(IFAC,1,1)=DUZRO+(DUZUX*COEF1)+(DUZUY*COEF2)+(DUZUZ
  469.      $           *COEF3)+(DUZP*COEF4)
  470.             UZUX.AM(IFAC,1,1)=(DUZUX/RC)-((UXC*GM1)*DUZP)
  471.             UZUY.AM(IFAC,1,1)=(DUZUY/RC)-((UYC*GM1)*DUZP)
  472.             UZUZ.AM(IFAC,1,1)=(DUZUZ/RC)-((UZC*GM1)*DUZP)
  473.             UZRET.AM(IFAC,1,1)=GM1*DUZP
  474. C
  475.             RETR.AM(IFAC,1,1)=DPRO+(DPUX*COEF1)+(DPUY*COEF2)+(DPUZ
  476.      $           *COEF3)+(DPP*COEF4)
  477.             RETUX.AM(IFAC,1,1)=(DPUX/RC)-((UXC*GM1)*DPP)
  478.             RETUY.AM(IFAC,1,1)=(DPUY/RC)-((UYC*GM1)*DPP)
  479.             RETUZ.AM(IFAC,1,1)=(DPUZ/RC)-((UZC*GM1)*DPP)
  480.             RETRET.AM(IFAC,1,1)=GM1*DPP
  481. C
  482.          ELSEIF(IJAC.EQ.2)THEN
  483.             RR.AM(IFAC,1,1)=DRORO
  484.             RUX.AM(IFAC,1,1)=DROUX
  485.             RUY.AM(IFAC,1,1)=DROUY
  486.             RUZ.AM(IFAC,1,1)=DROUZ
  487.             RRET.AM(IFAC,1,1)=DROP
  488. C
  489.             UXR.AM(IFAC,1,1)=DUXRO
  490.             UXUX.AM(IFAC,1,1)=DUXUX
  491.             UXUY.AM(IFAC,1,1)=DUXUY
  492.             UXUZ.AM(IFAC,1,1)=DUXUZ
  493.             UXRET.AM(IFAC,1,1)=DUXP
  494. C
  495.             UYR.AM(IFAC,1,1)=DUYRO
  496.             UYUX.AM(IFAC,1,1)=DUYUX
  497.             UYUY.AM(IFAC,1,1)=DUYUY
  498.             UYUZ.AM(IFAC,1,1)=DUYUZ
  499.             UYRET.AM(IFAC,1,1)=DUYP
  500. C
  501.             UZR.AM(IFAC,1,1)=DUZRO
  502.             UZUX.AM(IFAC,1,1)=DUZUX
  503.             UZUY.AM(IFAC,1,1)=DUZUY
  504.             UZUZ.AM(IFAC,1,1)=DUZUZ
  505.             UZRET.AM(IFAC,1,1)=DUZP
  506. C
  507.             RETR.AM(IFAC,1,1)=DPRO
  508.             RETUX.AM(IFAC,1,1)=DPUX
  509.             RETUY.AM(IFAC,1,1)=DPUY
  510.             RETUZ.AM(IFAC,1,1)=DPUZ
  511.             RETRET.AM(IFAC,1,1)=DPP
  512.          ENDIF
  513.  
  514.       ENDDO
  515. C
  516.       SEGDES MELEFC
  517. C
  518.       SEGSUP MLEMC
  519.       SEGSUP MLEMCB
  520.       SEGSUP MLEMF
  521. C
  522.       SEGDES MPNORM
  523.       SEGDES MPVOL
  524.       SEGDES MPSURF
  525.       SEGDES MPRC
  526.       SEGDES MPPC
  527.       SEGDES MPVC
  528.       SEGDES MPGAMC
  529.       SEGDES MPLIM
  530. C
  531.       SEGDES   RR  , RUX , RUY , RUZ ,RRET ,
  532.      &         UXR , UXUX , UXUY , UXUZ ,UXRET ,
  533.      &         UYR , UYUX , UYUY , UYUZ , UYRET ,
  534.      &         UZR , UZUX , UZUY , UZUZ , UZRET ,
  535.      &         RETR , RETUX , RETUY, RETUZ , RETRET
  536. C
  537.  9999 CONTINUE
  538.       RETURN
  539.       END
  540.  
  541.  
  542.  
  543.  
  544.  
  545.  
  546.  
  547.  
  548.  
  549.  
  550.  
  551.  
  552.  
  553.  

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