Comments (2)
Inconsistent in number of species appears to be caused by HAPS integration treatment. Depending on SMK_PROCESS_HAPS
setting (e.g., "NONE" or "ALL"), area_map_rwc_2018gg_18j.txt
may content map file for VOC along with VOC_INV, whereas area_rwc.ncf account for VOC.
from smoke.
The issue being isolated to the SETNONHAP subrountine ($SMOKE_HOME/src/smkinven/setnonhap.f)
The following is array of emitted species input to SETNONHAP when processing rwc source.
INVDNAME( 1) = FORMALD
INVDNAME( 5) = BENZENE
INVDNAME( 7) = ACETALD
INVDNAME( 9) = NAPHTH
INVDNAME( 11) = BUTADIE
INVDNAME( 12) = ACROLEI
INVDNAME( 22) = CO
INVDNAME( 79) = NH3
INVDNAME( 83) = NOX
INVDNAME( 88) = PM10
INVDNAME( 89) = PM2_5
INVDNAME(101) = SO2
INVDNAME(107) = VOC
INVDNAME(108) = VOC_INV
When SMK_PROCESS_HAPS
is set to "ALL" or "PARTIAL", a new pollutant NONHAPVOC is calculated as the remainder of subtracting HAPs emission from VOC or TOG emission. Around code line 614 - 620, VOC is renamed to NONHAPVOC, and INVSTAT (which keeps track of which pollutants from the INTABLE are valid from input emission source) is activated for NONHAPVOC.
C..................... Rename VOC to NONHAPVOC
C..................... increment a number of poll for NONHAP[VOC|TOG]
NONVPOS = INDEX1( NONHAPNAM( K ), MXIDAT, INVDNAM )
IPOSCOD ( VPOS ) = NONVPOS
TMPPOSCOD( NVPOS ) = NONVPOS
TMPPOLVAL( NVPOS,: )= POLVAL( VPOS,: )
INVSTAT ( NONVPOS )= 2
Output array of emitted species from SETNONHAP when SMK_PROCESS_HAPS = 'ALL'
INVDNAME( 1) = FORMALD
INVDNAME( 5) = BENZENE
INVDNAME( 7) = ACETALD
INVDNAME( 9) = NAPHTH
INVDNAME( 11) = BUTADIE
INVDNAME( 12) = ACROLEI
INVDNAME( 22) = CO
INVDNAME( 79) = NH3
INVDNAME( 82) = NONHAPVOC
INVDNAME( 83) = NOX
INVDNAME( 88) = PM10
INVDNAME( 89) = PM2_5
INVDNAME(101) = SO2
INVDNAME(108) = VOC_INV
Note how VOC was renamed to NONHAPVOC, and the total number of species remain unchanged at 14. However, INVSTAT for VOC is still activated. This resulted in the total number of tracked pollutant increased by one (15) instead of remain unchanged (14).
Potential solution: deactivate INVSTAT for VOC when SMK_PROCESS_HAPS = 'ALL'
C Add position index of VOC in INVSTAT
INTEGER VPOS_SAVE
then
VPOS_SAVE = IPOSCOD( VPOS ) ! UNC-IE, Feb 2024: Save location of VOC in INVSTAT
NONVPOS = INDEX1( NONHAPNAM( K ), MXIDAT, INVDNAM )
IPOSCOD ( VPOS ) = NONVPOS
TMPPOSCOD( NVPOS ) = NONVPOS
TMPPOLVAL( NVPOS,: )= POLVAL( VPOS,: )
INVSTAT ( NONVPOS )= 2
c UNC-IE: Feb 2024, as VOC is renamed to NONHAPVOC,
c INVSTAT(NVPOS) should be set to zero when SMK_PROCESS_HAPS = 'ALL'
IF ( PROC_HAPS == 'ALL' ) INVSTAT ( VPOS_SAVE ) = 0
from smoke.
Related Issues (20)
- potential updates to CHKGRID subroutine
- review potential updates to Tmpbeis4 code HOT 7
- SMOKE Array Out of Bound HOT 4
- make SPEED inventory optional HOT 2
- Missing information about the SPDIST speed distribution profile option HOT 2
- Incorrect format statement in subroutine RDMRCLIST
- notes on SMOKE Installation Instructions HOT 6
- Ambiguous definition of EFLAG variable in setoutdate HOT 1
- problems with Smkinven daily area source processing HOT 12
- Gentpro issue for areas without Daylight Saving Time HOT 2
- documentation updates for utility programs HOT 1
- Test inline2d and mrgpt program and see if updates are needed for SMOKE 5.0 User guide HOT 2
- SMOKE Utility Programs Geofac: Mismatch in number of variable issue; Documentation needs update HOT 1
- possible DST ending date bug in I/O API HOT 1
- EGU hourly FF10 processing: HAP calculation error HOT 7
- SMOKE (and IOAPI) compilations with various fortran compilers HOT 1
- miscalculated hourly emissions with SMOKE 5.1 HOT 3
- Assignment of meteorological values inconsistent with CMAQ methods HOT 4
- document that stack parameters <= 0 are treated as missing
- Making Movesmrg to apply adjustment factors by state
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from smoke.