EXTRACT emissions from EDGAR_HTAPv3 with limited geographical domain (South America)

Hi Alison, we resolved this issue using the cdo operators. So, we converted from longitudes (-180,180) from EDGAR-HTAPv3 to 0,360 longitudes. The line code in cdo is:

cdo sellonlatbox,0,360,-90,90 $input_file temp.nc

Another solution is with Python. I considered in the script slicing the EDGAR-HTAPv3 into two files:

sector = # add sector name with less than 10 character
range1 = slice(0,180)
range2 = slice(180,0)

v3_1 = v3.sel(lon=range1)
v3_2 = v3.sel(lon=range2)
v3_2['lon'] = v3_2.lon.values + 360
v3_new = xr.concat([v3_1[sector],v3_2[sector], dim='lon')

In SMOKE I run successfully with the year emission files. However, If I want to process monthly emission files in the SMOKE for the temporal step, it will stop with errors. I noted an interesting message in log files that could be related to the running issue:

SMKINVEN_MONTH not defined returning default 0

I understand that I have to set up in the SMK run the SMKINVEN_MONTH = 5 that represents my month (May), is that correct?

Best regards,
Alejandro

2 Likes