# SA setup for the PSAT run for EPA 2016v1 platform questions

**URL:** <https://forum.cmascenter.org/t/sa-setup-for-the-psat-run-for-epa-2016v1-platform-questions/2295>\
**Category:** Emissions Inventory\
**Tags:** admin\
**Created:** [April 27, 2021, 4:27pm UTC](https://forum.cmascenter.org/t/sa-setup-for-the-psat-run-for-epa-2016v1-platform-questions/2295 "2021-04-27T16:27:22Z")\
**Posts on this page:** 1\
**Page:** 1

<div class="post-metadata">

**Author:** ![Iqbal](https://avatars.discourse-cdn.com/v4/letter/i/bc8723/32.png) [@Iqbal](https://forum.cmascenter.org/u/Iqbal)\
**Post date:** [April 27, 2021, 4:27pm UTC](https://forum.cmascenter.org/t/sa-setup-for-the-psat-run-for-epa-2016v1-platform-questions/2295/1 "2021-04-27T16:27:22Z")

</div>

Hello,  
I would be grateful if anyone helps me with the following questions:

I have successfully simulated EPA 2016v1 runs by using CAMx. Now I want to run 2028 with PSAT. I am planning to tag 10 EGU facilities with SO4 and NO3 species tracer.

1. I am running two SMOKE runs: scenario-1 will create one file(Let’s say FILEA) for ptegu sector which is ( ALL\_EGU – 10 tag facilities); scenario-2 will create one file (FILEB) for ptegu sector which contains only ’ 10 tag facilities’
2. I am trying to generate a ’ CAMx control’ file for the SA PSAT.Please see below (_use of FILEA & FILEB in bold color, Is that how I should use these files?_):

!— Model options —

Diagnostic\_Error\_Check = .false., ! True = will stop after 1st timestep  
Flexi\_Nest = .false., ! True = expect flexi-nested inputs  
Advection\_Solver = ‘PPM’, ! (PPM,BOTT)  
Chemistry\_Solver = ‘EBI’, ! (CMC,IEH,LSODE)  
PiG\_Submodel = ‘NONE’, ! (None,GREASD,IRON)  
**Probing\_Tool = ‘SA’, ! (None,SA,DDM,HDDM,PA,IPR,IRR,RTRAC,RTCMC)**  
Chemistry = .true.,  
Drydep\_Model = ‘ZHANG03’, ! (NONE,WESELY89,ZHANG03)  
Bidi\_NH3\_Drydep = .false.,  
Wet\_Deposition = .true.,  
ACM2\_Diffusion = .false.,  
Surface\_Model = .false.,  
Inline\_Ix\_Emissions = .true.,  
Super\_Stepping = .true.,  
Gridded\_Emissions = .true.,  
Point\_Emissions = .true.,  
Ignore\_Emission\_Dates = .true.,

!— Output specifications —

Root\_Output\_Name = ‘$OUT/$today.$CGRID.$NZ.$CCASE.epa.camx’,  
Average\_Output\_3D = .false.,  
NetCDF\_Format\_Output = .false.,  
NetCDF\_Use\_Compression = .false.,  
Output\_3D\_Grid(1) = .false.,  
Output\_Species\_Names(1) = ‘ALL’,

PiG\_Sampling\_Grid = .false.,  
Sample\_BAckground = .false.,  
Number\_of\_Sampling\_Grids = 1,  
SG\_Beg\_I\_Index(1) = 17,  
SG\_End\_I\_Index(1) = 18,  
SG\_Beg\_J\_Index(1) = 21,  
SG\_End\_J\_Index(1) = 22,  
SG\_Mesh\_Factor(1) = 36.,

!— Input files —

Chemistry\_Parameters = ‘$CHEM/CAMx7.0.chemparam.CB6r4DMS\_CF\_SOAP\_ISORROPIA\_elements\_NH3RSCALE0’,  
Photolyis\_Rates = ‘RATE/rate.{YYYY}${MM}$DD.$CGRID.do.CB6.txt’,  
Ozone\_Column = ‘$OZCOL/ozcol.$YYYY$MM$DD.$CGRID.txt’,  
Initial\_Conditions = ‘$ICBC/init.$YYYY$MM$DD.CGRID.{NZ}L.$ICCASE.epa.camx’,  
Boundary\_Conditions = ‘$ICBC/bndr.$YYYY$MM$DD.CGRID.{NZ}L.$ICCASE.epa.camx’,  
Top\_Concentrations = ’ ',

**Point\_Sources = ‘$PTSRDIR/FILEA’,**  
Master\_Grid\_Restart = ‘$OUT/$YESTERDAY.$CGRID.$NZ.$CCASE.epa.camx.inst’,  
Nested\_Grid\_Restart = ‘$OUT/$YESTERDAY.$CGRID.$NZ.$CCASE.epa.camx.finst’,  
PiG\_Restart = ‘$OUT/$YESTERDAY.$CGRID.$NZ.$CCASE.epa.camx.pig’,  
Srfmod\_Grid(1) = ’ ',  
Srfmod\_Grid(2) = ’ ',

Surface\_Grid(1) = ‘$MET/landuse.151221.36US3.zhang.camx.salt’,  
Surface\_Grid(2) = ‘$MET2/landuse.151221.12US2b.zhang.camx.salt’,  
Met3D\_Grid(1) = ‘$MET/met3d.$YY$MM$DD.$CGRID.$NZ.wrf.camx’,  
Met3D\_Grid(2) = ‘$MET2/met3d.$YY$MM$DD.12US2b.$NZ.wrf.camx’,  
Met2D\_Grid(1) = ‘$MET/met2d.$YY$MM$DD.$CGRID.$NZ.wrf.camx’,  
Met2D\_Grid(2) = ‘$MET2/met2d.$YY$MM$DD.12US2b.$NZ.wrf.camx’,  
Vdiff\_Grid(1) = ‘$MET/vdif.$YY$MM$DD.$CGRID.$NZ.wrf.camx.0.1.YSU.kvpatch’,  
Vdiff\_Grid(2) = ‘$MET2/vdif.$YY$MM$DD.12US2b.$NZ.wrf.camx.0.1.YSU.kvpatch’,  
Cloud\_Grid(1) = ‘$MET/clra.$YY$MM$DD.$CGRID.$NZ.wrf.camx’,  
Cloud\_Grid(2) = ‘$MET2/clra.$YY$MM$DD.12US2b.$NZ.wrf.camx’,  
Emiss\_Grid(1,1) = ‘EMIS/emis2d.{YYYY}{MM}{DD}.{CGRID}.{CCASE2}+ss.camx’,  
Emiss\_Grid(2,1) = ‘EMIS2/emis2d.{YYYY}{MM}{DD}.12US2b.${CCASE2}+ss.camx’,

&

!-------------------------------------------------------------------------------  
&SA\_Control

SA\_File\_Root = ‘$OUT/$today.$CGRID.$NZ.$CCASE.epa.camx’,  
SA\_Summary\_Output = .true.,

SA\_Stratify\_Boundary = .false.,  
SA\_Deposition\_Output = .false.,  
**SA\_Number\_of\_Source\_Regions = 54,**  
**SA\_Number\_of\_Source\_Groups = 1,**  
Use\_Leftover\_Group = .false.,  
Number\_of\_Timing\_Releases = 0,  
SA\_Treat\_SULFATE\_Class = .true.,  
SA\_Treat\_NITRATE\_Class = .true.,  
SA\_Treat\_SOA\_Class = .false.,  
SA\_Treat\_PRIMARY\_Class = .false.,  
SA\_Treat\_MERCURY\_Class = .false.,  
SA\_Treat\_OZONE\_Class = .false.,  
SA\_Use\_APCA = .false.,  
SA\_Use\_APCA\_Ptoverride = .false.,  
SA\_PT\_Override = .false.,

SA\_Receptor\_Definitions = ‘’,  
**SA\_Source\_Area\_Map(1) = ‘$LAND/satmap.$CGRID.epa.txt’,**  
SA\_Source\_Area\_Map(2) = ‘’,  
SA\_Use\_Partial\_SourceMap = .false.,  
Partial\_Source\_Area\_Map(1,1) = ’ ', ! Map for SA group 1, grid 1  
Partial\_Source\_Area\_Map(2,1) = ’ ', ! Map for SA group 2, grid 1

SA\_Master\_Restart = ‘$OUT/$YESTERDAY.$CGRID.$NZ.$CCASE.epa.camx.sa.inst’,  
SA\_Nested\_Restart = ‘$OUT/$YESTERDAY.$CGRID.$NZ.$CCASE.epa.camx.sa.finst’,

**SA\_Points\_Group(1) = ‘$PTSRDIR/FILEB’,**  
SA\_Points\_Group(2) = ‘’,

SA\_Emiss\_Group\_Grid(1,1) = ‘’,  
SA\_Emiss\_Group\_Grid(1,2) = ‘’,

_My questions are:_

_1. What is the value / How do I determine the value for ’ SA\_Number\_of\_Source\_Regions =54’ for the EPA 2016v1 platform for this case?._  
_2. How do I generate this file ‘SA\_Source\_Area\_Map(1) = ‘$LAND/satmap.$CGRID.epa.txt’,’ ?_

Any help and suggestions would be greatly appreciated

Thank you