Home

ECOLE DES MINES DE DOUAI

1. CALGRID LILLE EPISODE ETE 100 x 100 4 km LILLE ETE Run title 3 lines CALGRID MODEL CONTROL FILE Additional user comments NO CLOUDS USES 1315 CALGRID DAT WITH SH MIXING HEIGHTS NO EMISSIONS INPUT GROUP 0 Input files name CALMET c Pollux simulateur ETUDES lille2 ete Resultats CALMET DAT ICONFL c Pollux simulateur ETUDES 1lille2 ete Donnees icon dat l BCONFL x TCONFL 5 CLDDAT c Pollux simulateur ETUDES lille2 ete Donnees CLOUD DAT PTECYC x PTEMRB c Pollux simulateur ETUDES lille2 ete Resultats Calgrid ptemarb dat 1 PTEMOB x AREMFL c Pollux simulateur ETUDES lille2 ete Resultats Calgrid arem dat VDEPFL c Pollux simulateur ETUDES lille2 ete Donnees VD DAT GRDLST c Pollux simulateur ETUDES lille2 ete Resultats Calgrid CALGRID LST 1 CNCDAT c Pollux simulateur ETUDES lille2 ete Resultats Calgrid CONC DAT RESTRT c Pollux simulateur ETUDES lille2 ete Resultats Calgrid RESTART DA 1 DEPDAT BNEST All file names will be converted to lower case if LCFILES T Otherwise if LCFILES F file names will be converted to UPPER CASE T lower case LCFILES T F UPPER CASE IEND INPUT GROUP 1 General run control parameters Starting date Y
2. NANTES CASE 80 x 90 2 km meteorological grid wind amp met model Met stations used 21 surface 20 upper air 0 precip Run title 3 lines CALMET MODEL CONTROL FILI INPUT GROUP 0 Input and Output File Names Subgroup a Default Name Type File Name GEO DAT input l GEODAT C pollux Simulateur ETUDES ExEtude ExSimu Donnees G SURF DAT input l SRFDAT C pollux Simulateur ETUDES ExEtude ExSimu Donnees S CLOUD DAT input CLDDAT Y PRECIP DAT input PRCDAT MM4 DAT input MM4DAT WIDAT DA input WIDAT i output y METLST C pollux Simulateur ETUDES ExEtude ExSimu Resul output y METDAT C pollux Simulateur ETUDES ExEtude ExSimu Resul PACOUT DA output PACDAT All file names will be converted to lower case if LCFIL Otherwise if LCFILES F file names will be converted to UPP T lower case LCFILES F F UPPER CASE NUMBER OF UPPER AIR amp OVERWATER STATIONS Number of upper air stations NUSTA No default Number of overwater met stations NOWSTA No default END Subgroup b Upper air files one per station Default Name Type File Name UPO1 DAT input me UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu UPO2 DAT input 2 A UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu UPO3 DAT input S A UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UPO4 DAT input a UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu UPO5 DAT input O a UPDAT C Pollux Simulateur ETUDES
3. Defaults 0 75 0 25 WTEM 0 75 0 25 NOTE NZEM values must be entered and must add up to 1 0 Height m of each emission layer face ZFEM nzem 1 Default 0 0 50 100 ZFEM 0 50 100 NOTE NZEM 1 values must be entered Default values distribute 75 of area source mass below 50 m and 25 between 50 100 m Ecole des Mines de Douai D partement Energ tique Industrielle 45 IEMIREC is the first guess for record length for PTECYC direct access fil TEMIREC END INPUT GROUP 2 Grid control parameters Horizontal grid definition No X grid cells NX No default NX 25 No Y grid cells NY No default NY 25 Grid spacing DGRIDKM km No default DGRIDKM 4 Reference UTM coordinates km of SOUTHWEST corner of grid point 1 1 X coordinate XORIGKM No default XORIGKM 454 Y coordinate YORIGKM No default YORIGKM 5558 UTM zone IUTMZN No default TUTMZN 31 Reference coordinates of CENTER of the domain used in the calculation of solar elevation angles Latitude deg XLAT No default XLAT 47 25 Longitude deg XLONG No default XLONG 1 5 Time zone XTZ No default XTZ 1 PST 8 MST 7 CST 6 EST 5 Vertical grid definition No of vertical layers in the CALGRID meteorological grid NZM No default NZM 5 Vertical CA
4. IFORMO Default 1 IFORMO 1 1 CALPUFF CALGRID type file CALMET DAT 2 MESOPUFF II type file PACOUT DAT LINE PRINTER OUTPUT OPTIONS Print met fields LPRINT Default F LPRINT F F Do not print T Print NOTE parameters below control which met variables are printed Print Interval IPRINF in hours Default 1 IPRINF 1 Meteorological fields are printed every 6 hours Specify which layers of U V wind component to print IUVOUT NZ NOTE NZ values must be entered 0 Do not print 1 Print used only if LPRINT T Defaults NZ 0 Y IUVOUT 0 0 0 0 0 0 0 0 0 O Specify which levels of the W wind component to print NOTE W defined at TOP cell face 6 values IWOUT NZ NOTE NZ values must be entered 0 Do not print 1 Print used only if LPRINT T amp LCALGRD T Defaults NZ 0 IWOUT 0 0 0 0 0 0 O 0 0 O Specify which levels of the 3 D temperature field to print ITOUT NZ NOTE NZ values must be entered 0 Do not print 1 Print used only if LPRINT T amp LCALGRD T Defaults NZ 0 M ITOUT 0 Op O 0 07 05 UD 0 00 4 Ecole des Mines de Douai D partement Energ tique Industrielle 34 Specify which meteorological fields to print used only if LPRINT T Variable Print 0 do not print 1 print STABILITY 0 PGT stability class USTAR 0 Friction velocity
5. INPUT GROUP 1 General run control parameters Starting date Year IBYR No default IBYR 98 Month IBMO No default IBMO 1 Day IBDY No default IBDY 19 Hour IBHR No default IBHR 1 Ecole des Mines de Douai D partement Energ tique Industrielle 32 Base time zone IBTZ No default IBTZ 1 PST 08 MST 07 CST 06 EST 05 Length of run hours IRLG No default IRLG 24 Run type IRTYPE Default 1 IRTYPE 1 0 Computes wind fields only 1 Computes wind fields and micrometeorological variables u w Dy zi ete IRTYPE must be 1 to run CALPUFF or CALGRID Compute special data fields required by CALGRID i e 3 D fields of W wind components and temperature in additional to regular Default T LCALGRD T fields LCALGRD LCALGRD must be T to run CALGRID Flag to stop run after SETUP phase ITEST Default 2 ITEST 2 Used to allow checking of the model inputs files etc ITEST 1 STOPS program after SETUP phase ITEST 2 Continues with execution of COMPUTATIONAL phase after SETUP IEND INPUT GROUP 2 Grid control parameters HORIZONTAL GRID DEFINITION No X grid cells NX No default NX 80 No Y grid cells NY No default NY 90 GRID SPACING DGRIDKM No default DGRIDKM 2
6. lot un comportement mixte Cette classification des lots selon les donn es INSEE nous semble tr s int ressante du point de vue de la repr sentativit physique de la mod lisation En outre ces donn es INSEE permettent galement de d finir un Indice de D marrage a Froid idf qui intervient dans le calcul des missions dues au trafic routier un moteur chaud pollue diff remment d un moteur froid Pour ce qui est de notre simulation il n a pas t possible de suivre cette m thodologie par manque de ces donn es INSEE Nous avons donc employ un zonage fictif pr sentant n anmoins un degr de repr sentativit physique voir figure 15 associ des valeurs d idf similaires a celles de la simulation de Nantes Ces zones sont les suivantes 1 mixte 2 et 3 conomiques 4 domicile et 5 zone reste en jaune Nous avons appliqu ce m me zonage pour les 5 jours de simulation Ecole des Mines de Douai D partement Energ tique Industrielle 23 Figure 15 Zonage Economique Zonage juridique Toujours afin de tester cette fonctionnalit nous avons appliqu un zonage juridique N ayant pas de donn es r elles concernant ce zonage nous l avons appliqu tr s localement sur quelques brins du r seau routier voir figure 16 afin de ne pas fausser les r sultats sur le r seau routier Figure 16 Zonage juridique Ecole des Mines de Douai D partement Energ tique Industrielle 24
7. 291 4 286 4 907 0 1000 290 8 287 5 881 0 1250 290 2 287 6 856 0 1500 289 3 286 8 700 0 3000 283 0 286 9 9999 2001 98 8 712 10 10 1017 0 20298 5 261 2 1013 0 50 297 1 262 2 1008 0 100 296 8 263 3 990 0 250 295 4 262 3 962 0 500 293 4 260 4 935 0 750 292 0 261 5 908 0 1000 291 0 263 6 881 0 1250 290 0 268 6 856 0 1500 289 2 276 7 700 0 3000 284 0 280 8 Ecole des Mines de Douai D partement Energ tique Industrielle 52 ANNEXE 5 Visualisation SIG par agr gation spatiale Emissions de NOx Heure de pointe matin du 07 08 1998 Metropole lilloise spe WKH d E A Agr gation Spatiale Parc Total f Emission Totale g km NOx 07 08 98 8h00 9h00 EN y E 9 000 to 37 600 15 M3600to 9000 12 BM 1e800to 3600 14 E 800to 1800 15 o Oto 800 14 AM TL ES L an ik Es ST a Al B o TR COS Si Sei sous pe Y d PEDRO Ya E aaa A FAR FER AT TE HE DER PERRET MIRA VEA ao LAN Spr er aH Pala ARS WAN AXE A AS fe mate VA 5 te as gt GA X AA aS q 7 A De an De is r Ecole des Mines de Douai D partement Energ tique Industrielle 53 ANNEXE 6 Visualisation SIG par agr gation spatiale Emissions de COV Heure de pointe matin du 07 08 1998 Metropole lilloise Agr gation Spatiale Parc Total Emission Totale g km COW 07 08 98 8h00 9h00 E 43 400 to 97 100 10 M 22700 to 43 400
8. Ecole des Mines de Douai D partement Energ tique Industrielle AVERTISSEMENT Le logiciel dont l valuation est pr sent e dans ce rapport est un produit en phase de d veloppement qui a volu et qui continue de l tre tout au long de cette tude et ce suite de nombreux changes avec la Soci t ACRI Ce logiciel initialement baptis Pollux a finalement t nomm SAMAA Par la suite nous ne ferons pas le distinguo entre ces deux noms Ce logiciel englobe en outre les logiciels Calmet et Calgrid pour les calculs respectifs de la m t orologie et de la chimie dispersion L valuation du code sur le calcul photochimique n a pas pu tre enti rement finalis e et certains points demandent tre tudi s plus finement Ecole des Mines de Douai D partement Energ tique Industrielle RESUME de l tude n 1 Etude r alis e par Jean Philippe VERMEULEN Tel 03 27 71 23 94 EVALUATION DU LOGICIEL POLLUX SAMAA Dans le cadre de Assistance maitrise d ouvrage en mod lisation nous avons proc d a l valuation du logiciel Pollux SAMAA Cette valuation a t men e suivant deux tapes compl mentaires La premi re fut de prendre en main le logiciel en effectuant une simulation existante concernant la r gion nantaise Nous pr sentons ainsi le fonctionnement global de ce logiciel au travers la pr sentation de cette phase La seconde fut d effectuer une simulation complete
9. MONIN 0 Monin Obukhov length MIXHT 0 Mixing height WSTAR 0 Convectiv PRECIP 0 Precipitation rate SENSHEAT 0 Sensible heat flux CONVZI 0 Convective mixing ht Testing and debug print options for micrometeorological module Print input meteorological data and internal variables LDB Default F F Do not print T print NOTE this option produces large amounts of output First time step for which debug data Defaults are printed NN1 Last time step are printed NN2 T for whic 0 Default 1 h debug data Default 1 Testing and debug print options for wind field module all of the following print options control output to wind field module s output files TEST PRT TEST OUT EST KIN TEST FRD and TEST SLP Control variable for writing the test debug wind fields to disk files IOUTD 0 Do not write l write Default 0 Number of levels starting at the surface to print NZPRN2 Default 1 Print the INTERPOLATED wind components IPRO O no l yes Default 0 Print the TERRAIN ADJUSTED SURFACE wind components IPR1 0 n0 l yes Default 0 l Print the SMOOTHED wind components and the INITIAL DIVERGENCE fields IPR2 0 n0 1 yes Default 0 l Print the FINAL wind speed And direction fields
10. Na a o E gt LOS UN aN ess y N Er JE VAR Oi A Etre A PoP Dee A A Ecole des Mines de Douai D partement Energ tique Industrielle 56 ANNEXE 9 Visualisation SIG l mentaire Emissions de NOx Heure de pointe matin du 07 08 1998 Metropole lilloise Emissions Parc Total Emission Totale g km NOx 07 08 98 8h00 9h00 3 300 to 23 700 1085 600to 3300 587 1900to 2600 893 o 1900 681 Oto 1500 666 o 1200 786 o 800 967 de PENSARA AS VA NZ ANT N NAS os X is ree wb EP PANG iy HS Ecole des Mines de Douai D partement Energ tique Industrielle 57 ANNEXE 10 Visualisation SIG l mentaire Emissions de NOx Heure de pointe soir du 07 08 1998 Metropole lilloise AS ANA aN Emissions Parc Total Emission Totale g km NOx 07 08 98 18h00 19h00 2 200 to 17 100 1261 Ps MS QA ne mek O Ors x RS RL LN A Bs 7 CAN DIO ESE A 4 AN a a NW j oy nig O gt lt NAS Y Ecole des Mines de Douai D partement Energ tique Industrielle 58 ANNEXE 11 Visualisation SIG par agr gation spatiale Emissions de NOx sources industrielles Metropole lilloise 0 am UD ina MA GA AL RU LA MAISON Agr
11. videmment et c est un moyen de v rifier la coh rence des calculs effectu s par Calmet que les champs sont uniformes puisque a partir d une station de surface nous en avons recr e deux autres comportant exactement les m mes donn es voir section 4 3 1 Les annexes 15 a 18 pr sentent respectivement les concentrations de O3 NO2 NO et xyl ne pour le 11 08 1998 entre 14 et 15 h pour une altitude de 10 m 4 4 Comparaison de r sultats A ce stade de l tude nous avons entrepris une comparaison entre les r sultats obtenus dans cette simulation et ceux pr sent s dans le rapport d activit n 3 de l Ecole des Mines de Douai de 1997 relatif aux inventaires d missions Cette comparaison porte sur les missions de NOx et de COV issus du trafic routier et concerne la maille associ e au centre de Lille Les annexes 19 et 20 pr sentent respectivement les missions annuelles de NOx et de COV en t an km pour 1996 Les annexes 21 et 22 pr sentent respectivement les valeurs de ces polluants calcul es par SAMAA Bien videmment il est quelque peu hasardeux de comparer ces deux jeux de r sultats dates diff rentes moyennes annuelles pour les donn es de 1996 et cumul horaire pour la pr sente simulation discr tisation spatiale diff rente mais il se d gage clairement les tendances suivantes Au point de vue qualitatif les r sultats de SAMAA montrent que l on identifie clairement les trois principaux axes de circu
12. 100 300 600 1000 1 IEND INPUT GROUP 3 Species list DRY SPECIES MODELED ADVECTED EMITTED DEPOSITED NAME O NO 1 YES O NO 1 YES O NO 1 YES O NO 1 COMPUTED GAS 2 COMPUTED PARTICLE 3 USER SPECIFIED BUILD UP SPECIES ACTIVE SPECIES NO 1 1 1 0 NO2 1 5 ky de VOSS O Li HONO a O 2 A HNO3 7 E MS PNA y FAO FOF N205 7 oe HOE VOL gf NO3 7 m Oye L A HO2 1 Og CH t IEC dr chee Ala 600s il FORM pel PE CEE ALD2 j e o Waya PAN P po 0 Ay L X02 E O O C203 7 r ORS TO al CRO 7 E MGLY 1 OES MO 4 PAR E cl de CO A ETH 7 dent 0 OLE ge ele lap AO TOL oa alee rae O XYL 7 ee On OPEN s 0 1 CRES Os 40 74 TO2 CO AO 1 ROR r Oe 08 H202 7 locr TE ISOP tr Dye Or MEOH po lige Tepe O 4 ETOH po ll Tee ON A L 502 ye ge SL B r 08 pe Sbe LE Oh UNR ge ee se PQ ut ISOPRD 1 1 0 O NTR pe er Gs SOP r 4E02 ew tp dee CI MOR S ep 00 0 0 I L O Hye 072 2042 0 1 L OLD Lax OR O O7 Ot 1 XOZN Lp 030 0 r H2O Lo 07 MO 40 M END INPUT GROUP 4 Chemical parameters for dry deposition of gases SPECIES DIFFUSIVITY ALPHA REACTIVITY MESOPHYLL HENRY S LAW COEFFICIENT NAME STAR RESISTANCE CM 2 S S CM DIMENSIONLESS Ecole des Mines de Douai D partement Energ tique Industr
13. DKHSTB 6 Defaults 224 9 62547 82g Dz Ds 0 USED ONLY IF KHMODE 0 1 or 3 LL DRASTE 224 96 7 32343 Dir 0 0 NOTE if KHMODE 1 or 3 these DKHSTB values are scaled by wind speed VERTICAL DIFFUSIVITY Minimum vertical diffusivity m 2 s DKZMIN Default 1 0 DKZMIN 1 Vertical diffusivity m 2 s at the model top DKZTOP Default 0 0 DKZTOP 0 IEND Ecole des Mines de Douai D partement Energ tique Industrielle 50 ANNEXE 3 Fichier SURF DAT D but du fichier simulation sur Lille pisode d t Le texte entre sert ici de commentaire et ne doit pas fiqurer dans le fichier En caract re gras les valeurs mesur es a Oh 6h En rouge les valeurs interpol es lin airement 98 219 0 98 224 18 1 3 ann e jour julien et heure du d but de la simulation ann e jour julien et heure de fin de simulation fuseau horaire nombre de stations 1001 1002 1003 Label des stations 98 219 0 lann e jour julien et heure relatifs aux donn es 2 20 232 20 90 1 288 300 88 1020 00 O donn es station n 1 2 20 232 20 9999 9999 288 300 88 1020 00 O donn es station n 2 2 20 232 20 9999 9999 288 300 88 1020 00 O donn es station n 3 98 219 1 2 07 2993304 90 1 288 167 88 1020 00 0 2 07 233 02 9999 9999 288 167 88 10200010 2 07 233 32 9999 9999 236 167 88 1020 00 0 98 219 2 1 93 234 43 90 1 288 033 88 1020 00 10 1 93 234 43 9999 9999 266 033 88 1020 00
14. Industrielle 22 port a 4 7 millions de v hicules par jour Par contre la longueur moyenne des trajets n a pas pu tre modifi e X En fait il est possible de changer cette valeur au niveau de l interface et de la valider mais lors d une nouvelle ouverture de Pollux cette valeur est automatiquement remise a la valeur par d faut qui est de 10 km Nous supposons ainsi que c est cette derni re valeur qui est prise en compte lors du calcul des mission ce bug a t mentionn a ACRI Le fichier de temp rature g n r par Caltemp fichier t dat a t utilis pour le d marrage a froid Onglet Donn es urbaines Zonage activit Afin de tester cette option nous avons cr e un zonage pour lequel l activit a t consid r e comme nulle entre 1 et 2 h du matin pour le premier jour de simulation Zonage conomique La m thodologie de cr ation du zonage conomique pr conis e par la soci t ACRI est bas e sur l emploi de donn es INSEE nombres d actifs r sidant dans un lot et nombre d actifs travaillant dans cet lot disponibles dans un fichier SIRENE au format lisible par MapInfo voir le document de sp cification Description des fichiers d entr e Le caract re conomique d un lot se fait par rapport a une valeur seuil par exemple 75 il en est de m me pour le caract re domicile d une zone Il se peut ainsi qu un lot se voit attribuer les deux caract res Dans ce cas cet
15. Units km REFERENCE COORDINATES of SOUTHWEST corner of grid cell 1 1 X coordinate XORIGKM No default XORIGKM 534 56 Y coordinate YORIGKM No default YORIGKM 5143 93 Units km Latitude XLATO No default XLATO 46 Longitude XLONO No default XLONO 2 UTM ZONE IUTMZN Default 0 TUTMZN 30 LAMBERT CONFORMAL PARAMETERS Rotate input winds from true north to map north using a Lambert conformal projection LLCONF Default F LLCONF F Latitude of 1st standard parallel Default 30 XLAT1 35 Latitude of 2nd standard parallel Default 60 XLAT2 45 XLAT1 and XLAT2 in NH in SH Longitude RLONO Default 90 RLONO 75 used only if LLCONF T Positive W Hemisphere Negative E Hemisphere Origin Latitude RLATO Default 40 RLATO 40 used only if IPROG gt 2 Positive N Hemisphere Ecole des Mines de Douai D partement Energ tique Industrielle 33 Negative S Hemisphere Vertical grid definition No of vertical layers NZ No default NZ 10 Cell face heights in arbitrary vertical grid ZFACE NZ 1 No defaults Units m ZFACE 0 0 20 0 100 300 600 1000 1500 2000 3000 4000 5000 IEND INPUT GROUP 3 Output Options DISK OUTPUT OPTION Save met fields in an unformatted output file LSAVE Default T LSAV F Do not save T Save BI Il H Type of unformatted output file
16. d marrage zones r sidentielles et les zones d arriv e des v hicules zones d activit conomique Ceci permet de tenir compte des missions li es au d marrage froid des v hicules Y Le zonage l gislatif permet de d finir une zone pour laquelle il est possible de modifier la r partition du parc automobile par ex bloquer la circulation PL dans un centre ville Donn es de sp ciation s lection des profils de sp ciation pour les COVNM et NOx A Initialement la r solution de l cran 800x600 ne permettait pas l affichage complet de l interface de saisie des param tres voir figure 6 la partie inf rieure de l interface n tant pas visible notamment le champ des donn es de sp ciation ainsi que les boutons de validation et d annulation La r solution a t ensuite augment e en 1024x768 afin de r soudre ce probleme voir figure 7 Ecole des Mines de Douai D partement Energ tique Industrielle 12 SS eee Figure 7 Affichage en 1024x768 4 2 4 2 Chauffage R sidentiel Les param tres d entr e concernant le chauffage r sidentiel sont group s en 6 cat gories correspondant aux 6 onglets de la fen tre de l interface Habitat Besoins en chauffage Facteurs d missions Profils de Sp ciation Eau chaude sanitaire et Fichiers Ecole des Mines de Douai D partement Energ tique Industrielle 13 Onglet Habitat Surfaces moyennes habitables par type de logement et par p riod
17. de Calmet s est ensuite termin sans aucun probleme temps de calcul de 113 secondes pour 132h de simulation soit 5 jours 2 Ecole des Mines de Douai D partement Energ tique Industrielle 21 4 3 2 Donn es g n rales Les donn es g n rales sont ensuite saisies sous Pollux D but et fin de simulation en jj mm aa et non pas en jj mm hh Il est noter que m me si Calmet permet de faire des simulations sur un nombre non entier de jours par exemple 5 jours 1 2 Pollux fera la simulation sur le nombre entier imm diatement inf rieur ou gal au nombre de jours d cimal entr sous Calmet dans la limite de 15 jours impos e par Pollux Dans notre cas la simulation a donc t r duite a 5 jours Polluants consid r s CO NOx COV PM SO2 CO2 et CH4 4 3 3 Calcul des missions 4 3 3 1 Trafic routier Onglet G n ral Routier A ce niveau on sp cifie apr s cr ation l emplacement de trois fichiers Le fichier Davis txt contenant pour chaque brin du r seau routier la longueur la vitesse a vide la capacit horaire le type de route le d bit et la vitesse des v hicules La version actuelle de Pollux permet maintenant de g rer des capacit s et d bits nuls Le fichier d intersection Davis Maillage en brin contenant les informations relatives l intersection entre le maillage utilis pour l agr gation et les brins du r seau routier voir m thode d obtention de ce fichi
18. gation Spatiale TGAP Industries g km NOx 07 08 98 8h00 9h00 14 M 700to 2050 2 M1 360 to 1700 1 M1 020 to 1 360 1 M 680to 1020 2 O 340to 680 3 Oto 340 9 ae az ete eva UN ae LO a LP Di POP MO eck A re TA TR o eae COSI Says soy hes ICE RCA ANA gt PLA NT ST RUE RTS RAR sen ars E SS es x pe o ate KEIN SPRY 5 3 iS 0 h I T ae 4 LA N h f E E rt p d q oe ae Ecole des Mines de Douai D partement Energ tique Industrielle 59 ANNEXE 12 Visualisation SIG par agr gation spatiale Emissions de SO2 sources industrielles Metropole lilloise Agr gation Spatiale TGAP Industries g km 502 07 08 98 8h00 9h00 M1 560 to 6870 2 M 1 300 to 1 560 3 M 230 to 1 300 4 E 40to 290 3 Oto 40 3 LA RATES wi AM E LOLA ai q Sy i Lia x le 1 FRET Li NA pets 4d SS APA BAS 2 Mist Al EE BE EC Les RER HE LEE ES OR Re Se ete DITES a Toig es REC an AS A ok Y Ta es aM ates A A K Am K B K Ecole des Mines de Douai D partement Energ tique Industrielle 60 ANNEXE 13 Visualisation ACRPlot du champ de temp rature M tropole lilloise Direction Y km E 02 5578 5588 5598 5608 5618 5628 5638 5648 5658 I I T I 5568 LE 02 LE 02 LE 02 LE 02 LE 02 E 02
19. mobiles sources d missions surfaciques AREM DAT vitesse de d position VD DAT Fichiers de sortie Fichier de v rification CALGRID LST fichier de sortie binaire CONC DAT Param tres g n raux Titre date de d but et dur e de la simulation esp ces chimiques mises en jeu dans la simulation nombre total esp ces advect es d pos es mises Types de fichiers d entr e Types de conditions limites concernant les variables m t o et les esp ces advect es fichiers relatifs aux sources d missions stationnaires ponctuelles constantes ou variables ponctuelles mobiles surfaciques fonction de distribution pour les sources surfaciques Maillage Projection D finition du maillage utilis e grille horizontale et grille verticale Esp ces chimiques Tableau contenant le nom des esp ces chimiques consid r es si elles sont mod lis es 0 oui 1 non advect es 0 oui 1 non mises 0 oui 1 non et d pos es diff rents choix possibles concernant le calcul D position des gaz Param tres chimiques pour la d position s che des gaz pour diff rentes esp ces chimiques D position des particules Param tres de taille pour la d position s cher des particules pas d esp ces chimiques de type particules dans cette mod lisation D position s che M thode d interpolation pour le calcul Options de sortie D finition des variables a enregistrer ai
20. sur l agglom ration lilloise Cette seconde phase fut l occasion de localiser et de r soudre certains bugs et d approfondir la phase pr c dente Nous pr sentons ainsi le cheminement que nous avons suivi afin de mener bien cette simulation Certains r sultats sont ensuite pr sent s et compar s afin de se prononcer sur leur repr sentativit physique Nous pr sentons galement certaines recommandations relatives la stabilit du logiciel qui reste encore un produit jeune ainsi que ses aspects positifs Ecole des Mines de Douai D partement Energ tique Industrielle 1 INTRODUCTION L objectif du travail pr sent dans ce rapport est de proc der a une valuation rapide du logiciel Pollux dans sa version disponible version 1 0 Dans un premier temps cette valuation a t faite a partir d une simulation effectu e par la soci t ACRI sur la r gion nantaise puis sur un pisode d t concernant la r gion lilloise 2 CONFIGURATION MATERIELLE UTILISEE Le logiciel a t install sur un PC Pentium III cadenc a 550 MHz et dot de 128 Mo de m moire vive de 10 Go de disque dur et fonctionnant sous Windows 98 La r solution de l cran initialement en 800x600 a t augment e a 1024x768 cf 4 2 4 3 PRESENTATION SUCCINCTE DE POLLUX 1 0 Le logiciel Pollux est principalement compos de trois entit s 1 Le logiciel Calmet qui calcule les param tres m t orologiques 2 L environnement
21. 0 1 93 234443 9999 9999 286 032 88 1020 00 0 98 219 3 1 80 235 55 90 1 2872900 89 1020 00 0 1 80 235 55 9999 9999 267 900 89 1020 00 0 1 80 235 905 9999 9999 287 900 89 1020 00 0 98 219 4 1 67 236 07 90 1 2814167 89 1020 00 10 1 67 236 67 9999 9999 267 167 69 1020 00 0 1 07 236 67 9999 9999 287 161 89 1020 00 0 98 219 5 1 53 297478 90 1 2872633 89 1020 00 0 1 53 237 718 9999 9999 267 633 69 1020 00 0 1 59 297 716 9999 9999 287 653 89 1080 00 0 98 219 6 1 40 238 90 90 1 287 500 90 1020 00 1 40 238 90 9999 9999 287 500 90 1020 00 1 40 238 90 9999 9999 287 500 90 1020 00 Be 219 7 1 47 242 63 90 1 289 500 82 1019 83 1 47 242 63 9999 9999 289 500 62 1019 83 1 47 242 63 9999 9999 289 500 82 1019 83 0 ooo oo Ecole des Mines de Douai D partement Energ tique Industrielle 51 ANNEXE 4 Fichier UP01 DAT D but du fichier simulation sur Lille pisode d t En gras valeurs extrapol es pour z 3000 m 98 219 0 98 224 12 700 T T 9999 2001 98 8 7 0 10 10 1018 0 20 289 5 232 3 1014 0 50 290 0 240 4 1008 0 100 290 2 250 4 991 0 250 291 9 267 4 962 0 500 292 8 283 5 935 0 750 292 1 292 5 908 0 1000 291 2 293 6 882 0 1250 290 3 289 6 856 0 1500 289 5 287 6 700 0 3000 285 0 287 7 9999 2001 98 8 7 6 10 10 1018 0 20 291 0 239 2 1014 0 50 292 1 248 3 1008 0 100 292 1 259 3 991 0 250 291 9 273 3 962 0 500 291 9 284 4 934 0 750
22. 11 17 900 to 32 700 11 Ml 7600to17 900 12 El 4000to 7600 11 El 2000to 4000 10 D 500t 2000 12 o Oto 500 12 AU Mba o de Y ANIMZERA EN D RI Le SS y IA OA TERT k RE PES be E EE LEE AO SANTE ng T lo HA SAS aem Ecole des Mines de Douai D partement Energ tique Industrielle 54 ANNEXE 7 Visualisation SIG par agr gation spatiale Emissions de NOx Heure de pointe soir du 07 08 1998 Metropole lilloise Agr gation Spatiale Parc Total Emission Totale g km NOx 07 08 98 16h00 19h00 M6 500 to 27 100 14 E1 100to 2000 12 El 800to 1100 5 O s00to 800 7 O Oto 500 14 LR amavon menmam eo gorros omar E av wan ma 2 ELA WS UN E ss 7 Pa ADA ie a ay o PA NU rch 2 EK TaN TY IRA ICL RATE aceon CR RE A EE AT H ES S F En ps 5 y ES E S naa RS AND E Re a Ecole des Mines de Douai D partement Energ tique Industrielle 55 ANNEXE 8 Visualisation SIG par agr gation spatiale Emissions de COV Heure de pointe soir du 07 08 1998 Metropole lilloise ner ur MA usna aranan uana eu end RER Agr gation Spatiale Parc Total i Emission Totale g km COW 07 08 98 18h00 19h00 M 36 200 to 70 800 11 Wi 24 700 to 36 200 13 E 12 200 to 24 700 13 fs a DE At OS BTR SETA US RA Pun ALL Pt is RATAS
23. ALD2Z SIRO ED 0 m0 A UBAN 55 0 BORG Oo 0 1 L EXOZ AO SOK GOO 0 I Y E203 2X05 DADA 0 O L CRO S 5 0 7 AO On Oye MGLY 5 0 5 0 0 0 LPAR Ss AO a0 5 A O r ETH 5207 BOK gt O6 0 l L OLE t5407 0 00 0 1 1 TOD r 5805 AO a 90 0 1 XYL omy toe re 0 OPEN 5 0 S 0 0 lt Q 2 CREST DOF OD Org O 1 OZ OOF AO ge Oo 0 1 r ROR 5x0 30 Oo 0 1 LOH202 DEL SrL apa 0 7 O 1 L TSOP s 50 gt ARMA o OA O t MEOH 5 0 5 0 0 O L ETOR 550 DEDO Ory oO 1 Y 502 GAL SFE Or O5 2 Ve SOS SEL CENT OCs 0 r UNR 507 300 0 0 1 t TSOPRD OS DAD SOKO GS 405 20 NTR 5 0 5 0 desa OK TI USEOZ SALE CA e 0 LI END INPUT GROUP 8 Vertical and Horizontal Diffusivity Parameters HORIZONTAL DIFFUSIVITY Method flag KHMODE Default 3 KHMODE 3 KHMODE 0 PBL diffusivities are determined based on stability class see DKHSTB array KHMODE 1 Same as above except diffusivities Ecole des Mines de Douai D partement Energ tique Industrielle 49 are adjusted for wind speed 2 Uses Smagorinsky method 3 Adds the results of methods 1 and 2 KHMOD KHMOD EH m Horizontal Diffusivity m 2 s above the DIFFBREAK height DKHUP Default 0 0 DKHUP USED ONLY IF KHMODE 0 1 or 3 Il o Fl ct Horizontal Diffusivity m 2 s below the DIFFBREAK height
24. COD 1 O1 amp wh RADIUS OF INFLUENCE PARAMETERS Use varying radius of influence Default F LVARY T if no stations are found within RMAX1 RMAX2 or RMAX3 then the closest station will be used aximum radius of influence over land in the surface layer RMAX1 No default RMAX1 100 Units km aximum radius of influence over land aloft RMAX2 No default RMAX2 100 Units km aximum radius of influence over water RMAX3 No default RMAX3 100 Units km OTHER WIND FIELD INPUT PARAMETERS inimum radius of influence used in the wind field interpolation RMIN Default 0 1 RMIN 2 Units km Radius of influence of terrain features TERRAD No default TERRAD 10 Units km Relative weighting of the first Ecole des Mines de Douai D partement Energ tique Industrielle guess field and observations SURFACE layer R1 R1 is the distance from an observational station at which the observation and first guess field are equally weighted Relative 37 weighting of the first guess field and observations layers ALOFT R2 R2 is applied in the upper layers Units km in the same manner as Rl is used in the surface layer in the No default R1 20 Units km in the No default R2 20 Relative weighting parameter of the prognostic wind field data RPROG No default RPROG 90 Used only
25. E 02 E 02 E 02 E 02 54 5558 464 474 1 1 1 484 494 504 514 524 534 544 554 Direction X km Champ de Temperature K le 07 08 98 a 9h niveau 10 m Ecole des Mines de Douai D partement Energ tique Industrielle 61 ANNEXE 14 Visualisation ACRPlot du champ de vitesse M tropole lilloise Direction Y km sa N 8 ie jx o L g MM 8 3 5 E par we a pe J 3 gt a Ze BL 1 MMS E m gt 8 1 gr LA e a m O ao sl J wo oa SH 5 ee Ped a we a ao BF 4 wo ao ar wo ed 7 ae we we ao BF wo ao BL 4 re wo sa a we al a a l i i l l 1454 464 474 484 494 504 514 524 534 544 554 Direction X km Vecteurs vitesse m s 07 08 98 a 9h niveau 10 m Ecole des Mines de Douai D partement Energ tique Industrielle 62 ANNEXE 15 Concentrations en O3 le 11 08 14 15 h en microgrammes m niveau 10 m Metropole lilloise Direction Y km 96 429 96 365 96 302 96 238 98 175 96 111 96 048 95 984 95 92 M 95 857 95 793 95 73 95 666 95 603 95 539 Mos 476 95 412 55 349 95 285 95 221 95 158 5658 5618 5628 5638 5648 T 5608 5598 T 5588 5578 T 5568 jo Ng gg S E E E E E 54 464 474 484 494 504 514 524 534 544 554 5558 Direction X km Concentration O3 Le 11 08 14h 15h Ecole des Mines de Doua
26. ECOLE DES MINES DE DOUAI DEPARTEMENT ENERGETIQUE INDUSTRIELLE ETUDE N 1 Partie n 2 EVALUATION DU LOGICIEL POLLUX SAMAA Jean Philippe VERMEULEN Decembre 2000 Table des Mati res Avertissement 3 R sun 4 1 Introduction 5 2 Configuration mat rielle utilis e 5 3 Pr sentation succincte de Pollux 1 0 5 3 1 Interface Calmet 5 3 2 Interface Pollux 5 3 3 Interface Calgrid 5 4 Evaluation de Pollux 1 0 6 4 1 Installation de Pollux 1 0 6 4 2 Etude existante Nantes 6 4 2 1 Saisie des donn es pour le calcul des conditions m t orologiques 6 4 2 2 Ex cution de Calmet 9 4 2 3 Saisie des donn es g n rales param tres g n raux 10 4 2 4 Saisie des param tres de calcul des missions 10 4 2 5 Calcul des missions 13 4 2 6 Ex cution du pr processeur de sp ciation Pr cal 15 4 2 7 Saisie des param tres de calcul de la chimie et de la dispersion Interface Calgrid 15 4 2 8 Ex cution de la chimie dispersion Calgrid 17 4 2 9 Visualisations des diff rentes sorties m t o missions dispersion __ 18 4 3 Nouvelle Etude Lille 19 4 3 1 Calcul de la m t o 20 4 3 2 Donn es g n rales 21 4 3 3 Calcul des missions 21 4 3 4 Ex cution du pr processeur de sp ciation 25 4 3 5 Ex cution de la chimie dispersion 25 4 3 6 Visualisation des r sultats m t o et chimie dispersion 26 4 4 Comparaison de r sultats 26 5 Recommandations et points forts 27 6 Conclusions 28 Annexes IIA 0
27. ExEtude ExSimu UPO6 DAT input oo UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu UPO7 DAT input Al UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu UPO8 DAT input or at UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Ecole des Mines de Douai D partement Energ tique Industrielle E E E 31 UPO9 DAT input ya l UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP09 DAT IEND UP10 DAT input LO UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP10 DAT IEND UP11 DAT input TA E UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP11 DAT END UP12 DAT input IZ 1 UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP12 DAT IEND UP13 DAT input 13 4 UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP13 DAT IEND UP14 DAT input 14 UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP14 DAT IEND UP15 DAT input Vo 1 UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP15 DAT IEND UP16 DAT input Lo 1 UPDAT C Pollux Sim
28. ION ICLDATT DAT FILE Default 1 ITICON 1 for a formatted text file containing one concentration value per layer per species for an unformatted file containing a full 3 D set of concentrations for each species NX NY NZ NSA values ITCLOD 0 U EFAULT 0 ITCLOD 1 NO CLOUD INFORMATION ID USED CLEAR SKY USER SHOULD MAKE SURE CALGRID WAS RUN WITH CLEAR SKY ITCLOD 1 A CLOUD ATTENUATION FILE IS TO BE OPENED AND READ Side Boundary Concentration BCON DAT file ITBCON 1 ITBCON 2 ITBCON 3 Top Boundary Default 1 ITBCON 3 for a formatted text file containing boundary types and time independent boundary conditions as a function of height for each advected species for an unformatted file containing a full set of boundary conditions for each advected species if BCON file is not used Time independent boundary conditions are taken from the initial concentration file Concentration TCON DAT file ITTCON 1 ITTCON 2 ITTCON 3 Default 1 ITTCON 3 for a formatted text file containing time independent top boundary conditions for each advected species for an unformatted file containing an full set of time and space dependent top boundary conditions for each advected species if TCON file is not used The top layer of the Ecole des Mines de Douai D partement Energ tique Ind
29. IPR3 0 n0 1 yes Default 0 Print the FINAL DIVERGENCE fields IPR4 0 n0 l yes Default 0 l Print the winds after kinematic effects are added IPR5 0 n0 1 yes Default 0 Print the winds after the Froude Number adjustment is made IPR6 O no l yes Default 0 Ecole des Mines de Douai D partement Energ tique Industrielle LD NN1 NN IO NZ B 2 UTD PRN2 PRO PR1 PR2 PR3 PR4 PR5 PR6 all variables 0 velocity scale 0 Print the w are added IPR7 0O no inds l yes after 35 slope FLOWS Default 0 Print the FINAL wind field components IPR8 0 n0 END l yes Default 0 IPR7 IPR8 H PUT GROUP 4 Meteorological data options NUMBER OF SURFACE amp PRECIP METEOROLOGICAL Number of surfa Number of precipita CLOUD DATA OPTIONS Griddid cloud f 0 Gr 1 Gr 2 Gr ICLOUD ICLOUD ICLOUD file FORMATS Surface EL 2 unformatte formatted Precipitation data file format 2 unformatte formatted Cloud data file 1 2 formatted END ce stations tion stations ields idded c d e g louds not used idded CLOUD DAT idded CLOUD DAT NSSTA No de NPSTA No de ICLOUD generated as read as INPUT meteorological data file format Default IFORMS SMERGE ou
30. LGRID grid type IVGTYP Default 2 IVGTYP 2 IVGTYP 0 uniform thickness above and below DIFFBREAK IVGTYP 1 dynamically varying layers IVGTYP 2 for arbitrary fixed grid Enter values for the following variables based on IVGTYP IVGTYP REQUIRED VARIABLES 0 NZ NZL ZMINB ZMAXB ZMINA ZMAXA ZTOP 1 NZ NZL DZMIN ZTOP 2 NZ ZFACE nz 1 IVGTYP 0 172 NZ No of vertical layers o default NZ 5 0 1 NZL No layers below DIFFBREAK o default NZL 0 ZMINB Minimum depth m of cells Default 20 ZMINB below DIFFBREAK amp above layer 1 0 ZMAXB Maximum depth m of cells Default 5000 ZMAXB below DIFFBREAK amp above layer 1 0 ZMINA Minimum depth m of cells Default 20 ZMINA above DIFFBREAK 0 ZMAXA Maximum depth m of cells Default 5000 ZMAXA above DIFFBREAK 0 1 ZTOP Top of domain m o default ZTOP 1 DZMIN Minimum cell depth m for Default 20 DZMIN x layers above layer 1 2 ZFACE nz 1 Cell face heights m o defaults in arbitrary vertical grid Ecole des Mines de Douai D partement Energ tique Industrielle 46 ZFACE 0 0 20 0
31. N Default 2400 CONSTN 2400 Overwater mixing ht equation CONSTW Default 0 16 CONSTW 0 Absolute value of Coriolis parameter FCORIOL Default 1 E 4 FCORIOL 0 Units 1 s SPATIAL AVERAGING OF MIXING HEIGHTS Conduct spatial averaging IAVEZI O no l yes Default 1 TAVEZI 1 ax search radius in averaging process MNMDAV Default 1 MNMDAV 5 Units Grid cells Half angle of upwind looking cone for averaging HAFANG Default 30 HAFANG 30 Units deg Layer of winds used in upwind averaging ILEVZI Default 1 ILEVZI 2 must be between 1 and NZ OTHER MIXING HEIGHT VARIABLES inimum potential temperature laps rate in the stable layer above th current convective mixing ht Default 0 001 DPTMIN 0 DPTMIN Units deg K m Depth of layer above current conv mixing height through which lapse Default 200 DZZI 200 rate is computed DZZI Units meters inimum overland mixing height Default 50 ZIMIN 50 ZIMIN Units meters aximum overland mixing height Default 3000 ZIMAX 3000 ZIMAX Units meters inimum overwater mixing height Default 50 ZIMINW 50 ZIMINW Not used if observed Units meters overwater mixing hts are used aximum overwater mixing height Default 3000 ZIMAXW 3000 ZIMAXW Not used if observed Units meters overwater mixing hts are used temperature Parameters Interpolation type 1 1 R 2 1 R 2 Default 1 TRAD 1 Radius
32. OV_bd prn relatif l activit des sources d missions industrielles en mettant jour la table de correspondance fichier Correspondance_COV txt CATS SN A SARRIA l LP x J AS ORS DOK JN GS F K i KK RP See ee TEE ae sal lees as SE Figure 17 uv K L aG 4 3 3 4 Chauffage r sidentiel N ayant pas notre disposition toutes les donn es n cessaires et au bon format nous n avons pas calcul les missions dues au chauffage r sidentiel Travaillant sur un pisode d t ces sources d missions n ont pas t consid r es comme tant significatives 4 3 4 Ex cution du pr processeur de sp ciation L ex cution de Pr cal a ensuite t effectu e sans probl mes majeurs a l exception d un bug portant sur la revalidation n cessaire des profils de sp ciation du module calculant les missions dues au trafic routier 4 3 5 Ex cution de la chimie dispersion Calgrid a ensuite t ex cut en tenant compte des sp cificit s de notre simulation nombre de jours maillage fichier de concentrations initiales Icon dat En ce qui concerne notre simulation le temps de calcul f t approximativement de 2 heures Ecole des Mines de Douai D partement Energ tique Industrielle 26 4 3 6 Visualisation des r sultats m t o et chimie dispersion Les annexes 13 et 14 pr sentent respectivement les champs de temp ratures et de vitesses pour le 07 08 98 a 9h On remarque bien
33. Onglet Donn es de sp ciation Les profils de sp ciation sont les m mes que ceux utilis s pour la simulation concernant la r gion nantaise Le calcul des missions routi res est ensuite ex cut Dans un premier temps Pollux proc de la cr ation d un fichier d agr gation n cessaire a la visualisation par agr gation SIG des missions Ensuite Pollux calcule les emissions pour les diff rents jours de la simulation et compacte ensuite au format gzip les fichiers r sultats Pour notre simulation le calcul des missions routi res dur une dizaine de minutes Les annexes 5 10 pr sentent les visualisations SIG par agr gation spatiale heures de pointe du matin et du soir des missions NOx et COV ainsi que les missions l mentaires pour les NOx heures de pointe du matin et du soir A Initialement Pex cutable calculant les missions dues au trafic routier routier exe ne pouvait pas g rer plus de 2000 brins ce qui provoquait une erreur d ex cution de ce module En outre en regardant plus pr cis ment le fichier Davis txt nous nous sommes apercu que certains brins de reseau comportait des capacit s et des d bits nuls Nous avons ensuite recu une version de ce module g rant les 5665 brins de notre r seau routier ainsi que des capacit s et d bits nuls De plus lors des premi res ex cutions de ce module le fichier d agr gation n a pas t cr e et d autres erreurs fatales entrainaient un planta
34. Pollux compos d interfaces permettant la saisie de certaines donn es ainsi que la pr sentation des r sultats sous forme d images interface Calmet pour les donn es m t o interface Pollux pour les donn es concernant les missions de polluants interface Calgrid pour les param tres n cessaires au calcul de la chimie dispersion interface de visualisations des r sultats 3 Le logiciel Calgrid qui calcule les param tres de la chimie dispersion des polluants 3 1 Interface Calmet A partir de fichiers de donn es m t orologiques mesur es Calmet calcule par interpolation les champs de vitesses et les variables microm t orologiques dans un domaine de calcul spatial et temporel d termin par l utilisateur 3 2 Interface Pollux Dans la version 1 0 de Pollux seuls deux types de sources d missions de polluants sont actifs Routier et Chauffage R sidentiel les autres Production Services Ferroviaire A rien et Milieu Naturel n tant pas op rationnels dans la version 1 0 du simulateur mais pr vus dans les versions ult rieures 3 3 Interface Calgrid Connaissant les champs de vitesses et les sources d missions Calgrid calcule la chimie r active des diff rents polluants ainsi que leurs concentrations en tout point du domaine de calcul temporel et spatial Ecole des Mines de Douai D partement Energ tique Industrielle 4 EVALUATION DE POLLUX 1 0 4 1 Installation de Pollux 1 0 La phase
35. R TOL RYL FORM ALD2 ETH MEOH ETOH ISOP 1 2 3 4 5 6 7 8 9 Figure 11 Fen tre DOS du pr processeur de sp ciation L ex cution de pr cal g n re un fichier arem dat dans Resultats Calgrid qui sera ensuite utilis par Calgrid Accessoirement le fichier pre_cal err dans Donnees contient des indications sur d ventuels probl mes ou erreurs d ex cution Le temps de calcul sur le cas de Nantes est de 5 mn approximativement pour cette phase de calcul 4 2 7 Saisie des param tres de calcul de la chimie et de la dispersion Interface Calgrid Ecole des Mines de Douai D partement Energ tique Industrielle 16 L interface Calgrid est lanc e a partir du menu Saisie Des Donn es Donn es Dispersion Param tres de l environnement Pollux Cette interface est similaire a celle de Calmet Les parametres relatifs a la chimie dispersion sont entr s au travers de diff rentes rubriques voir figure 12 Les param tres saisis au travers de cette interface sont lisibles dans le fichier texte calgrd inp voir Annexe 2 dans le sous r pertoire Donnees utile dans un souci de v rification des donn es entr es sous cette interface Dossiers et Fichiers Fichiers d entr e Fichiers propres a la m t o CALMET DAT CLOUD DAT de concentrations initiales ICON DAT et fichiers relatifs aux missions emissions ponctuelles cycliques ou constantes missions variables arbitrairement sources ponctuelles
36. aisable de v rifier si cet ordre est bien respect Nous sommes tout a fait conscients que cette gestion des erreurs est tr s d licate au point de vue de la programmation en effet il faut penser la majeure partie des erreurs que l utilisateur peut faire et qu elle engendrera un alourdissement important des ex cutables Par contre nous jugeons que ce d veloppement est primordial si l on d sire que le simulateur soit employ par l utilisateur en effet on h site utiliser un logiciel qui plante de mani re intempestive et dont les messages ne sont pas explicites Dans le m me registre il serait souhaitable dans la mesure du possible qu une aide en ligne au niveau du simulateur soit disponible au moment ou Pon aborde les phases les plus d licates d une simulation bien entendu cette aide devenant inutile pour un utilisateur confirm Finalement il serait galement souhaitable que la cr ation des fichiers d intersections ilots maillage brins maillage se fasse d une mani re automatique et non pas manuelle cette fonctionnalit tant pr vue dans les versions ult rieures de SAMAA Ecole des Mines de Douai D partement Energ tique Industrielle 28 En ce qui concerne les points forts du logiciel nous retenons les points suivants class s par ordre d importance La prise en compte de tr s nombreux param tres physiques repr sentatifs de la r alit pente de certains brins du r seau ro
37. as above except layer 1 data at upper air stations are ignored Extrapolate surface winds even if calm ICALM Default 0 ICALM 0 0 NO 1 YES Layer dependent biases modifying the weights of surface and upper air stations BIAS NZ 1 lt BIAS lt 1 egative BIAS reduces the weight of upper air stations e g BIAS 0 1 reduces the weight of upper air stations by 10 BIAS 1 reduces their weight by 100 Positive BIAS reduces the weight of surface stations e g BIAS 0 2 reduces the weight of surface stations by 20 BIAS 1 reduces their weight by 100 Zero BIAS leaves weights unchanged 1 R 2 interpolation Default NZ 0 BIAS 0 0 0 0 0 0 0 0 0 0 Minimum distance from nearest upper air station to surface station for which extrapolation of surface winds at surface station will be allowed RMIN2 Set to 1 for IEXTRP 4 or other situations where all surface stations should be extrapolated Default 4 RMIN2 1 Use gridded prognostic wind field model output fields as input to the diagnostic wind field model IPROG Default 0 IPROG 0 0 No IWFCOD 0 or 1 1 Yes use CSUMM prog winds as Step 1 field IWFCOD 0 Yes use CSUMM prog winds as initial guess field IWFCOD 1 Yes usr MM4 prog winds as Step 1 field IWFCOD 0 Yes use MM4 prog winds as initial guess field IWFCOD 1 Yes use MM4 prog winds as observations IWF
38. avons r alis nous avons obtenu les donn es horaires voir Annexe 3 UPnn DAT nn tant le nombre de stations m t orologiques d altitude Ce fichier contient en plus de diff rents ent tes voir pages 4 143 a 4 146 du guide de l utilisateur Calmet des donn es pour diff rentes altitudes voir Annexe 4 Initialement nous de disposions de donn es que pour une altitude maximale de 1500 m Lors de l ex cution de Calmet une erreur est survenue pour la 18 i me heure de simulation L ouverture du fichier calmet lst nous a permis de localiser facilement la provenance de cette erreur Lors du calcul de certains param tres Calmet a besoin de donn es mesur es pour une altitude sup rieure a HTOLD DZZI voir guide utilisateur Calmet Si cette quantit est sup rieure a l altitude maximale des donn es contenues dans le fichier UPnn DAT Calmet g n re une erreur d ex cution Pour contourner ce manque de donn es nous avons ajout un point pour l altitude 3000 m en extrapolant les donn es de celles a 1500 m Voir Annexe 4 loi log pour la vitesse et extrapolation lin aire pour les autres quantit s SEAmm DAT mm tant le nombre de stations m t orologiques maritimes Dans notre cas nous n avons pas de stations maritimes Initialement l absence de donn es relatives a des stations maritimes provoquait une erreur d ex cution de l interface Calmet Ce bug a ensuite t corrig par la soci t ACRI L ex cution
39. c routier et aux activit s de production et services Entre la premi re et la seconde phase de l valuation de Pollux nous avons recu une nouvelle version du logiciel permettant de traiter les missions dues aux activit s de production SAMAA version 1 1 4 Certains bugs de la premi re version y ont galement t corrig s Pour cette simulation nous avons utilis un maillage assez grossier domaine de calcul de 100 km par 100 km avec des mailles de 4 km par 4 km suivant x et y Ce choix de maillage a t guid par la volont d estimer plus la repr sentativit physique du logiciel ad quation entre les ordres de grandeur des r sultats de la simulation et certaines valeurs mesur es sur site que la repr sentativit num rique convergence en maillage En outre l utilisation d un maillage fin aurait augment consid rablement les temps de calculs Suivant l altitude z nous avons employ 5 mailles dont les faces sont situ es en 0 20 100 300 600 et 1000 m Le r seau routier consid r comporte 5665 brins et s tend sur une surface de 36 km sur 36 km voir figure 14 mulet LES bre ete Loir Toe are all i CAEN T rir eet Cure Er RL ENT mn CTI rae A dl vdd PRET FU a Raman SN AAA EC Petersen foe TEST Figure 14 R seau routier et maillage surfacique consid r s Ecole des Mines de Douai D partement Energ tique Industrielle 20 Dans cette section nous allons reprendre d
40. d installation du logiciel s est d roul e normalement La s quence de d sinstallation a galement t essay e avec succ s Comparativement l annexe A du manuel d installation du logiciel Pollux Doc N QDA MA ACR ST 001 du 31 05 2000 les diff rences suivantes ont t observ es dans la rubrique Fichiers ex cutables install s dans le r pertoire d installation Taille annonc e Octects Taille r elle Ko REGSVR32 EXE 3469312 RegTypLib exe 23552 En outre le fichier MDAC_TYP EXE est introuvable dans la version 1 0 4 2 Etude existante Nantes Le logiciel a tout d abord t test en reprenant la simulation existante effectu e sur la r gion nantaise Dans les paragraphes suivants nous allons exposer la d marche que nous avons suivie afin d effectuer la simulation en mentionnant les bugs que nous avons d tect s ainsi que quelques informations utiles Dans cette section nous ne pr senterons pas la d marche de mise en forme des fichiers contenant les donn es n cessaires a Pollux L utilisation de pollux est s quentielle c est a dire que l ordre d ex cution des diff rents programmes interfaces et logiciels est impos a savoir 1 Saisie des donn es pour le calcul des conditions m t orologiques Interface Calmet 2 Ex cution de la m t orologie logiciel Calmet 2Bis Visualisations des sorties m t orologiques facultatif a ce stade 3 Saisie des donn es g n rales environnement P
41. e d ach vement de construction Onglet Besoins en Chauffage Dates de d but et de fin de p riode de chauffe Besoins en chauffage kWh m et temp rature int rieure des logements Onglet Facteurs d missions Pour diff rents combustibles et pour les deux types de logements individuels et collectifs valeurs des facteurs d missions des polluants SO2 NOx CO COVNM CO2 et PS en g kWh Onglet Profil de Sp ciation Profils de sp ciation pour les COVNM par type de combustible et profil de sp ciation des NOx Onglet Eau chaude sanitaire D finition des besoins en ECS par type de logement kWh logement an Onglet Fichiers Emplacements et noms de diff rents fichiers de donn es d entr e n cessaires a l ex cution du pr processeur de sp ciation Pr cal X Dans l onglet Profil de sp ciation l activation du bouton r pertoire pour l eau chaude sanitaire entraine l ouverture d un fichier vide sous Wordpad Pour les autres combustibles le fichier prn est correctement ouvert k Pantie de ohaullage Po501_co_bd pr Foro _cov pd pm Poll _ cor bd pm Fra ow ed pri Se oo Dhiaga individual af colksctd oc races bd Figure 8 Interface Parametres de chauffage Onglet Profils de Sp ciation 4 2 5 Calcul des missions Les missions dues au trafic routier et au chauffage collectif sont ensuite calcul es par les deux ex cutables correspondants respective
42. e des Mines de Douai D partement Energ tique Industrielle
43. e mani re succincte les diff rentes tapes n cessaires a la simulation voir chapitre 4 2 en pr sentant les problemes auxquels nous avons fait face ainsi que les solutions correspondantes 4 3 1 Calcul de la m t o Le calcul de la m t o par le logiciel Calmet n cessite en outre quatre fichiers avec des formats sp cifiques qui il imp ratif de respecter Calmet ainsi que Calgrid sont des programmes crits en fortran GEO DAT Ce fichier comporte en plus des donn es relatives au maillage voir p 4 135 du guide de l utilisateur Calmet deux tableaux de donn es qui doivent obligatoirement tre entr es sous forme de tableaux de dimensions Nx Ny Nx et Ny tant respectivement le nombre de cellules du maillage suivant les directions X et Y SURF DAT Ce fichier au format libre contient les observations horaires en surface pour le nombre de stations prises en compte Ces donn es sont la vitesse du vent m s sa direction deg la hauteur de plafond centaine de pieds l opacit de la couverture nuageuse en dixi me la temp rature K l humidit relative la pression statique mb et un code de pr cipitations Dans notre cas nous avons affect les valeurs de la seule station a notre disposition a deux autres stations fictives donnant ainsi 3 stations de surface A partir d un fichier contenant les observations mesur es toutes les six heures et en utilisant un programme d interpolation lin aire que nous
44. ear IBYR No default IBYR 98 Month IBMO No default TBMO 8 Day IBDY No default IBDY 7 Hour IBHR No default IBHR 1 Length of run hours IRLG No default TRLG 120 Number of time steps per hour NSUBTS Default 3 NSUBTS 5 Time step sec 3600 NSUBTS Total number of chemical species advected steady state species NSPEC Default 50 NSPEC 41 umber of chemical species to be advected NSA Default 40 NSA 36 umber of chemical species to be deposited NSDD Default 15 NSDD 10 umber of chemical species to be emitted NSE Default 13 NSE 17 Ecole des Mines de Douai D partement Energ tique Industrielle 43 TYPE OF CHEMICAL MECHANISM IN USE ICHEM SAPRC 90 2 CB4 m COLUMNAR OZONE IN TO BE USED DURING DEFAULT 1 ICHEM 2 DOBSON UNITS RUN DEFAULT 295 TOZDU 320 Method flag for integration of chemical mechanism METHINT Default 2 METHINT 1 1 Hybrid method 2 QSSA method INPUT FILE TYPES METEOROLOGICAL VARIABLES CALGRID DAT FILE U EFAULT 1 ITIMET 1 a STANDARD HOUR ENDING CALGRID DAT FILE 2 NEW HOUR SPANNING METEOROLOGICAL DATA FILE Initial Concentration ICON DAT file ITICON Il m ITICON Il N SOLAR ATTENUAT
45. elles en NOx dues au trafic automobile en 1996 source Rapport d Activit n 3 EMD 1997 Ecole des Mines de Douai D partement Energ tique Industrielle 67 ANNEXE 20 Emissions annuelles en COV dues au trafic automobile en 1996 source Rapport d Activit n 3 EMD 1997 lt 5 10 25 50 75 100 gt 150 en an kra Ecole des Mines de Douai D partement Energ tique Industrielle 68 ANNEXE 21 Emissions en NOx dues au trafic automobile Cumul journalier du 07 08 1998 Parc Total Emissions Totales de NOx CANE 1 EE TY E NE on NEN RE Agr gation Mixte Transport Routier g km YL Essence NOx 07 08 98 1h 00 pas 24 M318 000 to 424000 1 m 212 000to318 000 2 106 000 to 212000 9 O Oto 106 000 58 PRA Te i de Al Y PAP CIN EEN NY LE JA TA Al EPP REN Ea EA AAN FA FAR E A AN E des LAINAN j At oy qu Y a ga PO ve RE ES jy cA oar le E Ecole des Mines de Douai D partement Energ tique Industrielle 69 ANNEXE 22 Emissions en COV dues au trafic automobile Cumul journalier du 07 08 1998 Parc Total Emissions Totales de COV Bar Agr gation Mixte Transport Routier g km YL Essence COY 07 08 98 1h 00 pas 23 Ml 572 000t0 716000 1 E 429 000to 572 000 7 E 286 000 to 429 000 12 E 143 000 to 286 000 23 Oto 143 000 46 IL Chie HS AT marae fH AER cies m Danit T Ecol
46. enue par SAMAA peut s expliquer par le fait que les emissions de COV sont tr s d pendantes de l indice de d marrage froid de l vaporation et du zonage conomique que nous avons estim d une mani re assez grossi re voir paragraphe 4 3 3 1 relatif au zonage conomique utilis En outre les missions de COV d pendent galement tr s fortement de la temp rature La division par 365 des donn es de 1996 n est ainsi pas forc ment tr s repr sentative des mois d t et elle int gre les faibles missions lors des p riodes hivernales 5 RECOMMANDATIONS ET POINTS FORTS Dans cette partie nous formulons certaines recommandations relatives l utilisation du logiciel SAMAA Nous pr sentons galement une liste de points qualifi s de points forts qui ont retenu notre attention Dans un premier temps il est tout a fait indispensable que le simulateur soit dot d une proc dure de gestion des erreurs afin de limiter de mani re importante le nombre de plantages intempestifs En outre les messages d erreurs sont insuffisamment explicites exemple erreur 7 Indice en dehors de la plage et ne permettent pas de rep rer facilement l origine des erreurs d ex cution Il est tout a fait envisageable par exemple de verifier si le format d un fichier est bien celui qui est attendu par le simulateur par exemple pour un fichier qui doit contenir des donn es class es suivant un certain ordre il est tout a fait f
47. er dans le fascicule Document de sp cification description des fichiers d entr e section Visualisation des missions du trafic routier Le fichier d intersection Davis Maillage Ce fichier est similaire au pr c dent sauf que les donn es sont class es par ordre croissant du num ro de brin pour le fichier pr c dent les donn es sont class es par ordre de nom de maille voir le fascicule fascicule Document de sp cification description des fichiers d entr e section Donn es d mission Trafic routier Fichier de positionnement des brins sur le maillage Calmet Les types de jours sont ceux pr d finis aucun changement de type de jours n a t op r La liste des polluants est celle initialement valid e dans les donn es g n rales Onglet Parc Automobile Nous avons utilis les donn es par d faut relatives aux carburants et au de Super Par contre nous avons proc d a un ajustement manuel concernant les valeurs relatives la r partition du parc automobile tout en v rifiant que la somme des de v hicules l gers essence et diesel est gale 100 et que la somme des de v hicules poids lourds et bus est galement de 100 Onglet Donn es du trafic Les champs relatifs aux variations temporelles poids lourds et Activit brin n ont pas t modifi s Le nombre total de v hicules pr sents sur le r seau routier a t Ecole des Mines de Douai D partement Energ tique
48. eur grossi re Une utilisation assez simple des logiciels Calmet et Calgrid la seule contrainte faible tant le format impos de certains fichiers d entr es avec des fichiers permettant de savoir si toutes les donn es ont t bien interpr t es fichier calmet Ist et galgrid st Une utilisation tr s simple de la visualisation SIG avec un maniement ais des couches Mapinfo La possibilit d utiliser d autres profils de sp ciation nous avons tr s facilement greff notre profil de sp ciation COV pour les missions industrielles Citons galement un point qui est actuellement en phase de d veloppement dans le Soci t ACRi et qui concerne l impl mentation de routines de cr ation de fichiers d entr e pour Calmet a partir des donn es ARPEGE et MEL Un dernier point non n gligeable et que l on peut inclure dans la rubrique des aspects positifs du logiciel est la grande disponibilit de l quipe de d veloppement de SAMAA de la Soci t ACRI 6 CONCLUSIONS Ecole des Mines de Douai D partement Energ tique Industrielle 29 Le logiciel SAMAA dans sa forme actuelle est un produit tres prometteur mais qu il est absolument n cessaire de stabiliser Le point le plus important qui nous a s duit est la prise en compte d un nombre important de param tres physiques dans le but de fournir des simulations tr s d taill es Bien entendu on accro t ainsi d une mani re inh rente le deg
49. fluxes IFPRT Default 0 IFPRT 0 Print top fluxes ITPRT Default 0 ITPRT 0 Print deposition vel IVDPRT Default 0 IVDPRT 0 0 Do not print 1 Print No layers of gridded total areatpoint missions to print IEPRT Default 0 TEPRT 0 IEPRT must be lt NZ Concentration print interval ICFRQ in hours Default 1 ICFRQ 1 Ecole des Mines de Douai D partement Energ tique Industrielle 48 Dry flux print interval IFFRQ in hours Default 1 TFFRO 1 Top flux print interval ITFRQ in hours Default 1 TTFRO 1 Deposition vel print interval IVDFRQ in hours Default 1 IVDFRQ 1 Emissions print interval IEFRQ in hours Default 1 TEFRO 1 SPECIES LIST FOR OUTPUT OPTIONS SITAS CONCENTRATIONS DRY PLUMES gt SPECIES PRINTE SAVED ON DISK PRINTED SAVED ON DISK NAME 1 VALUE LAYER 1 VALUE LAYER O NO 1 YES O NO 1 YES 0 NO 1 YES O NO 1 YES ANO Sat SET Dz wr 1 NO2 gt DAL DAL hy Oe 0 1 VOSS bed bed OO r HONO DF Oy gt OD 0 4 O 2 1 HNO3 IO 5 0 0 41 LPNA 5 0 SKOi 05 0 r oN205 D405 5 0 o 10 5 O 1 r NOB DEL D E O O 2 HOZ orl SAT O LEO A o O7 0 od FORM 5 0 5 0 0 0 l VS
50. ge de Pollux Suite a de fructueux changes t l phoniques avec la soci t ACRI nous avons identifi les causes de ces problemes nom et emplacement des fichier de d marrage a froid ainsi que du fichier du maillage Il est noter que ces bugs ont t engendr s parts gales par une mauvaise utilisation du logiciel et par des ex cutables d fectueux qui ont t corrig s par la suite 4 3 3 2 Trafic a rien N ayant a notre disposition aucune donn e concernant le trafic a rien nombre d avions r partition du parc suivant le type d appareils et r partition temporelle du trafic nous avons dans un premier temps utilis les valeurs par d faut du logiciel Les missions dues au trafic a rien ont ensuite t ex cut es sans probl me puis revalid es une seconde fois avec un nombre d appareils insignifiant 10 appareils pas an 4 3 3 3 Production et services En ce qui concerne le calcul des missions dues aux activit s de production et services NOx et SO2 nous n avons consid r que des sources d missions ponctuelles pas de stations service et pas de chauffage tertiaire Nous avons travaill sur 38 sources industrielles suivant une approche TGAP voir figure 17 pour la localisation et annexes 11 et 12 pour les visualisations SIG des missions Nous Ecole des Mines de Douai D partement Energ tique Industrielle 25 avons utilis notre propre profil de sp ciation COV fichier Prof701_C
51. gure 5 l utilisateur d finit le d but et la fin de la simulation soit en les saisissant directement dans les deux champs pr vus soit par l interm diaire du calendrier ainsi que les polluants consid r s dans le champ Choix polluant A ll est conseill de v rifier ces dates avant de lancer les calculs missions chimie Pollux ayant tendance r actualiser les champs avec la date courante ce qui entraine une erreur fatale d ex cution lors des calculs A Fasnia ginan aus Date de depa reida panse D bui de drid ahh jam cr 1998 Bireriar 1338 han El Finde dresin arar 19 han 26 11 12 11 14 18 18 20 21 5 26 Er 26 Pi Di Soop re Bec Pa pp i a us 14 ri T Chom polluard CO il DO PM 500 Oa CHA E ue Figure 5 Saisie des donn es g n rales 4 2 4 Saisie des param tres de calcul des missions 4 2 4 1 Trafic Routier Sans entrer plus dans les d tails ce que fait tr s bien le manuel d utilisation du logiciel Pollux les donn es n cessaires au calcul des missions dues au trafic routier sont de plusieurs types voir figure 7 G n ral Routier regroupe plusieurs donn es type de polluants type de jour il est possible de faire le distinguo entre un jour ouvr et un jour f ri de Ecole des Mines de Douai D partement Energ tique Industrielle 11 simuler des d parts en vacances fichier de donn es de trafic fichier d inter
52. i D partement Energ tique Industrielle 63 ANNEXE 16 Concentrations en NO2 le 11 08 14 15 h en microgrammes m niveau 10 m Metropole lilloise 5658 Direction Y km 5578 5588 5598 5608 5618 5628 5638 5648 T 5568 loto la la to la ll E 464 474 484 494 504 514 524 534 544 554 5558 g Direction X km Concentration NO2 Le 11 08 14h 15h niveau 10m Ecole des Mines de Douai D partement Energ tique Industrielle 64 ANNEXE 17 Concentrations en NO le 11 08 14 15 h en microgrammes m niveau 10 m Metropole lilloise 5658 Direction Y km 5588 5598 5608 5618 5628 5638 5648 T 1 L 1 L N 5578 5568 loto lala to lalo ly to 464 474 484 494 504 514 524 534 544 554 5558 g Direction X km Concentration NO Le 11 08 14h 15h niveau 10m Ecole des Mines de Douai D partement Energ tique Industrielle 65 ANNEXE 18 Concentrations en xyl ne le 11 08 14 15 h en microgrammes m niveau 10 m M tropole lilloise a 14 Petes 3 831 m5 4 m 5254 L 7 8 2 756 E J 2 189 i A 2 1 042 3 E 5 1 095 0 548 CR o 453 i 3 906 E d 5 8 J co 5 la a ta toto tot po 54 464 474 484 494 504 514 524 534 544 554 Direction X km Concentration Xylene Le 11 08 14h 15h niveau 10m Ecole des Mines de Douai D partement Energ tique Industrielle 66 ANNEXE 19 Emissions annu
53. ielle S N 1NO2 0 1656 1 00 8 07 520 4 E 2 LOS 019945 10500 Lry Los A eHaZ I HONO 0 1100 1 00 4 0 2 0 4 E 2 r NOS 0616567 1 00 820 lt 5 0 y AE 2 1 1 HO2 0 2402 1 00 12 0 0 0 4 E 2 FORM 0 2336 1 00 4 0 0 0 4 E 2 PAN 0 1050 1 00 4 0 1 0 4 E 2 H202 0 2402 1 00 12 0 0 0 4 E 2 IEND INPUT GROUP 5 Size parameters for dry deposition of particles SPECIES GEOMETRIC MASS MEAN GEOMETRIC STANDARD AME DIAMETER DEVIATION microns microns IEND INPUT GROUP 6 Miscellaneous dry deposition parameters REFERENCE CUTICLE RESISTANCE RCUTR S CM RCUTR 17 REFERENCE GROUND RESISTANCE RGR S CM RGR REFERENCE POLLUTANT REACTIVITY REACTR EACTR VEGETATION STATE IN UNIRRIGATED AREAS IVEG IVEG 1 IVEG 1 FOR ACTIVE AND UNSTRESSED VEGETATION IVEG 2 FOR ACTIVE AND STRESSED VEGETATION IVEG 3 FOR INACTIVE VEGETATION IEND INPUT GROUP 7 Output Options FILE DEFAULT VALUE VALUE THIS RUN Concentrations ICON T ICON 1 Fluxes IDRY lt IDRY 0 0 Do not create file 1 create fil LINE PRINTER OUTPUT OPTIONS Print concentrations ICPRT Default 0 ICPRT 1 Print dry
54. if IPROG 1 Units km aximum acceptable divergence in th divergence minimization procedure DIVLIM Default 5 E 6 DIVLIM O aximum number of iterations in the divergence min procedure NITER Default 50 NITER 50 Procedure OTE NZ Defa NS aximum n each laye data toa OTE NZ N Critical ALPHA to upper umber of passes in the smoothing NSMTH NZ values must be entered ult 2 mxnz 1 4 MIH 3y 087 8 12 1 25 125 A SE E Ay umber of stations used in r for the interpolation of grid point NINTR2 NZ values must be entered Default 99 TNTR2 99 99 9 97 99 99 99 7 995 995 1994 99 1 Froude number CRITFN Default 1 0 CRITFN 1 Empirical factor controlling the influence of kinematic effects ultiplicative scaling factor for extrapolation of surface observations layers FEXTR2 NZ FEXTR2 0 Oer O 0 Osy Used onl barrier Info Default 0 1 ALPHA 0 Default NZ 0 0 Org Dar Dijo Ob OF y if IEXTRP 3 or 3 rmation Number of barriers to interpolation of the wind fields NBAR Default 0 NBAR 0 THE FOLLOWING 4 VARIABLES ARE INCLUDED ONLY IF NBAR gt 0 NOTE NBAR values must be entered No defaults for each variable Units km X coordinate of BEGINNING of each barrier XBBAR NBAR XBBAR 0 Y coordinate of BEGINNING of each ba
55. ipitations Stations de surface Pour chaque station saisie du n du nom codifi de son Id coordonn es suivant X et Y Fuseau horaire et hauteur an mom trique d altitude Pour chaque station saisie du n du nom codifi de son Id coordonn es suivant X et Y et du fuseau horaire de pr cipitation Pour chaque station saisie du n du nom codifi de son Id et de ses coordonn es suivant X et Y De plus amples d tails concernant ces diff rents param tres sont disponibles dans le guide d utilisation du logiciel Calmet format des donn es contenues par exemple dans les fichiers UP DAT mM Nous avons d tect un probl me au niveau de l interface Calmet dans le menu Hauteurs de m lange voir figure 3 Les valeurs des hauteurs de m lange minimales et maximales sont invers es lorsque l on relance une seconde fois l interface probl me de relecture du fichier calmet inp ou probl me d interpr tation de ces donn es par l interface Ce bug galement mentionn a ACRI peut n anmoins tre contourn en ditant le fichier calmet inp voir Annexe 1 et en corrigeant manuellement ces valeurs dans la section INPUT GROUP 6 Mixing Height Temperature and Precipitation Parameters la sous section OTHER MIXING HEIGHT VARIABLES Ce fichier est tres bien comment et a la limite il est possible de ne pas passer par l interface pour rentrer les donn es m t orologiques du moins pour un
56. l O NO 1N02 103 1502 minimum 6 60886 4 61671 0 01878 0 00349 maximum 0 00570 0 01790 0 03553 6 68547 average 6 66642 4 61269 6 82261 6 886485 HEHE COMP HOUR 2 finished reading lateral boundary conditions finished reading top boundary conditions HORIZONTAL ADVECTION ist HALF TIME STEP NHR CENTRAL OPERATOR NHR 2 HORIZONTAL ADVECTION 2nd HALF TIME STEP NHR Figure 13 Fen tre DOS de l ex cution de Calgrid 4 2 9 Visualisations des diff rentes sorties m t o missions dispersion Deux types de visualisations sont propos es suivant l origine des fichiers de donn es visualiser La premi re est la visualisation SIG qui concerne uniquement les missions la seconde ACRPlot est relative a la visualisation des donn es issues de Calmet m t o et Galgrid chimie dispersion La visualisation SIG permet de visualiser les missions l mentaires par agr gation spatiale et ou temporelle d une mani re tr s simple voir par exemple les annexes 5 11 Le second mode de visualisation est plus complexe d utilisation dans la mesure ou la proc dure de cr ation de ces images est plus lourde et quelque peu fig e Dans un premier temps il est n cessaire de d finir des parametres n cessaires a ACRplot au travers du menu Pr paration des donn es fichier d entr e i e calmet dat ou calgrid dat param tres de simulation param tres de sortie Dans un
57. lation qui sont axe Sud centre de Lille axe centre de Lille m tropole Roubaix Tourcoing direction Nord Est axe centre de Lille Nord Ouest direction autoroute Dunkerque Le centre de l agglom ration lilloise est aussi tr s clairement identifi Au point de vue quantitatif et pour le centre de Lille nous avons le tableau suivant Tableau 1 Comparatif des missions routi res de NOx et COV Donn es Donn es R sultats Donn es Donn es R sultats 1996 NOx 1996 NOx SAMAA NOx 1996 COV 1996 COV SAMAA t an km kg j km kg j km Van km kg j km COV Kg km 100 274 318 424 75 205 572 716 Les donn es 1996 exprim es en kg j km ont t obtenues en divisant les donn es exprim es en t an km2 par 365 Les donn es SAMAA concernant les NOx ont t obtenues en faisant un cumul sur les 24 premi res heures de simulation ce qui correspond au premier jour de simulation voir annexes 21 et 22 Ecole des Mines de Douai D partement Energ tique Industrielle 27 Attention sur ces figures les l gendes VL essence NOx et VL essence COV sont fausses et les valeurs repr sent es correspondent au parc total de v hicules pour le cumul journalier du 07 08 1998 Ainsi d un point de vue quantitatif on observe que l ordre de grandeur des missions est respect les NOx tant mieux estim s que les COV La sur estimation des valeurs obt
58. m US1 UAA 2001 538 37 5149 51 1 1 US2 UAB 2002 576 73 5149 87 Sne US3 UAC 2003 615 10 5150 48 1 1 US4 UAD 2004 653 47 5151 33 sles US5 UAE 2005 691 83 5152 42 1 41 US6 UAF 2006 538 01 5205 07 At US7 UAG 2007 576 03 5205 43 1 1 US8 UAH 2008 614 04 5206 04 oe US9 UAI 2009 652 05 5206 89 1 4 US10 UAJ 2010 690 07 5207 98 sl US11 UAK 2011 537 66 5260 63 a1 4 US12 UAL 2012 575 32 5261 00 sleet US13 UAM 2013 612 97 5261 60 e US14 UAN 2014 650 63 5262 45 sheet r US15 UAO 2015 688 28 5263 54 1 4 US16 UAP 2016 537 30 5316 20 ls US17 UAQ 2017 574 60 5316 57 i US18 UAR 2018 611 90 5317 17 Le r US19 UAS 2019 649 19 5318 02 A US20 UAT 2020 686 49 5319 11 1 1 1 Four character string for station name MUST START IN COLUMN 9 2 Five digit integer for station ID END INPUT GROUP 9 Precipitation station parameters PRECIPITATION STATION VARIABLES One record per station 0 records in all NOT INCLUDED IF NPSTA 0 a 2 NAME Station X Coord Y Coords Code km km 1 Four character string for station name MUST START IN COLUMN 9 2 Six digit station code composed of state code first 2 digits and station ID last 4 digits IEND Ecole des Mines de Douai D partement Energ tique Industrielle 42 ANNEXE 2 Exemple de fichier calgrd inp Simulation sur Lille pisode t
59. ment routier exe et chauffage exe Ici encore le progression des calculs en terme d tapes est affich e l cran voir figures 9 et 10 Ecole des Mines de Douai D partement Energ tique Industrielle MS ere amp routier 14 procedure lecture lecture du fichier entree_interface procedure lectDavis procedure Debit vitesse procedure calc_gamma debut du calc des emissions calcul calcul calcul calcul calcul calcul calcul calcul fin du calcul des CO des s NOx des COU des PM des e des des des calcul des emissions rfaciques Soak et diurne phase d ecriture des resultats Figure 9 Fen tre DOS du calcul des missions dues au trafic routier lecture des donnees generales lecture du fichier Meteo debut du calcul des emissions ecriture des resultats ecriture des fichiers binaires pour Calgrid Figure 10 Fen tre DOS du calcul des missions dues au chauffage r sidentiel L ex cution du calcul des missions dues au trafic routier g n re l criture de fichiers texte du type Rout _ _ _0 txt dans le r pertoire Resultats fichiers lisibles sous Wordpad ainsi que des fichiers binaires du type Rout _ _ _u txt dans le r pertoire Resultats Calgrid La d nomination de ces fichiers est explicite par exemple le fichier RoutVE MC CO O txt contient les missions de CO dues aux V hicules Essence lors de D marrages Froid Ces missions sont donn es sous forme de
60. nsi que leurs formats fr quence d enregistrement Diffusivit D finition des diffusivit s horizontale et verticale Ecole des Mines de Douai D partement Energ tique Industrielle 17 Figure 12 Interface Calgrid 4 2 8 Ex cution de la chimie dispersion Calgrid L ex cution de Calgrid a partir de l environnement Pollux entra ne l ouverture d une fen tre DOS voir figure 13 avec affichage l cran de le progression des calculs Apr s ex cution Calgrid cr e en outre le fichier texte calgrid Ist contenant divers renseignements param tres entr s sous l interface progression des calculs temps de calcul A titre d exemple le temps de calcul sur la simulation de Nantes est de 3h Calgrid doit tre lanc a partir de l environnement Pollux et non pas a partir de l interface Calgrid Ecole des Mines de Douai D partement Energ tique Industrielle 18 HORIZONTAL ADUECTION 2nd HALF TIME STEP NHR Hr index cputime step elapsed time secs HORIZONTAL ADUECTION ist HALF TIME STEP NHR CENTRAL OPERATOR NHR 1 HORIZONTAL ADUECTION 2nd HALF TIME STEP NHR Hr index cputime step elapsed time secs HORIZONTAL ADUECTION ist HALF TIME STEP NHR CENTRAL OPERATOR NHR 1 HORIZONTAL ADUECTION 2nd HALF TIME STEP NHR Hr index cputime step amp elapsed time secs 2666066066666666 Hourly Chemica
61. of influence for temperatur interpolation TRADKM Default 500 TRADKM 500 Units km Ecole des Mines de Douai D partement Energ tique Industrielle 40 Maximum Number of stations to include in temperature interpolation NUMTS Default 5 NUMTS 5 Conduct spatial averaging of temp eratures IAVET O no l yes Default 1 TAVET 1 will use mixing ht MNMDAV HAFANG so make sure they are correct Il o Default temperature gradient Default 0098 TGDEFB below the mixing height over water K m TGDEFB Il o Default temperature gradient Default 0045 TGDEFA above the mixing height over water K m TGDEFA Beginning JWAT1 And ending JWAT2 land use categories for temperature t JWAT1 55 interpolation over water Mak JWAT2 55 bigger than largest land use to disable PRECIP INTERPOLATION PARAMETERS Method of interpolation NFLAGP Default 2 NFLAGP 3 1 1 R 2 1 R 2 3 EXP R 2 Radius of Influence km SIGMAP Default 0 0 gt use half dist btwn nearest stns w amp w out precip when NFLAGP 3 Minimum Precip Rate Cutoff mm hr Default 0 01 CUTP 1 values lt CUTP 0 0 mm hr 100 0 SIGMAP 1 END INPUT GROUP 7 Surface meteorological station parameters SURFACE STATION VARIABLES One record per station 12
62. ollux 4 Saisie des parametres de calcul des missions Routier et Chauffage R sidentiel 5 Ex cution du calcul des missions Airemis 5Bis Visualisations des sorties d missions facultatif ce stade 6 Ex cution du pr processeur de sp ciation Pr cal 7 Saisie des param tres de calcul de la chimie et de la dispersion 8 Ex cution de la chimie dispersion Calgrid 9 Visualisations des diff rentes sorties m t o missions dispersion 4 2 1 Saisie des donn es pour le calcul des conditions m t orologiques Ecole des Mines de Douai D partement Energ tique Industrielle La s lection de la rubrique saisie des donn es m t o active l interface Calmet L ouverture du fichier calmet inp entra ne l affichage d un graphique pr sentant quelques param tres de la simulation cf figure 1 Figure 1 Interface Calmet Les param tres d entr e de Calmet sont ensuite saisis via diff rents menus voir figure 2 et paragraphes suivants Figure 2 Param tres d entr e de Calmet Dossiers et Fichiers Fichiers d entr e G ographie GEO DAT observation horaire en surface SURF DAT n bulosit et pr cipitations observ es donn e MM4 station WT stations lev es fichier UP DAT stations sur l eau SEA DAT donn es diagnostiques et pronostiques Fichiers de sortie Fichiers ASCII et binaires Test et d bogages Impression ou non de variables Ecole des Mines de Douai D pa
63. r de complexit du logiciel ce qui implique donc un investissement personnel tr s important au niveau de son utilisation Citons en outre quelques pr requis qui sont relatifs l utilisation de Mapinfo et dans une moindre mesure Access et des notions de bases de Fortran pour comprendre le format d entr e de certains fichiers utilis s par Calmet et Calgrid La repr sentativit physique des r sultats de SAMAA est en ce qui concerne les emissions de NOx et COV coh rente Notons pour conclure cette valuation que SAMAA est a post riori relativement ais d utilisation pour le peu que l on s en donne les moyens puisque nous avons r ussi mener une simulation assez compl te sur une agglom ration cons quente Finalement le fait de r aliser une simulation tout fait diff rente de celle effectu e sur la r gion nantaise en termes de nombre de brins du r seau routier nombre de jours a permis de localiser un certain nombres de bugs qui initialement n ont pas t d tect s mais qui ont t corrig s par l quipe de d veloppement de la Soci t ACRI donnant ainsi un logiciel plus stable Ecole des Mines de Douai D partement Energ tique Industrielle 30 ANNEXE 1 Exemple de fichier calmet inp Simulation sur Nantes 3 overwater
64. records in all 1 2 NAME ID X coord Y coord Time Anem km km zone Ht m SSI AAA 1001 615 10 5150 48 1 10 1 SS2 AAB 1002 653 47 5151 33 I TO SSS AAC 1003 691 83 5152 42 i 10 SS4 AAD 1004 576 03 5205 43 ll TO SS5 AAE 1005 614 04 5206 04 1 10 SS6 AAF 1006 652 05 5206 89 1 TO 1 SS7 AAG 1007 690 07 5207 98 10 SS8 AAH 1008 537 66 5260 63 ll TO SS9 AAI 1009 575 32 5261 00 1 10 SS10 AAJ 1010 612 97 5261 60 1 10 1 SS11 AAK 1011 650 63 5262 45 ll 10 1 SS12 AAL 1012 688 28 5263 54 i TO SS13 AAM 1013 537 30 5316 20 a 10 SS14 AAN 1014 574 60 5316 57 ml TO 1 1 SS15 AAO 1015 69053 ee 17 i 10 SS16 AAP 1016 649 19 5318 02 ll TO SS17 AAQ 1017 686 49 5319 11 a 10 SS18 AAR 1018 605 93 5225 22 I 10 L SS19 AAS 1019 608 85 5232 48 ll 10 SS20 AAT 1020 584 61 5239 53 ll 10 SS21 AAU 1021 562 87 5243 14 ell 10 Four character string for station name MUST START IN COLUMN 9 Ecole des Mines de Douai D partement Energ tique Industrielle 41 2 Five digit integer for station ID END INPUT GROUP 8 Upper air meteorological station parameters UPPER AIR STATION VARIABLES One record per station 3 records in all 1 2 AME ID X coord Y coord Time zone km k
65. rrier YBBAR NBAR YBBAR 0 X coordinate of ENDING of each barrier XEBAR NBAR XEBAR 0 Y coordinate of ENDING of each barrier YEBAR NBAR YEBAR 0 Ecole des Mines de Douai D partement Energ tique Industrielle DIAGNO Sur Dom rat Dom ID Obs for STIC MODULE DATA INPUT OPTIONS face temperature IDIOPTI O Compute internally from hourly surface observation 1 Read preprocessed values f a data file DIAG DAT Surface met station to use fo the surface temperature ISURF Must be a value from 1 to NSS Used only if IDIOPT1 0 ain averaged temperature laps e IDIOPT2 0 Compute internally from twice daily upper air obse 38 1 Read hourly preprocessed v from a data file DIAG DAT Upper air station to use for the domain scale lapse rate I Must be a value from 1 to NUS Used only if IDIOPT2 0 Depth through which the domain lapse rate is computed ZUPT Used only if IDIOPT2 0 ain averaged wind components IOPT3 0 Compute internally from twice daily upper air obse 1 Read hourly preprocessed v a data file DIAG DAT Upper air station to use for the domain scale winds IUPWND Must be a value from 1 to NU Used only if IDIOPT3 0 Bottom and top of layer throug which the domain scale winds are computed ZUPWND 1 ZUPWND 2 Used only if IDIOPT3 0 erved s
66. rtement Energ tique Industrielle Param tres G n raux Titre de simulation date de d but et dur e de la simulation type de variables calcul es champ de vents et ou variables microm t orologiques Maillage Projection UTM ou Lambert nombre de cellules discr tisation spatiale suivant X Y et Z Maillage uniforme suivant X et Y maillage uniforme ou non suivant Z avec d finition du nombre et de la position des cellules Origine du maillage point de r f rence et espacement pas d int gration spatiale suivant X et Y Options de sortie Binaires et ou ASCII avec diff rents param tres imprimer Tests et d bogages Options du modules m t orologique impression ou non des variables m t orologiques dans un fichier Ist fichier ASCII Options du module vents impression ou non des donn es de test d bogage dans un fichier TEST Donn es M t o Nombre de stations m t orologiques discr tisation de la n bulosit format de fichier Champs de vents mod les de vent barri res module diagnostique brises rayon d influence autres param tres Hauteurs de m lange Constantes empiriques moyenne spatiale autres variables associ es hauteurs de m lange minimales et maximales Temp ratures Type d interpolation utilis e pour le calcul et gradient de temp rature sur l eau par d faut Pr cipitations Param tres des pr cipitations Type d interpolation utilis e pour le calcul des pr c
67. second temps on d finit les attributs des images que l on souhaite cr er taille de l image titres sous titres tailles des caract res axes variable a imprimer mode d affichage chelles On peut galement d finir des attributs diff rents pour plusieurs images Les images au format gif sont ensuite cr es par ACRPlot Ces attributs sont enregistr s dans un fichier plt que l on peut r ouvrir ult rieurement afin de cr er une nouvelle image Par contre il est impossible par l interface graphique de changer un des attributs d une image d j cr e si l on souhaite modifier une image existante il est n cessaire de la red finir comme une nouvelle image Par contre le fichier plt tant un fichier ASCII il est possible de le modifier manuellement sous un diteur de texte Wordpad condition de comprendre les commandes qui y sont contenues Un des points p nible de ce mode de visualisation est que si un des attributs n a pas t correctement renseign Pollux se plante sans avertissement Ecole des Mines de Douai D partement Energ tique Industrielle 19 4 3 Nouvelle Etude Lille Afin de compl ter l valuation de Pollux nous avons effectu une nouvelle simulation sur la m tropole lilloise Cette simulation concerne un pisode de pollution par l ozone en t d une dur e de cing jours du 07 08 98 au 11 08 98 Pour cet pisode nous avons pris en compte les missions dues au trafi
68. section entre les brins du trafic et le maillage Calmet Calgrid Parc Automobile R partition des parcs VL essence et diesel et PL et Bus super pourcentages de v hicules roulant au Sans Plomb de v hicules dot s d un systeme d injection plut t que de carburateurs de v hicules quip s d un systeme de contr le d vaporation Carburants caract ristiques des carburants Donn es de trafic Cette partie permet de calculer les d bits de v hicules VL et PL pour chaque brin du r seau routier et pour chaque heure de simulation D marrage Froid Cette section permet de tenir compte du fait qu un moteur chaud pollue diff remment qu un moteur froid On entre ainsi la longueur moyenne des trajets ainsi que le champ de temp ratures calcul pr c demment par Caltemp Donn es urbaines Ce sont des zonages qui couvrent la zone de simulation et qui permettent de d crire plus finement le trafic urbain et ou interurbain Y Le zonage Activit permet de s lectionner une zone g ographique et d appliquer l ensemble des brins la constituant une activit facteur multiplicatif compris entre 0 et 1 et ce pour des heures pr cises de simulation L int r t est de simuler une route barr e ou une interdiction de circuler dans cette zone Malheureusement le trafic n est pas redistribu sur les autres brins lorsque le trafic est coup sur une zone Y Le zonage Economique permet de localiser les zones de
69. tableau avec en ligne l heure de simulation et en colonne le nombre de brins de trafic consid r s Ecole des Mines de Douai D partement Energ tique Industrielle 15 L ex cution du calcul des missions dues au chauffage r sidentiel g n re l criture de fichiers du type ChaufRes_ _ _0 txt et ChaufECS_ _ _0 txt dans le r pertoire Resultats fichiers lisibles sous Wordpad ainsi que des fichiers binaires du type ChaufRes_ _u txt et ChaufECS_ _ _u txt dans le r pertoire Resultats Calgrid La d nomination de ces fichiers est ici encore tr s explicite par exemple le fichier ChaufRes_BO_NOx_0 txt contient les missions de NOx dues a l utilisation de BOis comme combustible pour le chauffage r sidentiel Ces missions sont donn es sous forme de tableau avec en ligne l heure de simulation et en colonne le nombre de zones consid r es Le calcul des missions est rapide de l ordre d une vingtaine de secondes pour l ensemble des deux types d missions avec 1434 brins de trafic routier et 287 zones de chauffage r sidentiel 4 2 6 Ex cution du pr processeur de sp ciation Pr cal Le pr processeur de sp ciation voir figure 11 est lanc partir du menu Saisie Des Donn es Donn es Dispersion Pr processeur sp ciation de l environnement Pollux 2 pre_cal2000 40 deal al SES Al NUMBER OF COMPOUNDS READ FROM COMPOUND DATA BASE SPECIES LIST FROM COMPOUND DATA BASE OLE PA
70. tput free formatted user input d e g IFORMP PMERGE output free formatted user input format unformatted CALM IFORMC ET unformatted output Default Default Default STATIONS fault fault 0 OUTPUT 2 NSSTA NPSTA ICLOUD IFORMS IFORMP IFORMC free formatted CALMET output or user input 21 H PUT GROUP WIND FIELD MODEL OPTIONS Model selection variable 5 Wind Field Options and Parameters IWFCOD 0 Objective analysis only 1 Diagnostic wind module Compute Froud number adjustment effects 0 NO 1 YE IFRADJ S Compute kinemat 0 NO 1 YE use O Brien procedure for adjustment of the vertical velocity NO 1 YES 0 Comput flow ic S ffects ffects slop 0 NO Il K ES IKINE IOBR ISLOPE Ecole des Mines de Douai D partement Energ tique Industrielle Default Default Default Default Default IWFCO IFRAD IKIN El IOBR ISLOP D J ti 36 Extrapolate surface wind observations to upper layers IEXTRP Default 4 TEXTRP 1 l 1 no extrapolation is done 2 power law extrapolation used 3 user input multiplicative factors for layers 2 Z used s FEXTRP array 4 similarity theory used ty y 3 4 same
71. ulateur ETUDES ExEtude ExSimu Donnees UP16 DAT IEND UP17 DAT input ey UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP17 DAT IEND UP18 DAT input 18 UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP18 DAT IEND UP19 DAT input L9 11 UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP19 DAT IEND UP20 DAT input 20 1 UPDAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees UP20 DAT IEND Subgroup c Overwater station files one per station Default Name Type File Name SEA01 DA input dt SEADAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees SEA01 DAT IEND SEA02 DA input Zi l SEADAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees SEA02 DAT IEND SEA03 DA input 3s l SEADAT C Pollux Simulateur ETUDES ExEtude ExSimu Donnees SEA03 DAT IEND Subgroup d Other file names Default Name Type File Name DIAG DAT input DIADAT PROG DAT input PRGDAT EST PR output TSTPRT x EST OU output TSTOUT E EST KIN output TSTKIN Es EST FRD output TSTFRD x EST SLP output TSTSLP El NOTES 1 File path names can be up to 70 characters in length 2 Subgroups a and d must have ONE END surround by delimiters at the end of the group 3 Subgroups b and c must have an delimiters at the end of EACH LINE E ND surround by END
72. urface wind components wind field module IDIOPT4 O Read WS WD from a surface Data file SURF DAT 1 Read hourly preprocessed U a data file DIAG DAT Observed upper air wind components for wind field module IDIOPT5 0 Read WS WD from an upper air data file UP1 DAT UP 1 Read hourly preprocessed U a data file DIAG DAT E BREEZE INFORMATION Default 0 IDIOPTI 0 s rom r T No default TSURFT 9 TA Default 0 IDIOPT2 Qu el rvations alues UPT No default TUPT TO TA scale Default 200 ZUPT 200 Units meters Default 0 IDIOPT3 Qu al rvations alues Default 1 TUPWND 41 STA h Defaults 1000 ZUPWND 1 1000 Units meters Default 0 IDIOPT4 V from Default 0 IDIOPT5 2 DAT etc V from Use Lake Breeze Modul LLBRE ZE Ecole des Mines de Douai D partement Energ tique Industrielle 39 Default F LLBREZE 0 Number of lake breeze regions NBOX NBOX 0 END INPUT GROUP 6 Mixing Height Temperature and Precipitation Parameters EMPIRICAL MIXING HEIGHT CONSTANTS Neutral mechanical equation CONSTB Default 1 41 CONSTB 1 Convective mixing ht equation CONSTE Default 0 15 CONSTE 0 Stable mixing ht equation CONST
73. ustrielle 44 initial concentration distribution is used to obtain time independent top boundary conditions Stationary Point Source Emissions File with cyclical or constant missions PTECYC DAT Default 1 ITEM1 2 ITEM1 1 if the unformatted direct access PTECYC input file is used 2 if the PTECYC file is not used H E Il Stationary Point Source Emissions File with arbitrarily varying missions PTEMRB DAT Default 1 ITEM2 2 1 ITEM2 1 if the unformatted PTEMRB input file is used ITEM2 2 if the PTEMRB file is not used Mobile Point Source Emissions File PTEMOB DAT constant missions Default 1 ITEM3 2 ITEM3 1 if the unformatted PTEMOB input file is used ITEM3 2 if the PTEMOB file is not used Area Source Emissions File AREM DAT Default 1 ITEM4 1 ITEM4 1 if the unformatted AREM input file is used ITEM4 2 if the AREM file is not used Distribution Function for Area Source Emissions Area source missions can be distributed among several layers according user specified weighting factors Number of missions layers NZEM Default 2 NZEM 2 NOTE Up to MXNZ layers are allowed wher XNZ is defined in the parameter fil Fraction of mass distributed into each user defined missions layer WTEM nzem
74. utier r partition horaire et journali re du trafic routier v hicules munis de syst me de contre vaporation diff rents types de combustibles pour le calcul des missions dues au chauffage r sidentiel r partition du type de logements r partition du trafic a rien suivant le type d appareils qui permettront a condition de conna tre ces param tres d obtenir des simulations tr s repr sentatives de la r alit Au premier abord les valeurs des missions NOx et COV nous semblent tout fait bien valu es pas SAMAA Des fonctionnalit s int ressantes concernant l emploi des diff rents types de zonages pour les missions dues au trafic routier m thodologie de d finition des zones conomiques et r sidentielles bas e sur les donn es INSEE possibilit au travers du zonage juridique de n autoriser que le roulage de certains types de v hicules Un s quen age rationnel des diff rentes phases d ex cution des modules tout d abord calculs m t o d finition et validation des missions ex cution des missions et finalement ex cution de la chimie dispersion avec la possibilit de visualiser les r sultats interm diaires Cette phase de visualisation interm diaire permet de valider ou non la coh rence des r sultats par exemple m t o et missions avant de lancer le calcul de la chimie dispersion ceci permettant d viter de perdre du temps dans le cas ou l utilisateur a commis une err
75. utilisateur confirm Il est m me conseill de lire ce fichier afin de v rifier si toutes les donn es ont t correctement interpr t es par l interface Ecole des Mines de Douai D partement Energ tique Industrielle Figure 3 Probleme de relecture du fichier calmet inp 4 2 2 Ex cution de Calmet L ex cution de Calmet entra ne l ouverture d une fen tre DOS avec affichage de la progression des calculs voir figure 4 Calmet engendre galement un fichier texte de sortie calmet Ist que l on peut diter sous Wordpad par exemple pour v rifier si les donn es ont t bien interpr t es par l ex cutable Cual ee eee Tisan Smile Lec eae Miche Pape J DARIE HE en E a IEE Nm a Figure 4 Fen tre DOS ouverte par Calmet et progression des calculs A la fin de l ex cution de Calmet Caltemp module de calcul des champs de temp rature est ensuite lanc automatiquement Ecole des Mines de Douai D partement Energ tique Industrielle 10 En ce qui concerne le temps de calcul celui ci est tout a fait acceptable Le calcul des conditions m t orologiques sur le cas de Nantes a dur 739 secondes dur e donn e a la fin du fichier calmet lst soit environ 12 mn pour un maillage X Y Z de 80 90 10 pour une simulation sur 24 heures avec 20 stations d altitude et 3 stations maritimes 4 2 3 Saisie des donn es g n rales param tres g n raux Dans cette fen tre voir fi

ECOLE DES MINES DE DOUAI

Contents

Download Pdf Manuals

Related Search

Related Contents