Télécharger elfres.eso

Retour à la liste

Numérotation des lignes :

  1. C ELFRES SOURCE BP208322 15/06/22 21:17:42 8543
  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. -INC CCOPTIO
  16. -INC CCHAMP
  17. -INC SMELEME
  18. -INC SMCHPOI
  19. -INC SMSOLUT
  20. -INC SMLREEL
  21. C
  22. SEGMENT MANBN
  23. POINTEUR KAB(NSGA).ANBN
  24. ENDSEGMENT
  25. C
  26. C NSGA : NOMBRE DE BLOCS ANBN
  27. C
  28. SEGMENT ANBN
  29. REAL*8 AB(NTANBN,LANBN)
  30. ENDSEGMENT
  31. C
  32. C AB(I,K) : TERME I DE LA MATRICE A OU B D'UN ELEMENT AU TEMPS K
  33. C NTANBN : NOMBRE DE TERMES DES MATRICES A ET B DE TOUS LES ELEMENTS
  34. C LANBN : NOMBRE DE PAS DE TEMPS STOCKES DANS UN BLOC ANBN
  35. C
  36. SEGMENT MDNCN
  37. POINTEUR KDC(NSGD).DNCN
  38. ENDSEGMENT
  39. C
  40. C NSGD : NOMBRE DE BLOCS DNCN
  41. C
  42. SEGMENT DNCN
  43. REAL*8 DC(NIDNCN,LDNCN)
  44. ENDSEGMENT
  45. C
  46. C DC(I,K) : DDL I AU TEMPS K
  47. C NIDNCN : NOMBRE TOTAL D'INCONNUES
  48. C LDNCN : NOMBRE DE PAS DE TEMPS STOCKES DANS UN BLOC DNCN
  49. C
  50. SEGMENT MNREFE
  51. INTEGER NREFE(8,NSTR)
  52. INTEGER NTANBN
  53. INTEGER NIDNCN
  54. INTEGER NTVN
  55. POINTEUR NREPA.MPASS
  56. POINTEUR NRECA.MCARA
  57. POINTEUR NRENO.MNORM
  58. POINTEUR NRECPR.ICPR
  59. POINTEUR NREMEL.MELEME
  60. POINTEUR NREDEN.MDEN
  61. ENDSEGMENT
  62. C
  63. C NSTR : NOMBRE D'ELEMENTS
  64. C NREFE(1,I) : MELEME
  65. C NREFE(2,I) : MSOSTU
  66. C NREFE(3,I) : TYPE DE L'ELEMENT
  67. C NREFE(4,I) : NOMBRE DE POINTS DU MELEME
  68. C NREFE(5,I) : NOMBRE DE DDL PAR POINT
  69. C NREFE(6,I)=IVN :LE 1ER DDL DE L'ELEMENT EST LE IVN+1 IEME DE VN
  70. C NREFE(7,I)=IAN :LE 1ER TERME DE LA MATRICE A EST LE IAN IEME DE ANBN
  71. C NREFE(8,I)= 1 :LE IEME ELEMENT EST RIGIDE (OU PARTIELLEMENT) SINON 0
  72. C NTANBN : NOMBRE DE TERMES DES MATRICES A ET B POUR TOUS LES ELEMENTS
  73. C NIDNCN : NOMBRE TOTAL D'INCONNUES DE DNCN
  74. C NTVN : LONGUEUR DU TABLEAU VN
  75. C
  76. SEGMENT MPASS
  77. REAL*8 PASS(9*NSTR)
  78. ENDSEGMENT
  79. C
  80. C PASS(9) : MATRICE DE CHANGEMENT DE REPERE POUR L'ELEMENT I
  81. C
  82. SEGMENT VN(NTVN)
  83. C
  84. SEGMENT MACHAR
  85. INTEGER LACHAL(NCHAR)
  86. POINTEUR LACHAT(NCHAR).ICHATE
  87. POINTEUR LACHAF(NCHAR).ICHAFO
  88. REAL*8 BCHAR1(NV1)
  89. REAL*8 BCHAR2(NV1)
  90. ENDSEGMENT
  91. C
  92. C LACHAL(I) : NOMBRE DE TERMES RELATIFS AU CHARGEMENT I DANS BCHAR1
  93. C LACHAT(I) : POINTEUR SUR LA LISTE DES TEMPS DU CHARGEMENT I
  94. C LACHAF(I) : POINTEUR SUR LA LISTE DES F(T) DU CHARGEMENT I
  95. C BCHAR1 : VECTEUR FORME PAR L'EMPILEMENT DE LA PARTIE CHARGEMENT
  96. C DES SECOND-MEMBRES BLIAI, POUR TOUTES LES LIAISONS CONCERNEES
  97. C PAR UN SOUS-CHARGEMENT ET POUR TOUS LES SOUS-CHARGEMENTS
  98. C BCHAR2 : IDEM MODULE PAR LES FONCTIONS TEMPORELLES
  99. C
  100. SEGMENT MOMALI
  101. POINTEUR NOMALI(NBLIPE).MALIAI
  102. POINTEUR KCPR2.ICPR2
  103. ENDSEGMENT
  104. C
  105. SEGMENT MALIAI
  106. REAL*8 ALIA1(NALI)
  107. REAL*8 ALIA2(NBLI)
  108. REAL*8 ALIA3(NL/2,NL/2)
  109. REAL*8 ALIA4(NL/2,NL/2)
  110. REAL*8 BLIAI(NL)
  111. REAL*8 XLIAI(NL)
  112. INTEGER NLBLI(2,NBSTLI)
  113. INTEGER NDCLIA(NL)
  114. INTEGER NVNLIA(NL/2)
  115. INTEGER IBCHA(NCHAL)
  116. POINTEUR KWLIAI.MWLIAI
  117. POINTEUR KLIMAS.MLIMAS
  118. ENDSEGMENT
  119. C
  120. C MATRICE DE LIAISON PUIS SON INVERSE EN 4 BLOCS:
  121. C ALIA1 : MATRICE BANDE I-A0
  122. C -1
  123. C ALIA2 : MATRICE BANDE -B0
  124. C ALIA3 : MATRICE DE LIAISON TERMES EN DEPLACEMENT
  125. C ALIA4 : MATRICE DE LIAISON TERMES EN CONTRAINTE
  126. C NLBLI : TABLEAU DE CORRESPONDANCE ENTRE LES BLOCS
  127. C BLIAI : VECTEUR SECOND MEMBRE -> VN ET CHARGEMENT EVENTUEL
  128. C XLIAI : VECTEUR PREMIER MEMBRE -> DNCN
  129. C NDCLIA(NJC)=IKID : LA NJC-IEME INCONNUE DE XLIAI EST LE IKID-IEME
  130. C DDL DE DNCN
  131. C DNCN(IKID) <- XLIAI(NJC)
  132. C NVNLIA(NJL)=IKIV : LA NJL-IEME INCONNUE DE BLIAI EST LE IKIV-IEME
  133. C DDL DE VN
  134. C BLIAI(NJL) <- VN(IKIV)
  135. C NCHAL : NOMBRE DE SOUS-CHARGEMENTS AGISSANT SUR LA LIAISON
  136. C IBCHA(I)=K : LA PARTIE CHARGEMENT DE BLIAI POUR LE I-IEME SOUS-
  137. C CHARGEMENT SE TROUVE A PARTIR DE LA K+1-IEME PLACE DANS BCHAR2
  138. C BLIAI(NL/2+NJL) <- BCHAR2(K+NJL)
  139. C
  140. SEGMENT MLIMAS
  141. REAL*8 BLIMAS(2,NJON)
  142. REAL*8 ALIMAS(LIM,NJON)
  143. INTEGER NLIMAS(LIM)
  144. ENDSEGMENT
  145. C
  146. C SEGMENT POUR LE CALCUL DE L'ACCELERATION DU CENTRE DE GRAVITE DE
  147. C L'ELEMENT RIGIDE DE LIAISON
  148. C BLIMAS(1,I) <--- DC(NDCLIA(I),N-1)
  149. C BLIMAS(2,I) <--- DC(NDCLIA(I),N-2)
  150. C NLIMAS(J)=K : BLIAI(NL+K)=-2*C(I,K)*BLIMAS(1,I)+C(I,K)*BLIMAS(2,I
  151. C
  152. SEGMENT MNORM
  153. REAL*8 DNORM(LNORM)
  154. ENDSEGMENT
  155. C
  156. C DNORM : VECTEUR DE NORMALISATION
  157. C
  158. C
  159. SEGMENT MCARA
  160. REAL*8 CARA(LCAR*NSTR)
  161. ENDSEGMENT
  162. C
  163. C LCAR : NOMBRE DE CARACTERISTIQUES DE L'ELEMENT
  164. C
  165. SEGMENT ICPR(XCOOR(/1)/(IDIM+1))
  166. C
  167. C IKID=ICPR(NUM(I,J))+K
  168. C LA COMPOSANTE NOMD(K) DU POINT NUM EST LE IKID-IEME DDL DE DNCN( , ,1
  169. C
  170. SEGMENT ICPR2(NIDNCN)
  171. C
  172. C NJ=ICPR2(IKID) : LE IKID-IEME DDL DE DNCN EST LE NJ-IEME DANS
  173. C L'ENSEMBLE DES XLIAI MIS BOUT-A-BOUT
  174. C
  175. WRITE(IOIMP,*) ' DEBUT DE ELFRES '
  176. MNREFE=KNREFE
  177. MNORM=NRENO
  178. MPASS=NREPA
  179. C
  180. C 0 - CREATION DE L'OBJET SOLUTION DYNAMIQUE
  181. C
  182. N=NPAS+1
  183. NIPO=10
  184. SEGINI MSOLUT
  185. ITYSOL='DYNAMIQU'
  186. SEGINI MSOLRE,MSOLEN
  187. MSOLIS(1)=MSOLRE
  188. MSOLIT(1)=0
  189. MSOLIS(5)=MSOLEN
  190. MSOLIT(5)=2
  191. SEGDES MSOLUT
  192. C
  193. C
  194. C - CREATION DU CHAMPOINT SOLUTION INITIALE
  195. C
  196. NSOUPO=1
  197. NAT=1
  198. SEGINI MCHPOI
  199. IFOPOI=IFOMOD
  200. C Dans les solutions il n'y que des chpo diffus
  201. JATTRI(1) = 1
  202. NDDL=NREFE(5,1)
  203. NC=2*NDDL
  204. SEGINI MSOUPO
  205. IPCHP(1)=MSOUPO
  206. IGEOC=NREMEL
  207. DO 5 I=1,NDDL
  208. NOCOMP(I)= NOMDD(I)
  209. NOCOMP(I+NDDL) = NOMDU(I)
  210. 5 CONTINUE
  211. MELEME=NREMEL
  212. N=NUM(/2)
  213. SEGDES MELEME
  214. SEGINI MPOVAL
  215. CALL ZERO(VPOCHA,N,NC)
  216. IPOVAL=MPOVAL
  217. SEGDES MPOVAL,MSOUPO,MCHPOI
  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=IFOMOD
  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. SEGDES MPOVAL,MSOUPO,MCHPOI
  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
  481. SEGDES MLREE2
  482. 100 CONTINUE
  483. SEGSUP MACHAR
  484. C
  485. DO 110 NBLI=1,NBLIPE
  486. MALIAI=NOMALI(NBLI)
  487. IF(KLIMAS.NE.0) THEN
  488. MLIMAS=KLIMAS
  489. SEGSUP MLIMAS
  490. ENDIF
  491. SEGSUP MALIAI
  492. 110 CONTINUE
  493. ICPR2=KCPR2
  494. SEGSUP ICPR2
  495. SEGSUP MOMALI
  496. C
  497. SEGSUP MNORM
  498. SEGSUP MPASS
  499. ICPR=NRECPR
  500. MCARA=NRECA
  501. SEGSUP ICPR
  502. SEGSUP MCARA
  503. SEGSUP MNREFE
  504. KSOLUT=MSOLUT
  505. WRITE (IOIMP,*) 'FIN DE ELFRES '
  506. RETURN
  507. END
  508.  
  509.  
  510.  
  511.  
  512.  

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