So here we go. The code can be found here It's a linux version but should run on a PC OK. Make sure your map dataset "maps.states" is uncompressed. I run it from the directory $EMS_BIN/cleaner/ptclean. You'll need a ems_run.ptemis data set and the file ems_grd/project_xy.pt ascii file. These datasets should be created with a succesful run of the point source model. If you do not have one try running uni_point.mej as discussed in the What's new in ems-2001 document.
To run it make sure all your environment variables are set(like "source uni_980719.mej). The environment variables you want to make sure to set are $EMSDAY, $CASE, $EMS_RUN, and $EMS_GRD. To run the code type the following command:
sas -sysparm "NOX" circle.sas
Where NOX is the pollutant you want to map. You can also specific either PS for postscript or GIF for gif image types. If your printing to a postscript printer PS is the best option. This is done by using the command :
sas -sysparm "NOX;PS" circle.sas
or
sas -sysparm "SO2;GIF" circle.sas
IF you do not put the second variable the program will assume GIF.
I have documented the code with lines that begin with "PERSONALIZE" where there are options the user might want to changed for their area. The most common that you might want to change include.
On line 6 and 7 you can change the location on the map of the legend, This is based on a lattitude and longitude of the point.
On line 11-16 you can change the scale of the legend.
On line 20-30 you can add new device types like JPG or screen if your feeling brave
On line 36 you can change fonts and colors of the output plots.
On line 225 you can change the list of states that will be included in the map. here is also where you might add canada or mexico.
On line 265 you can change all the title and foot note information.
Circle Plots have been used for some time to show the largest point sources. The problem with these older plotting techniques is that they have relied on photochemical modelers to generate the plots with complex Fortran code that had little chance of user modificiation to local needs. The techniques used had limited flexibility in generating specific variations in the plots like different countries, projections, or color and shading techniques. One of the biggest drawbacks of past methods was the dynamic re-calculation of circle size and legend cutpoints relative to the largest source in the domain. When control strategy modeling used these older techniques and the legend would change and it became difficult to compare base year runs with future year runs. Given these past limitiation LADCO wanted a new method with the following attributes:
1. Written in SAS so it would be flexible
2. Works for different pollutants(NOX, SOX, ROG, ETC...)
3. Scale is both dynamic and definable.
The results are: