Errors when running Spatial Allocator

Spatial Allocator
1

Dear Sir or Madam,

I would like to include sea-salt emissions in my CMAQ model. According to the official documentation, this requires generating an OCEAN file using the Spatial Allocator software. However, I have encountered a series of issues while running the program.

First, the input Shapefile I used is a global ocean shapefile (in .shp format) to which I added a SURF_ZONE field. I will also provide the information from the corresponding .prj file. The GRIDDESC.txt file was generated by MCIP and contains the grid information for CMAQ. The two files have different spatial resolutions.

When I run the program, it crashes and stops automatically. I will attach the run script and other related files. Could you please advise me on how to fix these problems? In addition, I would like to know if there are other ways to include sea-salt emission data in CMAQ.

Thank you very much for your help.

Best regards,

PROJCS[“LCC”,GEOGCS[“GCS_WGS_1984”,DATUM[“D_WGS_1984”,SPHEROID[“WGS_1984”,6378137.0,298.257223563]],PRIMEM[“Greenwich”,0.0],UNIT[“Degree”,0.0174532925199433]],PROJECTION[“Lambert_Conformal_Conic”],PARAMETER[“False_Easting”,0.0],PARAMETER[“False_Northing”,0.0],PARAMETER[“Central_Meridian”,101.738],PARAMETER[“Standard_Parallel_1”,30.0],PARAMETER[“Standard_Parallel_2”,60.0],PARAMETER[“Scale_Factor”,1.0],PARAMETER[“Latitude_Of_Origin”,36.493],UNIT[“Meter”,1.0]]

alloc_surf_zone_to_oceanfile.csh (1.5 KB)

GRIDDESC.txt (197 Bytes)

e7aa2e1bbb9506c2566e93086158bc9a
e7aa2e1bbb9506c2566e93086158bc9a1500×987 232 KB

e2ed82f59e22f612141f0107171da89a
e2ed82f59e22f612141f0107171da89a1494×978 237 KB

e99e709c30a3d90c610d6d7bd37a1686
e99e709c30a3d90c610d6d7bd37a16861665×1053 44.6 KB

2

Your warning specifically points to the ALLOC_ATTR_TYPE not being set prior to running the Spatial Allocator. Please verify that this gets set (e.g. echo back).

I can’t tell from the image if the grid cells are part of your shapefile. The surf zone shapefile should only include land area polygons that extend to the ocean shores. Are you able to run the Spatial Allocator with the provided example surf zone shapefile described in the CMAQ documentation (surfzone_poly_st.shp)?

3

Sir, thank you very much for your reply.

After reviewing the documentation again, I realized that the shapefile does not need to contain all ocean areas; it only needs to include the near-shore coastal zone. I re-generated a new shapefile in ArcGIS, created a field named SURF_ZONE, and assigned all values to 1. I also modified the script accordingly, and the output file was successfully generated.

I would like to ask you to help me interpret the run log.
In addition, I would like to know what characteristics the output file should have so that I can confirm it is correct.

Finally, I would like to ask: if I do not plan to simulate DMS and only need sea-salt emissions, is this output file already sufficient to be used by CMAQ?

Thank you, and I look forward to your guidance.

alloc_surf_zone_to_oceanfile.csh (1.4 KB)

log.txt (5.9 KB)

4

When I input this file into CMAQ, it tells me that the variable OPEN is missing. How should I generate the OPEN variable? Thank you very much.

 "OCEAN_1" opened as OLD:READ-ONLY   
 File name "/home/guoxh/SpatialAllocatorMaster/Spatial-Allocator-master/output//ocean_file_202301.ncf"
 File type GRDDED3 
 Execution ID "????????????????"
 Grid name "202301"
 Dimensions: 161 rows, 183 cols, 1 lays, 1 vbles
 NetCDF ID:   1048576  opened as READONLY            
 Time-independent data.
 
 >>--->> WARNING in subroutine XTRACT3
 Requested variable OPENnot available in file OCEAN_1which contains variables:
 
 SURF_ZONE_1


 *** ERROR ABORT in subroutine retrieve_ocean_d on PE 021
 Could not read OPEN from OCEAN_1
5

Could you possibly tell me the location and names of the shapefiles required to generate the OCEAN file in the official documentation? Thank you very much.