Télécharger c2dp.eso

Retour à la liste

Numérotation des lignes :

  1. C C2DP SOURCE CHAT 05/01/12 21:44:59 5004
  2. SUBROUTINE C2DP(SIG,DSIG,YOUN,ANU,RT1,RT2,RDP,ADP,HDP,
  3. 1 XLAM1,XLAM2,XLAM3,IDAM,ANG,EPPLDP,ALPHA,KERRE)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6.  
  7. -INC PPARAM
  8. -INC CCOPTIO
  9. DIMENSION SIG(3),DSIG(3),SGG(3),DSGG(3),SFG(3),DSFG(3),
  10. 1 A(3),B(3),IDAM(3),EPC(3),SSI(3),SI(3),EPC1(3),
  11. 1 EPDP(3),EC(3),EPPLDP(3),SG(3),AT(4),BT(2),TENS(3)
  12. ITER=0
  13. Y=YOUN/(1.D0-ANU*ANU)
  14. SIREF=1.E-6*YOUN
  15. DO 10 ITYP=1,3
  16. EPDP(ITYP)=0.D0
  17. 10 EPPLDP(ITYP)=0.D0
  18. DL3=0.D0
  19. XT=0.D0
  20. IF(IIMPI.EQ.9) WRITE(IOIMP,9999)
  21. 9999 FORMAT(1X,'C2DP COUPLAGE COMP 2 DRUCKER',/)
  22. 9998 FORMAT(1X,'C2DP ',I4,'ITERATIONS INTERNES',/)
  23. C
  24. C------------------------------------------------
  25. C COUPLAGE COMPRESSION 2 DRUCKER
  26. C------------------------------------------------
  27. C----------------------------------------------
  28. C DEFINITION D UNE CONTRAINTE DE REFERENCE
  29. C SI RMAX DU DRUCKER INFERIEURE A CETTE VALEUR
  30. C ALORS SIGMA=0
  31. C----------------------------------------------
  32. C
  33. RMAX=MAX((RDP/1.73),(RDP/(1.D0-2.D0*ADP)))
  34. IF(SIREF.GT.RMAX) THEN
  35. DO 20 ITYP=1,3
  36. 20 SIG(ITYP)=SIG(ITYP)+DSIG(ITYP)
  37. RDP=0.D0
  38. HDP=0.D0
  39. ADP=0.D0
  40. CALL CPHOIN(SIG,EPPLDP,YOUN,ANU)
  41. CALL NORME(SIG,DL3)
  42. DL3=DL3/YOUN*10.D0
  43. XLAM3=XLAM3+DL3
  44. DO 30 ITYP=1,3
  45. IDAM(ITYP)=0
  46. DSIG(ITYP)=0.D0
  47. 30 SIG(ITYP)=0.D0
  48. IF(IIMPI.EQ.9) WRITE(IOIMP,9998) ITER
  49. RETURN
  50. ENDIF
  51. C-----------------------------------------------
  52. C ON SEPLACE DANS LE REPERE DE FISSURATION
  53. C-----------------------------------------------
  54. C
  55. CALL CHREP(SIG,SFG,ANG)
  56. CALL CHREP(DSIG,DSFG,ANG)
  57. C
  58. C-----------------------------------------------------
  59. C ESTIMATION DU PAS D'INCREMENT DE CONTRAINTES
  60. C-----------------------------------------------------
  61. C
  62. CALL SIJ(SFG,SI,SIETOI)
  63. CALL EPSDP(SI,SIETOI,ADP,EPC)
  64. CALL CPHOMO(EPC,SSI,YOUN,ANU,ALPHA)
  65. CALL SCAL2(EPC,SSI,VAL)
  66. VAL=VAL-HDP
  67. IF(VAL.LT.0.D0) THEN
  68. KERRE=459
  69. RETURN
  70. ENDIF
  71. CALL SCAL(SSI,DSFG,VAL)
  72. CALL NORME(DSFG,VA1)
  73. CALL NORME(SSI,VA2)
  74. IF(VA1.EQ.0.D0) THEN
  75. IDAM(1)=0
  76. IDAM(2)=0
  77. IDAM(3)=0
  78. RETURN
  79. ENDIF
  80. RMIN=MIN((RDP/1.73),(RDP/(1.D0+2.D0*ADP)))
  81. X=VAL/VA1/VA2
  82. IF(X.GE.1.D0) THEN
  83. X=1.D0
  84. GOTO 41
  85. ENDIF
  86. X=1.D0/SQRT(1.0001D0-X*X)*RMIN/VA1/8.D0
  87. IF(X.GE.1.D0) X=1.D0
  88. C
  89. C----------------------------------------------------------
  90. C ON ECOULE ET ON REGARDE LA VARIATION DE NORMALE
  91. C----------------------------------------------------------
  92. C
  93. C
  94. C---------------------------------------
  95. C RESOLUTION DU SYSTEME EN DL3
  96. C---------------------------------------
  97. 41 X=X/2.D0
  98. EC(1)=0.D0
  99. EC(2)=1.D0
  100. EC(3)=0.D0
  101. CALL CPHOOB(EC,SG,YOUN,ANU)
  102. 100 X=X*2.D0
  103. CALL SIJ(SFG,SI,SIETOI)
  104. CALL EPSDP(SI,SIETOI,ADP,EPC)
  105. CALL CPHOMO(EPC,SSI,YOUN,ANU,ALPHA)
  106. CALL SCAL2(EPC,SSI,VAL)
  107. VAL=VAL-HDP
  108. IF(VAL.LT.0.D0) THEN
  109. KERRE=459
  110. RETURN
  111. ENDIF
  112. IF((XT+X).GT.1.D0) X=1.D0-XT
  113. 110 RMAX=MAX((RDP/1.73),(RDP/(1.D0-2.D0*ADP)))
  114. C
  115. C --------------------------------------------------
  116. C CAS OU LE RAYON MAXIMAL DU DRUCKER-PRAGER
  117. C INFERIEUR A LA CONTRAINTE DE REFERENCE
  118. C --------------------------------------------------
  119. C
  120. IF(SIREF.LT.RMAX) GOTO 131
  121. DO 120 ITYP=1,3
  122. 120 SFG(ITYP)=SFG(ITYP)+(1.D0-XT)*DSFG(ITYP)
  123. RDP=0.D0
  124. ADP=0.D0
  125. HDP=0.D0
  126. CALL NORME(SFG,DL3)
  127. DL3=DL3/YOUN*10.D0
  128. XLAM3=XLAM3+DL3
  129. CALL CPHOIN(SFG,EPPLDP,YOUN,ANU)
  130. DO 130 ITYP=1,3
  131. IDAM(ITYP)=0
  132. EPDP(ITYP)=EPDP(ITYP)+EPPLDP(ITYP)
  133. DSIG(ITYP)=0.D0
  134. 130 SIG(ITYP)=0.D0
  135. CALL CHREP(EPDP,EPPLDP,-ANG)
  136. IF(IIMPI.EQ.9) WRITE(IOIMP,9998) ITER
  137. RETURN
  138. 131 ITER=ITER+1
  139. IF(ITER.GT.200) THEN
  140. KERRE=460
  141. RETURN
  142. ENDIF
  143. C
  144. C---------------------------------
  145. C ON REALISE LE COUPLAGE
  146. C---------------------------------
  147. C
  148. E=-(SFG(2)+X*DSFG(2))/SG(2)
  149. F=SSI(2)/SG(2)
  150. DO 140 ITYP=1,3
  151. DSGG(ITYP)=X*DSFG(ITYP)
  152. A(ITYP)=SFG(ITYP)+DSGG(ITYP)+E*SG(ITYP)
  153. 140 B(ITYP)=-SSI(ITYP)+F*SG(ITYP)
  154. CALL DLAMDP(A,B,DL3,RDP,ADP,HDP,ITEST)
  155. IF(ITEST.EQ.1) THEN
  156. X=X/2.D0
  157. GOTO 131
  158. ENDIF
  159. DL2=-E-F*DL3
  160. DO 145 ITYP=1,3
  161. 145 SGG(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  162. C
  163. C----------------------------------------------
  164. C SI UN DES DLAMDA < 0 FAUX COUPLAGE
  165. C----------------------------------------------
  166. C
  167. IF(DL2.GT.0.D0) IDAM(2)=0
  168. IF(DL3.LT.0.D0) IDAM(3)=0
  169. IF(IDAM(2).EQ.0.OR.IDAM(3).EQ.0) THEN
  170. DL2=0.D0
  171. DL3=0.D0
  172. X=0.D0
  173. GOTO 3000
  174. ENDIF
  175. C
  176. C-----------------------------------------------------
  177. C ON VERIFIE QUE LA NORMALE SUIVANT LE DRUCKER
  178. C PRAGER TOURNE PEU
  179. C-----------------------------------------------------
  180. C
  181. CALL SIJ(SGG,SI,SIETOI)
  182. CALL EPSDP(SI,SIETOI,ADP,EPC1)
  183. CALL SCAL(EPC,EPC1,VAL)
  184. CALL NORME(EPC,VA1)
  185. CALL NORME(EPC1,VA2)
  186. CO=VAL/VA1/VA2
  187. IF(CO.LT.0.99) THEN
  188. X=X/2.D0
  189. GOTO 131
  190. ENDIF
  191. IF((DL2+XLAM2).GE.0.D0) GOTO 200
  192. C
  193. C------------------------------------------------
  194. C CAS OU LA FISSURE EST COMPLETEMENT FERMEE
  195. C------------------------------------------------
  196. C
  197. IDAM(2)=0
  198. DL2=-XLAM2
  199. IF(DSFG(2).EQ.0.D0) GOTO 160
  200. E=(DL2*SG(2)-SFG(2))/DSFG(2)
  201. F=SSI(2)/DSFG(2)
  202. DO 150 ITYP=1,3
  203. A(ITYP)=SFG(ITYP)-DL2*SG(ITYP)+E*DSFG(ITYP)
  204. 150 B(ITYP)=F*DSFG(ITYP)-SSI(ITYP)
  205. CALL DLAMD(A,B,DL31,DL32,RDP,ADP,HDP,ITEST)
  206. DO 151 ITYP=1,3
  207. SGG(ITYP)=A(ITYP)+DL31*B(ITYP)
  208. 151 DSGG(ITYP)=A(ITYP)+DL32*B(ITYP)
  209. CALL SIJ(SGG,SI,SIETOI)
  210. CALL EPSDP(SI,SIETOI,ADP,EPC1)
  211. CALL SCAL(EPC,EPC1,VAL)
  212. CALL NORME(EPC,VA1)
  213. CALL NORME(EPC1,VA2)
  214. CO1=VAL/VA1/VA2
  215. CALL SIJ(DSGG,SI,SIETOI)
  216. CALL EPSDP(SI,SIETOI,ADP,EPC1)
  217. CALL SCAL(EPC,EPC1,VAL)
  218. CALL NORME(EPC,VA1)
  219. CALL NORME(EPC1,VA2)
  220. CO2=VAL/VA1/VA2
  221. XIN1=E+F*DL31
  222. XIN2=E+F*DL32
  223. IF(XIN1.GT.-1.D-10.AND.XIN1.LE.X.AND.CO1.GT.0.9) THEN
  224. DL3=DL31
  225. X=XIN1
  226. GOTO 200
  227. ENDIF
  228. IF(XIN2.GT.-1.D-10.AND.XIN2.LE.X.AND.CO2.GT.0.9) THEN
  229. DL3=DL32
  230. X=XIN2
  231. GOTO 200
  232. ENDIF
  233. IF(IIMPI.EQ.9) WRITE(IOIMP,10101)
  234. 10101 FORMAT(1X,'ERREUR DANS C2DP FERM ',/)
  235. DL3=DL31
  236. X=XIN1
  237. GOTO 200
  238. 160 DL3=(SFG(2)-DL2*SG(2))/SSI(2)
  239. RD=RDP-HDP*DL3
  240. IF(RD.LT.0.D0) THEN
  241. RDP=0.D0
  242. GOTO 110
  243. ENDIF
  244. DO 170 ITYP=1,3
  245. A(ITYP)=SFG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  246. 170 B(ITYP)=DSFG(ITYP)
  247. CALL XDP(A,B,RD,ADP,X,ITEST)
  248. 200 CONTINUE
  249. DO 205 ITYP=1,3
  250. DSGG(ITYP)=X*DSFG(ITYP)
  251. 205 SGG(ITYP)=SFG(ITYP)+DSGG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  252. C
  253. C-----------------------------------------------
  254. C ON REGARDE SI ON ENDOMMAGE PAS LES AUTRES
  255. C CRITERES
  256. C-----------------------------------------------
  257. C-------------------------------------
  258. C ON ENDOMMAGE LA TRACTION 1
  259. C-------------------------------------
  260. C
  261. CALL CTRAF(SGG(1),RT1,VCTR)
  262. IF(VCTR.LE.0.D0) GOTO 1000
  263. IDAM(1)=1
  264. IDAM(2)=-1
  265. IDAM(3)=1
  266. AT(1)=-DSFG(1)
  267. AT(4)=Y
  268. AT(2)=Y*ANU
  269. AT(3)=-DSFG(2)
  270. BT(1)=SFG(1)-RT1
  271. BT(2)=SFG(2)
  272. CALL SYLIN2(AT,BT,X1,DL21)
  273. BT(1)=-SSI(1)
  274. BT(2)=-SSI(2)
  275. CALL SYLIN2(AT,BT,X2,DL22)
  276. DO 230 ITYP=1,3
  277. A(ITYP)=SFG(ITYP)-DL21*SG(ITYP)+X1*DSFG(ITYP)
  278. 230 B(ITYP)=-DL22*SG(ITYP)+X2*DSFG(ITYP)-SSI(ITYP)
  279. CALL DLAMD(A,B,DL31,DL32,RDP,ADP,HDP,ITEST)
  280. DO 231 ITYP=1,3
  281. SGG(ITYP)=A(ITYP)+DL31*B(ITYP)
  282. 231 DSGG(ITYP)=A(ITYP)+DL32*B(ITYP)
  283. CALL SIJ(SGG,SI,SIETOI)
  284. CALL EPSDP(SI,SIETOI,ADP,EPC1)
  285. CALL SCAL(EPC,EPC1,VAL)
  286. CALL NORME(EPC,VA1)
  287. CALL NORME(EPC1,VA2)
  288. CO1=VAL/VA1/VA2
  289. CALL SIJ(DSGG,SI,SIETOI)
  290. CALL EPSDP(SI,SIETOI,ADP,EPC1)
  291. CALL SCAL(EPC,EPC1,VAL)
  292. CALL NORME(EPC,VA1)
  293. CALL NORME(EPC1,VA2)
  294. CO2=VAL/VA1/VA2
  295. DLIN1=DL21+DL31*DL22
  296. DLIN2=DL21+DL32*DL22
  297. XIN1=X1+DL31*X2
  298. XIN2=X1+DL32*X2
  299. IF(XIN1.GT.-1.D-10.AND.DLIN1.LT.1.E-10.AND.XIN1.LE.X
  300. 1 .AND.CO1.GT.0.9) THEN
  301. DL3=DL31
  302. X=XIN1
  303. DL2=DLIN1
  304. GOTO 3000
  305. ENDIF
  306. IF(XIN2.GT.-1.D-10.AND.DLIN2.LT.1.E-10.AND.XIN2.LE.X
  307. 1 .AND.CO2.GT.0.9) THEN
  308. DL3=DL32
  309. X=XIN2
  310. DL2=DLIN2
  311. GOTO 3000
  312. ENDIF
  313. IF(IIMPI.EQ.9) WRITE(IOIMP,20202)
  314. 20202 FORMAT(1X,'ERREUR DANS C2DP TRAC ',/)
  315. DL3=DL31
  316. DL2=DLIN1
  317. X=XIN1
  318. GOTO 3000
  319. C
  320. C-------------------------------------
  321. C ON ENDOMMAGE LA COMPRESSION 1
  322. C-------------------------------------
  323. C
  324. 1000 CALL CCOAF(SGG(1),XLAM1,VCCO1)
  325. IF(VCCO1.LE.0.D0) GOTO 2000
  326. IDAM(1)=-1
  327. IDAM(3)=1
  328. IDAM(2)=-1
  329. AT(1)=-DSFG(1)
  330. AT(4)=Y
  331. AT(2)=Y*ANU
  332. AT(3)=-DSFG(2)
  333. BT(1)=SFG(1)
  334. BT(2)=SFG(2)
  335. CALL SYLIN2(AT,BT,X1,DL21)
  336. BT(1)=-SSI(1)
  337. BT(2)=-SSI(2)
  338. CALL SYLIN2(AT,BT,X2,DL22)
  339. DO 1010 ITYP=1,3
  340. A(ITYP)=SFG(ITYP)-DL21*SG(ITYP)+X1*DSFG(ITYP)
  341. 1010 B(ITYP)=-DL22*SG(ITYP)+X2*DSFG(ITYP)-SSI(ITYP)
  342. CALL DLAMD(A,B,DL31,DL32,RDP,ADP,HDP,ITEST)
  343. DO 1011 ITYP=1,3
  344. SGG(ITYP)=A(ITYP)+DL31*B(ITYP)
  345. 1011 DSGG(ITYP)=A(ITYP)+DL32*B(ITYP)
  346. CALL SIJ(SGG,SI,SIETOI)
  347. CALL EPSDP(SI,SIETOI,ADP,EPC1)
  348. CALL SCAL(EPC,EPC1,VAL)
  349. CALL NORME(EPC,VA1)
  350. CALL NORME(EPC1,VA2)
  351. CO1=VAL/VA1/VA2
  352. CALL SIJ(DSGG,SI,SIETOI)
  353. CALL EPSDP(SI,SIETOI,ADP,EPC1)
  354. CALL SCAL(EPC,EPC1,VAL)
  355. CALL NORME(EPC,VA1)
  356. CALL NORME(EPC1,VA2)
  357. CO2=VAL/VA1/VA2
  358. * DLIN1=DL11+DL31*DL12 CORRECTION PV ???
  359. * DLIN2=DL11+DL32*DL12 SUITE ERREUR APPOLO
  360. DLIN1=DL21+DL31*DL22
  361. DLIN2=DL21+DL32*DL22
  362. XIN1=X1+DL31*X2
  363. XIN2=X1+DL32*X2
  364. IF(XIN1.GT.-1.D-10.AND.DLIN1.LT.1.E-10.AND.XIN1.LE.X
  365. 1 .AND.CO1.GT.0.9) THEN
  366. DL3=DL31
  367. X=XIN1
  368. DL2=DLIN1
  369. GOTO 3000
  370. ENDIF
  371. IF(XIN2.GT.-1.D-10.AND.DLIN2.LT.1.E-10.AND.XIN2.LE.X
  372. 1 .AND.CO2.GT.0.9) THEN
  373. DL3=DL32
  374. X=XIN2
  375. DL2=DLIN2
  376. GOTO 3000
  377. ENDIF
  378. IF(IIMPI.EQ.9) WRITE(IOIMP,30303)
  379. 30303 FORMAT(1X,'ERREUR DANS C2DP COMP ',/)
  380. DL3=DL31
  381. DL2=DLIN1
  382. X=XIN1
  383. GOTO 3000
  384. 2000 CONTINUE
  385. C
  386. C------------------------------------------------
  387. C CAS OU AUCUN AUTRE CRITERE EST ENDOMMAGE
  388. C------------------------------------------------
  389. C-----------------------------------
  390. C REMISE A JOUR DES VARIABLES
  391. C PUIS NOUVEAU PAS SI TOUT
  392. C L INCREMENT N EST PAS ECOULE
  393. C-----------------------------------
  394. C
  395. XT=XT+X
  396. RDP=RDP-DL3*HDP
  397. RMAX=MAX((RDP/1.73),(RDP/(1.D0-2.D0*ADP)))
  398. IF(SIREF.GT.RMAX) THEN
  399. RDP=0.D0
  400. GOTO 110
  401. ENDIF
  402. DO 2010 ITYP=1,3
  403. SFG(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  404. 2010 EPDP(ITYP)=DL3*EPC(ITYP)+EPDP(ITYP)
  405. XLAM2=XLAM2+DL2
  406. XLAM3=XLAM3+DL3
  407. IF(XLAM2.LT.1.E-10) XLAM2=0.D0
  408. IF(XLAM2.LT.1.E-10) IDAM(2)=0
  409. C=1.D0-1.D-10
  410. IF(XT.LT.C) GOTO 2020
  411. DO 2015 ITYP=1,3
  412. IDAM(ITYP)=0
  413. 2015 DSIG(ITYP)=0.D0
  414. CALL CHREP(SFG,SIG,-ANG)
  415. CALL CHREP(EPDP,EPPLDP,-ANG)
  416. IF(IIMPI.EQ.9) WRITE(IOIMP,9998) ITER
  417. RETURN
  418. 2020 IF(IDAM(2).EQ.-1) GOTO 100
  419. XLAM2=0.D0
  420. GAMDP=10.D0
  421. GAMTR1=10.D0
  422. GAMTR2=10.D0
  423. GAMCO1=10.D0
  424. DO 2030 ITYP=1,3
  425. DSFG(ITYP)=DSFG(ITYP)*(1.D0-XT)
  426. 2030 TENS(ITYP)=SFG(ITYP)+DSFG(ITYP)
  427. CALL CTRAF(TENS(1),RT1,VCTR1)
  428. CALL CTRAF(TENS(2),RT2,VCTR2)
  429. CALL CCOAF(TENS(1),XLAM1,VCCO1)
  430. CALL CDP(TENS,ADP,RDP,VCDP)
  431. IF(VCTR1.GT.0.D0) CALL GAMTAF(SFG(1),DSFG(1),RT1,GAMTR1)
  432. IF(VCTR2.GT.0.D0) CALL GAMTAF(SFG(2),DSFG(2),RT2,GAMTR2)
  433. IF(VCCO1.GT.0.D0) CALL GAMCAF(SFG(1),DSFG(1),GAMCO1)
  434. IF(VCDP.GT.0.D0) CALL GDP(SFG,DSFG,RDP,ADP,GAMDP)
  435. IDAM(1)=0
  436. IDAM(2)=0
  437. IDAM(3)=0
  438. GAM=MIN(GAMTR1,GAMTR2,GAMCO1,GAMDP)
  439. IF(GAM.GE.1.D0) THEN
  440. DO 2040 ITYP=1,3
  441. SFG(ITYP)=SFG(ITYP)+DSFG(ITYP)
  442. 2040 DSIG(ITYP)=0.D0
  443. CALL CHREP(SFG,SIG,-ANG)
  444. RETURN
  445. ENDIF
  446. IF(ABS(GAM-GAMTR1).LE.1.E-10) IDAM(1)=1
  447. IF(ABS(GAM-GAMTR2).LE.1.E-10) IDAM(2)=1
  448. IF(ABS(GAM-GAMCO1).LE.1.E-10) IDAM(1)=-1
  449. IF(ABS(GAM-GAMDP).LE.1.E-10) IDAM(3)=1
  450. DO 2050 ITYP=1,3
  451. SFG(ITYP)=SFG(ITYP)-GAM*DSFG(ITYP)
  452. 2050 DSFG(ITYP)=(1.D0-GAM)*DSFG(ITYP)
  453. CALL CHREP(SFG,SIG,-ANG)
  454. CALL CHREP(DSFG,DSIG,-ANG)
  455. CALL CHREP(EPDP,EPPLDP,-ANG)
  456. RETURN
  457. 3000 RDP=RDP-DL3*HDP
  458. RMAX=MAX((RDP/1.73),(RDP/(1.D0-2.D0*ADP)))
  459. IF(SIREF.GT.RMAX) THEN
  460. RDP=0.D0
  461. GOTO 110
  462. ENDIF
  463. XT=XT+X
  464. DO 3010 ITYP=1,3
  465. SFG(ITYP)=SFG(ITYP)+X*DSFG(ITYP)-DL2*SG(ITYP)-DL3*SSI(ITYP)
  466. DSFG(ITYP)=(1.D0-XT)*DSFG(ITYP)
  467. 3010 EPDP(ITYP)=DL3*EPC(ITYP)+EPDP(ITYP)
  468. CALL CHREP(SFG,SIG,-ANG)
  469. CALL CHREP(DSFG,DSIG,-ANG)
  470. CALL CHREP(EPDP,EPPLDP,-ANG)
  471. XLAM2=XLAM2+DL2
  472. XLAM3=XLAM3+DL3
  473. IF(IIMPI.EQ.9) WRITE(IOIMP,9998) ITER
  474. RETURN
  475. END
  476.  
  477.  

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