version 2020r1

Training courseMarch 8, 9, 10, 11 2021

PM10 surface concentrations in the Alps with 1km resolution.

2

Training course program

3

1 List of participants

# Name Institute City/Country e-mail

1 Farouk Lemmouchi PhD, LISA Paris, France flemmouchi@lisa.ipsl.fr

2 Sergio Ibarra Espinosa Post Doc, USP Sao Paulo, Brazil sergio.ibarra@usp.br

3 Ludovico di Antonio PhD, LISA Paris, France ludovico.diantonio@lisa.ipsl.fr

4 Alexis Squarcioni PhD, CEREA Paris, France squarcionialexis@yahoo.fr

5 Marko Ravina Post Doc, Turin Polytechnic Turin, Italy marco.ravina@polito.it

6 Taos Benoussasd PhD, LISA Paris, France taos.benoussaid@lisa.ipsl.fr

7 Krisztina LÃazÃar Meteorologist, Hungarian Meteorological Service Hungary lazar.k@met.hu,toth.a@met.hu

8 Anita TÃsth Meteorologist, Hungarian Meteorological Service Hungary toth.a@met.hu

9 Zita Ferenczi Meteorologist, Hungarian Meteorological Service Hungary ferenczi.z@met.hu

10 Quentin Marsal Ms, IPGP, Paris, France Quentin.Marsal98@outlook.fr

11 Maxence Descheemaecker IE, Atmo-Occitanie Toulouse, France Maxence.DESCHEEMAECKER@atmo-occitanie.org

12 Zeinab Salah Mahmoud Researcher, Egyptian Meteorological Authority Egypt zeinabsalah@gmail.com

13 Andy Delcloo Researcher at the Royal Meteorological Institute in Belgium Belgium Andy.Delcloo@meteo.be

14 Karen De Causmaecker Researcher at the Royal Meteorological Institute in Belgium Belgium karendc@meteo.be

15 Prashant Singh PhD, University of Mumbai Mumbai, India prashant.singh@physics.mu.ac.in

16 Sunny Kumar PhD, LISA Paris, France skumar@lisa.ipsl.fr

17 Harry Dupont Engineer, Atmo Auvergne RhÃt’ne Alpes France hdupont@atmo-aura.fr

18 Virginie HutsemÃl’kers IR, irCELine Belgium hutsemekers@irceline.be

19 Elke Trimpeneers IR, irCELine Belgium trimpeneers@irceline.be

20 Nishi Srivastava Birla Institute of Technology, Mesra Mesra, India nsrivastava@bitmesra.ac.in

21 Nyaga Ezekiel Waiguru PhD, LISA Paris, France ezekielwaiguru24@gmail.com

22 Carla Gama Researcher, University of Aveiro Aveiro, Portugal carlagama@ua.pt

23 CÃatia GonÃgalves Researcher, University of Aveiro Aveiro, Portugal catialavinia@ua.pt

24 Serafim Kontos PhD, Aristotle University of Thessaloniki Thessaloniki, Greece serkontos@gmail.com

25 Neeldip Barman PhD, Indian Institute of Technology Guwahati Guwahati, India neeldip@iitg.ac.in

26 Juan Henao Researcher, Universidad de Antioquia Antioquia, Colombia jjose.henao@udea.edu.co

27 Victor Almanza Researcher, UNAM Mexico City, Mexico victor.almanza@atmosfera.unam.mx

28 Agustin Garcia Researcher, UNAM Mexico City, Mexico agustin@atmosfera.unam.mx

29 Alex Resovsky Engineer, LSCE, Paris, France aresovsk@gmail.com

30 Bo Zheng Assistant Professor, Tsinghua University China bozheng@sz.tsinghua.edu.cn

31 Carlos Ross Engineer, DFM Consultores Ambientales Chile cross@dfmconsultores.cl

32 Zhizhao WANG PhD, CEREA Paris, France zhizhao.wang@enpc.fr

33 Eduardo Illueca FernÃandez PhD, University of Murcia Spain eduardo.illueca@um.es

34 Alejandro Martà nez MÃl’ndez PhD, University of Murcia Spain alejandro.martinez@hopu.eu

35 Phong Nguyen Hoang Research Assistant, Digital Control and System Engineering Vietnam nhphong@dcselab.edu.vn

36 Damien Piga Engineer, AtmoSud France damien.piga@atmosud.org

37 Sonia Oppo Engineer, AtmoSud France sonia.oppo@atmosud.org

38 Fernanda Garcia CONAE Argentina fgarciaferreyra@conae.gov.ar

39 Arthur Elessa Postdoc, LISA Paris, France aelessa@lisa.ipsl.fr

40 Sandrine Guhirwa Forecasting officer, Meteo Rwanda Rwanda s.guhirwa@meteorwanda.gov.rw

4

Contents

1 List of participants 3

2 Install and run the CHIMERE model 52.1 Prepare for the compilation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52.2 Exercise 1 : Compilation of OASIS and CHIMERE . . . . . . . . . . . . . . . . . . . . . . . 62.3 Exercise 2 : Compilation of WRF and WPS . . . . . . . . . . . . . . . . . . . . . . . . . . . 62.4 Exercise 3: Run the offline test-case . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

3 The output files 83.1 Exercise 4 : Check all files are there . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 83.2 Exercise 5 : Check files contain what they should . . . . . . . . . . . . . . . . . . . . . . . . 9

4 The coupled WRF-CHIMERE system 104.1 Exercise 6 : Run the online test-case . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 104.2 Exercise 7: Aerosol direct and indirect effects seen with the WRF-CHIMERE system . . . . . 10

5 The CHIMERE anthropogenic emissions preprocessor 125.1 Exercise 8 : Install and run emiSURF . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 125.2 Exercise 9 : The emiSURF output files . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13

6 Horizontal domain and nesting 156.1 Exercise 10 : Create your own domain . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 156.2 Exercise 11 : Create a nest inside your domain . . . . . . . . . . . . . . . . . . . . . . . . . . 16

5

2 Install and run the CHIMERE model

Each participant has a personal login (pollution1, pollution2...). The first step for all tutorials is to connect tothe lmd system. The command you need to type in your terminal is:

ssh -Y login@calcul.lmd.polytechnique.fr

Once connected you are at your home directory: /home/loginAt your home directory you will find two directories : the afternoon_session and morning_session. Go to thedirectory corresponding to your session.

cd /home/login/morning_sessionorcd /home/login/afternoon_sessionls

For this training and on each account the required files are already downloaded.You will find two directories :

• chimere_v2020r1 : This is the CHIMERE directory. It includes three models: WRF, WPS (its pre-processor)and CHIMERE, as well as all necessary scripts to compile and run the WRF-CHIMERE model.

• emiSURF2020 : This is the anthropogenic emissions program

2.1 Prepare for the compilation

The first tutorial consists of compiling the WRF-CHIMERE system and running a first preconfigured test run.

cd /home/login/mysession/chimere_v2020r1

This directory will be called CHIMERE directory for the rest of the tutorials.You need to compile the four components of the system :

1. OASIS3-MCT

2. CHIMERE

3. WRF

4. WPS

The obtained executables will be :

1. chimere.e : the CHIMERE model

2. wrf_n0.exe : the WRF model ( used in online mode )

3. the WRF Preprocessing System : geogrid.exe, metgrid.exe, ungrib.exe. geogrid is needed both inoffline and online mode to build the domain grid.

4. the oasis3-mct library, needed for online simulation coupled with WRF.

Everything has to be compiled the first time :

1. oasis3-mct is strongly embedded into CHIMERE and needs to be compiled even if the user uses only theoffline mode.

6

2. WPS is used for the preparation of the CHIMERE domains even in offline mode. Therefore, WPSneeds to be compiled even for offline simulations. For WPS to compile properly you need to have WRFcompiled. So WRF too must be compiled even if the user uses only the offline mode.

Compilers and external libraries path as well as compiler options are controlled using an architecture configurefile in ./mychimere/mychimere-’my_arch’, where ’my_arch’ is a name the user decides to give to its computingarchitecture. This name format for the achitecture configure is mandatory, ’my_arch’ suffix is left free.

For this Training we will be using the file CHIMERE/mychimere/mychimere-chns.gfortran configuration filefor the compilation. Therefore, here ’my_arch’ = chns.gfortran.Two compilation scripts take care of those tasks : build-chimere.sh and build-wrf.sh. General help messagesare displayed this way :

./build-chimere.sh --help

./build-wrf.sh --help

2.2 Exercise 1 : Compilation of OASIS and CHIMERE

After filling out the mychimere/mychimere-my_arch file, to compile OASIS and CHIMERE, the user has toexecute the following :

./build-chimere.sh --arch chns.gfortran

Normally, last lines of the output should be :

Compiling OASIS...Compiling CHIMERE...Compilation OK. Saving compiled code to CHIMERE/exe_PROD...

Log files of the compilations will be found in : compilogs/make.oasis.???.log and compilogs/make.???.log.

2.3 Exercise 2 : Compilation of WRF and WPS

To compile WRF and WPS, the user has to execute the following :

./build-wrf.sh --arch chns.gfortran

Normally, last lines of the output should be :

Compiling WRF1. Compiling CHIMERE/WRFV3_cpl371/main/real_nO.exe without OASIS2. Compiling CHIMERE/WRFV3_cpl371/main/wrf_O.exe with OASIS

Compiling WPS...Compilation Success

Log files of the compilations will be found in : compilogs/make.em_real_without_oasis.????.log, compilogs/-make.wrf_with_oasis.????.log and make.wps.????.log.

2.4 Exercise 3: Run the offline test-case

The first thing you need to do is to check that all executables have been created.Now you may run a simple offline test-case.There is a single file you need to configure in order to run this first case. It is theCHIMERE/chimere.par.test.offline file. Open the file and search for the _ _FILL_ _ signs. You will haveto replace _ _FILL_ _ with the corresponding directories.

7

[simuldir] Simulation output directory : __FILL__[bigfilesdir] Big files (AEROMIN.bin, etc) : __FILL__[emissdir] Anthropogenic emission Dir (stage 1) : __FILL__[$meteo_DIR] Where are your meteo files : __FILL__ # Path to file wrfout[$wrf_geog] If WRF, where is geogfile : __FILL__# Path to file geo_em.d01.nc

Note that when running a new simulation, if you don’t wish to lose the previous simulation you need to changethe [simuldir] flag

• [simuldir] : where outputs will be written. Here you can simply put them in /home/login/mysession/CHI-MOUT

• [bigfilesdir] : where the BIGFILES directory is on the disc. Here they are : /data/-PLT8/pub/mvalari/BIGFILES2020

• [emissdir] : where are the anthropogenic emission files for the simulation. Here they are : /data/-PLT8/pub/mvalari/TestCase2020/INPUTS/EMI

• [meteo_DIR] : full path to the directory where the wrfout file is Here it is: /data/-PLT8/pub/mvalari/TestCase2020/INPUTS/WRFtest

• [wrf_geog] : full path to the geog file. This file is obtained by the WRF model preprocessor (WPS)and contains the geographical coordinates of the CHIMERE simulation domain. Here it is: /data/-PLT8/pub/mvalari/TestCase2020/INPUTS/WRFtest/geo_em.d01.nc

Once these paths are set, the program can be run with the following command:

./chimere.sh chimere.par.test.offline f 2017010100 24

8

3 The output files

You have hopefully run the offline test-case simulation.In this section you are guided on how to check that the simulation run as expected.

3.1 Exercise 4 : Check all files are there

Check the simuldir directory and make sure all output files are creted : /home/myaccount/mysession/CHI-MOUT.ls /home/login/mysession/CHIMOUT

Normally you should find there three netcdf files:

• out.2017010100_24_test.offline.nc• dep.2017010100_24_test.offline.nc• end.2017010100_24_test.offline.nc

CHIMERE has preprocessors to estimate anthropogenic and natural emission fluxes.All output of these programs are stored in the /home/login/mysession/CHIMOUT/data_$lab

• AEMISSIONS : Anthropogenic emissions• BEMISSIONS : Biogenic emissions• SEMISSIONS : Sea-dalt emissions• DEMISSIONS : Dust emissions

Figure 1: Off-line test-case simulation

Boundary and initial condion files are stored in the /home/login/mysession/CHIMOUT/IBC directory :

• BOUN_CONCS.2017010100_24_test.offline.list• BOUN_CONCS.2017010100_24_test.offline.nc-aer• BOUN_CONCS.2017010100_24_test.offline.nc-dust• BOUN_CONCS.2017010100_24_test.offline.nc-gas• INI_CONCS.2017010100_24_test.offline.list• INI_CONCS.2017010100_24_test.offline.nc-aer• INI_CONCS.2017010100_24_test.offline.nc-dust• INI_CONCS.2017010100_24_test.offline.nc-gas

9

3.2 Exercise 5 : Check files contain what they should

• You may use the ncdump -h name_of_the_netcdf_file command to check that all time-slots are written in theoutput file

• You may use the ncview name_of_the_netcdf_file to plot maps of variables of interest• Click on point in the map to see the diurnal cycle is realistic• Plot vertical profiles to see if the variable evolves reasonable

What are the dimensions of the AEMISSIONS, DEMISSIONS, SEMISSIONS files?What is the difference between pH2SO4 and bH2SO4?How does the AOD change depending on the region of the electormagnetic spectrum?

10

4 The coupled WRF-CHIMERE system

In this section you are guided to run an online simulation. We also propose to investigate the direct and indirectaerosol effects on atmospheric dynamics and chemical composition as seen by the coupled WRF-CHIMEREsystem.

4.1 Exercise 6 : Run the online test-case

CHIMERE and WRF will run together exchanging information at $dtphys frequency.This is the cpl_case=1.The user has to modify chimere.par.test.online to run this test case. The lines that need to be modified areindicated with the __FILL__ sign. Some modifications are common in both offline and online cases (seeprevious section) :

• [simuldir] : where outputs will be written. Here you can simply put them in /home/login/mysession/CHI-MOUT

• [bigfilesdir] : where the BIGFILES directory is on the disc. Here they are : /data/-PLT8/pub/mvalari/BIGFILES2020

• [emissdir] : where are the anthropogenic emission files for the simulation. Here they are : /data/-PLT8/pub/mvalari/TestCase2020/INPUTS/EMI

In addition, for the online run you also need to specify the path to the necessary input meteorological data forthe WRF simulation (boundary and initial conditions):

[dgrb] Dir for input meteo grib files : /data/PLT8/pub/mvalari/TestCase2020/INPUTS/AVN

Once these paths are set, the program can be run with :

./chimere.sh chimere.par.test.online f 2017010100 24

Additionally, the online run, this simulation has also produced two new files :

• wrfout file : wrfout_d01_2017-01-01_00:00:00

• wrf restart file for the next simulation : wrfrst_d01_2017-01-02_00:00:00

4.2 Exercise 7: Aerosol direct and indirect effects seen with the WRF-CHIMERE system

In this section we invite you to explore the output of a simulation that was ran with the four coupling optionsprovided with the WRF-CHIMERE system.In all cases WRF sends meteorological fields to CHIMERE every $dtpphys minutes.

• cpl_case = 1 : WRF recieves no feedback from CHIMERE• cpl_case = 2 : CHIMERE sends aerosol concentration and optical properties to WRF every $dtpphys minutes.

Aerosol direct effect on SWR and LWR radiation is accounted for.• cpl_case = 3 : CHIMERE sends aerosol concentration and optical properties to WRF every $dtpphys minutes.

Aerosol indirect effect on cloud formation by increase in CCN is accounted for.• cpl_case = 4 : Both direct and indirect effects are accounted for.

The output of these runs is accessible here : /data/PLT8/pub/mvalari/AFR50_COUPLAGE

11

• out.2010022700_24_AFR50-cpl1.nc• out.2010022700_24_AFR50-cpl2.nc• out.2010022700_24_AFR50-cpl3.nc• out.2010022700_24_AFR50-cpl4.nc

In the same directory you will find three netcdf files that contain differences between each coupled run and thereference simulation cpl_case=1

• diffout.2010022700_24_AFR50_cpl2-cpl1.nc => Differences in the variables in this file highlights the directeffect of aerosols in atmospheric dynamics and chemistry

• diffout.2010022700_24_AFR50_cpl3-cpl1.nc => Differences in the variables in this file highlights the indi-rect effect of aerosols in atmospheric dynamics and chemistry

• diffout.2010022700_24_AFR50_cpl4-cpl1.nc => Differences in the variables in this file highlights both di-rect and indirect effects of aerosols in atmospheric dynamics and chemistry

Figure 2: On-line simulations. Investigation of the direct and indirect effect of aersols seen by the WRF-CHIMERE system.

What variables in the CHIMERE output should me most affected by each type of coupling?In what regions of the domain should one expect to see the largest effect?Which effect is stronger the direct or indirect?NB : You can use the ncview utility to answer these questions in a qualitative manner.

12

5 The CHIMERE anthropogenic emissions preprocessor

The emiSURF model is the pre-processor of the anthropogenic emission fluxes used in atmospheric compositionsimulations with the CHIMERE chemistry transport model. The emiSURF model can be used as a stand-alonemodel but in this case, the user needs to be aware of all the modeling stages before adapting the model outputto other uses. The model takes as input gridded monthly or annual emission totals from several availableinventories. The main operations applied by the emiSURF model on the inventory data are listed bellow:

• horizontal regridding of emissions from the inventory grid to the CHIMERE simulation grid• distribution of emissions to vertical layers• chemical speciation profiles to map emissions from inventory pollutants to ’model species’.• temporal profiles to obtain hourly fluxes

This top-down approach results in a set of output files: one file per emitted "model species" and month. Monthlytotals in these files are distributed across the seven days of the week and the 24 hours of the day through the"type_day" and "Time" dimensions respectively.

Figure 3: Top-down processing of anthropogenic emissions with the emiSURF program.

5.1 Exercise 8 : Install and run emiSURF

The emiSURF is a program available for download ih the CHIMERE webpage. Here it is already downloadedin you account.Go to the emiSURF directory :

cd /home/login/mysession/

NB: To compile emiSURF you need to add in the emisurf.par namelist the path to your CHIMERE directoryThen you may compile emiSURF with the following command:

./build-emisurf.sh --par emisurf.par --arch {\bf my\_arch}

where my_arch = chns.gfortran here

[$chimeredir] : __FILL_[$outdir] : __FILL__[$geogfile] : __FILL__[$bigfiles] : __FILL__[$camsdir] : __FILL__

13

A test case for the emiSURF code is available at the CHIMERE model download page. Here all necessary dataare downloaded. You may fill in the emisurf.par.test file as follows:

• [chimeredir] : /home/login/mysession/chimere_v2020r1

• [outdir] : /home/login/mysession/EMISURFOUT

• [geogfile] : /data/PLT8/pub/mvalari/TestCaseEMISURF2020/INPUTS/geog/geo_em.d01.nc

• [bigfiles] : /data/PLT8/pub/mvalari/BIGFILES2020

• [camsdir] : /data/PLT8/pub/mvalari/TestCaseEMISURF2020/INPUTS/CAMS

Once you have finished configuring the emisurf.par file you may execute the program with the command :

./emisurf.sh emisurf.par.test

5.2 Exercise 9 : The emiSURF output files

The emiSURF program provides one file for each model species and each month. Monthly totals for a givenpollutant are further distributed in hourly fluxes in two steps:

• Weekly profiles : This is the "type_day" dimension.• Daily profiles : This is the "Time" dimension

Most anthropogenic emissions are emitted on the surface layer. However some sectors may emit at highervertical layers. The emiSURF program applies vertical profiles on the emission of some sectors. The verticallayers are the 11 EMEP vertical layers. These data are in the file inputdata/misc/HEIGHT_EMEP_update.Look the structure oh the emiSURF program with the following command:

ncdump -h /home/login/mysession/EMISURFOUT/EMI_2017/EMIS.IDF5.01.NO2.s.nc

netcdf EMIS.IDF5.01.NO2.s {dimensions:Time = 24 ;west_east = 45 ;south_north = 48 ;bottom_top = 7 ;type_day = 7 ;SpStrLen = 23 ;nlentzone = 50 ;nmaxdate = 62 ;nmaxzone = 1 ;dlen = 19 ;variables:

....float NO2(Time, type_day, bottom_top, south_north, west_east) ;NO2:units = "molecule/cm2/s" ;NO2:long_name = "NO2 Emission" ;....

Make maps, time-series and vertical profiles for different emitted pollutants.

14

Figure 4: emiSURF test-case simulation

15

6 Horizontal domain and nesting

For the following two exercises you will use the chimere.par.test.nest

6.1 Exercise 10 : Create your own domain

Use the WrfDomainWizard found in your /home/login/afternoon_session directory to to create a domain whereyou would like to run CHIMERE.The command is :cd /home/login/afternoon\_session/WrfDomainWizardjava -Xmx675m -jar WRFDomainWizard.jar

You need to provide some information :

• WPS Programs : /home/login/mysession/chimere_v2020r1/WPS4.1• Geography : /data/PLT8/pub/mvalari/BIGFILES2020/geog_data_v4.1

and the you may create the new domain

• Create the new domain mydomain• Chose the Projection• Update the map• Define the resolution (Grid points distance)• Adjust your map by playing with the horizontal dimensions X and Y• exit

You should now find the new directory /home/login/afternoon_session/WrfDomainWizard/mydomain.The namelist.wps file there contains all the necessary information for the domainlist.nml

• e_we• e_sn• dx• dy• ref_lat• ref_lon• Proj

NB: In some cases, depending on the latitude, you may need to adjust some hard-coded values or add newvariables in the namelists_wrf/namelist.wps.sed file.Now to run CHIMERE with this new domain you will use the chimere.par.test.nest.You will need to search the __FILL__ signs and replace it with the values corresponding to your case.

• provide the new domain name : mydomain to the chimere.par.test.nest file• provide a prefix in the simuldir variable letting the postfix lab(/home/mylogin/mysession/CHIMOUTNEST/lab)• If not activated, activate the istopdom flag in chimere.par.test.nest

For more accuracy you may also want to change truelat1 and truelat2 in the namelists_wrf/namelist.wps.sedfile to the values of the namelist.wps provided with WrfDomainWizard.truelat1 = 30.0,truelat2 = 60.0,

Once all this done you can :

• Run CHIMERE• Check the geog file in domains/mydomain

16

6.2 Exercise 11 : Create a nest inside your domain

We will now create a nest inside the domain you created in Exercise 10.We will assume that the parent domain is called "myparent" and has :

• nx1 cells in the west-east dimension• ny1 cells in the south-north dimension• dx1 grid-spacing in the west-east dimension• dy1 grid-spacing in the south-north dimension

You can now try and add a nest inside this domain so that the following conditions are satisfied:

• the nest starts at the 1/2 of the coarse domain at each dimension (west-east and south-north)• the nest stops at the 2/3 of the coarse domain at each dimension (west-east and south-north)• the grid-spacing is 1/3 of the grid-spacing in the coarse domain

We will call this nest "mynest". You need to add two lines in the domains/domainlist.nml for the new nesteddomain mynest :In the first section of the domains/domainlist.nml you have to add a line with the number of grid-cells andgrid-cpacing in both dimensions for this nest "mynest".

# Each domain definition:# Domain e_we e_sn dx dy geores dt ref_lat ref_lon Projmynest nx2 nx2 dx2 dx2 same_as_parent ....

where nx2, nx1, dx1 and dx2 are the number of grid-cells and grid-spacing in the west-east and south-northdimensions of the mynest domainBased on the conditions given above you may estimate:

• iparst = nx1 / 2• jparst = ny1 /2• iparend = (2/3) * nx1• jparend = (2/3) * ny1• dx2 = dx1 / 3• dy2 = dy2 / 3• nx2 = (iparend - iparst) * 3• ny2 = (jparend - jparst) * 3

In the second section of the domains/domainlist.nml you need to provide iparst and jparst.Of course nx, ny, iparst, jparst being indexes in the arrays they are integers

# Nested domains (parent i,j are given by DomainWizard or manually chosen)# Dom1 Dom2 grid_ratio iparst jparstNEST mynest myparent 3 iparst jparst

Now you can go back to the chimere.par.test.nest and do the following:

runs : 2[$dom] CHIMERE domain : myparent,mynest[$clab] Coarse domain label : none,test.nest.myparent[$simuldir] Simulation output directory : /home/mylogin/mysession/CHIMOUT_TRAINING/${lab}[$istopdom] Make domains and exit : 1

Once all this done you can :

• Run CHIMERE

17

• Check the geog files in domains/mydomain

To see if these files correspond to what you expect you can use the python script check_nest.py as follows :

python check_nest.py chimeredir prentdomname nestdomname

With

• chimeredir : the chimere directory• prentdomname : the name of your parent domain• nestdomname : the name of the nest

You should obtain something like this :

Figure 5: Parent and nest domains

import sysimport matplotlib.pyplot as pltimport xarray as xr

chimeredir=sys.argv[1]parentdom=sys.argv[2]nestdom=sys.argv[3]

d1=xr.open_dataset(chimeredir+’/domains/’+parentdom+’/geog_USGS_’+parentdom+’.nc’)d2=xr.open_dataset(chimeredir+’/domains/’+nestdom+’/geog_USGS_’+nestdom+’.nc’)

18

plt.plot(d1.XLONG_M[0,:,:],d1.XLAT_M[0,:,:],’k.’)plt.plot(d2.XLONG_M[0,:,:],d2.XLAT_M[0,:,:],’r.’)

plt.show()