Télécharger la plaquette Cast3M


fr en

Télécharger tubturb.dgibi

Retour à la liste

Numérotation des lignes : 

  1. *****************************************************
  2. ************************************************************************
  3. ************************************************************************
  4. * fichier :  tubturb.dgibi                          *
  5. ** modifie le 15/06/2014 passage EQPR -> EQEX       *
  6. *****************************************************
  7. COMPLET = FAUX ;
  8. GRAPH = 'N' ;
  9.  
  10. KPRESS='CENTRE';
  11. DISCR='LINE';
  12. BETA=1.;
  13.  
  14. SI ( COMPLET ) ;
  15. along = 3 ; n = 24 ; nlong = 6 ;
  16. SINON ;
  17. along = 3 ; n = 24 ; nlong = 15 ;
  18. FINSI ;
  19.  
  20. ******************************************************************
  21. * TESTE DE STABILITÉ DU MODELE K-E ET DES FONCTIONS DE PAROI EN  * 
  22. * CONVECTION FORCE - ÉCOULEMENT EN CONDUITE LISSE CYLINDRIQUE    *     
  23. * LES PROFILS RADIAUX DES VARIABLES V,K,E OBTENUS À 3 DIAMÈTRES  *
  24. * DE L'ENTRÉE SONT COMPARÉS AU PROFIL OBTENU LORS DE L'ÉTABLISSE * 
  25. * DE LA FICHE DE VALIDATION CORRESPONDANTE                       * 
  26. ******************************************************************
  27.  
  28. 'OPTI' 'DIME' 2 'ELEM' 'QUA4'  'MODE' 'AXIS'  ;
  29. 'OPTI' 'ISOV' 'SURF'  ;
  30.  
  31. ** dimension
  32. Rh = 0.05 ; DH = 2. * RH ; LH = along*DH ;
  33.  
  34. ** maillage
  35. A = 0.0 0. ; B = RH  0.  ; C = RH  LH  ; D = 0.  LH  ;
  36. bas = a DROIT n b ;
  37. paroi = b DROIT (2*n*along/nlong) c ;
  38. haut = c droit n d ;
  39. axe = d droit  (2*n*along/nlong) a ;
  40. mgeo = dalle bas paroi haut axe plan ;
  41. mgeo = orien mgeo ;
  42. cgeo = contour mgeo ;
  43.  
  44. *-----------------------------------------------------------------
  45. **** donnée hydraulique de l écoulement 
  46. NU = 6.53e-4 / 992.2D0 ;
  47. CNU = 0.09 ;
  48. AKER = 0.41 ;
  49. REY = 100000. ;
  50. UI = REY*NU/DH ; 
  51.  
  52. *-----------------------------------------------------------------
  53. **** calcul de la perte de charge 
  54.  
  55. LAM = LOG REY / (LOG 10.) * 1.8 - 1.64 ** (-2.) ;
  56. mess 'Lambda = ' LAM ;
  57. REPE BL1 5 ;
  58.   lam = lam ** 0.5 * REY log / (log 10) * 2. - 0.8 ** (-2.) ;
  59.   mess 'Lambda = ' LAM ;
  60. FIN bl1 ;
  61.  
  62. DPSRDZ = LAM*0.5/DH*UI*UI ;
  63.  
  64. *-----------------------------------------------------------------
  65. **** calcul de la vitesse de glissement à partir des pertes de charge
  66. TO0     = DPSRDZ * Rh * 0.5 ;
  67. TO0_IMP = (-1.*to0) 0. ;
  68. UETH    = TO0 ** 0.5 ; 
  69.  
  70. *-----------------------------------------------------------------
  71. **** restriction du domaine à la zone logarithmique (Y+>40)
  72. R = MGEO COOR 1 ;
  73. Y = RH - R ;
  74. YCL= RH * 0.05 ; 
  75. DEPLA MGEO AFFI (RH - YCL/RH) A B ;
  76. R = MGEO COOR 1 ;
  77. Y = RH - R ;
  78.  
  79. *-----------------------------------------------------------------
  80. **** construction des objets modeles
  81. Mmgeo = chan mgeo QUAF ;
  82. Mbas  = chan bas  QUAF ;
  83. ELIM Mmgeo Mbas 0.0001 ;
  84. $DOMT = 'MODE' Mmgeo 'NAVIER_STOKES' DISCR ;
  85. $DBAS = 'MODE' Mbas  'NAVIER_STOKES' DISCR ;
  86. mgeo=doma $DOMT 'MAILLAGE';
  87. bas=doma $DBAS 'MAILLAGE';
  88. DOMA $DOMT 'IMPR' ;
  89. VTBAS = PI * RH * RH ; 
  90. Mhaut  = chan haut  QUAF ;
  91. nbp = nbel paroi ;
  92. paroi = elem paroi (lect 2 pas 1 nbp) ;
  93. Mparoi = chan paroi QUAF ;
  94. ELIM (Mmgeo et Mbas et Mhaut et Mparoi) 0.0001 ;
  95. $DHAU = 'MODE' Mhaut 'NAVIER_STOKES' DISCR ;
  96. $DPAR = 'MODE' Mparoi 'NAVIER_STOKES' DISCR ;
  97. haut= doma $DHAU 'MAILLAGE';
  98. paroi= doma $DPAR 'MAILLAGE';
  99. YPLUS = Y * UETH / NU ;
  100.  
  101. *-----------------------------------------------------------------
  102. **** estimation du profil de vitesse (hypothèse de couche limite log.)
  103. VZ =  YPLUS LOG / AKER + 5.5 * UETH ;
  104. VZ = NOMC VZ 'UY' ;
  105. VZ1 = KCHT $DOMT VECT SOMMET VZ ;
  106. DZ = DBIT $DBAS (VZ1 + (VZ1 mini * PI * YCL* RH)) ;
  107. UZ = DZ / VTBAS ;
  108. **** correction du profil de vitesse pour le respect du REYNOLDS  
  109. MESS 'FACTEUR DE CORRECTION DES VITESSES ' (UI/UZ) ;
  110. VZ = KCHT $DOMT VECT SOMMET (VZ1/UZ*UI) ;
  111.  
  112. *-----------------------------------------------------------------
  113. **** estimation de la viscosité turbulente  
  114. NUT = AKER * Y * UETH * (R / RH) + (Y/RH*nu) ;
  115. NUTET = NUT/UETH/RH*0.75+(Y/RH*0.015) ;
  116. NI = kcht $domt scal sommet (NUTET*UETH*RH) ;
  117.  
  118. *-----------------------------------------------------------------
  119. **** calcul de la longueur de mélange
  120. R2 = R/RH ** 2 ;
  121. LP = R2 * 0.06 + 0.08 * R2 * (-1.) + 0.14 ;
  122. LP = LP * RH ;
  123.  
  124. *-----------------------------------------------------------------
  125. **** estimation de l'énergie cinétique turbulente à partir
  126. **** de la longueur de mélange 
  127. KI = NI / LP ** 2. / (CNU**0.5) ;
  128.  
  129. *-----------------------------------------------------------------
  130. **** estimation de la dissipation visqueuse
  131. EI = CNU * KI * KI / NI ;
  132.  
  133. *-----------------------------------------------------------------
  134. **** conditions en entrée 
  135. VEN = exco (REDU VZ BAS) 'UY'   '2UN'  ; 
  136. KEN = exco (REDU KI BAS) 'SCAL' 'KN'   ;
  137. EEN = exco (REDU EI BAS) 'SCAL' 'EN'   ; 
  138.  
  139. MT  = table ;
  140. MT.1 = 'MARQ TRIA REGU TITR SORTIE' ;
  141. MT.2 = 'TIRR TITR ENTREE' ; 
  142. HAUT  = INVE HAUT ;
  143. ZC = MGEO COOR 2 ; ZC = (MAXI ZC) - ZC ; 
  144.  
  145. RV  = EQEX $DOMT 'ITMA' 160 'ALFA' 0.9 
  146. 'OPTI' 'SUPG' 
  147.  'ZONE' $DOMT 'OPER' 'FILTREKE' UI Dh      'INCO' 'KN' 'EN'
  148.  'ZONE' $DOMT 'OPER' 'NSKE'   NU   'NUT' 'INCO' 'UN' 'KN' 'EN'
  149.  'ZONE' $DPAR 'OPER' 'FPU' NU 'UET' YCL 'INCO' 'UN' 'KN' 'EN' ;              
  150.  
  151. RV = EQEX RV 'OPTI' 'CENTREE'
  152.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'UN' 'DELTAT' INCO 'UN'
  153.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'KN' 'DELTAT' INCO 'KN'
  154.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'EN' 'DELTAT' INCO 'EN'
  155. ;
  156.  
  157. RV = EQEX RV  
  158.    'CLIM'  'UN' 'VIMP'  BAS  VEN  
  159.    'CLIM'  'KN' 'TIMP'  BAS  KEN  
  160.    'CLIM'  'EN' 'TIMP'  BAS  EEN  
  161.    'CLIM'  'UN' 'UIMP'  CGEO 0.   ;
  162.  
  163. RVP = EQEX 'OPTI' 'EF' KPRESS
  164.  'ZONE' $DOMT OPER  KBBT  -1. beta  INCO 'UN' 'PRES'
  165. ;
  166.  
  167.     rvp.'METHINV'.TYPINV=1 ;
  168.     rvp.'METHINV'.IMPINV=0 ;
  169.     rvp.'METHINV'.NITMAX=300;
  170.     rvp.'METHINV'.PRECOND=3 ;
  171.     rvp.'METHINV'.RESID  =1.e-8 ;
  172.     rvp.'METHINV' . 'FCPRECT'=100 ;
  173.     rvp.'METHINV' . 'FCPRECI'=100 ;
  174.  
  175.   RV.'PROJ'     =RVP ;
  176.  
  177. RV.'INCO' = table 'INCO' ;
  178. RV.INCO.'UN' =  VZ ;
  179. RV.INCO.'PRES' = KCHT $DOMT 'SCAL' KPRESS 0.;
  180. RV.INCO.'NUT' = NOEL $DOMT NI ;
  181. RV.INCO.'KN' =  KCHT $DOMT SCAL SOMMET KI ;
  182. RV.INCO.'EN' =  KCHT $DOMT SCAL SOMMET EI ;
  183. RV.'INCO'.'UET' = KCHT $DPAR SCAL CENTRE UETH ;
  184.  
  185. *-----------------------------------------------EXEC
  186.   EXEC RV ; 
  187. *-----------------------------------------------EXEC
  188.  
  189. SI ( COMPLET ) ;
  190. ** On recherche le régime permanent :
  191.  
  192.  
  193. repe bl2 15 ;
  194.  
  195. K = RV.'INCO'.'KN' ;
  196. E = RV.'INCO'.'EN' ;
  197. N =  RV.'INCO'.'NUT' ;
  198. U =  kops (RV.'INCO'.'UN') '*' (-1.)  ;
  199. UT = RV.'INCO'.'UET' ;
  200.  
  201. RV  = EQEX $DOMT 'ITMA' 400 'ALFA' 0.8 
  202. 'OPTI' 'SUPG' 
  203.  'ZONE' $DOMT 'OPER' 'FILTREKE' UI Dh      'INCO' 'KN' 'EN'
  204.  'ZONE' $DOMT 'OPER' 'NSKE'  NU  'NUT' 'INCO' 'UN' 'KN' 'EN'
  205.  'ZONE' $DPAR 'OPER' 'FPU' NU 'UET' YCL 'INCO' 'UN' 'KN' 'EN' ;
  206.  
  207. RV = EQEX RV 'OPTI' 'CENTREE'
  208.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'UN' 'DELTAT' INCO 'UN'
  209.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'KN' 'DELTAT' INCO 'KN'
  210.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'EN' 'DELTAT' INCO 'EN'
  211. ;
  212.  
  213. RVP = EQEX 'OPTI' 'EF' KPRESS
  214.  'ZONE' $DOMT OPER  KBBT  -1. beta  INCO 'UN' 'PRES'
  215. ;
  216.  
  217.     rvp.'METHINV'.TYPINV=1 ;
  218.     rvp.'METHINV'.IMPINV=0 ;
  219.     rvp.'METHINV'.NITMAX=300;
  220.     rvp.'METHINV'.PRECOND=3 ;
  221.     rvp.'METHINV'.RESID  =1.e-8 ;
  222.     rvp.'METHINV' . 'FCPRECT'=100 ;
  223.     rvp.'METHINV' . 'FCPRECI'=100 ;
  224.  
  225.   RV.'PROJ'     =RVP ;
  226.  
  227. RV.'INCO' = table 'INCO' ;
  228. RV.INCO.'UN' =  U ;
  229. RV.INCO.'PRES' = KCHT $DOMT 'SCAL' KPRESS 0.;
  230. RV.INCO.'NUT' = N ;
  231. RV.INCO.'KN' =  K ;
  232. RV.INCO.'EN' =  E ;
  233. RV.'INCO'.'UET' = UT ;
  234.  
  235. VEN        = exco (redu HAUT U) 'UY' 'SCAL' ; 
  236. KEN        =  redu HAUT K ; 
  237. EEN        =  redu HAUT E ; 
  238.  
  239.   RV.CLIM = 0 ;
  240.   RV = EQEX RV  
  241.    'CLIM'  'UN' 'VIMP'  HAUT VEN  
  242.    'CLIM'  'KN' 'TIMP'  HAUT KEN  
  243.    'CLIM'  'EN' 'TIMP'  HAUT EEN  
  244.    'CLIM'  'UN' 'UIMP'  CGEO 0.   ;
  245.  
  246. *-----------------------------------------------EXEC aller
  247.   EXEC RV ; 
  248. *-----------------------------------------------EXEC aller
  249.  
  250. K2 =  copi RV.'INCO'.'KN' ;
  251. K = RV.'INCO'.'KN' ;
  252. E = RV.'INCO'.'EN' ;
  253. N =  RV.'INCO'.'NUT' ;
  254. U =  kops (RV.'INCO'.'UN') '*' (-1.)  ;
  255. UT = RV.'INCO'.'UET' ;
  256.  
  257. RV  = EQEX $DOMT 'ITMA' 400 'ALFA' 0.8 
  258. 'OPTI' 'SUPG' 
  259.  'ZONE' $DOMT 'OPER' 'FILTREKE' UI Dh    'INCO' 'KN' 'EN'
  260.  'ZONE' $DOMT 'OPER' 'NSKE'   NU   'NUT' 'INCO' 'UN' 'KN' 'EN'
  261.  'ZONE' $DPAR 'OPER' 'FPU' NU 'UET' YCL 'INCO' 'UN' 'KN' 'EN' ;
  262.  
  263. RV = EQEX RV 'OPTI' 'CENTREE'
  264.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'UN' 'DELTAT' INCO 'UN'
  265.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'KN' 'DELTAT' INCO 'KN'
  266.  'ZONE' $DOMT 'OPER' 'DFDT' 1. 'EN' 'DELTAT' INCO 'EN'
  267. ;
  268.  
  269. RVP = EQEX 'OPTI' 'EF' KPRESS
  270.  'ZONE' $DOMT OPER  KBBT  -1. beta  INCO 'UN' 'PRES'
  271. ;
  272.  
  273.     rvp.'METHINV'.TYPINV=1 ;
  274.     rvp.'METHINV'.IMPINV=0 ;
  275.     rvp.'METHINV'.NITMAX=300;
  276.     rvp.'METHINV'.PRECOND=3 ;
  277.     rvp.'METHINV'.RESID  =1.e-8 ;
  278.     rvp.'METHINV' . 'FCPRECT'=100 ;
  279.     rvp.'METHINV' . 'FCPRECI'=100 ;
  280.  
  281.   RV.'PROJ'     =RVP ;
  282.  
  283. RV.'INCO' = table 'INCO' ;
  284. RV.INCO.'UN' =  U ;
  285. RV.INCO.'PRES' = KCHT $DOMT 'SCAL' KPRESS 0.;
  286. RV.INCO.'NUT' = N ;
  287. RV.INCO.'KN' =  K ;
  288. RV.INCO.'EN' =  E ;
  289. RV.'INCO'.'UET' = UT ;
  290.  
  291. VEN        = exco (redu BAS U) 'UY' 'SCAL' ; 
  292. KEN        =  redu BAS K ; 
  293. EEN        =  redu BAS E ; 
  294.  
  295.   RV.CLIM = 0 ;
  296.   RV = EQEX RV  
  297.    'CLIM'  'UN' 'VIMP'  BAS VEN  
  298.    'CLIM'  'KN' 'TIMP'  BAS KEN  
  299.    'CLIM'  'EN' 'TIMP'  BAS EEN  
  300.    'CLIM'  'UN' 'UIMP'  CGEO 0.   ;
  301.  
  302. *-----------------------------------------------EXEC retour
  303.   EXEC RV ; 
  304. *-----------------------------------------------EXEC retour
  305.  
  306. erm = k - k2 / (k maxi) abs maxi ;
  307. mess 'maxi (k - k2 / (k maxi)) : ' erm ;
  308. si ( erm < 2.e-5) ;  
  309.   quitter bl2 ;
  310. finsi ;
  311. menage ;
  312. FIN bl2 ;
  313.  
  314. *********************************************
  315. *                                           *
  316. *      POSTRAITEMENT DES RÉSULTATS          * 
  317. *                                           *
  318. *********************************************
  319.  
  320. SI ( 'EGA' graph 'O') ;
  321.  
  322. EVV = EVOL CHPO  (RV.INCO.UN/UETH) 'UY' (DOMA $DBAS MAILLAGE)  ;
  323. EVV0 = EVOL VERT CHPO (VZ/UETH) 'UY' (DOMA $DBAS MAILLAGE) ;
  324. RAY = EXTR evv 'ABSC' ; RAY = RAY / RH ;
  325. YY = EXTR evv 'ORDO' ; YY0 = EXTR evv0 'ORDO' ;
  326. EVV  = EVOL MANU  RAY  YY 'UY' ; 
  327. EVV0 = EVOL MANU  RAY  YY0 'UY' ; 
  328.  
  329. RV.INCO.'NNT' = ELNO $DOMT RV.INCO.'NUT' ;
  330. EVN = EVOL CHPO  (RV.INCO.'NNT'/(UETH*RH)) 'SCAL' (DOMA $DBAS MAILLAGE);
  331. EVN0 = EVOL VERT CHPO  NUTET 'SCAL' (DOMA $DBAS MAILLAGE)  ;
  332. YY = EXTR evN 'ORDO' ; YY0 = EXTR evN0 'ORDO' ;
  333. EVN  = EVOL MANU  RAY  YY 'SCAL' ; 
  334. EVN0 = EVOL MANU  RAY  YY0 'SCAL' ; 
  335.  
  336. EVK = EVOL CHPO (RV.INCO.'KN'/(UETH*UETH)) 'SCAL' (DOMA $DBAS MAILLAGE);
  337. EVK0 = EVOL VERT CHPO  (KI/UETH/UETH) 'SCAL' (DOMA $DBAS MAILLAGE)  ;
  338. YY = EXTR evK 'ORDO' ; YY0 = EXTR evK0 'ORDO' ;
  339. EVK  = EVOL MANU  RAY  YY 'SCAL' ; 
  340. EVK0 = EVOL MANU  RAY  YY0 'SCAL' ; 
  341.  
  342. EVE = EVOL CHPO (RV.INCO.'EN'/(UETH*UETH*UETH/RH))
  343.  'SCAL'  (DOMA $DBAS MAILLAGE)  ;
  344. EVE0 = EVOL VERT CHPO (EI/(UETH*UETH*UETH/RH)) 
  345. 'SCAL' (DOMA $DBAS MAILLAGE)  ;
  346. YY = EXTR evE 'ORDO' ; YY0 = EXTR evE0 'ORDO' ;
  347. EVE  = EVOL MANU  RAY  YY  'SCAL' ; 
  348. EVE0 = EVOL MANU  RAY  YY0 'SCAL' ; 
  349.  
  350. RV.INCO.'LN'  = KOPS (0.09**(-0.25)) '*'
  351. (KOPS RV.INCO.'NNT' '/' (KOPS RV.INCO.'KN' '**' 0.5)) ;
  352. EVL = EVOL CHPO ((RV.INCO.'LN')/RH) 'SCAL'  (DOMA $DBAS MAILLAGE)  ;
  353. EVL0 = EVOL VERT CHPO (LP/RH) 'SCAL' (DOMA $DBAS MAILLAGE)  ;
  354. YY = EXTR evL 'ORDO' ; YY0 = EXTR evL0 'ORDO' ;
  355. EVL  = EVOL MANU  RAY  YY  'SCAL' ; 
  356. EVL0 = EVOL MANU  RAY  YY0 'SCAL' ; 
  357.  
  358. nn = nbel dbas.maillage ;
  359. EVT0 =  EVOL MANU  RAY  RAY 'SCAL' ; 
  360. nn = nbel dPAR.maillage ;
  361. X = (dbas.maillage COOR 1) MAXI / RH ;
  362. TOC = SOMT RV.INCO.'UET' / nn / UETH ** 2 ;
  363. EVT =  EVOL MANU  (PROG X)  (PROG TOC)  'SCAL' ;
  364.  
  365. MT = TABLE ;
  366. MT.1 = 'MARQ CROI REGU TIRR' ; 
  367. MT.2 = 'MARQ CARR  ' ; 
  368.  
  369. DESS MT (EVT0  ET EVT)
  370. titr 'PROFIL RADIAL DU CISAILLEMENT ADIMENSIONNELLE'
  371. TITX 'RAYON/RH' TITY  'TO/UETH/UETH'
  372.  
  373. DESS (EVN ET EVN0) 
  374. titr 'PROFIL RADIAL DE LA VISCOSITE ADIMENSIONNELLE'
  375. TITX 'RAYON/RH' TITY  'NUT/UETH/RH';
  376.  
  377. DESS (EVE ET EVE0) 
  378. titr 'PROFIL RADIAL DU TAUX DE DISS. ADIMENSIONNEL'
  379. TITX 'RAYON/RH' TITY  'E/UETH*RH';
  380.  
  381. DESS (EVV ET EVV0) 
  382. titr 'PROFIL RADIAL DE LA VITESSE ADIMENSIONNELLE'
  383. TITX 'RAYON/RH' TITY  'V/UETH';
  384.  
  385. DESS (EVK ET EVK0) titr 
  386. 'PROFIL RADIAL DE L ENERGIE TURBULENTE ADIMENSIONNELLE'
  387. TITX 'RAYON/RH' TITY  'K/UETH/UETH';
  388.  
  389. DESS (EVL ET EVL0) titr 
  390. 'PROFIL RADIAL DE LA LONGUEUR DE MELANGE ADIMENSIONNELLE'
  391. TITX 'RAYON/RH' TITY  'L/RH';
  392.  
  393. FINSI ;
  394.  
  395. SINON ;
  396. EVV = EVOL CHPO  (RV.INCO.UN/UETH) 'UY' (DOMA $DBAS MAILLAGE)  ;
  397. v = EXTR evv 'ORDO' ;
  398. list v ;
  399. v0 = prog 
  400. 26.642 26.535 26.423 26.306 26.184 26.056 25.921 25.780 25.630
  401. 25.472 25.303 25.123 24.930 24.721 24.495 24.248 23.976 23.672  
  402. 23.330 22.936 22.474 21.914 21.204 20.231 18.682 ;
  403. VDIF0 = V - V0 / V0 abs maxi ; list VDIF0 ;
  404. SI ( VDIF0  > 5.E-5 ) ;
  405.             ERREUR 5 ;
  406. FINSI ;
  407.  
  408. EVK = EVOL CHPO (RV.INCO.'KN'/(UETH*UETH)) 'SCAL' (DOMA $DBAS MAILLAGE);
  409. K = EXTR evK 'ORDO' ;
  410. list k ;
  411. K0 = PROG
  412. 3.98961E-02 .11874 .22847 .36068 .50826 .66518 .82637 .98756 1.1453
  413. 1.2966 1.4394 1.5718 1.6926 1.8010 1.8965 1.9789 2.0483  2.1050  
  414. 2.1495 2.1822 2.2041 2.2158 2.2181 2.2120 2.1982 ; 
  415. KDIF0 = K - K0 / K0 ABS MAXI ; list KDIF0 ;
  416. SI (KDIF0  > 5.E-5 ) ;
  417.             ERREUR 5 ;
  418. FINSI ;
  419.  
  420. EVE = EVOL CHPO (RV.INCO.'EN'/(UETH*UETH*UETH/RH))
  421.  'SCAL' (DOMA $DBAS MAILLAGE)  ;
  422. E = EXTR eve 'ORDO' ;
  423. list e ;
  424. E0 = PROG
  425. 9.35296E-03 4.80687E-02 .12863 .25632 .43158 .65176 .91234 1.2079      
  426. 1.5329 1.8824 2.2521 2.6394 3.0434 3.4653 3.9095 4.3841 4.9025 
  427. 5.4863 6.1700 7.0111 8.1107 9.6651 12.110 16.648 28.290 ;
  428. EDIF0 = E - E0 / E0 ABS MAXI ; list EDIF0 ;
  429. SI ( EDIF0  > 5.E-5 ) ;
  430.             ERREUR 5 ;
  431. FINSI ;
  432.  
  433. FINSI ;
  434.  
  435. FIN;
  436.  
  437.  
  438.  
  439.  
  440.  
  441.  
  442.  
  443.  
  444.  
  445.  
  446.  
  447.  
  448.  
  449.  
  450.  

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