Télécharger la plaquette Cast3M


fr en

Télécharger konv_ther_cons2.dgibi

Retour à la liste

Numérotation des lignes : 

  1. * fichier :  konv_ther_cons2.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***********************************************************
  5. **** APPROCHE VF "Cell-Centred Formulation" pour la    ****
  6. **** solution des                                      **** 
  7. **** Equations d'Euler pour un gaz parfait             ****
  8. **** OPERATEURS PRET, KONV                             ****
  9. ****                                                   ****
  10. **** Cas gaz multi-especes  "thermally perfect"        ****
  11. **** Consistence                                       ****
  12. **** Splitting des scalaires passifs                   ****
  13. ****                                                   ****
  14. **** Methodes: VLH, Colella                            ****
  15. ****                                                   ****
  16. **** A. BECCANTINI DRN/DMT/SEMT/LTMF    FEVRIER 2000   ****
  17. ***********************************************************
  18.  
  19. 'OPTION'  'DIME' 2 ;
  20. 'OPTION'  'ELEM' QUA4 ;
  21. 'OPTION'  'ECHO'  0 ;
  22. 'OPTION'  'TRAC' 'X' ;
  23.  
  24. *
  25. *** GRAPH
  26. *
  27.  
  28.  
  29. GRAPH = FAUX ;
  30. * GRAPH = VRAI ;
  31.  
  32. ***********************
  33. **** LA TABLE PGAZ ****
  34. ***********************
  35.  
  36. PGAZ = 'TABLE' ;
  37. PGAZ2 = 'TABLE' ;
  38.  
  39. *
  40. **** Ordre des polynoms cv_i
  41. *
  42.  
  43. PGAZ . 'NORD' = 1 ;
  44. PGAZ2 . 'NORD' = 1 ;
  45.  
  46. *
  47. **** Especes qui sont dans les equations d'Euler
  48. *
  49.  
  50. PGAZ . 'ESPEULE' = 'MOTS' 'H2  ' 'O2  ' 'H2O ' ;
  51. PGAZ2 . 'ESPEULE' = 'MOTS' 'H2  ' 'O2  ' 'H2O ' ;
  52.  
  53. *
  54. **** Espece qui n'y est pas
  55. *
  56.  
  57.  
  58. PGAZ . 'ESPNEULE' = 'N2  ';
  59. PGAZ2 . 'ESPNEULE' = 'N2  ';
  60.  
  61. *
  62.  
  63. PGAZ .  'H2  ' = 'TABLE'  ;
  64. PGAZ .  'H2O ' = 'TABLE'  ;
  65. PGAZ .  'N2  ' = 'TABLE'  ;
  66. PGAZ .  'O2  ' = 'TABLE'  ;
  67. PGAZ2 .  'H2  ' = 'TABLE'  ;
  68. PGAZ2 .  'H2O ' = 'TABLE'  ;
  69. PGAZ2 .  'N2  ' = 'TABLE'  ;
  70. PGAZ2 .  'O2  ' = 'TABLE'  ;
  71.  
  72. *
  73. **** R (J/Kg/K)
  74. *
  75.  
  76. PGAZ .  'H2  ' . 'R' = 4130.0 ;
  77. PGAZ .  'H2O ' . 'R' = 461.4 ;
  78. PGAZ .  'N2  ' . 'R' = 296.8 ;
  79. PGAZ .  'O2  ' . 'R' = 259.8 ;
  80. PGAZ2 .  'H2  ' . 'R' = 4130.0 ;
  81. PGAZ2 .  'H2O ' . 'R' = 461.4 ;
  82. PGAZ2 .  'N2  ' . 'R' = 296.8 ;
  83. PGAZ2 .  'O2  ' . 'R' = 259.8 ;
  84.  
  85.  
  86. *
  87. **** Regressions polynomials
  88. *
  89.  
  90. PGAZ .  'H2  ' . 'A' = 'PROG'  9834.91866 0.54273926 ;
  91. PGAZ .  'H2O ' . 'A' = 'PROG' 1155.95625 0.768331151 ;
  92. PGAZ .  'N2  ' . 'A' = 'PROG' 652.940766 0.288239099 ;
  93. PGAZ .  'O2  ' . 'A' = 'PROG' 575.012333  0.350522002 ;
  94. PGAZ2 .  'H2  ' . 'A' = 'PROG'  9834.91866 0.54273926 ;
  95. PGAZ2 .  'H2O ' . 'A' = 'PROG' 1155.95625 0.768331151 ;
  96. PGAZ2 .  'N2  ' . 'A' = 'PROG' 652.940766 0.288239099 ;
  97. PGAZ2 .  'O2  ' . 'A' = 'PROG' 575.012333  0.350522002 ;
  98.  
  99.  
  100. *
  101. **** "Enthalpies" (ou energies) de formations a OK (J/Kg)
  102. *     Note: ce sont des entites numeriques
  103.  
  104. PGAZ .  'H2  ' . 'H0K' = -4.195D6 ;
  105. PGAZ .  'H2O ' . 'H0K' = -1.395D7 ;
  106. PGAZ .  'N2  ' . 'H0K' = -2.953D5 ;
  107. PGAZ .  'O2  ' . 'H0K' = -2.634D5 ;
  108. PGAZ2 .  'H2  ' . 'H0K' = -4.195D6 ;
  109. PGAZ2 .  'H2O ' . 'H0K' = -1.395D7 ;
  110. PGAZ2 .  'N2  ' . 'H0K' = -2.953D5 ;
  111. PGAZ2 .  'O2  ' . 'H0K' = -2.634D5 ;
  112.  
  113.  
  114. * PGAZ2: transport de scalaires passifs
  115.  
  116. PGAZ2 . 'SCALPASS' = 'MOTS' 'TUTU' 'TATA' ;
  117.  
  118.  
  119. *************************** 
  120. ***** DOMAINE SPATIAL  ****
  121. ***************************
  122.  
  123.  
  124. A1 = 0.0D0 0.0D0;
  125. A2 = 1.0D0 0.0D0;
  126. A3 = 2.0D0 0.0D0;
  127. A4 = 2.0D0 1.0D0;
  128. A5 = 1.0D0 1.0D0;
  129. A6 = 0.0D0 1.0D0;
  130.  
  131. L12 = A1 'DROIT' 1 A2;
  132. L23 = A2 'DROIT' 1 A3;
  133. L34 = A3 'DROIT' 1 A4;
  134. L45 = A4 'DROIT' 1 A5;
  135. L56 = A5 'DROIT' 1 A6;
  136. L61 = A6 'DROIT' 1 A1;
  137. L25 = A2 'DROIT' 1 A5;
  138.  
  139.  
  140. DOM10  = 'DALL' L12 L25 L56 L61
  141.         'PLANE';
  142. DOM20  = 'DALL' L23 L34 L45 ('INVERSE' L25)
  143.         'PLANE';
  144. *
  145. *** Point ou on controlle la consistence
  146. *
  147.  
  148. P10 = 1.0 0.5;
  149.  
  150. *
  151. *** Etats gauche et droite 
  152. *
  153.  
  154. rogd = 1.412 ;
  155. pgd = 1.17D6 ;
  156. ungd  = 8.34 ;
  157. utgd  = 18.2 ;
  158. yh2gd = 0.1 ;
  159. yh2Ogd = 0.2 ;
  160. yN2gd = 0.15 ;
  161. yO2gd = 1 '-' (yh2gd '+' yh2Ogd '+' yn2gd) ;
  162. tutugd = 14.9 ;
  163. tatagd = 23.8 ;
  164.  
  165. *
  166. **** tg
  167. *
  168.  
  169. tgd = pgd '/' (rogd *
  170.       ((yh2gd  * (PGAZ .  'H2  ' . 'R')) '+'
  171.        (yh2Ogd * (PGAZ .  'H2O ' . 'R')) '+'
  172.        (yn2gd  * (PGAZ .  'N2  ' . 'R')) '+'
  173.        (yo2gd  * (PGAZ .  'O2  ' . 'R')))) ;
  174.  
  175.  
  176. *
  177. **** gamg, gamd , htg , htd
  178. *
  179. *
  180.  
  181. A0H2 = 'EXTRAIRE' (PGAZ . 'H2  ' . 'A') 1 ;
  182. A1H2 = 'EXTRAIRE' (PGAZ . 'H2  ' . 'A') 2 ;
  183.  
  184. A0H2O= 'EXTRAIRE' (PGAZ . 'H2O ' . 'A') 1 ;
  185. A1H2O= 'EXTRAIRE' (PGAZ . 'H2O ' . 'A') 2 ;
  186.  
  187. A0N2 = 'EXTRAIRE' (PGAZ . 'N2  ' . 'A') 1 ;
  188. A1N2 = 'EXTRAIRE' (PGAZ . 'N2  ' . 'A') 2 ;
  189.  
  190. A0O2 = 'EXTRAIRE' (PGAZ . 'O2  ' . 'A') 1 ;
  191. A1O2 = 'EXTRAIRE' (PGAZ . 'O2  ' . 'A') 2 ;
  192.  
  193. T = tgd ;
  194. T2 = T * T ;
  195.  
  196.  
  197.  ETGD = (yh2gd *
  198.         ((A0H2 * T) '+' (A1H2 * T2 '/' 2)));
  199.  
  200.  ETGD = ETGD '+' (yH2Ogd *
  201.         ((A0H2O * T) '+' (A1H2O * T2 '/' 2)));
  202.  
  203.  ETGD = ETGD '+' (yN2gd *
  204.         ((A0N2 * T) '+' (A1N2 * T2 '/' 2)));
  205.  
  206.  ETGD = ETGD '+' (yO2gd *
  207.         ((A0O2 * T) '+' (A1O2 * T2 '/' 2))) ;
  208.  
  209.  ecingd = 0.5D0 '*'  ((ungd '*' ungd) '+' (utgd '*' utgd));
  210.  
  211.  retgd = rogd '*' (ETGd '+' ecingd) ;
  212.  
  213. *
  214. ***  flux en (n,t)
  215. *
  216.  
  217. f1gd = ungd  '*'  rogd ;
  218. f2gd = (rogd '*' (ungd '*' ungd)) '+' pgd ;
  219. f3gd = rogd '*' ungd '*' utgd ;
  220. f4gd = ungd '*' (retgd '+' pgd);
  221. f5gd = ungd '*' rogd '*' yh2gd ;
  222. f6gd = ungd '*' rogd '*' yo2gd ;
  223. f7gd = ungd '*' rogd '*' yh2ogd ;
  224. ftutugd =  ungd  '*'  rogd '*' tutugd ;
  225. ftatagd =  ungd  '*'  rogd '*' tatagd ;
  226.  
  227.  
  228. ****************************************************
  229. ****************************************************
  230. ******** Boucle sur les angles             *********
  231. ****************************************************
  232. ****************************************************
  233.  
  234. DANGLE = 360 '/' 7.15;
  235. ANGLE = 11.3 ;
  236.  
  237. 'REPETER' BLMET 2 ;
  238. SI (&BLMET 'EGA' 1) ;
  239. METO = 'VLH' ;
  240. 'FINSI' ;
  241. SI (&BLMET 'EGA' 2) ;
  242. METO = 'SS' ;
  243. 'FINSI' ;
  244.  
  245.  
  246. 'REPETER' BLOC 8;
  247.  
  248. *
  249. *** Rotation
  250. *
  251.  
  252.  
  253. ANGLE = ANGLE '+' DANGLE;
  254. ORIG = 0.0D0 0.0D0;
  255.  
  256. MESSAGE;
  257. MESSAGE (CHAIN 'Angle de rotation= ' ANGLE);
  258. MESSAGE;
  259.  
  260. DOM1   = DOM10   'TOURNER' ANGLE ORIG;
  261. DOM2   = DOM20   'TOURNER' ANGLE ORIG;
  262. P1     = P10     'TOURNER' ANGLE ORIG;
  263.  
  264.     'ELIMINATION' (DOM1 ET DOM2) 1D-6;
  265.     DOMTOT = DOM1 ET DOM2;
  266.  
  267.  $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  268.  $DOM1 = 'MODELISER'  DOM1 'EULER';
  269.  $DOM2 = 'MODELISER'  DOM2 'EULER';
  270.  TDOMTOT = 'DOMA'  $DOMTOT 'VF';
  271.  TDOM1 = 'DOMA'  $DOM1 'VF';
  272.  TDOM2 = 'DOMA'  $DOM2 'VF';
  273.  MDOM1 = TDOM1 . 'QUAF' ;
  274.  MDOM2 = TDOM2 . 'QUAF' ;
  275.  MDOMTOT = TDOMTOT . 'QUAF' ;
  276.  
  277.  'ELIM' (MDOMTOT 'ET' MDOM1 'ET' MDOM2) 1.E-6 ;
  278.  
  279.     'SI' GRAPH;
  280.       'TRACER' (('DOMA' $DOMTOT 'MAILLAGE' ) 'ET' 
  281.       ('DOMA' $DOMTOT 'FACEL') 'ET' P1) 'TITRE' 'Domaine et FACEL';
  282.     'FINSI' ;
  283. *
  284. **** Redefinition de P1 dans $DOMTOT 'FACE'
  285. *
  286.  
  287. P1 = ('DOMA' $DOMTOT 'FACE') 'POIN' 'PROC' P1;
  288.  
  289. ***********************
  290. **** Les CHPOINTs  ****
  291. ***********************
  292.  
  293. uxgd = (ungd '*' ('COS' ANGLE)) '-' (utgd '*' ('SIN' ANGLE));
  294. uygd = (ungd '*' ('SIN' ANGLE)) '+' (utgd '*' ('COS' ANGLE));
  295.  
  296.  
  297. RN1   = 'KCHT' $DOM1   'SCAL' 'CENTRE' rogd;
  298. RN2   = 'KCHT' $DOM2   'SCAL' 'CENTRE' rogd;;
  299. RN    = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' (RN1 'ET' RN2);
  300.  
  301. VN1   = 'KCHT' $DOM1   'VECT' 'CENTRE' (uxgd uygd);
  302. VN2   = 'KCHT' $DOM2   'VECT' 'CENTRE' (uxgd uygd);
  303. VN    = 'KCHT' $DOMTOT 'VECT' 'CENTRE' (VN1 'ET' VN2);
  304.  
  305. PN1   = 'KCHT' $DOM1   'SCAL' 'CENTRE' pgd;
  306. PN2   = 'KCHT' $DOM2   'SCAL' 'CENTRE' pgd;
  307. PN    = 'KCHT' $DOMTOT 'SCAL' 'CENTRE' (PN1 'ET' PN2);
  308.  
  309. * attention: PGAZ . 'ESPEULE' = 'MOTS' 'H2  ' 'O2  ' 'H2O ' ;
  310.  
  311. YH2 = ('KCHT' $DOM1   'SCAL' 'CENTRE' 'COMP' 'H2  ' yh2gd) 'ET'
  312.       ('KCHT' $DOM2   'SCAL' 'CENTRE' 'COMP' 'H2  ' yh2gd) ;
  313.  
  314. YO2 = ('KCHT' $DOM1   'SCAL' 'CENTRE' 'COMP' 'O2  ' yo2gd) 'ET'
  315.       ('KCHT' $DOM2   'SCAL' 'CENTRE' 'COMP' 'O2  ' yo2gd) ;
  316.  
  317. YH2O = ('KCHT' $DOM1   'SCAL' 'CENTRE' 'COMP' 'H2O ' yh2ogd) 'ET'
  318.        ('KCHT' $DOM2   'SCAL' 'CENTRE' 'COMP' 'H2O ' yh2ogd) ;
  319.  
  320. YN = YH2 'ET' YO2 'ET' YH2O ;
  321.  
  322. *
  323.  
  324.  
  325. SCAN1   = 'KCHT' $DOM1   'VECT' 'CENTRE'
  326.           'COMP' (PGAZ2 . 'SCALPASS') (tutugd tatagd) ;
  327. SCAN2   = 'KCHT' $DOM2   'VECT' 'CENTRE'
  328.           'COMP' (PGAZ2 . 'SCALPASS') (tutugd tatagd) ;
  329. SCAN    = 'KCHT' $DOMTOT 'VECT' 'CENTRE'
  330.           'COMP' (PGAZ2 . 'SCALPASS') (SCAN1 'ET' SCAN2);
  331.  
  332. ***************************
  333. ****  L'operateur  PRET****
  334. ***************************
  335.  
  336. ORDESP = 1;
  337. ORDTEM = 1;
  338.  
  339.  
  340. ROF VITF PF YF =  'PRET' 'PERFTEMP'  ORDESP ORDTEM
  341.     $DOMTOT PGAZ RN VN PN YN ;
  342.  
  343.  
  344. ROF2 VITF2 PF2 YF2 SCAF  =  'PRET' 'PERFTEMP'  ORDESP ORDTEM
  345.     $DOMTOT PGAZ2 RN VN PN YN SCAN ;
  346.  
  347. GEOP1    = ('DOMA' $DOMTOT 'FACEL') 'ELEM' 'APPUYE' 'LARGEMENT' P1;
  348. GEOP2    = ('DOMA' $DOMTOT 'FACE')  'ELEM' 'APPUYE' 'LARGEMENT' P1;
  349.  
  350.  
  351. REFGEOP1 = 'REDU' VITF GEOP2;
  352.  
  353. *
  354. **** Orientation de la normal n de castem par raport a la
  355. *    notre; t est par consequence
  356. *
  357.  
  358. NCOS = 'EXTRAIRE' REFGEOP1 'NX' 1 1 1;
  359. NSIN = 'EXTRAIRE' REFGEOP1 'NY' 1 1 1;
  360.  
  361. 'SI' (('ABS' NCOS) > ('ABS' NSIN));
  362.    ORIENT = ('COS' ANGLE) '/' NCOS;
  363. 'SINON'; 
  364.    ORIENT = ('SIN' ANGLE) '/' NSIN;
  365. 'FINSI' ;
  366.  
  367. ORIENT = 'SIGN' ORIENT;   
  368.  
  369.  
  370. **************************
  371. **************************
  372. **** 'KONV', TEST VLH ****
  373. **************************
  374. **************************
  375.  
  376. LISTINCO = 'MOTS' 'RN' 'RUX' 'RUY' 'RETN' 'H2' 'O2' 'H2O' ;
  377.  
  378.  
  379. TFLUX DELTAT =  'KONV' 'VF' 'PERFTEMP' 'FLUX' METO
  380.          $DOMTOT PGAZ LISTINCO  ROF VITF PF YF ;
  381.  
  382. *
  383. **** Les flux aux interfaces sont dans de
  384. *    CHPOINT FACE
  385. *
  386.  
  387. FLUX1 = 'EXCO' 'RN' TFLUX ;
  388. FLUX2X = 'EXCO' 'RUX' TFLUX ;
  389. FLUX2Y = 'EXCO' 'RUY' TFLUX ;
  390. FLUX3 = 'EXCO' 'RETN' TFLUX ;
  391. FLUX4 = 'EXCO' ('MOTS' 'H2' 'O2' 'H2O') TFLUX 
  392.  ('MOTS' 'H2' 'O2' 'H2O') ;
  393.  
  394. FLUX2N = (FLUX2X '*' ('COS' ANGLE)) '+' (FLUX2Y * ('SIN' ANGLE));
  395. FLUX2T = (FLUX2Y '*' ('COS' ANGLE)) '-' (FLUX2X * ('SIN' ANGLE));
  396.  
  397. *
  398. **** Test 'KONV'
  399. *
  400.  
  401. f1 = 'EXTRAIRE' FLUX1 'SCAL' P1;
  402. ERRO = 'ABS' (1D-12 '*' f1gd);
  403. LOGI1 = ('ABS' ((f1 '*' ORIENT) '-' f1gd)) < ERRO;
  404.  
  405. f2 = 'EXTRAIRE' FLUX2N 'SCAL' P1;
  406. ERRO = 'ABS' (1D-12 '*' f2gd);
  407. LOGI2 = ('ABS' ((f2 '*' ORIENT) '-' f2gd)) < ERRO;
  408. LOGI1 = LOGI1 'ET' LOGI2;
  409.  
  410. f3 = 'EXTRAIRE' FLUX2T 'SCAL' P1;
  411. ERRO = 'ABS' (1D-11 '*' f3gd) ;
  412. LOGI2 = ('ABS' ((f3 '*' ORIENT) '-' f3gd)) < ERRO;
  413. LOGI1 = LOGI1 'ET' LOGI2;
  414.  
  415. f4 = 'EXTRAIRE' FLUX3 'SCAL' P1;
  416. ERRO = 'ABS' (1D-12 '*' f4gd) ;
  417. LOGI2 = ('ABS' ((f4 '*' ORIENT) '-' f4gd)) < ERRO;
  418. LOGI1 = LOGI1 'ET' LOGI2;
  419.  
  420. f5 = 'EXTRAIRE' FLUX4 'H2  ' P1;
  421. ERRO = 'ABS' (1D-12 '*' f5gd) ;
  422. LOGI2 = ('ABS' ((f5 '*' ORIENT) '-' f5gd)) < ERRO;
  423. LOGI1 = LOGI1 'ET' LOGI2;
  424.  
  425. f6 = 'EXTRAIRE' FLUX4 'O2  ' P1;
  426. ERRO = 'ABS' (1D-12 '*' f6gd) ;
  427. LOGI2 = ('ABS' ((f6 '*' ORIENT) '-' f6gd)) < ERRO;
  428. LOGI1 = LOGI1 'ET' LOGI2;
  429.  
  430. f7 = 'EXTRAIRE' FLUX4 'H2O ' P1;
  431. ERRO = 'ABS' (1D-12 '*' f7gd) ;
  432. LOGI2 = ('ABS' ((f7 '*' ORIENT) '-' f7gd)) < ERRO;
  433. LOGI1 = LOGI1 'ET' LOGI2;
  434.  
  435. 'SI' ('NON' LOGI1);
  436.     'MESSAGE' ;
  437.     'MESSAGE' 'OPERATEUR KONV';
  438.     'MESSAGE' ('CHAINE' METO);
  439.     'MESSAGE' ;
  440.     'MESSAGE' ('CHAINE' 'df1 ' ('ABS' (f1gd '-'  (ORIENT '*'f1)))
  441.                         ' f1 ' f1);                       
  442.     'MESSAGE' ('CHAINE' 'df2 ' ('ABS' (f2gd '-'  (ORIENT '*'f2)))                       
  443.                         ' f2 ' f2);                       
  444.     'MESSAGE' ('CHAINE' 'df3 ' ('ABS' (f3gd '-'  (ORIENT '*'f3)))                       
  445.                         ' f3 ' f3);                       
  446.     'MESSAGE' ('CHAINE' 'df4 ' ('ABS' (f4gd '-'  (ORIENT '*'f4)))                           
  447.                         ' f4 ' f4);
  448.     'MESSAGE' ('CHAINE' 'df5 ' ('ABS' (f5gd '-'  (ORIENT '*'f5)))                           
  449.                         ' f5 ' f5);
  450.     'MESSAGE' ('CHAINE' 'df6 ' ('ABS' (f6gd '-'  (ORIENT '*'f6)))                           
  451.                         ' f6 ' f6);
  452.     'MESSAGE' ('CHAINE' 'df7 ' ('ABS' (f7gd '-'  (ORIENT '*'f7)))                           
  453.                         ' f7 ' f7);
  454.     'ERREUR' 5;
  455. 'FINSI' ;
  456.  
  457. ************************************************
  458. ***** Splitting des scalaires passifs **********
  459. ************************************************
  460. *
  461. LISTINCO = 'MOTS' 'RN' 'RUX' 'RUY' 'RETN' 'H2' 'O2' 'H2O' 'TUTU' 'TATA'
  462. ;
  463.  
  464. TFLUX2 DELTAT =  'KONV' 'VF' 'PERFTEMP' 'FLUX' METO
  465.          $DOMTOT PGAZ2 LISTINCO  ROF2 VITF2 PF2 YF2 SCAF ;
  466.  
  467. TFLUX3 = 'EXCO' ( 'MOTS' 'RN' 'RUX' 'RUY' 'RETN' 'H2' 'O2' 'H2O' )
  468.   TFLUX2 ( 'MOTS' 'RN' 'RUX' 'RUY' 'RETN' 'H2' 'O2' 'H2O' ) ;
  469.  
  470. ERRO = 'MAXIMUM' (TFLUX3 '-' TFLUX) 'ABS' ;
  471.  
  472. 'SI' (ERRO > 1.0D-15);
  473.     'MESSAGE' 'Probleme dans les transport de scalaires passifs' ;
  474.     'ERREUR' 5;
  475. 'FINSI' ;
  476.  
  477. *
  478. **** Test 'KONV' 
  479. *
  480.  
  481. fceltutu = 'EXTRAIRE' TFLUX2  'TUTU' P1;
  482. LOGI1 = ('ABS' ((fceltutu '*' ORIENT) '-' ftutugd)) <
  483.         ('ABS' (1.0D-12 * fceltutu)) ;
  484.  
  485. fceltata = 'EXTRAIRE' (TFLUX2)  'TATA' P1;
  486. LOGI2 = ('ABS' ((fceltata '*' ORIENT) '-' ftatagd)) <
  487.         ('ABS' (1.0D-12 * fceltata)) ;
  488. LOGI1 = LOGI1 'ET' LOGI2;
  489.  
  490. 'SI' ('NON' LOGI1);
  491.     'MESSAGE' ;
  492.     'MESSAGE' 'OPERATEUR KONV';
  493.     'MESSAGE' 'Transport de scalaires passifs' ;
  494.     'ERREUR' 5;
  495. 'FINSI' ;
  496.  
  497.  
  498. ****************************************************
  499. ****************************************************
  500. ******** Fin boucle sur les angles         *********
  501. ****************************************************
  502. ****************************************************
  503.  
  504. 'FIN' BLOC;
  505.  
  506. 'FIN' ;
  507.  
  508.  
  509.  
  510.  
  511.  
  512.  
  513.  
  514.  
  515.  
  516.  
  517.  
  518.  

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