Incorrect PM25_NA calculation in the ELMO module

CMAQ
2

Thanks for pointing this out Bonyoung. You are correct. Try modifying the content of the ID_PM25_NA case in ELMO_PROC.F to:

CASE ( ID_PM25_NA )
     N_VARS = 4
     CALL INIT_ELMO_COEFFS( IDG, N_VARS, 'PM25' )
     CALL SET_ELMO_COEFF_MAP( IDG, 1, 'AERO', 'ANA', 1.0 )
     CALL SET_ELMO_COEFF_MAP( IDG, 2, 'AERO', 'ASEACAT', 0.8373 )
     CALL SET_ELMO_COEFF_MAP( IDG, 3, 'AERO', 'ASOIL', 0.0626 )
     CALL SET_ELMO_COEFF_MAP( IDG, 4, 'AERO', 'ACORS', 0.0023 )

Note that the analogous update should be made for MG, K and Ca. And you are right that PM1 and PM25TO10 should be updated, although PM1 should not be affected too much of course. I apologize for the oversight. This will be corrected in the next release. In the meantime, I will issue an update to the public repo.

Best wishes,
Ben

6 Likes