Télécharger la plaquette Cast3M


fr en

Télécharger cli282.eso

Retour à la liste

Numérotation des lignes : 

cli282

  1. C CLI282    SOURCE    OF166741  24/12/13    21:15:42     12097          
  2.       SUBROUTINE CLI282(NSP,MELEMF,MELEMC,MELECB,MELEFC,MELRES,INORM,
  3.      &     ICHPVO,ICHPSU,LRECP,LRECV,
  4.      &     IROC,IVITC,IPC,IYC,ICHLIM,ILIINC,ILIINP,IJAC,IJACO)
  5. C************************************************************************
  6. C
  7. C PROJET            :  CASTEM 2000
  8. C
  9. C NOM               :  CLI282
  10. C
  11. C DESCRIPTION       :  Subroutine appellée par CLIM22
  12. C                      OPTION: 'INOU' -- Jacobian
  13. C
  14. C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  15. C
  16. C AUTEUR            :  S. Kudriakov, DEN/DM2S/SFME/LTMF
  17. C
  18. C************************************************************************
  19. C
  20. C APPELES (Calcul)  :
  21. C
  22. C************************************************************************
  23. C
  24. C HISTORIQUE (Anomalies et modifications éventuelles)
  25. C
  26. C HISTORIQUE :
  27. C
  28. C************************************************************************
  29. C
  30.       IMPLICIT INTEGER(I-N)
  31. C
  32.       INTEGER MELEMF,MELEMC,MELECB,INORM,ICHPVO,ICHPSU, IROC,IVITC,IPC
  33.      &     ,IGAMC,ICHLIM,ICEL,NFAC,IFAC,MELRES,IJACO
  34.      &     ,NGF,NGC,NLF,NLC,NLCB
  35.      &     ,ILIINC,ILIINP,IJAC,II,JJ,K
  36.      &     ,MP, NBEL, NBME, NBSOUS, NKID, NKMT, NMATRI, NP, NRIGE
  37.      &     ,NSP,I, IYC,J, LRECP,LRECV,KV
  38.       REAL*8 VOLU,SURF,RC,PC,UXC,UYC,GAMC,CNX,CNY,CTX,CTY
  39.      &     ,PSRF,COEF
  40.      &     ,BR1,BOT,TOP,DUNDUY,GAMF,UNC
  41.      &     ,COEF5,DPSRUN,DUXDUN,DRDUN,DUYDUN,DUNDUX,DPDUN
  42.      &     ,GM1,USGM1,RF,UNF,UTF,UXF,UYF,HTF,PF,SF,ECIN
  43.      &     ,TVECT(2),NVECT(2),WVEC_L(4),WVEC_R(4)
  44.       CHARACTER*(8) TYPE
  45.  
  46. -INC PPARAM
  47. -INC CCOPTIO
  48. -INC SMLMOTS
  49. -INC SMELEME
  50.       POINTEUR MELEFC.MELEME
  51. -INC SMLENTI
  52.       POINTEUR MLEMC.MLENTI, MLEMCB.MLENTI,MLEMF.MLENTI
  53. -INC SMCHPOI
  54.       POINTEUR MPNORM.MPOVAL, MPVOL.MPOVAL, MPSURF.MPOVAL, MPRC.MPOVAL,
  55.      &     MPVC.MPOVAL, MPPC.MPOVAL, MPLIM.MPOVAL, MPYC.MPOVAL
  57.       POINTEUR CELL.IZAFM
  58. C-------------------------------------------------------
  59. -INC SMLREEL
  60.       POINTEUR MLRECP.MLREEL, MLRECV.MLREEL
  61. C-------------------------------------------------------
  62. C*********  Les Jacobians ******************************
  63. C-------------------------------------------------------
  64.       SEGMENT JACEL
  65.          REAL*8 JAC(3+NSP,3+NSP)
  66.       ENDSEGMENT
  67.       POINTEUR JLL.JACEL, WL.JACEL, JRR.JACEL, JTT.JACEL, JT2.JACEL,
  68.      &         JFL.JACEL, JFR.JACEL, JR.JACEL, JPL.JACEL, JTL.JACEL
  69. C-------------------------------------------------------------
  70. C******* Les fractionines massiques **************************
  71. C-------------------------------------------------------------
  72.        SEGMENT FRAMAS
  73.          REAL*8 YET(NSP)
  74.       ENDSEGMENT
  75.       POINTEUR YC.FRAMAS, YF.FRAMAS
  76. C-------------------------------------------------------
  77. C**********  Les CP's and CV's   ***********************
  78. C-------------------------------------------------------
  79.       SEGMENT GCONST
  80.          REAL*8 GC(NSP)
  81.       ENDSEGMENT
  82.       POINTEUR CP.GCONST, CV.GCONST
  83. C-------------------------------------------------------------
  84. C**********  Segments for the vectors  ***********************
  85. C-------------------------------------------------------------
  86.        SEGMENT VECEL
  87.          REAL*8 VV(NSP)
  88.        ENDSEGMENT
  89.        POINTEUR DGDYC.VECEL
  90. C--------------------------------------------------
  91.       SEGINI JTT
  92.       SEGINI JTL
  93.       SEGINI JLL
  94.       SEGINI JPL
  95.       SEGINI JRR
  96.       SEGINI JFL
  97.       SEGINI JFR
  98.       SEGINI WL
  99.       SEGINI JR
  100.       SEGINI JT2
  101. C----------------------------------------------------
  102. C**** KRIPAD pour la correspondance global/local
  103. C----------------------------------------------------
  104.       CALL KRIPAD(MELEMC,MLEMC)
  105.       CALL KRIPAD(MELECB,MLEMCB)
  106.       CALL KRIPAD(MELEMF,MLEMF)
  107. C----------------------------------------------------
  108. C**** CHPOINTs de la table DOMAINE
  109. C----------------------------------------------------
  110.       CALL LICHT(INORM,MPNORM,TYPE,ICEL)
  111.       CALL LICHT(ICHPVO,MPVOL,TYPE,ICEL)
  112.       CALL LICHT(ICHPSU,MPSURF,TYPE,ICEL)
  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(IYC,MPYC,TYPE,ICEL)
  120.       CALL LICHT(ICHLIM,MPLIM,TYPE,ICEL)
  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 = (3+NSP)*(3+NSP)
  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.       DO 1 J=1,(NSP+3)
  156.         KV=(J-1)*(3+NSP)
  157.         IMATRI.LISPRI(KV+1) = MLMOTS.MOTS(1)
  158.         IMATRI.LISPRI(KV+2) = MLMOTS.MOTS(2)
  159.         IMATRI.LISPRI(KV+3) = MLMOTS.MOTS(3)
  160.         IMATRI.LISPRI(KV+4) = MLMOTS.MOTS(4)
  161.         DO 2 I=1,(NSP-1)
  162.           IMATRI.LISPRI(KV+4+I) = MLMOTS.MOTS(4+I)
  163.  2      CONTINUE
  164.  1    CONTINUE
  165. C-----------------------------------------------
  166.       SEGDES MLMOTS
  167.       MLMOTS=ILIINC
  168.       SEGACT MLMOTS
  169. C-----------------------------------------------
  170.       DO 3 J=1,(NSP+3)
  171.         KV=(J-1)*(3+NSP)
  172.         IMATRI.LISDUA(KV+1) = MLMOTS.MOTS(j)
  173.         IMATRI.LISDUA(KV+2) = MLMOTS.MOTS(j)
  174.         IMATRI.LISDUA(KV+3) = MLMOTS.MOTS(j)
  175.         IMATRI.LISDUA(KV+4) = MLMOTS.MOTS(j)
  176.         DO 4 I=1,(NSP-1)
  177.           IMATRI.LISDUA(KV+4+I) = MLMOTS.MOTS(j)
  178.  4     CONTINUE
  179.  3    CONTINUE
  180. C-----------------------------------------------
  181. C-----------------------------------------------
  182.       SEGDES MLMOTS
  183.       NBEL = NFAC
  184.       NBSOUS = 1
  185.       NP = 1
  186.       MP = 1
  187. C-----------------------------------------------------------
  188. C-----------------------------------------------------------
  189.       DO 5 I=1,NBME
  190.         SEGINI CELL
  191.         IMATRI.LIZAFM(1,I) = CELL
  192.  5    CONTINUE
  193. C---------------------------------
  194. C---------------------------------
  195. C**** Fin definition MATRIK
  196. C---------------------------------
  197.       DO IFAC=1,NFAC,1
  198.          NGF=MELEFC.NUM(1,IFAC)
  199.          NGC=MELEFC.NUM(2,IFAC)
  200.          NLF=MLEMF.LECT(NGF)
  201.          NLC=MLEMC.LECT(NGC)
  202.          NLCB=MLEMCB.LECT(NGF)
  203.          VOLU=MPVOL.VPOCHA(NLC,1)
  204.          SURF=MPSURF.VPOCHA(NLF,1)
  205. C        In CASTEM les normales sont sortantes
  206.          CNX=MPNORM.VPOCHA(NLF,1)
  207.          CNY=MPNORM.VPOCHA(NLF,2)
  208.          CTX=-1.0D0*CNY
  209.          CTY=CNX
  210. C----------------------------------------------
  211.          SEGINI CP, CV
  212.          MLRECP = LRECP
  213.          MLRECV = LRECV
  214.          SEGACT MLRECP, MLRECV
  215.          DO 10 I=1,(NSP-1)
  216.            CP.GC(I)=MLRECP.PROG(I)
  217.            CV.GC(I)=MLRECV.PROG(I)
  218.  10      CONTINUE
  219.          CP.GC(NSP)=MLRECP.PROG(NSP)
  220.          CV.GC(NSP)=MLRECV.PROG(NSP)
  221. C---------------------------------
  222. C        Variables au centre
  223. C---------------------------------
  224.          RC=MPRC.VPOCHA(NLC,1)
  225.          PC=MPPC.VPOCHA(NLC,1)
  226.          UXC=MPVC.VPOCHA(NLC,1)
  227.          UYC=MPVC.VPOCHA(NLC,2)
  228.          SEGINI YC
  229.          SEGACT MPYC
  230.          DO 100 I=1,(NSP-1)
  231.            YC.YET(I)=MPYC.VPOCHA(NLC,I)
  232.  100     CONTINUE
  233. C---------------------------------
  234. C        Variables à la face
  235. C---------------------------------
  236.          SEGACT MPLIM
  237.          HTF=MPLIM.VPOCHA(NLCB,1)
  238.          SF=MPLIM.VPOCHA(NLCB,2)
  239.          SEGINI YF
  240.          DO 101 I=1,(NSP-1)
  241.            YF.YET(I)=MPLIM.VPOCHA(NLCB,2+I)
  242.  101     CONTINUE
  243.          PF=MPLIM.VPOCHA(NLCB,NSP+2)
  244. c-------------------------------------------------------------
  245. c  Computing GAMMA at the face
  246. c-------------------------------------------------------------
  247.          top=0.0D0
  248.          bot=0.0D0
  249.          do 102  i=1,(nsp-1)
  250.           top=top+yf.yet(i)*(cp.gc(i)-cp.gc(nsp))
  251.           bot=bot+yf.yet(i)*(cv.gc(i)-cv.gc(nsp))
  252.  102     continue
  253.          top=cp.gc(nsp)+top
  254.          bot=cv.gc(nsp)+bot
  255.          GAMF=top/bot
  256.          GM1=GAMF-1.0D0
  257.          USGM1=1.0D0/GM1
  258. c-------------------------------------------------------------
  259. c  Computing GAMMA at the cell-centre
  260. c-------------------------------------------------------------
  261.          top=0.0D0
  262.          bot=0.0D0
  263.          do 103  i=1,(nsp-1)
  264.           top=top+yc.yet(i)*(cp.gc(i)-cp.gc(nsp))
  265.           bot=bot+yc.yet(i)*(cv.gc(i)-cv.gc(nsp))
  266.  103     continue
  267.          top=cp.gc(nsp)+top
  268.          bot=cv.gc(nsp)+bot
  269.          GAMC=top/bot
  270. C-------------------------------------------------------------
  271.          SEGINI DGDYC
  272.          do 41 i=1,(nsp-1)
  273.            dgdyc.vv(i)=(cp.gc(i)-cp.gc(nsp)-
  274.      &       GAMC*(cv.gc(i)-cv.gc(nsp)))/bot
  275.  41    continue
  276. c-------------------------------------------------------------
  277. c  matrix wl(i,j) represents the derivative of the i-component
  278. c  of the vector of primitive variables of the left state with
  279. c  respect to the j-component of the vector of the conservative
  280. c  variables of the left state.
  281. c
  282. c  Here: (rho, ux, uy, p, Y_1,...,Y_(nsp-1)) -
  283. c                               vector of primitive variables;
  284. c  (rho, rho ux, rho uy, rho e, rho Y_1,..., rho Y_(nsp-1)) -
  285. c                            vector of conservative variables.
  286. c-------------------------------------------------------------
  287.        wl.jac(1,1)=1.0d0
  288.        wl.jac(1,2)=0.0d0
  289.        wl.jac(1,3)=0.0d0
  290.        wl.jac(1,4)=0.0d0
  291.        do 83 i=1,(nsp-1)
  292.          wl.jac(1,4+i)=0.0d0
  293.  83    continue
  294. c------------------------------
  295.        wl.jac(2,1)=-UXC/RC
  296.        wl.jac(2,2)=1.0d0/RC
  297.        wl.jac(2,3)=0.0d0
  298.        wl.jac(2,4)=0.0d0
  299.        do 84 i=1,(nsp-1)
  300.          wl.jac(2,4+i)=0.0d0
  301.  84    continue
  302. c------------------------------
  303.        wl.jac(3,1)=-UYC/RC
  304.        wl.jac(3,2)=0.0d0
  305.        wl.jac(3,3)=1.0d0/RC
  306.        wl.jac(3,4)=0.0d0
  307.        do 85 i=1,(nsp-1)
  308.          wl.jac(3,4+i)=0.0d0
  309.  85    continue
  310. c------------------------------
  311.        br1=0.0d0
  312.        do 86 i=1,(nsp-1)
  313.          br1=br1+dgdyc.vv(i)*yc.yet(i)
  314.  86    continue
  315.        br1=br1*PC/(RC*(GAMC-1.0D0))
  316.        wl.jac(4,1)=(GAMC-1.0D0)*(UXC*UXC+UYC*UYC)/2.0d0-br1
  317.        wl.jac(4,2)=-UXC*(GAMC-1.0D0)
  318.        wl.jac(4,3)=-UYC*(GAMC-1.0D0)
  319.        wl.jac(4,4)=(GAMC-1.0D0)
  320.        do 87 i=1,(nsp-1)
  321.          wl.jac(4,4+i)=dgdyc.vv(i)*PC/(RC*(GAMC-1.0D0))
  322.  87    continue
  323. c------------------------------
  324.        do 88 i=1,(nsp-1)
  325.         do 89 j=1,4
  326.           wl.jac(4+i,j)=0.0d0
  327.           if(j.eq.1) wl.jac(4+i,j)=-yc.yet(i)/RC
  328.  89     continue
  329. c------------
  330.         do 890 j=5,(4+nsp-1)
  331.           wl.jac(4+i,j)=0.0d0
  332.           if(4+i.eq.j)  then
  333.             wl.jac(4+i,j)=1.0d0/RC
  334.           endif
  335.  890    continue
  336.  88    continue
  337. c------------------------------------------------
  338.        UNC=(UXC*CNX)+(UYC*CNY)
  339.        IF(UNC .LT. 0.0D0) THEN
  340. c       write(*,*) 'impl_insiiiiii', UNC, NLCB
  341. C-----------------------------------------
  342. C   Option 'INSU' - inlet BC
  343. C-----------------------------------------
  344.          CNX=-1.0D0*CNX
  345.          CNY=-1.0D0*CNY
  346.          CTX=-1.0D0*CNY
  347.          CTY=CNX
  348.          UNC=(UXC*CNX)+(UYC*CNY)
  349.          UNF=UNC
  350.          UTF=0.0D0
  351.          UXF=UNF*CNX+UTF*CTX
  352.          UYF=UNF*CNY+UTF*CTY
  353. C------------------------------------------------
  354. C******* Densite, vitesse, pression sur le bord
  355. C------------------------------------------------
  356.          ECIN=0.5D0*((UXF*UXF)+(UYF*UYF))
  357.          PSRF=(GM1/GAMF)*(HTF-ECIN)
  358.          RF=PSRF/SF
  359.          RF=RF**(1.0D0/GM1)
  360.          PF=SF*(RF**GAMF)
  361. C------------------------------
  362. C******* Derivatives w.r.t. PC
  363. C------------------------------
  364.          DPSRUN=-1*(GM1/GAMF)*UNC
  365.          DRDUN=USGM1*RF/PSRF*DPSRUN
  366.          DPDUN=GAMF*PSRF*DRDUN
  367.          DUXDUN=CNX
  368.          DUYDUN=CNY
  369. C-------------------------------
  370.          DUNDUX=CNX
  371.          DUNDUY=CNY
  372. C--------------------------------------------------------------
  373. C------------------------------
  374. C******* Derivatives
  375. C------------------------------
  376.         wvec_l(1)=RF
  377.         wvec_l(2)=UXF
  378.         wvec_l(3)=UYF
  379.         wvec_l(4)=PF
  380. C--------------------------
  381.         wvec_r(1)=RC
  382.         wvec_r(2)=UXC
  383.         wvec_r(3)=UYC
  384.         wvec_r(4)=PC
  385. C--------------------------
  386.         nvect(1)=CNX
  387.         nvect(2)=CNY
  388.         tvect(1)=CTX
  389.         tvect(2)=CTY
  390.         call copms2(nsp,jfl,jfr,wvec_l,wvec_r,nvect,tvect,
  391.      &                 yf,mpyc,lrecp,lrecv,nlc)
  392. C--------------------------------------------------------------
  393.         COEF=SURF/VOLU
  394.         JLL=JFL
  395.         JRR=JFR
  396.         SEGACT JLL
  397.         SEGACT JRR
  398. C-------------------------------------------------------------
  399. C   Taking in to account the dependance of the variables
  400. C   at the face on the variables at the centre (UNC)
  401. C-------------------------------------------------------------
  402.        DO 105 I=1,(3+NSP)
  403.        DO 1050 J=1,(3+NSP)
  404.         IF(J .EQ. 2) THEN
  405.           COEF5=(JLL.JAC(I,1)*DRDUN)+(JLL.JAC(I,2)*DUXDUN)+
  406.      &        (JLL.JAC(I,3)*DUYDUN)+(JLL.JAC(I,4)*DPDUN)
  407.           JR.JAC(I,J)=(COEF5*DUNDUX+JRR.JAC(I,J))*COEF
  408.         ELSEIF(J .EQ. 3) THEN
  409.           COEF5=(JLL.JAC(I,1)*DRDUN)+(JLL.JAC(I,2)*DUXDUN)+
  410.      &        (JLL.JAC(I,3)*DUYDUN)+(JLL.JAC(I,4)*DPDUN)
  411.           JR.JAC(I,J)=(COEF5*DUNDUY+JRR.JAC(I,J))*COEF
  412.         ELSE
  413.           JR.JAC(I,J)=JRR.JAC(I,J)*COEF
  414.         ENDIF
  415.  1050  CONTINUE
  416.  105   CONTINUE
  417. C-------------------------------------------------------------
  418.        do 114 i=1,(3+nsp)
  419.         do 115 j=1,(3+nsp)
  420.          jt2.jac(i,j)=0.0d0
  421.          do 116 k=1,(3+nsp)
  422.           jt2.jac(i,j)=jt2.jac(i,j)+
  423.      &      jr.jac(i,k)*wl.jac(k,j)/COEF
  424.  116     continue
  425.  115    continue
  426.  114  continue
  427. C----------------------------------------------------------
  428. C******* Jacobian with respect to conservative variables
  429. C----------------------------------------------------------------
  430.          IF(IJAC.EQ.1)THEN
  431.            DO 9 II = 1,(3+NSP)
  432.               DO 15 JJ = 1,(3+NSP)
  433.                 KV = (II-1)*(3+NSP)
  434. C----------------------------------
  435.                 CELL = IMATRI.LIZAFM(1,KV+JJ)
  436.                 CELL.AM(IFAC,1,1) = JT2.JAC(II,JJ)*COEF
  437.  15           CONTINUE
  438.  9         CONTINUE
  439.          ELSEIF(IJAC.EQ.2)THEN
  440.            DO 20 II = 1,(3+NSP)
  441.               DO 25 JJ = 1,(3+NSP)
  442.                 KV = (II-1)*(3+NSP)
  443. C----------------------------------
  444.                 CELL = IMATRI.LIZAFM(1,KV+JJ)
  445.                 CELL.AM(IFAC,1,1) = JR.JAC(II,JJ)
  446.  25           CONTINUE
  447.  20         CONTINUE
  448.          ENDIF
  449. c--------------------------------------------------
  450.       ELSE
  451. C---------------------------------------------
  452. C    Option --- 'OUTP'  - outlet BC
  453. C---------------------------------------------
  454. C------------------------------
  455. C******* Derivatives
  456. C------------------------------
  457.         wvec_l(1)=RC
  458.         wvec_l(2)=UXC
  459.         wvec_l(3)=UYC
  460.         wvec_l(4)=PC
  461. C--------------------------
  462.         wvec_r(1)=RC
  463.         wvec_r(2)=UXC
  464.         wvec_r(3)=UYC
  465.         wvec_r(4)=PF
  466. C--------------------------
  467.         nvect(1)=CNX
  468.         nvect(2)=CNY
  469.         tvect(1)=CTX
  470.         tvect(2)=CTY
  471.         call copmsp(nsp,jpl,jll,wvec_l,wvec_r,nvect,tvect,
  472.      &                   mpyc,lrecp,lrecv,nlc,nlc)
  473. C----------------------------------------
  474.          COEF=-SURF/VOLU
  475. C----------------------------------------
  476.          JTT=JLL
  477.          JTL=JPL
  478.          SEGACT JTT
  479.          SEGACT JTL
  480. C----------------------------------------------------------------
  481. C******* Jacobian with respect to conservative variables
  482. C----------------------------------------------------------------
  483.          IF(IJAC.EQ.1)THEN
  484.            DO 99 II = 1,(3+NSP)
  485.               DO 155 JJ = 1,(3+NSP)
  486.                 KV = (II-1)*(3+NSP)
  487. C----------------------------------
  488.                 CELL = IMATRI.LIZAFM(1,KV+JJ)
  489.                 CELL.AM(IFAC,1,1) = JTT.JAC(II,JJ)*COEF
  490.  155          CONTINUE
  491.  99        CONTINUE
  492.          ELSEIF(IJAC.EQ.2)THEN
  493.            DO 200 II = 1,(3+NSP)
  494.               DO 255 JJ = 1,(3+NSP)
  495.                 KV = (II-1)*(3+NSP)
  496. C----------------------------------
  497.                 CELL = IMATRI.LIZAFM(1,KV+JJ)
  498.                 CELL.AM(IFAC,1,1) = JTL.JAC(II,JJ)*COEF
  499.  255          CONTINUE
  500.  200        CONTINUE
  501.          ENDIF
  502. C---------------------------------------------------------------
  503.         ENDIF
  504. C----------------------------------------------------------------
  505.       ENDDO
  506. C
  507.       SEGDES MELEFC
  508. C
  509.       SEGDES MLEMC
  510.       SEGDES MLEMCB
  511.       SEGDES MLEMF
  512. C
  513.       SEGDES MPNORM
  514.       SEGDES MPVOL
  515.       SEGDES MPSURF
  516.       SEGDES MPRC
  517.       SEGDES MPPC
  518.       SEGDES MPVC
  519.       SEGDES MPYC
  520.       SEGDES MPLIM
  521.       SEGDES YC
  522.       SEGDES YF
  523.       SEGDES CP
  524.       SEGDES CV
  525.       SEGDES JLL
  526.       SEGDES JRR
  527.       SEGDES JTT
  528.       SEGDES JTL
  529.       SEGDES JFL
  530.       SEGDES JFR
  531.       SEGDES JR
  532.       SEGDES WL
  533.       SEGDES  DGDYC
  534.       SEGDES MATRIK
  535.       DO 80 II=1,NBME
  536.        CELL = IMATRI.LIZAFM(1,II)
  537.        SEGDES CELL
  538.  80   CONTINUE
  539.       SEGDES IMATRI
  540. C---------------------------------------------
  541.  9999 CONTINUE
  542.       RETURN
  543.       END
  544.  
  545.  
  546.  
  547.  
  548.  
  549.  
  550.  
  551.  
  552.  
  553.  
  554.  
  555.  
  556.  
  557.  
  558.  
  559.  
  560.  

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