Dear all,
I just downloaded CAMS data for the year 2017 for a simulation with WRF-Chem.
I was wondering if you have any advice on how to prepare the data to be ingested by WRF-Chem.
Hope to have some feedback soon
Kind Regards
Joana
Dear all,
I just downloaded CAMS data for the year 2017 for a simulation with WRF-Chem.
I was wondering if you have any advice on how to prepare the data to be ingested by WRF-Chem.
Hope to have some feedback soon
Kind Regards
Joana
30 Comments
Xiaobo Yang
Dear Joana,
I cannot promise you anything at the moment. But I want to let you know that we are looking into this for you and hopefully something will come up in the next couple of weeks/months.
Kind regards,
Xiaobo
Joana Medeiros
Dear Xiaobo Yang
I see that Bojan Kasic posted an answer for this How to run WRF-Chem using CAMS initial and boundary conditions , on the confluence documentation : [confluence] User Documentation > How to run WRF-Chem using CAMS initial and boundary conditions
Although the answer provided is exactly what I need , I can't access the folder/scripts to solve the question, so it's a bit pointless having the answer but being unable to realise it .
Could you help please ?
Many thanks
Kind Regards
Joana
Xiaobo Yang
Dear Joana,
Well spotted. It is a document in working progress. My colleague doing a lot of tests for the time being. Once we feel it is ready, we will certainly share the information with you. So stay tuned.
Kind regards,
Xiaobo
Joana Medeiros
Thanks Xiaobo Yang, I only spotted it because we got an e-mail notification regarding this and when I try to follow the link I couldn't access it.
Will wait for it to be available then
Many Thanks
Regards
J
Joana Medeiros
Dear Xiaobo Yangthe page is now public (How to run WRF-Chem using CAMS initial and boundary conditions on ECMWF's Cray HPC (CCA[B])), but is a bit pointless if you don't have access to the scripts. I posted this question many months ago, and have been waiting for an answer which Bojan Kasicreplied , but I can't solve it because I don't have access to CAMS_to_WRF_utility/ folder. Is there anyway I can access the scripts ? Hope I can have some sort of help.
Many Thanks
Kind Regards
Joana
Bojan Kasic
Dear Joana Medeiros,
Thank you for your interest in this topic. As you have already noted, the solution is only available to users with access to our HPC system. We are in the process of providing the solution for general public here and we are waiting for permission to share the peace of code developed externally. Once we have it, we will share it here.
Thanks for your patience,
Bojan
Joana Medeiros
That would be great,
I will hold for any updates on this...hopefully won't take as many months again
Many Thanks
Kind Regards
Joana
Bojan Kasic
Hi Joana Medeiros ,
Here it is:
How to run the WRF-Chem model using CAMS data as initial and boundary conditions (BC)?
Sorry you had to wait so much.
Let us know if you have any problems or questions.
Bojan
Joana Medeiros
Oh wow,
Thank you so much Bojan Kasic
Will try it as soon as possible
Kind Regards
Joana
Joana Medeiros
Dear Bojan KasicI tried to click the link for the public_box.tar.gz but it comes empty, can't download it :
Step 3: Using mozbc to interpolate CAMS data for WRF
With a set of additional utilities, it is possible to use the mozbc package provided by the Atmospheric Chemistry Observations and Modelling Lab (ACOM) of NCAR to pre-process CAMS data with some minor modifications to the original code. The set of additional utilities has been developed as part of the AQMEII project and kindly provided by A. Lupascu (IASS-Posdam) and C. Knote (NCAR): public_box.tar.gz. The package contains:
When I click , it sends me back to the top of the page
Is it available publicly ?
Thanks
J
Joana Medeiros
Dear Xiaobo,
Thanks for the reply,
Will hold for some feedback
Many Thanks
J
Michela Giusti
Could you please refresh the article page and try again?
Let me know if it works.
Thanks.
Michela
Joana Medeiros
Hi, yes it does, just downloaded it with no issues
Many Thanks
Cheers
J
Liz Coleman
Hi all, sorry i'm late to this party, but I have also prepared some scripts to get and process the operational near real time CAMS data and process for use in mozbc. THe species mapping would have to be tailored to your specific wrf-chem set up. I use mostly python. If there's any interest in these, let me know.
Liz
Joana Medeiros
Hi,
Oh That would be very helpful yes, all the help is more than welcoming,
If you don't mind sharing the scripts it would be much appreciated
Thank you very much Liz Coleman ,
Joana
Liz Coleman
get_merge_proc_cams_generic.pyHi , so i'll try talk through the code here, it's not extensively commented.
I use python and cdo system command, called through the python script.
Function 1: i get the data from cams via wget
function 2: resample each file to a smaller domain to speed up the processing. this task is farmed out over available processors. the lon lat box of the new datafile is defined here : domain_red=('-46,46,30,70')
function 3: we merge all the datafiles to one file per timestep. again, this is performed in parallel.
function 4: create a mozart-type file from the merged, reduced domain CAMS dataset. For this, i got the hyam and hybm variable from this site: https://www.ecmwf.int/en/forecasts/documentation-and-support/137-model-levels and wrote them to a text file to be read in by the script.
I changed the timestep to the way it is understood by mozbc.
o_date=YYMMDD and o_datesec= second of the day.
Wrote surpface pressure (PS), taken from lnsp in cams file.
finally, wrote the species to the files, looped, converting gas units to mol/mol .
This file can then be read in by mozbc. Please let me know if any questions, and if you do use this script, an acknowledgement would be much appreciated!
Liz Coleman, C-CAPS, NUI Galway, Ireland
Xiaobo Yang
Many thanks for sharing this, Liz!
Joana Medeiros
This is great
Many thanks for the sharing Liz Coleman
J
dasa damaskova
Hello all,
please, can anybody explain me why we need temperature? I understand surface pressure as we need to calculate model level pressure from P0,PS and hyam hybm. But I am not sure why we need temperature. I went through the codes and I just cant see where it is used. thank you
Bojan Kasic
Hello dasa damaskova ,
I can confirm that in the example given on page: How to run the WRF-Chem model using CAMS data as initial and boundary conditions (BC)?
MACC_BC2MOZART.ncl is not processing temperature on model levels from AER_20100101_20100102_eac4_EU.nc file even if the field is present.
Temperature is needed in meteorological input and boundary conditions.
Libo Zhang
Hello all,
In the example: How to run the WRF-Chem model using CAMS data as initial and boundary conditions (BC)?
I don't know why levels are cut out from level 7 to 60 (from top) in MACC_BC2MOZART.ncl, and it may not consider to convert kg/kg to mol/mol for gases.
Many thanks for any suggestion.
Thank you.
Angel Vela
Hello there,
I am trying to download CAMS data using the GRG python script but I got the following error message
ecmwfapi.api.APIException: 'ecmwf.API error 1: Resource not found: datasets/cams_reanalysis'
Did the name of the dataset change to a new one?
Thanks
Maria Francisca Cardell
Dear all Xiaobo Yang Bojan Kasic Liz Coleman ,
I downloaded CAMS data from global atmospheric composition forecast for a simulation with WRF-Chem following the instructions How to run the WRF-Chem model using CAMS data as initial and boundary conditions (BC)?#Step4Runningthemodel
The process of downloading the grib data with a python request, converting it to netcdf using the script of section 2.2 of the instructions worked well for 2021-08-01 data. So, I could follow next instructions to finally run mozbc and wrf.exe incorporating successfully BC from CAMS to WRF-Chem.
Get_CAMS_fc_grg_2023.py
However, when I downloaded same kind of data for 2022 or 2023 (python request attached before) I got problems converting grib data to netcdf. Specifically running the command cdo -f nc -b 32 -invertlat GRG_XXXX.grib GRG_XXXX.netcdf.
I wonder how could I convert grib data to netcdf for NO2, NO, CO, SO2, pan, CH2O CAMS forecast composition chemical species (for 2022 and 2023), obtaining all the parameters with 4 dimensions (time, lev, lat, lon) and the hybrid coordinates hyai, hybi, hyam, hybm in the same output nc file.
Many thanks,
Kind regards,
Maria
Bojan Kasic
Dear Maria Francisca Cardell ,
Thank you for your detail report. The issue with CDO conversion:
of GRG parameters from latest CAMS is happening because your CDO is built with one of the older versions of ecCodes libraries which isn't compatible with new parameters definition. This can be solved by re-installing CDO against ecCodes 2.28.0.
I have tested with cdo/2.1.1 build with ecCodes/2.28.0 and the conversion to netcdf worked for both, old and new CAMS GRG data.
I hope this helps.
Kind regards,
Bojan
Maria Francisca Cardell
Thank you very much for your help and you promt response Bojan Kasic !
I could finally fix the error following your instructions,
Kind regards,
Maria
Maria Francisca Cardell
Dear Bojan Kasic Xiaobo Yang Liz Coleman
I'm presently using the CAMS global forecast compostion to update the boundary conditions of gaseous species in WRF-Chem using mozbc utility, and found it to be very useful. Now is working using the MOZART option.
I also want to update the aerosol species from CAMS to use the mozcart option to simulate dust storm episodes. Is there a mapping of mozbc utility to convert the aerosol species from CAMS aerosol module to MOZCART aerosol schemes of WRF-Chem?
Here is the list of aerosol species available from CAMS.
Short name Long name
aermr01 Sea Salt Aerosol (0.03 - 0.5 um) Mixing Ratio
aermr02 Sea Salt Aerosol (0.5 - 5 um) Mixing Ratio
aermr03 Sea Salt Aerosol (5 - 20 um) Mixing Ratio
aermr04 Dust Aerosol (0.03 - 0.55 um) Mixing Ratio
aermr05 Dust Aerosol (0.55 - 0.9 um) Mixing Ratio
aermr06 Dust Aerosol (0.9 - 20 um) Mixing Ratio
aermr07 Hydrophobic Organic Matter Aerosol Mixing Ratio
aermr08 Hydrophilic Organic Matter Aerosol Mixing Ratio
aermr09 Hydrophobic Black Carbon Aerosol Mixing Ratio
aermr10 Hydrophilic Black Carbon Aerosol Mixing Ratio
aermr11 Sulphate Aerosol Mixing Ratio
I've seen that there is a mozbc namelist to convert aerosol data from MOZART global model to use it in WRF-Chem with MOZCART option. Find it attached
However, I do not know how to compute these mapping factors since CAMS only has 3 dust bins and 3 sea salt bins while MOZCART WRF-Chem has 5 dust bins and 4 sea salt bins.
MOZCART.inp
Many thanks,
Best regards,
Maria
Xiaobo Yang
Sorry for the reply, Maria. I've asked our scientist to have a look.
Maria Francisca Cardell
Hi Xiaobo Yang
Aren't there any updates about my last question ?
The mapping of mozbc utility to convert the aerosol species from CAMS aerosol module to MOZCART aerosol schemes of WRF-Chem ?
Many thanks,
Maria
Xiaobo Yang
Hi Maria Francisca Cardell,
Unfortunately I haven't got any rely. I sent another message today. I will keep you updated.
Sorry about the delay.
Regards,
Xiaobo
Johannes Flemming
Dear Maria,
Please note, the aerosol and chemistry modelling has been upgraded in June 2023 (48r1). We have now a more detailed representation of aerosols including nitrates, ammonium and secondary aerosols. Please have a look into the newly written documentation (Part 8) . https://www.ecmwf.int/en/publications/ifs-documentation
I can not advice you on the mapping on the dust bin from CAMS to MOZCART. But, all the required information about the size distribution for the CAMS dust and sea salt is given in the documentation. Please note, that sea salt in CAMS is at RH=80%. To calculate dry seal salt you have to convert radius and mass mixing ratio.
Johannes