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

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