Télécharger la plaquette Cast3M


fr en

Télécharger cli123.eso

Retour à la liste

Numérotation des lignes : 

cli123

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

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