Bugs for EPA_ANTHRO_EMIS

[Zongrun Li]

Dear all,

I recently compared the EPA_ANTHRO_EMIS processed inline stack emission with SMOKE inputs. I found there are some potential issues:

1. In anthro_emis.f90, instead of overwriting the wrk_emis with different stack emission files, should it be accumulation?

Before:
      wrk_emis(:,:,:) = 0.
      do stk = 1,data_file%stack%nStk
        if( data_file%stack%dataMask(stk) ) then
          il = data_file%stack%mdl_i(stk) ; jl = data_file%stack%mdl_j(stk)
          k = data_file%stack%mdl_k(stk)
          wrk_emis(il,jl,k) = data_file%stack%src_data(stk)   ! overwrites!
        endif
      enddo

Fixed:
      wrk_emis(:,:,:) = 0.
      do stk = 1,data_file%stack%nStk
        if( data_file%stack%dataMask(stk) ) then
          il = data_file%stack%mdl_i(stk) ; jl = data_file%stack%mdl_j(stk)
          k = data_file%stack%mdl_k(stk)
          wrk_emis(il,jl,k) = wrk_emis(il,jl,k) + data_file%stack%src_data(stk)   ! accumulate!
        endif
      enddo

2. In stack.f90 file, stkFileInit subroutine. The utility uses the stack height from SMOKE outputs. However, in CMAQ, Briggs is applied to calculate the effective plume heights, which considers not only stack height but also exit velocity, temperature, etc.

Thanks,

Zongrun