Télécharger la plaquette Cast3M


fr en

Télécharger cli182.eso

Retour à la liste

Numérotation des lignes : 

cli182

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

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