Possibility of adding layers to previously defined layers?

@cjcoats @bbaek

Hi,

I already have processed emission files, but with different vertical grid definitions, and I need to add layers to the predefined one. ( Processed data has 11 layer, and I need to change it to 30 by adding extra layers, both has same vertical height). Is there any way I can do this with IOAPI ?

Thank you,

Not quite sure what you are exactly looking for but there is a tool that I developed called Layalloc which is one of SMOKE utility tool that allows you to allocate the 2d emissions into vertical layers. You can define your own layer profiles using a simple CSV input file. You can find the detail from the link:

https://www.cmascenter.org/smoke/documentation/4.7/html/ch05s03s09.html

Hope this meets your needs.

Hi,

I really do appreciate your quick response and your help. I am sorry that I did not mentioned that clear.
I have SMOKE output files with vertical grid definition as follows:
VGTOP = 5000.f ;
VGLVLS = 1.f, 0.9958f, 0.9907f, 0.9846f, 0.9774f, 0.9688f, 0.9585f, 0.9463f, 0.9319f, 0.9148f, 0.8946f, 0.f

And my MCIP vertical definition is as follows:
:VGTOP = 5000.f

VGLVLS = 1.f, 0.9958f, 0.9907f, 0.9846f, 0.9774f, 0.9688f, 0.9585f, 0.9463f, 0.9319f, 0.9148f, 0.8946f, 0.8709f, 0.8431f, 0.8107f, 0.7733f, 0.7301f, 0.6809f, 0.6254f, 0.5635f, 0.4954f, 0.4218f, 0.3536f, 0.293f, 0.2394f, 0.1917f, 0.1495f, 0.112f, 0.0788f, 0.0493f, 0.0232f, 0.f

And I need to modify my SMOKE output (not processed by me) in a way that I can use it with my MCIP output. Is there any way to do this? ( I think extracting 30 layers out of 11 layers is possible, is that right?)

Thank you again,

There’s nothing “canned” that will do what you ask. You can get a reasonable approximation by converting the 11-layer SMOKE units from mass per 3-D grid cell per unit time to mass per cubic meter per unit time, describing both grids in pressure-coordinates interpolating that to the new grid in pressure-coordinates, and then re-converting back to mass per (now, smaller) 3-D grid cell per unit time. Not a trivial piece of work…

1 Like

Alright. Thank you for your response.