InMAP not reading in emissions data

[Broderik Craig]

When I run inmap.exe run stead --Config="config.toml" I get
2024/09/03 13:57:53 Loading emissions shapefile: .../Data/InMAP/inputs/test_emissions.shp.
2024/09/03 13:57:54 Loading CTM data...
2024/09/03 13:57:54 Reading input data...
2024/09/03 13:57:55 Loading population and mortality rate data...
2024/09/03 13:57:55 Initializing model...
2024/09/03 13:57:55 Emission totals:
2024/09/03 13:57:55 gOrg, 0 μg/s
2024/09/03 13:57:55 pOrg, 0 μg/s
2024/09/03 13:57:55 PM2_5, 0 μg/s
2024/09/03 13:57:55 gNH, 0 μg/s
2024/09/03 13:57:55 pNH, 0 μg/s
2024/09/03 13:57:55 gS, 0 μg/s
2024/09/03 13:57:55 pS, 0 μg/s
2024/09/03 13:57:55 gNO, 0 μg/s
2024/09/03 13:57:55 pNO, 0 μg/s

My emissions file has a PM25 column and a geometry column like this
        PM25                 geometry
0 98.009382 MULTILINESTRING ((-1243902.104 179806.229, -12...
1 98.031927 MULTILINESTRING ((-1246668.739 188252.389, -12...
2 101.078069 MULTILINESTRING ((-1243554.297 179874.291, -12...
3 100.145103 MULTILINESTRING ((-1244113.632 178168.202, -12...
4 100.145103 MULTILINESTRING ((-1244113.632 178168.202, -12...
... ... ...
326 72.305537 MULTILINESTRING ((-1241056.073 188074.166, -12...
327 96.410829 MULTILINESTRING ((-1241056.073 188074.166, -12...
328 72.306666 MULTILINESTRING ((-1241056.073 188074.166, -12...
329 64.271337 MULTILINESTRING ((-1241056.073 188074.166, -12...
330 56.504580 MULTILINESTRING ((-1238926.044 186312.086, -12...

I'm stumped because running a test model from the evaldata returns meaningful results. Any help is appreciated.

[Chris Tessum]

Hi Broderick,

Usually when this happens it is because the emissions are not located within the spatial domain that InMAP is running on. You can check this by opening your emissions input file and the InMAP output file in a GIS program and seeing if they overlap. Also, your attribute column should be named "PM2_5" instead of "PM25" as described here: https://inmap.run/docs/emissions/#emissions-shapefiles

Best regards,
Chris

[Broderik Craig]

Ok, I've changed the name of PM_25 to match the expected format. It looks like the spatial domain isn't overlapping, do you know why and what I can do to fix it?
Thanks!

Broderik

[Chris Tessum]

I think it should be PM2_5

[Broderik Craig]

That's what I meant, that was a typo. How do I specify the spatial domain? is that VariableGridXo, VariableGridYo, VariableGridDx, and VariableGridDy? Do those need to be in the same projection as the .ncf file?

[Chris Tessum]

Yes. But you can only change those if you're creating a new grid to run the model on. If you're using a previously created grid, then you're stuck with the values used to create that grid.

[Broderik Craig]

Thank you so much for your help so far, Chris!

Does InMAP have prebuilt atmospheric data? If so, is it in the evaldata or somewhere else, and if not, how do you recommend obtaining the baseline atmospheric data?

Thank you again!

[Chris Tessum]

Here is the default input data:

https://github.com/spatialmodel/inmap/releases/tag/v1.9.0

For the US: https://zenodo.org/records/3403934