Télécharger la plaquette Cast3M


fr en

Télécharger elfres.eso

Retour à la liste

Numérotation des lignes : 

elfres

  1. C ELFRES    SOURCE    CB215821  25/04/23    21:15:17     12247          
  2.       SUBROUTINE ELFRES(KANBN,KDNCN,KNREFE,KOMALI,KACHAR,DELTAT,M,NPAS,
  3.      *NFOIS,KSOLUT)
  4.       IMPLICIT INTEGER(I-N)
  5.       IMPLICIT REAL*8(A-H,O-Z)
  6. C
  7. C  ====================================================================
  8. C    SOUS-PROGRAMME EFFECTUANT LA RESOLUTION TEMPORELLE
  9. C    SOUS-PROGRAMME APPELE PAR ELFE
  10. C    N'APPELLE QUE DES SOUS-PROGRAMMES FORTRAN
  11. C    CREATION    :  3/11/86      MODIFICATION LE  22/01/88
  12. C    PROGRAMMEUR : GUILBAUD                  PAR LIONEL VIVAN
  13. C  ====================================================================
  14. C
  15.  
  16. -INC PPARAM
  17. -INC CCOPTIO
  18. -INC CCHAMP
  19. -INC SMELEME
  20. -INC SMCHPOI
  21. -INC SMCOORD
  22. -INC SMSOLUT
  23. -INC SMLREEL
  24. C
  25.       SEGMENT MANBN
  26.         POINTEUR KAB(NSGA).ANBN
  27.       ENDSEGMENT
  28. C
  29. C   NSGA  : NOMBRE DE BLOCS ANBN
  30. C
  31.       SEGMENT ANBN
  32.         REAL*8 AB(NTANBN,LANBN)
  33.       ENDSEGMENT
  34. C
  35. C   AB(I,K) : TERME I DE LA MATRICE A OU B D'UN ELEMENT AU TEMPS K
  36. C   NTANBN : NOMBRE DE TERMES DES MATRICES A ET B DE TOUS LES ELEMENTS
  37. C   LANBN  : NOMBRE DE PAS DE TEMPS STOCKES DANS UN BLOC ANBN
  38. C
  39.       SEGMENT MDNCN
  40.         POINTEUR KDC(NSGD).DNCN
  41.       ENDSEGMENT
  42. C
  43. C   NSGD  : NOMBRE DE BLOCS DNCN
  44. C
  45.       SEGMENT DNCN
  46.         REAL*8 DC(NIDNCN,LDNCN)
  47.       ENDSEGMENT
  48. C
  49. C   DC(I,K) : DDL I AU TEMPS K
  50. C   NIDNCN : NOMBRE TOTAL D'INCONNUES
  51. C   LDNCN  : NOMBRE DE PAS DE TEMPS STOCKES DANS UN BLOC DNCN
  52. C
  53.       SEGMENT MNREFE
  54.         INTEGER NREFE(8,NSTR)
  55.         INTEGER NTANBN
  56.         INTEGER NIDNCN
  57.         INTEGER NTVN
  58.         POINTEUR NREPA.MPASS
  59.         POINTEUR NRECA.MCARA
  60.         POINTEUR NRENO.MNORM
  61.         POINTEUR NRECPR.ICPR
  62.         POINTEUR NREMEL.MELEME
  63.         POINTEUR NREDEN.MDEN
  64.       ENDSEGMENT
  65. C
  66. C   NSTR       : NOMBRE D'ELEMENTS
  67. C   NREFE(1,I) : MELEME
  68. C   NREFE(2,I) : MSOSTU
  69. C   NREFE(3,I) : TYPE DE L'ELEMENT
  70. C   NREFE(4,I) : NOMBRE DE POINTS DU MELEME
  71. C   NREFE(5,I) : NOMBRE DE DDL PAR POINT
  72. C   NREFE(6,I)=IVN :LE 1ER DDL DE L'ELEMENT EST LE IVN+1 IEME DE VN
  73. C   NREFE(7,I)=IAN :LE 1ER TERME DE LA MATRICE A EST LE IAN IEME DE ANBN
  74. C   NREFE(8,I)= 1 :LE IEME ELEMENT EST RIGIDE (OU PARTIELLEMENT) SINON 0
  75. C   NTANBN : NOMBRE DE TERMES DES MATRICES A ET B POUR TOUS LES ELEMENTS
  76. C   NIDNCN : NOMBRE TOTAL D'INCONNUES DE DNCN
  77. C   NTVN   : LONGUEUR DU TABLEAU VN
  78. C
  79.       SEGMENT MPASS
  80.         REAL*8 PASS(9*NSTR)
  81.       ENDSEGMENT
  82. C
  83. C   PASS(9) : MATRICE DE CHANGEMENT DE REPERE POUR L'ELEMENT I
  84. C
  85.       SEGMENT VN(NTVN)
  86. C
  87.       SEGMENT MACHAR
  88.         INTEGER  LACHAL(NCHAR)
  89.         POINTEUR LACHAT(NCHAR).ICHATE
  90.         POINTEUR LACHAF(NCHAR).ICHAFO
  91.         REAL*8   BCHAR1(NV1)
  92.         REAL*8   BCHAR2(NV1)
  93.       ENDSEGMENT
  94. C
  95. C     LACHAL(I) : NOMBRE DE TERMES RELATIFS AU CHARGEMENT I DANS BCHAR1
  96. C     LACHAT(I) : POINTEUR SUR LA LISTE DES TEMPS DU CHARGEMENT I
  97. C     LACHAF(I) : POINTEUR SUR LA LISTE DES F(T)  DU CHARGEMENT I
  98. C     BCHAR1 : VECTEUR FORME PAR L'EMPILEMENT DE LA PARTIE CHARGEMENT
  99. C          DES SECOND-MEMBRES BLIAI, POUR TOUTES LES LIAISONS CONCERNEES
  100. C          PAR UN SOUS-CHARGEMENT ET POUR TOUS LES SOUS-CHARGEMENTS
  101. C     BCHAR2 : IDEM MODULE PAR LES FONCTIONS TEMPORELLES
  102. C
  103.       SEGMENT MOMALI
  104.         POINTEUR NOMALI(NBLIPE).MALIAI
  105.         POINTEUR KCPR2.ICPR2
  106.       ENDSEGMENT
  107. C
  108.       SEGMENT MALIAI
  109.         REAL*8  ALIA1(NALI)
  110.         REAL*8  ALIA2(NBLI)
  111.         REAL*8  ALIA3(NL/2,NL/2)
  112.         REAL*8  ALIA4(NL/2,NL/2)
  113.         REAL*8  BLIAI(NL)
  114.         REAL*8  XLIAI(NL)
  115.         INTEGER NLBLI(2,NBSTLI)
  116.         INTEGER NDCLIA(NL)
  117.         INTEGER NVNLIA(NL/2)
  118.         INTEGER IBCHA(NCHAL)
  119.         POINTEUR KWLIAI.MWLIAI
  120.         POINTEUR KLIMAS.MLIMAS
  121.       ENDSEGMENT
  122. C
  123. C     MATRICE DE LIAISON PUIS SON INVERSE EN 4 BLOCS:
  124. C     ALIA1 : MATRICE BANDE  I-A0
  125. C                               -1
  126. C     ALIA2 : MATRICE BANDE  -B0
  127. C     ALIA3 : MATRICE DE LIAISON  TERMES EN DEPLACEMENT
  128. C     ALIA4 : MATRICE DE LIAISON  TERMES EN CONTRAINTE
  129. C     NLBLI : TABLEAU DE CORRESPONDANCE ENTRE LES BLOCS
  130. C     BLIAI : VECTEUR SECOND MEMBRE -> VN ET CHARGEMENT EVENTUEL
  131. C     XLIAI : VECTEUR PREMIER MEMBRE -> DNCN
  132. C     NDCLIA(NJC)=IKID : LA NJC-IEME INCONNUE DE XLIAI EST LE IKID-IEME
  133. C     DDL DE DNCN
  134. C               DNCN(IKID)      <-  XLIAI(NJC)
  135. C     NVNLIA(NJL)=IKIV : LA NJL-IEME INCONNUE DE BLIAI EST LE IKIV-IEME
  136. C     DDL DE VN
  137. C               BLIAI(NJL)      <-  VN(IKIV)
  138. C     NCHAL : NOMBRE DE SOUS-CHARGEMENTS AGISSANT SUR LA LIAISON
  139. C     IBCHA(I)=K : LA PARTIE CHARGEMENT DE BLIAI POUR LE I-IEME SOUS-
  140. C     CHARGEMENT SE TROUVE A PARTIR DE LA K+1-IEME PLACE DANS BCHAR2
  141. C               BLIAI(NL/2+NJL) <-  BCHAR2(K+NJL)
  142. C
  143.       SEGMENT MLIMAS
  144.         REAL*8  BLIMAS(2,NJON)
  145.         REAL*8  ALIMAS(LIM,NJON)
  146.         INTEGER NLIMAS(LIM)
  147.       ENDSEGMENT
  148. C
  149. C     SEGMENT POUR LE CALCUL DE L'ACCELERATION DU CENTRE DE GRAVITE DE
  150. C     L'ELEMENT RIGIDE DE LIAISON
  151. C     BLIMAS(1,I) <--- DC(NDCLIA(I),N-1)
  152. C     BLIMAS(2,I) <--- DC(NDCLIA(I),N-2)
  153. C     NLIMAS(J)=K  : BLIAI(NL+K)=-2*C(I,K)*BLIMAS(1,I)+C(I,K)*BLIMAS(2,I
  154. C
  155.       SEGMENT MNORM
  156.         REAL*8 DNORM(LNORM)
  157.       ENDSEGMENT
  158. C
  159. C   DNORM : VECTEUR DE NORMALISATION
  160. C
  161. C
  162.       SEGMENT MCARA
  163.         REAL*8 CARA(LCAR*NSTR)
  164.       ENDSEGMENT
  165. C
  166. C   LCAR    : NOMBRE DE CARACTERISTIQUES DE L'ELEMENT
  167. C
  168.       SEGMENT ICPR(nbpts)
  169. C
  170. C                 IKID=ICPR(NUM(I,J))+K
  171. C  LA COMPOSANTE NOMD(K) DU POINT NUM EST LE IKID-IEME DDL DE DNCN( , ,1
  172. C
  173.       SEGMENT ICPR2(NIDNCN)
  174. C
  175. C     NJ=ICPR2(IKID) : LE IKID-IEME DDL DE DNCN EST LE NJ-IEME DANS
  176. C     L'ENSEMBLE DES XLIAI MIS BOUT-A-BOUT
  177. C
  178.       WRITE(IOIMP,*) ' DEBUT DE ELFRES '
  179.       MNREFE=KNREFE
  180.       MNORM=NRENO
  181.       MPASS=NREPA
  182. C
  183. C   0 -  CREATION DE L'OBJET SOLUTION DYNAMIQUE
  184. C
  185.       N=NPAS+1
  186.       NIPO=10
  187.       SEGINI MSOLUT
  188.       ITYSOL='DYNAMIQU'
  189.       SEGINI MSOLRE,MSOLEN
  190.       MSOLIS(1)=MSOLRE
  191.       MSOLIT(1)=0
  192.       MSOLIS(5)=MSOLEN
  193.       MSOLIT(5)=2
  194.       SEGDES MSOLUT
  195. C
  196. C
  197. C     -  CREATION DU CHAMPOINT SOLUTION INITIALE
  198. C
  199.       NSOUPO=1
  200.       NAT=1
  201.       SEGINI MCHPOI
  202.       IFOPOI=IFOUR
  203. C Dans les solutions il n'y que des chpo  diffus
  204.       JATTRI(1) = 1
  205.       NDDL=NREFE(5,1)
  206.       NC=2*NDDL
  207.       SEGINI MSOUPO
  208.       IPCHP(1)=MSOUPO
  209.       IGEOC=NREMEL
  210.       DO 5 I=1,NDDL
  211.       NOCOMP(I)= NOMDD(I)
  212.       NOCOMP(I+NDDL) = NOMDU(I)
  213.     5 CONTINUE
  214.       MELEME=NREMEL
  215.       N=NUM(/2)
  216.       SEGINI MPOVAL
  217.       CALL ZERO(VPOCHA,N,NC)
  218.       IPOVAL=MPOVAL
  219.       ISOLEN(1)=MCHPOI
  220. C     CALL ECCHPO(MCHPOI)
  221.       SOLRE(1)=0.D0
  222. C
  223.       SEGINI VN
  224.       NSTR=NREFE(/2)
  225.       LREF=NREFE(/1)
  226.       MANBN=KANBN
  227.       ANBN=KAB(1)
  228.       SEGACT ANBN
  229.       LANBN=AB(/2)
  230.       MDNCN=KDNCN
  231.       DNCN=KDC(1)
  232.       SEGACT DNCN
  233.       LDNCN=DC(/2)
  234.       MOMALI=KOMALI
  235.       NBLIPE=NOMALI(/1)
  236.       MACHAR=KACHAR
  237.       NCHAR=LACHAL(/1)
  238. C
  239. C    BOUCLE SUR LES PAS EN TEMPS
  240. C
  241.       NFF=0
  242.       NFS=0
  243.       NPAS1=NPAS+1
  244.       DO 70 NN=2,NPAS1
  245.       NFF=NFF+1
  246.       MNN=MIN(NN,M)
  247.       IF (IIMPI.EQ.1) THEN
  248.          WRITE(IOIMP,*) 'NN=',NN,'MNN=',MNN
  249.       END IF
  250.       DO 10 I=1,NTVN
  251.       VN(I)=0.D0
  252.   10  CONTINUE
  253. C
  254. C
  255. C   1 -  CALCUL DE LA CONVOLUTION, REMPLISSAGE DE VN
  256. C
  257. C   LA CONVOLUTION SE FAIT DANS L ORDRE CROISSANT DES INDICES POUR DNCN
  258. C   MNN-1 : LONGUEUR SUR LAQUELLE SE FAIT LA CONVOLUTION
  259. C   IDEP : INDICE DE DEPART POUR LE CALCUL DE LA CONVOLUTION
  260. C   IARR : INDICE D ARRIVEE POUR LE CALCUL DE LA CONVOLUTION
  261. C   INDICES DE DEPART ET D'ARRIVEE A L INTERIEUR DES BLOCS:
  262. C   POUR ANBN : DEPART NDANBN, ARRIVEE NAANBN
  263. C   POUR DNCN : DEPART NDDNCN, ARRIVEE NADNCN
  264. C   LAB : NUMERO DU BLOC ANBN
  265. C   LDC : NUMERO DU BLOC DNCN
  266. C
  267.       IDEP=NN-(MNN-1)
  268.       IARR=NN-1
  269. C     WRITE (IOIMP,*) 'IDEP=',IDEP,'IARR=',IARR
  270.       J=MNN
  271.       LAB=((J-1)/LANBN)+1
  272. C     LDC=((IDEP-1)/LDNCN)+1
  273. C     DECALAGE DES BLOCS DNCN LORSQUE LE PREMIER NE SERT PLUS
  274.       IF(IDEP.NE.1.AND.MOD(IDEP,LDNCN).EQ.1) THEN
  275.         KKK=KDC(1)
  276.         NDCL=KDC(/1)-1
  277.         DO 1 LDC1=1,NDCL
  278.         KDC(LDC1)=KDC(LDC1+1)
  279.     1   CONTINUE
  280.         KDC(NDCL+1)=KKK
  281. C     WRITE (IOIMP,*) 'DECALAGE IDEP ',IDEP
  282.       ENDIF
  283.       LDC=1
  284. C     WRITE (IOIMP,*) 'LAB=',LAB,'LDC=',LDC
  285.       ANBN=KAB(LAB)
  286.       DNCN=KDC(LDC)
  287.       SEGACT ANBN
  288.       SEGACT DNCN
  289. C
  290. C  BOUCLE SUR LES TERMES DU PRODUIT DE CONVOLUTION
  291. C
  292.       NDDNCN=MOD(IDEP,LDNCN)
  293.       IF(NDDNCN.EQ.0) NDDNCN=LDNCN
  294.       NDANBN=MOD(J,LANBN)
  295.       IF(NDANBN.EQ.0) NDANBN=LANBN
  296.       DO 40 I=IDEP,IARR
  297.       NADNCN=MOD(I,LDNCN)
  298.       NAANBN=MOD(J,LANBN)
  299. C
  300. C  PRODUIT DE CONVOLUTION PARTIEL UTILISANT LES DEUX BLOCS EN MEMOIRE
  301. C
  302.       IF(NADNCN.EQ.0.OR.NAANBN.EQ.1.OR.I.EQ.IARR) THEN
  303.         IF(NADNCN.EQ.0) NADNCN=LDNCN
  304.         IF(NAANBN.EQ.0) NAANBN=LANBN
  305. C       WRITE (IOIMP,*) 'NDANBN=',NDANBN,'NAANBN=',NAANBN
  306. C       WRITE (IOIMP,*) 'NDDNCN=',NDDNCN,'NADNCN=',NADNCN
  307. C       WRITE (IOIMP,*) I,LAB,ANBN,LDC,DNCN
  308.       CALL ELFRE1(VN(1),AB,DC,NREFE,NSTR,NTANBN,NIDNCN,NDANBN,NDDNCN,
  309.      *NADNCN,LREF)
  310. C       WRITE(IOIMP,*)'VN'
  311. C       WRITE(IOIMP,*)(VN(IJ),IJ=1,VN(/1))
  312. C
  313.         NDDNCN=MOD(I+1,LDNCN)
  314.         IF(NDDNCN.EQ.0) NDDNCN=LDNCN
  315.         NDANBN=MOD(J-1,LANBN)
  316.         IF(NDANBN.EQ.0) NDANBN=LANBN
  317. C
  318. C  FIN DE BLOC DNCN
  319. C
  320.         IF(NADNCN.EQ.LDNCN) THEN
  321.           IF(I.NE.IARR) THEN
  322.             SEGDES DNCN
  323.           ELSE
  324.             SEGDES DNCN
  325.           ENDIF
  326.           IF(LDC.LT.KDC(/1)) THEN
  327.             LDC=LDC+1
  328.             DNCN=KDC(LDC)
  329. C     WRITE(IOIMP,*) KDC(/1),LDC,DNCN
  330.             SEGACT DNCN
  331.           ENDIF
  332.         ENDIF
  333. C
  334. C  FIN DE BLOC ANBN
  335. C
  336.         IF(NAANBN.EQ.1) THEN
  337.           SEGDES ANBN
  338.           IF(I.NE.IARR) THEN
  339.             LAB=LAB-1
  340.             ANBN=KAB(LAB)
  341.             SEGACT ANBN
  342.           ENDIF
  343.         ENDIF
  344.       ENDIF
  345.       J=J-1
  346.    40 CONTINUE
  347. C
  348. C
  349. C   2 -  CALCUL DU CHARGEMENT - BOUCLE SUR LES SOUS-CHARGEMENTS
  350. C
  351.       T=(NN-1)*DELTAT
  352.       IDCH=1
  353. C     WRITE(IOIMP,*) ' NCHAR =',NCHAR
  354.       DO 50 NCH=1,NCHAR
  355.       MLREE1=LACHAT(NCH)
  356.       MLREE2=LACHAF(NCH)
  357.       LON=MLREE1.PROG(/1)
  358.       CALL INTE22(MLREE1.PROG(1),MLREE2.PROG(1),T,FT,LON)
  359. C     CALL INTER1(LACHAT(NCH),LACHAF(NCH),T,FT)
  360.       LCH=LACHAL(NCH)
  361. C     WRITE(IOIMP,*) ' LCH ',LCH,' NCH ',NCH
  362.       CALL MULVEC(BCHAR2(IDCH),BCHAR1(IDCH),FT,LCH)
  363.       IDCH=IDCH+LCH
  364.    50 CONTINUE
  365. C     WRITE(IOIMP,*) ' BCHAR2 '
  366. C     WRITE(IOIMP,1002) (BCHAR2(IJ),IJ=1,BCHAR2(/1))
  367. C
  368. C   3 -  CALCUL DE DNCN  - BOUCLE SUR LES LIAISONS
  369. C
  370.       DO 60 NBLI=1,NBLIPE
  371.       MALIAI=NOMALI(NBLI)
  372.       NL=NDCLIA(/1)
  373.       NLS2=NL/2
  374. C
  375. C  RANGEMENT DES SECONDS MEMBRES VN DANS LES NL PREMIERES
  376. C  COMPOSANTES DU VECTEUR BLIAI DE LA NBLI- EME LIAISON
  377. C
  378.       CALL ELFBL1(BLIAI,VN(1),NVNLIA,NLS2)
  379. C     WRITE(IOIMP,*) ' BLIAI 1 '
  380. C     WRITE(IOIMP,1002) (BLIAI(IJ),IJ=1,NLS2)
  381. C
  382. C  CALCUL DE LA RESULTANTE DES NCHAL CHARGEMENTS CONTENUS DANS BCHAR2
  383. C  QUI AGISSENT SUR LA NBLI-IEME LIAISON,
  384. C  ET RANGEMENT DANS LA DEUXIEME MOITIE DU VECTEUR SECOND MEMBRE BLIAI
  385. C  A L'AIDE DU TABLEAU IBCHA
  386. C
  387.       NCHAL=IBCHA(/1)
  388.       NLS21=NLS2+1
  389.       IF(NCHAL.NE.0) CALL ELFBL2(BLIAI(NLS21),BCHAR2,IBCHA,NCHAL,NLS2)
  390.       IF(KLIMAS.NE.0) THEN
  391. C
  392. C  CALCUL DE L'ACCELERATION DU CENTRE DE GRAVITE DE L'ELEMENT DE LIAISON
  393. C
  394.         MLIMAS=KLIMAS
  395.         LIM=NLIMAS(/1)
  396.         CALL ELFBL6(BLIAI(NLS21),BLIMAS,ALIMAS,NLIMAS,NLS2,LIM)
  397.       ENDIF
  398. C     WRITE(IOIMP,*) ' BLIAI 2 '
  399. C     WRITE(IOIMP,1002) (BLIAI(IJ),IJ=NLS21,NL)
  400. C
  401. C  MULTIPLICATION ALIAI*BLIAI=XLIAI
  402. C
  403.       NBSTLI=NLBLI(/2)
  404.       CALL MULMAV(ALIA1,ALIA2,ALIA3,ALIA4,BLIAI,XLIAI,NLS2,NCHAL,NLBLI,
  405.      *NBSTLI)
  406. C     WRITE(IOIMP,*) ' XLIAI '
  407. C     WRITE(IOIMP,1002) (XLIAI(IJ),IJ=1,NL)
  408.       IF(KLIMAS.NE.0) THEN
  409. C
  410. C  RANGEMENT DE XLIAI DANS BLIMAS
  411. C
  412.         CALL ELFBL5(BLIMAS,XLIAI,NLS2)
  413.       ENDIF
  414. C
  415. C  RANGEMENT DE XLIAI DANS DNCN
  416. C
  417.       CALL ELFBL3(DC(1,NDDNCN),XLIAI,NDCLIA,NL)
  418. C     WRITE(IOIMP,*) ' DNCN  '
  419. C     WRITE(IOIMP,1002) (DC(IJ,NDDNCN),IJ=1,NIDNCN)
  420. C
  421.       IF(IIMPI.EQ.1) THEN
  422.         WRITE(IOIMP,1001) NBLI
  423.  1001   FORMAT(//I5,' IEME SECOND MEMBRE BLIAI'//)
  424.         WRITE (IOIMP,1002) (BLIAI(L),L=1,NL)
  425.  1002   FORMAT(1X,10(1PE12.5,1X))
  426.         WRITE(IOIMP,1003)
  427.  1003   FORMAT(//' SOLUTION XLIAI'//)
  428.         WRITE (IOIMP,1002) (XLIAI(L),L=1,NL)
  429.       ENDIF
  430.    60 CONTINUE
  431. C
  432. C   4 -  CREATION DU CHAMPOINT SOLUTION A L'INSTANT T
  433. C
  434.       IF(NFF.EQ.NFOIS) THEN
  435.         NFF=0
  436.         NFS=NFS+1
  437.         NAT=1
  438.         SEGINI MCHPOI
  439.         IFOPOI=IFOUR
  440.         JATTRI(1) = 1
  441.         SEGINI MSOUPO
  442.         IPCHP(1)=MSOUPO
  443.         IGEOC=NREMEL
  444.         DO 65 I=1,NDDL
  445.         NOCOMP(I)= NOMDD(I)
  446.         NOCOMP(I + NDDL)= NOMDU(I)
  447.    65   CONTINUE
  448.         SEGINI MPOVAL
  449.         IPOVAL=MPOVAL
  450.         CALL ELFBL4(N,DC(1,NDDNCN),DNORM,PASS,VPOCHA)
  451.         ISOLEN(NFS)=MCHPOI
  452. C       CALL ECCHPO(MCHPOI)
  453.         SOLRE(NFS)=T
  454.       ENDIF
  455. C
  456.    70 CONTINUE
  457. C
  458. C   5 -  SUPPRESSION DES SEGMENTS DE TRAVAIL ET DESACTIVATION
  459. C
  460.       SEGDES MSOLRE,MSOLEN
  461. C
  462.       SEGSUP VN
  463. C
  464.       NSGA=KAB(/1)
  465.       DO 80 NSG=1,NSGA
  466.       ANBN=KAB(NSG)
  467.       SEGSUP ANBN
  468.    80 CONTINUE
  469.       NSGD=KDC(/1)
  470.       DO 90 NSG=1,NSGD
  471.       DNCN=KDC(NSG)
  472.       SEGSUP DNCN
  473.    90 CONTINUE
  474.       SEGSUP MANBN
  475.       SEGSUP MDNCN
  476. C
  477.       DO 100 NCH=1,NCHAR
  478.       MLREE1=LACHAT(NCH)
  479.       MLREE2=LACHAF(NCH)
  480.       SEGDES MLREE1,MLREE2
  481.   100 CONTINUE
  482.       SEGSUP MACHAR
  483. C
  484.       DO 110 NBLI=1,NBLIPE
  485.       MALIAI=NOMALI(NBLI)
  486.       IF(KLIMAS.NE.0) THEN
  487.         MLIMAS=KLIMAS
  488.         SEGSUP MLIMAS
  489.       ENDIF
  490.       SEGSUP MALIAI
  491.   110 CONTINUE
  492.       ICPR2=KCPR2
  493.       SEGSUP ICPR2
  494.       SEGSUP MOMALI
  495. C
  496.       SEGSUP MNORM
  497.       SEGSUP MPASS
  498.       ICPR=NRECPR
  499.       MCARA=NRECA
  500.       SEGSUP ICPR
  501.       SEGSUP MCARA
  502.       SEGSUP MNREFE
  503.       KSOLUT=MSOLUT
  504.       WRITE (IOIMP,*) 'FIN DE ELFRES '
  505.       RETURN
  506.       END
  507.  
  508.  
  509.  
  510.  
  511.  
  512.  
  513.  
  514.  
  515.  

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