# EXTRACT emissions from EDGAR\_HTAPv3 with limited geographical domain (South America) **URL:** \ **Category:** SMOKE\ **Tags:** emissions\ **Created:** [January 8, 2024, 3:26am UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610 "2024-01-08T03:26:02Z")\ **Posts on this page:** 20\ **Page:** 1
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [January 8, 2024, 3:26am UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/1 "2024-01-08T03:26:02Z")
Hi there, I was wondering if I can process with SMOKE my little domain (São Paulo) from EDGAR\_HTAPv3 emission files that covers only a big geographical domain (e.g. South America). I was tried to run but unfortunately I got the log message related to smkinven: ```auto Value for WKDAY_NORMALIZE: N returning FALSE NOTE: Setting inventory to use full-week normalizer for weekly profiles Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file Bad NDIMS for "Transport" in "TMPFILE" Could not read "Transport" from gridded inventory file ``` I don’t know if it is mandatory to use all dimensions of EDGAR-HTAPv3, I mean 1800 rows, 3600 cols (all the world and not only South American). Do you have any suggestions? Best regards, Alejandro ---
**Author:** ![eyth.alison](https://avatars.discourse-cdn.com/v4/letter/e/7ab992/32.png) [@eyth.alison](https://forum.cmascenter.org/u/eyth.alison)\ **Post date:** [January 10, 2024, 12:30am UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/2 "2024-01-10T00:30:11Z")
It may be that the gridded inventory list input (EMISLST) is not set up correctly. The error seems to indicate that SMOKE is not parsing the gridded inventories properly, as if there is an incorrect path, or the variable name in the EMISLST file does not match what is in the gridded inventory. The EMISLST input for the road transport sector should look something like this: #LIST GRID #SCC, Pollutant, Variable\_Name, Month, File\_location\_name “TRANSPORT”,“PEC”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_BC.nc” “TRANSPORT”,“CO”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_CO.nc” “TRANSPORT”,“NH3”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_NH3.nc” “TRANSPORT”,“VOC”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_NMVOC.nc” “TRANSPORT”,“NOX”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_NOx.nc” “TRANSPORT”,“POC”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_OC.nc” “TRANSPORT”,“PM10”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_PM10.nc” “TRANSPORT”,“PM2\_5”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_PM2.5.nc” “TRANSPORT”,“PM2\_5\_OTH”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_PM2.5\_OTHER.nc” “TRANSPORT”,“SO2”,“Transport\_Road”,6,“/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_SO2.nc” ---
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [January 10, 2024, 1:07am UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/3 "2024-01-10T01:07:44Z")
Hi Eyth. Thank you for your answer and support. I understand that the problem was in an additional dimension in my EDGAR-HTAPv3 files ( **time** , see below), because when I ran the SMOKE yesterday with two dimensions (lat, lon) and for all the Globe (1800, 3600) , I got “results” but with a wrong location of emissions. ```auto netcdf edgar_HTAPv3_2018_PM2.5 { dimensions: time = 12 ; lat = 230 ; lon = 350 ; variables: float Transport(time, lat, lon) ; Transport:_FillValue = 9.96921e+36f ; Transport:description = "Road_Transport" ; Transport:grid_mapping = "spatial_ref" ; Transport:long_name = "Road Transport" ; Transport:units = "kg m-2 s-1" ; float Industry(time, lat, lon) ; Industry:_FillValue = 9.96921e+36f ; Industry:description = "Industry" ; Industry:grid_mapping = "spatial_ref" ; Industry:long_name = "Industry" ; Industry:units = "kg m-2 s-1" ; ``` On the other hand, I was wondering if I can use a small domain of EDGAR-HTAP that covers only South America as an input to SMOKE or I must create a GRIDMASK\_EDGAR.ncf that covers South America as well as geocodes files? As I was saying, a problem that I’m experiencing is emission outputs from SMOKE. I saw that emission locations are not consistent with my domain (Southeastern Brazil, center in São Paulo). I don’t know if the GRIDDESC for São Paulo domain that was processed with MCIPv5.4 has errors or maybe the IOAPI\_ISPH=20 in my ASSIGN file. I attached my files if you want to see my two scripts (ASSIGN and sms\_edgar\_TRANSPORT\_SEBRA.csh). [ASSIGNS.EDGAR.TRANSPORT.SEBRA.csh](https://forum.cmascenter.org/uploads/short-url/vvUR19rsO1tK5PHOqmoRSnqxMzY.csh) (17.9 KB) ```auto ' ' 'LAM_28S55W' 2 -60.000 -30.000 -55.000 -55.000 -45.000 ' ' '2018_SEBRA' 'LAM_28S55W' -40499.000 1813357.500 9000.000 9000.000 247 177 1 ' ' ``` [smk\_edgar\_TRANSPORT\_SEBRA.csh](https://forum.cmascenter.org/uploads/short-url/xxXXnifMh1rg0Dq08eDUMVCIHjt.csh) (6.7 KB) Best regards, Alejandro ---
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [January 10, 2024, 12:03pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/4 "2024-01-10T12:03:22Z")
Hi Eyth, SMOKE is an interesting tool, but I thought that it could process pre-gridded emissions based on lat lon coordinates like the “anthro\_emis” tool for WRF-Chem. Yesterday, I ran the SMOKE model with the EDGAR-HTAPv2 (the base year 2010), and I got consistent locations of the emissions according to my domain, and that is because the map of EDGAR-HTAPv2 starts with a time zone 0 in Africa that is different compared with the EDGAR-HTAPv3. So, the problem was not the GRIDDESC. This is the EDGAR-HTAPv2: ![htapv2](https://canada1.discourse-cdn.com/flex027/uploads/cmas/original/2X/2/25ba89c51341e9628fc7a8ea8c1c09dde957e275.png) This is the EDGAR-HTAPv3: ![htapv3](https://canada1.discourse-cdn.com/flex027/uploads/cmas/original/2X/e/e61bc9fed44c15d54e67a5cd252177e2178ce6e9.png) BTW, is it possible that you could share us the geocodes and the GRIDMASK\_EDGAR.ncf adapted to the EDGAR-HTAPv3 order? Best regards, Alejandro ---
**Author:** ![eyth.alison](https://avatars.discourse-cdn.com/v4/letter/e/7ab992/32.png) [@eyth.alison](https://forum.cmascenter.org/u/eyth.alison)\ **Post date:** [January 12, 2024, 4:24pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/5 "2024-01-12T16:24:40Z")
Thank you for pointing this out. We will soon be starting some work with HTAPv3 but we have not yet. We will be providing some information on how to process HTAPv3 data in the next couple of months. If I learn anything new before then, I will let you know. ---
**Author:** ![eyth.alison](https://avatars.discourse-cdn.com/v4/letter/e/7ab992/32.png) [@eyth.alison](https://forum.cmascenter.org/u/eyth.alison)\ **Post date:** [January 12, 2024, 5:14pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/6 "2024-01-12T17:14:13Z")
Also, we note the variables are organized a bit differently in the HTAPv3 gridded inventories. Instead of a single variable called “Transport” there is now “Transport\_Road”, “Transport\_BWTW”, “Transport\_Other”, etc. You will need to update the EMISLST input accordingly. Please let us know if you are able to get it to work. ---
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [January 12, 2024, 7:20pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/7 "2024-01-12T19:20:33Z")
Alison, thank you for your answer and support. I understand that previously to list in the EMISLST, I must change variable names in each netCDF with no more than 16 character. I can do that with nco and cdo operators using a script. The only thing that I couldn’t do is to create a new GRIDMASK\_EDGAR\_htapv3.ncf. I’ve read about it how can I do. I will return with you in case I will create one to process EDGA-HTAPv3. Best regards, Alejandro ---
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [February 26, 2024, 2:36pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/8 "2024-02-26T14:36:52Z")
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: ```auto 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: ```auto 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 ---
**Author:** ![eyth.alison](https://avatars.discourse-cdn.com/v4/letter/e/7ab992/32.png) [@eyth.alison](https://forum.cmascenter.org/u/eyth.alison)\ **Post date:** [February 27, 2024, 1:49am UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/9 "2024-02-27T01:49:39Z")
SMKINVEN\_MONTH should be set to 5 if your EDGAR-HTAPv3 files are monthly. We aren’t sure about the starting point for your scripts. If you are using scripts from one of our modeling platforms, then the following applies: We have these two helper scripts which are called by the top-level SMOKE script (along with “3D” versions of each of these for sectors with plume rise): smk\_ar\_gridded\_annual\_emf.csh smk\_ar\_gridded\_monthly\_emf.csh When using the monthly\_emf.csh script, SMKINVEN\_MONTH is set to the appropriate month automatically. So if you are using scripts based on one of our EMPs, they should switch from the annual\_emf.csh script to the monthly\_emf.csh script (near the bottom of the top-level SMOKE script that they run directly) and that will take care of SMKINVEN\_MONTH. ---
**Author:** ![Heomin](https://avatars.discourse-cdn.com/v4/letter/h/9fc348/32.png) [@Heomin](https://forum.cmascenter.org/u/Heomin)\ **Post date:** [December 7, 2024, 11:18pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/10 "2024-12-07T23:18:05Z")
If I set `SMKINVEN_MONTH`, will SMOKE be able to directly process monthly data with dimensions (12, 1800, 3600)? ---
**Author:** ![eyth.alison](https://avatars.discourse-cdn.com/v4/letter/e/7ab992/32.png) [@eyth.alison](https://forum.cmascenter.org/u/eyth.alison)\ **Post date:** [December 9, 2024, 3:14pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/11 "2024-12-09T15:14:22Z")
The raw HTAPv3 gridded files are structured so that all months are in a single yearly file with 12 timesteps. However when processing this data in SMOKE we typically split the months into separate files, and then list each monthly file separately in the Gridded List Inventory input (EMISLST), with the 4th column pertaining to the month number. See the below example. We hope that helps. #LIST GRID #SCC, Pollutant, Variable\_Name, Month, File\_location\_name “ENERGY”,“PEC”,“Energy”,1,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-01\_BC.nc” “ENERGY”,“PEC”,“Energy”,2,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-02\_BC.nc” “ENERGY”,“PEC”,“Energy”,3,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-03\_BC.nc” “ENERGY”,“PEC”,“Energy”,4,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-04\_BC.nc” “ENERGY”,“PEC”,“Energy”,5,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-05\_BC.nc” “ENERGY”,“PEC”,“Energy”,6,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_BC.nc” “ENERGY”,“PEC”,“Energy”,7,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-07\_BC.nc” “ENERGY”,“PEC”,“Energy”,8,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-08\_BC.nc” “ENERGY”,“PEC”,“Energy”,9,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-09\_BC.nc” “ENERGY”,“PEC”,“Energy”,10,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-10\_BC.nc” “ENERGY”,“PEC”,“Energy”,11,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-11\_BC.nc” “ENERGY”,“PEC”,“Energy”,12,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-12\_BC.nc” “ENERGY”,“CO”,“Energy”,1,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-01\_CO.nc” “ENERGY”,“CO”,“Energy”,2,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-02\_CO.nc” “ENERGY”,“CO”,“Energy”,3,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-03\_CO.nc” “ENERGY”,“CO”,“Energy”,4,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-04\_CO.nc” “ENERGY”,“CO”,“Energy”,5,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-05\_CO.nc” “ENERGY”,“CO”,“Energy”,6,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-06\_CO.nc” “ENERGY”,“CO”,“Energy”,7,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-07\_CO.nc” “ENERGY”,“CO”,“Energy”,8,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-08\_CO.nc” “ENERGY”,“CO”,“Energy”,9,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-09\_CO.nc” “ENERGY”,“CO”,“Energy”,10,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-10\_CO.nc” “ENERGY”,“CO”,“Energy”,11,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-11\_CO.nc” “ENERGY”,“CO”,“Energy”,12,“/work/MOD3DEV/ktalgo/WR713/work/inputs/HTAPv3/data/fluxes/monthly/edgar\_HTAPv3\_2018-12\_CO.nc” ---
**Author:** ![Heomin](https://avatars.discourse-cdn.com/v4/letter/h/9fc348/32.png) [@Heomin](https://forum.cmascenter.org/u/Heomin)\ **Post date:** [December 20, 2024, 6:38am UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/12 "2024-12-20T06:38:02Z")
Hi Jano I have a question. The original EDGAR file has data with dimensions (1800, 3600). If I crop this data to match the modeling domain (e.g., (350, 550)), can it still be used to run SMOKE? ---
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [December 24, 2024, 1:54am UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/13 "2024-12-24T01:54:34Z")
Hi Heomin, I tried cropping the data, and the SMOKE couldn’t recognize it. So, the SMOKE only works with global dimensions (1800, 3600). Best regards, Jano ---
**Author:** ![yan2024](https://avatars.discourse-cdn.com/v4/letter/y/59ef9b/32.png) [@yan2024](https://forum.cmascenter.org/u/yan2024)\ **Post date:** [December 30, 2024, 3:40pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/14 "2024-12-30T15:40:54Z")
Hi Jano, Thank you for your guide and I followed what you advised by cdo to process HTAPv3 files but I still get some errors. I doubt if it is the reason for the different attribute of variables. Could you please see my latest thread and kindly give me some suggestions? ---
**Author:** ![eyth.alison](https://avatars.discourse-cdn.com/v4/letter/e/7ab992/32.png) [@eyth.alison](https://forum.cmascenter.org/u/eyth.alison)\ **Post date:** [December 30, 2024, 5:34pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/15 "2024-12-30T17:34:56Z")
Sorry, but we won’t be able to offer more support on this at this time – maybe someone else on the forum has suggestions. Sometime later this winter, we plan to post a package with all needed information to run HTAPv3 through SMOKE, but it likely won’t be for a couple of months. ---
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [January 13, 2025, 9:38pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/16 "2025-01-13T21:38:54Z")
Hi Yan, Sorry for my delay. I have a useful Python script for reducing source names and correcting east coordinates. Those converted files work when I process them in SMOKE v5.0. Please see the file attached. [HTAP\_processing.py](https://forum.cmascenter.org/uploads/short-url/mOHAXq43ByTw9JqvpskAYHhsko3.py) (6.3 KB) Best regards, Jano ---
**Author:** ![yan2024](https://avatars.discourse-cdn.com/v4/letter/y/59ef9b/32.png) [@yan2024](https://forum.cmascenter.org/u/yan2024)\ **Post date:** [January 14, 2025, 2:06am UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/17 "2025-01-14T02:06:33Z")
Hi Jano, Thank you for your kindly help, I used what you advised with cdo to resolve that problem. I found that it is the issue of netcdf format of HTAP-v3 files. So I transferred them to netcdf3 and it worked within smoke successfully. However, I encountered errors when I runned CMAQ 5.4 after processing htap-v3 files with CB6-AE7 MECHS gspro/gsref. I am wondering if you could provide suggesstions if you had runned successfully by CMAQ with htap-v3 inventory. Could you help me check my error logs that are in my latest thread? Best Regards, Yan ---
**Author:** ![Heomin](https://avatars.discourse-cdn.com/v4/letter/h/9fc348/32.png) [@Heomin](https://forum.cmascenter.org/u/Heomin)\ **Post date:** [January 16, 2025, 1:33pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/18 "2025-01-16T13:33:31Z")
Hi Jano, Thank you for the code you provided. However, I have a question. HTAPv3 files are provided on a monthly basis (12,1800,3600). In the ARINV file, you set all months to 0. Could you clarify if you merged all the monthly HTAP files into an annual dataset, or if you extracted only the values corresponding to each month and saved them in the LST file? Best Regards, Heomin ---
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [January 20, 2025, 7:02pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/19 "2025-01-20T19:02:52Z")
Hi Heomin, I downloaded HTAPv3 files on an annual basis (1800, 3600). After that, I also considered setting 0 in the ARINV file. During the temporal profile, I specify the sector variation that could come from monthly HTPAv3 files or other criteria (sources behavior). Best regards, Jano ---
**Author:** ![Jano](https://avatars.discourse-cdn.com/v4/letter/j/9dc877/32.png) [@Jano](https://forum.cmascenter.org/u/Jano)\ **Post date:** [January 20, 2025, 7:11pm UTC](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610/20 "2025-01-20T19:11:13Z")
Hi Yan, I also experienced those bad errors under SMOKE running. I mentioned that because gspro/gsref are related with that program. I fixed the error when I deleted output files and premerged files. It is important review if all species name are correct and the path files. I usually run HTAPv-3 in SMOKE by sector or I sum some sectors to create a new merge sector called “Others”. I processed by each sector for example Industry, Road exhaust, Road resuspension, and Others. After that, I run merge to create files for all time series based on the temporal profile chosen. Finally, depending on the period, I create other emission files related to sea spray, biogenic, biomass burning. Best regards, Jano [Next page](https://forum.cmascenter.org/t/extract-emissions-from-edgar-htapv3-with-limited-geographical-domain-south-america/4610.md?page=2)