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

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