Setting up a new domain in SMOKE

#1

First time SMOKE user so apologies for the basics - I’m using the EPA’s 2016 alpha modeling platform to set up emissions over a few domains and nests over the U.S. but I’m running into problems setting it up over my domain. I’m starting with my outer domain, which is a 12km over the US and parts of Canada and Mexico. The error I’m receiving is:

 ERROR: Columns, rows, x-origin, or y-origin for US12KM_474X336   in
      surrogates file are inconsistent with initialized values.

My GRIDDESC is:

’ ’
‘LAM_040N97W’
2 33.000 45.000 -97.000 -97.000 39.000
’ ’
‘US12_1’
‘LAM_040N97W’ -2615999.500 -1677427.500 12000.000 12000.000 469 309 1

My question is, do I need to generate my own spatial surrogates for each domain if it is not exactly the one defined in the surrogates provided? I’ve looked at the spatial allocator generate_surrogates tool but it does not seem consistent with the shapefiles that are provided with the EPA platform. Any help is appreciated - Happy to provide additional details. Thanks.

#2

First, the map projection of modeling domains have be same, and then if the surrogates cover a bigger region than your domain, SMOKE can zoom in for you. However, if your domain is not covered by these surrogates, then you need to generate your own surrogates using the spatial allocator.

#3

Thank you @bbaek - My domain uses the same projection as the surrogates (LAM_40N97W) so that should not be the issue. I’ve also subsetted my WRF domain in MCIP to match one of the domains in the EPA-platform provided GRIDDESC: the file is named: griddesc_lambertonly_08dec2017_v0.txt and my domain matches: ‘12US1_459X299’. However, I’m still getting this error in GRDMAT:

 ERROR: Columns, rows, x-origin, or y-origin for US12KM_474X336   in
      surrogates file are inconsistent with initialized values.
 ERROR: Grid parameters for US12KM_474X336   in
      surrogates file are inconsistent with initialized values.

What am I doing wrong?

#4

@bbaek - I should say, here is the grid defined in teh surrogates:

#GRID CONUS12k_444x336 -2736000.000000 -2088000.000000 12000.000000 12000.000000 444 336 1 LAMBERT meters 33.000000 45.000000 -97.000000 -97.000000 40.000000

And here is my GRIDDESC:

’ ’
‘LAM_40N97W’
2 33.000 45.000 -97.000 -97.000 40.000
’ ’
‘US12_1’
‘LAM_40N97W’ -2556000.000 -1728000.000 12000.000 12000.000 459 299 1
’ ’

#5

A smaller domain yields a similar error:

‘LAM_40N97W’, -2436000.000, -1608000.000, 12000.000, 12000.000, 439, 279, 1

#6

Your defined grid is larger in the X dimension than the surrogates: yours is 459 299 and the surrogates are 444x236

#7

@eyth.alison - thank you, just tried with the above domain (439x279) and I am getting the same error, though the XORIG and YORIG are different

#8

Everything has to align or the grids may not fully overlap