Télécharger konv_resi_dem_constant_state_11.dgibi

Retour à la liste

Numérotation des lignes :

  1. * fichier : konv_resi_dem_constant_state_11.dgibi
  2. ************************************************************************
  3. ************************************************************************
  4. ***********************************************************
  5. **** APPROCHE VF "Cell-Centred Formulation" pour la ****
  6. **** solution des ****
  7. **** Equations d'Euler pour un gaz parfait ****
  8. **** Approche DEM pour la combustion ****
  9. **** OPERATEURS PRET, KONV ****
  10. **** ****
  11. **** Consistency ****
  12. **** Left state = phase 1 ****
  13. **** Right state = phase 1 ****
  14. **** ****
  15. **** Methodes: VLH, SS, AUSMPUP ****
  16. **** ****
  17. **** A. BECCANTINI DM2S/SFME/LTMF DECEMBRE 2009 ****
  18. ***********************************************************
  19.  
  20.  
  21. 'OPTION' 'DIME' 2 ;
  22. 'OPTION' 'ELEM' QUA4 ;
  23. 'OPTION' 'ECHO' 1 ;
  24. 'OPTION' 'TRAC' 'X' ;
  25.  
  26. *
  27. *** GRAPH
  28. *
  29.  
  30. GRAPH = FAUX ;
  31. * GRAPH = VRAI ;
  32.  
  33. **** Cas homogene : la table proprieté de gaz
  34.  
  35. *************************************************
  36. **** The table for the properties of the gas ****
  37. *************************************************
  38. *
  39. PGAS = 'TABLE' ;
  40. *
  41. **** Order of the polynomial order for cv = cv(T)
  42. * For T > TMAX, cv(T) = cv(Tmax)
  43. *
  44. PGAS . 'TMAX' = 6000.0 ;
  45. PGAS . 'NORD' = 4 ;
  46. *
  47. **** Species involved in the mixture (before or after
  48. * the chemical reaction)
  49. *
  50. PGAS . 'SPECIES' = 'MOTS' 'H2 ' 'O2 ' 'H2O ' 'N2 ' ;
  51. *
  52. *
  53. **** Coefficient of the chemical reaction.
  54. * Note that for the first species this coefficient should be positive
  55. * Normal, we take it equal to 1.
  56. *
  57. * H2 '+' 0.5 O2 ---> H2O
  58. *
  59. PGAS . 'CHEMCOEF' = 'PROG' 1.0 0.5 -1.0 0.0 ;
  60. *
  61. **** Mass fraction of the first species before and after the combustion
  62. * Final mass fractions of the species with positive coefficients.
  63. * Final mass fractions of the species with non-positive coefficient.
  64. * The mass fraction of the last species is not given.
  65. *
  66. PGAS . 'MASSFRA' = 'PROG' 0.285219E-01 0.964039E-11 0.765104E-10
  67. 0.127442E-10 ;
  68. PGAS . 'MASSFRA' = 'PROG' 0.285219E-01 0.964039E-11 0.765104E-10
  69. 0.127442E-10 ;
  70. *
  71. **** Coef with the gas properties
  72. *
  73. PGAS . 'H2 ' = 'TABLE' ;
  74. PGAS . 'H2O ' = 'TABLE' ;
  75. PGAS . 'N2 ' = 'TABLE' ;
  76. PGAS . 'O2 ' = 'TABLE' ;
  77. *
  78. **** Runiv (J/mole/K)
  79. *
  80. PGAS . 'RUNIV' = 8.31441 ;
  81. *
  82. **** W (kg/mole). Gas constant (J/kg/K = Runiv/W)
  83. *
  84. PGAS . 'H2 ' . 'W' = 2.016E-3 ;
  85. PGAS . 'O2 ' . 'W' = 31.999E-3 ;
  86. PGAS . 'H2O ' . 'W' = 18.0155E-3 ;
  87. PGAS . 'N2 ' . 'W' = 28.013E-3 ;
  88. *
  89. **** Polynomial coefficients
  90. *
  91. PGAS . 'H2 ' . 'A' = 'PROG' 9834.91866 0.54273926 0.000862203836
  92. -2.37281455E-07 1.84701105E-11 ;
  93. PGAS . 'H2O ' . 'A' = 'PROG' 1155.95625 0.768331151 -5.73129958E-05
  94. -1.82753232E-08 2.44485692E-12 ;
  95. PGAS . 'N2 ' . 'A' = 'PROG' 652.940766 0.288239099 -7.80442298E-05
  96. 8.78233606E-09 -3.05514485E-13 ;
  97. PGAS . 'O2 ' . 'A' = 'PROG' 575.012333 0.350522002 -0.000128294865
  98. 2.33636971E-08 -1.53304905E-12;
  99. *
  100. **** Formation enthalpies (energies) at 0K (J/Kg)
  101. *
  102. PGAS . 'H2 ' . 'H0K' = -4.195D6 ;
  103. PGAS . 'H2O ' . 'H0K' = -1.395D7 ;
  104. PGAS . 'N2 ' . 'H0K' = -2.953D5 ;
  105. PGAS . 'O2 ' . 'H0K' = -2.634D5 ;
  106.  
  107.  
  108. ***************************
  109. ***** DOMAINE SPATIAL ****
  110. ***************************
  111.  
  112. A1 = 0.0D0 0.01D0;
  113. A2 = 1.011D0 0.21D0;
  114. A3 = 2.12D0 0.021D0;
  115. A4 = 2.11D0 1.09D0;
  116. A5 = 1.023D0 1.132D0;
  117. A6 = 0.098D0 1.199D0;
  118.  
  119. L12 = A1 'DROIT' 1 A2;
  120. L23 = A2 'DROIT' 1 A3;
  121. L34 = A3 'DROIT' 1 A4;
  122. L45 = A4 'DROIT' 1 A5;
  123. L56 = A5 'DROIT' 1 A6;
  124. L61 = A6 'DROIT' 1 A1;
  125. L25 = A2 'DROIT' 1 A5;
  126.  
  127. DOM10 = 'DALL' L12 L25 L56 L61
  128. 'PLANE';
  129. DOM20 = 'DALL' L23 L34 L45 ('INVERSE' L25)
  130. 'PLANE';
  131.  
  132. *
  133. *** Point ou on controlle la consistence
  134. *
  135.  
  136. P10 = 1.0 0.5;
  137.  
  138. *
  139. EPS = 1.0D-16 ;
  140. K0 = 10.0 ;
  141. *
  142. *******************************************
  143. *** Etats gauche et droite ****************
  144. *******************************************
  145. *
  146. * There are 3 fluxes to control :
  147. * 1 -> 2 can be checked by considering an exact reactive shock
  148. * 2 -> 1 can be checked by considering an exact reactive shock
  149. * 1 -> 1 can be checked via consistency (or via an exact non-reactive
  150. * or by comparing the results with a non-reactive one)
  151. * 2 -> 2 can be checked via consistency (or via an exact non-reactive
  152. * or by comparing the results with a non-reactive one)
  153. *
  154. * Case 1 -> 1
  155. * Constant state
  156. *
  157. tg1 = 293.16 ;
  158. td1 = 293.16 ;
  159. pg1 = 1.023e5 ;
  160. pd1 = 1.023e5 ;
  161. ung1 = 122.0 ;
  162. und1 = 122.0 ;
  163. utg1 = 101.0 ;
  164. utd1 = 101.0 ;
  165. * utg1 = 0.0 ;
  166. * utd1 = 0.0 ;
  167. tg2 = 2800.15 ;
  168. td2 = 2800.15 ;
  169. pg2 = 1.023e5 ;
  170. pd2 = 1.023e5 ;
  171. ung2 = 150.0 ;
  172. und2 = 150.0 ;
  173. utg2 = 160.0 ;
  174. utd2 = 160.0 ;
  175.  
  176. TN1 = 'PROG' tg1 td1 ;
  177. TN2 = 'PROG' tg2 td2 ; ;
  178. PN1 = 'PROG' pg1 pd1 ;
  179. PN2 = 'PROG' pg2 pd2 ;
  180. UN1 = 'PROG' ung1 und1 ;
  181. UN2 = 'PROG' ung2 und2 ;
  182. UT1 = 'PROG' utg1 utd1 ;
  183. UT2 = 'PROG' utg2 utd2 ;
  184. ALPHA1 = 'PROG' 1.0 1.0 ;
  185. ALPHA2 = ('PROG' ('DIME' ALPHA1) '*' 1.0) '-' ALPHA1 ;
  186.  
  187. *
  188. * ETHER = int_0^T cv(T') dT' T < TMAX
  189. * = int_0^TMAX cv(T') dT' '+'
  190. * cv(TMAX) T >= TMAX
  191. *
  192.  
  193. ESP1 = 'EXTRAIRE' (PGAS . 'SPECIES') 1 ;
  194. * DY1 = y_i - y_f for the species 1
  195. DY1 = (('EXTRAIRE' (PGAS . 'MASSFRA') 1) '-'
  196. ('EXTRAIRE' (PGAS . 'MASSFRA') 2)) ;
  197. COEF1 = ('EXTRAIRE' (PGAS . 'CHEMCOEF') 1) '*'
  198. (PGAS . ESP1 . 'W') ;
  199. YFINPH1 = 1.0 ;
  200. YFINPH2 = 1.0 ;
  201. 'SI' (COEF1 > 0) ;
  202. YPH2 = 'EXTRAIRE' (PGAS . 'MASSFRA') 2 ;
  203. YPH1 = YPH2 '+' DY1 ;
  204. 'SINON' ;
  205. YPH1 = 'EXTRAIRE' (PGAS . 'MASSFRA') 2 ;
  206. YPH2 = YPH1 '-' DY1 ;
  207. 'FINSI' ;
  208. YFINPH1 = YFINPH1 '-' YPH1 ;
  209. YFINPH2 = YFINPH2 '-' YPH2 ;
  210. PRYPH1 = 'PROG' YPH1 ;
  211. PRYPH2 = 'PROG' YPH2 ;
  212. 'REPETER' BLESP (('DIME' (PGAS . 'SPECIES')) '-' 2) ;
  213. ESP = 'EXTRAIRE' (PGAS . 'SPECIES') (&BLESP '+' 1) ;
  214. COEF = ('EXTRAIRE' (PGAS . 'CHEMCOEF') (&BLESP '+' 1))
  215. '*' (PGAS . ESP . 'W') ;
  216. DY = (DY1 * (COEF '/' COEF1)) ;
  217. 'SI' (COEF > 0) ;
  218. YPH2 = 'EXTRAIRE' (PGAS . 'MASSFRA') (&BLESP '+' 2) ;
  219. YPH1 = YPH2 '+' DY ;
  220. 'SINON' ;
  221. YPH1 = 'EXTRAIRE' (PGAS . 'MASSFRA') (&BLESP '+' 2) ;
  222. YPH2 = YPH1 '-' DY ;
  223. 'FINSI' ;
  224. PRYPH1 = PRYPH1 'ET' ('PROG' YPH1) ;
  225. PRYPH2 = PRYPH2 'ET' ('PROG' YPH2) ;
  226. YFINPH1 = YFINPH1 '-' YPH1 ;
  227. YFINPH2 = YFINPH2 '-' YPH2 ;
  228. 'FIN' BLESP ;
  229. PRYPH1 = PRYPH1 'ET' ('PROG' YFINPH1) ;
  230. PRYPH2 = PRYPH2 'ET' ('PROG' YFINPH2) ;
  231. 'LISTE' PRYPH1 ;
  232. 'LISTE' PRYPH2 ;
  233. *
  234. TMAX = ('PROG' ('DIME' ALPHA1) '*' (PGAS . 'TMAX')) ;
  235. * TCAL1 = MIN TN1, TMAX
  236. TCAL1 = 0.5D0 '*' ((TMAX '+' TN1) '-' ('ABS' (TN1 '-' TMAX))) ;
  237. DTN1 = TN1 '-' TCAL1 ;
  238. * TCAL1 = MIN TN1, TMAX
  239. TCAL2 = 0.5D0 '*' ((TMAX '+' TN2) '-' ('ABS' (TN2 '-' TMAX))) ;
  240. DTN2 = TN2 '-' TCAL2 ;
  241. *
  242. * Internal energy (J/kg in SI)
  243. *
  244. ETHER1 = ('PROG' ('DIME' ALPHA1) '*' 0.0) ;
  245. CV1 = ('PROG' ('DIME' ALPHA1) '*' 0.0) ;
  246. ETHER2 = ('PROG' ('DIME' ALPHA1) '*' 0.0) ;
  247. CV2 = ('PROG' ('DIME' ALPHA1) '*' 0.0) ;
  248. FUNTN1 = ('PROG' ('DIME' ALPHA1) '*' 1.0) ;
  249. FUNTN2 = ('PROG' ('DIME' ALPHA1) '*' 1.0) ;
  250. 'REPETER' BLPO ((PGAS . 'NORD') '+' 1) ;
  251. 'REPETER' BLESP ('DIME' (PGAS . 'SPECIES')) ;
  252. ESP = 'EXTRAIRE' (PGAS . 'SPECIES') &BLESP ;
  253. YCEL1 = 'EXTRAIRE' PRYPH1 &BLESP ;
  254. YCEL2 = 'EXTRAIRE' PRYPH2 &BLESP ;
  255. AA = 'EXTRAIRE' (PGAS . ESP . 'A') &BLPO ;
  256. DCV1 = (AA * YCEL1 * FUNTN1) ;
  257. DCV2 = (AA * YCEL2 * FUNTN2) ;
  258. CV1 = CV1 '+' DCV1 ;
  259. CV2 = CV2 '+' DCV2 ;
  260. ETHER1 = ETHER1 '+' (DCV1 * TCAL1 '/' (&BLPO)) ;
  261. ETHER2 = ETHER2 '+' (DCV2 * TCAL2 '/' (&BLPO)) ;
  262. 'FIN' BLESP ;
  263. FUNTN1 = FUNTN1 '*' TCAL1 ;
  264. FUNTN2 = FUNTN2 '*' TCAL2 ;
  265. 'FIN' BLPO ;
  266. ETHER1 = ETHER1 '+' (CV1 '*' DTN1) ;
  267. ETHER2 = ETHER2 '+' (CV2 '*' DTN2) ;
  268. *
  269. * Formation energy/enthalpy (J/kg in SI) and gas constant (J/kg/K)
  270. *
  271. EFORM1 = 0.0 ;
  272. EFORM2 = 0.0 ;
  273. RGAS1 = 0.0 ;
  274. RGAS2 = 0.0 ;
  275. 'REPETER' BLESP ('DIME' (PGAS . 'SPECIES')) ;
  276. ESP = 'EXTRAIRE' (PGAS . 'SPECIES') &BLESP ;
  277. YCEL1 = 'EXTRAIRE' PRYPH1 &BLESP ;
  278. YCEL2 = 'EXTRAIRE' PRYPH2 &BLESP ;
  279. EFORM1 = EFORM1 '+' (YCEL1 * (PGAS . ESP . 'H0K')) ;
  280. EFORM2 = EFORM2 '+' (YCEL2 * (PGAS . ESP . 'H0K')) ;
  281. RGAS1 = RGAS1 '+' (YCEL1 * (PGAS . 'RUNIV') '/'
  282. (PGAS . ESP . 'W')) ;
  283. RGAS2 = RGAS2 '+' (YCEL2 * (PGAS . 'RUNIV') '/'
  284. (PGAS . ESP . 'W')) ;
  285. 'FIN' BLESP ;
  286. EFORM1 = 'PROG' ('DIME' ETHER1) '*' EFORM1 ;
  287. EFORM2 = 'PROG' ('DIME' ETHER2) '*' EFORM2 ;
  288.  
  289. *
  290. * Computation of the conservative variables
  291. *
  292. RN1 = PN1 '/' (RGAS1 '*' TN1) ;
  293. RN2 = PN2 '/' (RGAS2 '*' TN2) ;
  294. GNX1 = RN1 * UN1 ;
  295. GNY1 = RN1 * UT1 ;
  296. GNX2 = RN2 * UN2 ;
  297. GNY2 = RN2 * UT2 ;
  298. ECIN1 = 0.5D0 '*' ('PROG' ((ung1 * ung1) '+' (utg1 * utg1))
  299. ((und1 * und1) '+' (utd1 * utd1)));
  300. ECIN2 = 0.5D0 '*' ('PROG' ((ung2 * ung2) '+' (utg2 * utg2))
  301. ((und2 * und2) '+' (utd2 * utd2)));
  302. RETN1 = RN1 '*' (ETHER1 '+' ECIN1 '+' EFORM1) ;
  303. RETN2 = RN2 '*' (ETHER2 '+' ECIN2 '+' EFORM2) ;
  304. *
  305. * Computation of the resi contribution in (n,t)
  306. *
  307. retg1 = MAXI RETN1 ;
  308. rhog1 = 'MAXIMUM' RN1 ;
  309. f1gd = ung1 '*' rhog1 ;
  310. f2gd = (f1gd '*' ung1) '+' pg1 ;
  311. f3gd = f1gd '*' utg1 ;
  312. f4gd = ung1 '*' (retg1 '+' pg1);
  313.  
  314. ****************************************************
  315. ****************************************************
  316. ******** Boucle sur les angles *********
  317. ****************************************************
  318. ****************************************************
  319.  
  320. DANGLE = 360 '/' 7.15;
  321. ANGLE = 11.3 ;
  322.  
  323. 'REPETER' BLOC 8;
  324.  
  325. *
  326. *** Rotation
  327. *
  328.  
  329. ANGLE = ANGLE '+' DANGLE;
  330. ORIG = 0.0D0 0.0D0;
  331.  
  332. 'MESSAGE' ;
  333. 'MESSAGE' (CHAIN 'Angle de rotation= ' ANGLE);
  334. 'MESSAGE' ;
  335.  
  336. DOM1 = DOM10 'TOURNER' ANGLE ORIG;
  337. DOM2 = DOM20 'TOURNER' ANGLE ORIG;
  338. P1FAC = P10 'TOURNER' ANGLE ORIG;
  339.  
  340. DOMTOT = DOM1 ET DOM2;
  341. 'ELIMINATION' DOMTOT 1D-6;
  342.  
  343. $DOMTOT = 'MODELISER' DOMTOT 'EULER';
  344. $DOM1 = 'MODELISER' DOM1 'EULER';
  345. $DOM2 = 'MODELISER' DOM2 'EULER';
  346. TDOMTOT = 'DOMA' $DOMTOT 'VF';
  347. TDOM1 = 'DOMA' $DOM1 'VF';
  348. TDOM2 = 'DOMA' $DOM2 'VF';
  349. MDOM1 = TDOM1 . 'QUAF' ;
  350. MDOM2 = TDOM2 . 'QUAF' ;
  351. MDOMTOT = TDOMTOT . 'QUAF' ;
  352. GRALP1 = 'NOMC' ('MOTS' 'UX' 'UY') ('DOMA' $DOMTOT 'XXNORMAF')
  353. ('MOTS' 'P1DX' 'P1DY') ;
  354.  
  355. 'ELIM' (MDOMTOT 'ET' MDOM1 'ET' MDOM2) 1.E-6 ;
  356.  
  357. *
  358. VINF = 'MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' 100. ;
  359. *
  360. **** Redefinition de P1FAC dans $DOMTOT 'FACE'
  361. *
  362.  
  363. P1FAC = ('DOMA' $DOMTOT 'FACE') 'POIN' 'PROC' P1FAC ;
  364. XSURF = 'EXTRAIRE' ('DOMA' $DOMTOT 'XXSURFAC') P1FAC 'SCAL' ;
  365. NX = 'EXTRAIRE' ('DOMA' $DOMTOT 'XXNORMAF') P1FAC 'UX' ;
  366. NY = 'EXTRAIRE' ('DOMA' $DOMTOT 'XXNORMAF') P1FAC 'UY' ;
  367. TX = -1 * NY ;
  368. TY = NX ;
  369.  
  370. GEOPC = ('DOMA' $DOMTOT 'FACEL') 'ELEM' 'APPUYE' 'LARGEMENT'
  371. P1FAC ;
  372. AA = 'CHANGER' GEOPC 'POI1' ;
  373. PC1 = 'POIN' 1 AA ;
  374. P2 = 'POIN' 1 ('DOMA' $DOMTOT 'CENTRE') ;
  375. PC2 = 'POIN' 3 AA ;
  376. * Si P3 n'existe pas, probleme en FACEL.
  377. 'SI' ('NEG' PC1 P2) ;
  378. ORIENT = -1 ;
  379. PCD = PC1 ;
  380. PCG = PC2 ;
  381. XVOLD = 'MAXIMUM' (TDOM1 . 'XXVOLUM') ;
  382. XVOLG = 'MAXIMUM' (TDOM2 . 'XXVOLUM') ;
  383. 'SINON' ;
  384. ORIENT = 1 ;
  385. PCD = PC2 ;
  386. PCG = PC1 ;
  387. XVOLG = 'MAXIMUM' (TDOM1 . 'XXVOLUM') ;
  388. XVOLD = 'MAXIMUM' (TDOM2 . 'XXVOLUM') ;
  389. 'FINSI' ;
  390.  
  391. uxg1 = ((ung1 '*' NX) '+' (utg1 '*' TX)) '*' ORIENT ;
  392. uyg1 = ((ung1 '*' NY) '+' (utg1 '*' TY)) '*' ORIENT ;
  393. uxg2 = ((ung2 '*' NX) '+' (utg2 '*' TX)) '*' ORIENT ;
  394. uyg2 = ((ung2 '*' NY) '+' (utg2 '*' TY)) '*' ORIENT ;
  395. uxd1 = ((und1 '*' NX) '+' (utd1 '*' TX)) '*' ORIENT ;
  396. uyd1 = ((und1 '*' NY) '+' (utd1 '*' TY)) '*' ORIENT ;
  397. uxd2 = ((und2 '*' NX) '+' (utd2 '*' TX)) '*' ORIENT ;
  398. uyd2 = ((und2 '*' NY) '+' (utd2 '*' TY)) '*' ORIENT ;
  399.  
  400. 'SI' GRAPH;
  401. 'TRACER' (('DOMA' $DOMTOT 'MAILLAGE') 'ET'
  402. ('COULEUR' ('DOMA' $DOMTOT 'FACEL') 'ROUGE') 'ET' P1FAC)
  403. 'TITRE' 'Domaine et FACEL';
  404. 'FINSI' ;
  405.  
  406. MAILIM = 'DIFF' ('DOMA' $DOMTOT 'FACE') ('MANUEL' 'POI1' P1FAC) ;
  407.  
  408. ***********************
  409. **** Les CHPOINTs ****
  410. ***********************
  411.  
  412. CHRN1 = ('MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' RN1) ;
  413. CHRN2 = ('MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' RN2) ;
  414.  
  415. CHVN1 = ('MANUEL' 'CHPO' (TDOM1 . 'CENTRE') 2 'UX' uxg1
  416. 'UY' uyg1) '+'
  417. ('MANUEL' 'CHPO' (TDOM2 . 'CENTRE') 2 'UX' uxd1
  418. 'UY' uyd1) ;
  419. CHVN2 = ('MANUEL' 'CHPO' (TDOM1 . 'CENTRE') 2 'UX' uxg2
  420. 'UY' uyg2) '+'
  421. ('MANUEL' 'CHPO' (TDOM2 . 'CENTRE') 2 'UX' uxd2
  422. 'UY' uyd2) ;
  423.  
  424. CHGN1 = CHRN1 '*' CHVN1 ;
  425. CHGN2 = CHRN2 '*' CHVN2 ;
  426.  
  427. CHRET1 = ('MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' RETN1) ;
  428. CHRET2 = ('MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' RETN2) ;
  429.  
  430. CHTN1 = ('MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' TN1) ;
  431. CHTN2 = ('MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' TN2) ;
  432.  
  433. CHAL1 = ('MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' ALPHA1) ;
  434. CHAL2 = ('MANUEL' 'CHPO' (TDOMTOT . 'CENTRE') 1 'SCAL' ALPHA2) ;
  435.  
  436. R1 R2 V1 V2 P1 P2 T1 T2 = 'PRIM' 'DEM' PGAS CHAL1 CHAL2
  437. (CHAL1 * CHRN1) (CHAL2 * CHRN2) (CHAL1 * CHGN1)
  438. (CHAL2 * CHGN2) (CHAL1 * CHRET1) (CHAL2 * CHRET2)
  439. CHTN1 CHTN2 EPS ;
  440.  
  441. GRADAL1 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 2
  442. 'P1DX' 0.0 'P1DY' 0.0 ;
  443. GRADAL2 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 2
  444. 'P1DX' 0.0 'P1DY' 0.0 ;
  445. LIMAL1 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 1
  446. 'P1' 0.0 ;
  447. LIMAL2 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 1
  448. 'P1' 0.0 ;
  449. GRADR1 = 'MANU/ 'CHPO' (TDOMTOT . 'CENTRE') 2
  450. 'P1DX' 0.0 'P1DY' 0.0 ;
  451. GRADR2 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 2
  452. 'P1DX' 0.0 'P1DY' 0.0 ;
  453. LIMR1 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 1
  454. 'P1' 0.0 ;
  455. LIMR2 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 1
  456. 'P1' 0.0 ;
  457. GRADV1 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 4
  458. 'P1DX' 0.0 'P1DY' 0.0 'P2DX' 0.0 'P2DY' 0.0 ;
  459. GRADV2 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 4
  460. 'P1DX' 0.0 'P1DY' 0.0 'P2DX' 0.0 'P2DY' 0.0 ;
  461. LIMV1 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 2
  462. 'P1' 0.0 'P2' 0.0 ;
  463. LIMV2 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 2
  464. 'P1' 0.0 'P2' 0.0 ;
  465. GRADP1 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 2
  466. 'P1DX' 0.0 'P1DY' 0.0 ;
  467. GRADP2 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 2
  468. 'P1DX' 0.0 'P1DY' 0.0 ;
  469. LIMP1 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 1
  470. 'P1' 0.0 ;
  471. LIMP2 = 'MANU' 'CHPO' (TDOMTOT . 'CENTRE') 1
  472. 'P1' 0.0 ;
  473.  
  474. *
  475. **** L'operateur 'PRET'
  476. *
  477.  
  478. CHFAL1 CHFAL2 CHFRN1 CHFRN2 CHFVN1 CHFVN2 CHFPN1 CHFPN2 =
  479. 'PRET' 'DEM' $DOMTOT
  480. CHAL1 (0.0 * GRADAL1) LIMAL1
  481. CHAL2 (0.0 * GRADAL2) LIMAL2
  482. R1 (0.0 * GRADR1) LIMR1
  483. R2 (0.0 * GRADR2) LIMR2
  484. V1 (0.0 * GRADV1) LIMV1
  485. V2 (0.0 * GRADV2) LIMV2
  486. P1 (0.0 * GRADP1) LIMP1
  487. P2 (0.0 * GRADP2) LIMP2 ;
  488.  
  489.  
  490. *
  491. **** L'operateur 'KONV'
  492. *
  493.  
  494. 'REPETER' BLMETO 3 ;
  495.  
  496. 'SI' ('EGA' &BLMETO 1) ;
  497. METO = 'SS' ;
  498. 'FINSI' ;
  499. 'SI' ('EGA' &BLMETO 2) ;
  500. METO = 'VLH' ;
  501. 'FINSI' ;
  502. 'SI' ('EGA' &BLMETO 3) ;
  503. METO = 'AUSMPUP' ;
  504. 'FINSI' ;
  505.  
  506. 'MESSAGE' ('CHAINE' 'METO = ' METO) ;
  507. LISTINC1 = ('MOTS' 'ALF1' 'RN1' 'RNX1' 'RUY1' 'RET1'
  508. 'ALF2' 'RN2' 'RUX2' 'RUY2' 'RET2') ;
  509.  
  510. SI ('EGA' METO 'AUSMPUP') ;
  511. CHPRES DT SURF = 'KONV' 'VF' 'DEM' 'RESI' METO 'CONS'
  512. $DOMTOT PGAS LISTINC1 CHAL1 CHAL2 CHFAL1 CHFAL2 CHFRN1 CHFRN2
  513. CHFVN1 CHFVN2 CHFPN1 CHFPN2 K0 GRALP1 EPS MAILIM VINF VINF ;
  514. 'SINON' ;
  515. CHPRES DT SURF = 'KONV' 'VF' 'DEM' 'RESI' METO 'CONS'
  516. $DOMTOT PGAS LISTINC1 CHAL1 CHAL2 CHFAL1 CHFAL2 CHFRN1 CHFRN2
  517. CHFVN1 CHFVN2 CHFPN1 CHFPN2 K0 GRALP1 EPS MAILIM ;
  518. 'FINSI' ;
  519.  
  520. * 'LISTE' ('CHAINE' 'FORMAT' '(E16.10)' ('EXTRAIRE' ETHER1 1)) ;
  521. * 'LISTE' ('CHAINE' 'FORMAT' '(E16.10)' ('EXTRAIRE' ECIN1 1)) ;
  522. * 'LISTE' ('CHAINE' 'FORMAT' '(E16.10)' ('EXTRAIRE' EFORM1 1)) ;
  523. * 'LISTE' ('CHAINE' 'FORMAT' '(E16.10)' (pg1)) ;
  524. * 'LISTE' ('CHAINE' 'FORMAT' '(E16.10)' ('EXTRAIRE' RN1 1)) ;
  525. * 'LISTE' ('CHAINE' 'FORMAT' '(E16.10)' (ung1)) ;
  526. * 'OPTION' DONN 5 ;
  527.  
  528. RESX1 = 'EXCO' ('EXTRAIRE' LISTINC1 2) CHPRES ;
  529. RESX2X = 'EXCO' ('EXTRAIRE' LISTINC1 3) CHPRES ;
  530. RESX2Y = 'EXCO' ('EXTRAIRE' LISTINC1 4) CHPRES ;
  531. RESX3 = 'EXCO' ('EXTRAIRE' LISTINC1 5) CHPRES ;
  532.  
  533. RESX2N = (RESX2X '*' NX) '+' (RESX2Y * NY);
  534. RESX2T = (RESX2X '*' TX) '+' (RESX2Y * TY);
  535.  
  536. f1 = 'EXTRAIRE' RESX1 'SCAL' PCD ;
  537. f1bis = 'EXTRAIRE' RESX1 'SCAL' PCG ;
  538. ERRO = 1D-8 '*' f1gd 'ABS' ;
  539. LOGI1 = ('ABS' ((f1 * XVOLD) + (f1bis * XVOLG))) < ERRO ;
  540. LOGI2 = ('ABS' ((f1gd * XSURF) + (f1bis * XVOLG))) < ERRO ;
  541. LOGI1 = LOGI1 'ET' LOGI2;
  542.  
  543. f2 = 'EXTRAIRE' RESX2N 'SCAL' PCD ;
  544. f2bis = 'EXTRAIRE' RESX2N 'SCAL' PCG ;
  545. ERRO = 1D-8 '*' f2gd 'ABS' ;
  546. LOGI2 = ('ABS' ((f2 * XVOLD) + (f2bis * XVOLG))) < ERRO ;
  547. LOGI1 = LOGI1 'ET' LOGI2;
  548. LOGI2 = ('ABS' ((f2gd * XSURF) + (f2bis * XVOLG))) < ERRO ;
  549. LOGI1 = LOGI1 'ET' LOGI2;
  550.  
  551. f3 = 'EXTRAIRE' RESX2T 'SCAL' PCD ;
  552. f3bis = 'EXTRAIRE' RESX2T 'SCAL' PCG ;
  553. ERRO = 1D-8 '*' f3gd 'ABS' ;
  554. LOGI2 = ('ABS' ((f3 * XVOLD) + (f3bis * XVOLG))) < ERRO ;
  555. LOGI1 = LOGI1 'ET' LOGI2;
  556. LOGI2 = ('ABS' ((f3gd * XSURF) + (f3bis * XVOLG))) < ERRO ;
  557. LOGI1 = LOGI1 'ET' LOGI2;
  558.  
  559. f4 = 'EXTRAIRE' RESX3 'SCAL' PCD ;
  560. f4bis = 'EXTRAIRE' RESX3 'SCAL' PCG ;
  561. ERRO = 5D-7 '*' f4gd 'ABS' ;
  562. LOGI2 = ('ABS' ((f4 * XVOLD) + (f4bis * XVOLG))) < ERRO ;
  563. LOGI1 = LOGI1 'ET' LOGI2;
  564. LOGI2 = ('ABS' ((f4gd * XSURF) + (f4bis * XVOLG))) < ERRO ;
  565. LOGI1 = LOGI1 'ET' LOGI2;
  566.  
  567. 'SI' ('NON' LOGI1);
  568. 'MESSAGE' ;
  569. 'MESSAGE' 'OPERATEUR KONV';
  570. 'MESSAGE' ('CHAINE' METO);
  571. 'MESSAGE' ;
  572. 'ERREUR' 5 ;
  573. 'FINSI' ;
  574.  
  575. 'FIN' BLMETO ;
  576.  
  577. ****************************************************
  578. ****************************************************
  579. ******** Fin boucle sur les angles *********
  580. ****************************************************
  581. ****************************************************
  582.  
  583. 'FIN' BLOC;
  584.  
  585. 'FIN' ;
  586.  
  587.  
  588.  
  589.  
  590.  
  591.  
  592.  
  593.  
  594.  
  595.  
  596.  
  597.  
  598.  
  599.  
  600.  
  601.  
  602.  
  603.  

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