# CMAQ 5.3.1 IPR mass closure problem **URL:** \ **Category:** CMAQ\ **Created:** [October 24, 2020, 7:24am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863 "2020-10-24T07:24:02Z")\ **Posts on this page:** 20\ **Page:** 1
**Author:** ![Zhizhua](https://avatars.discourse-cdn.com/v4/letter/z/90db22/32.png) [@Zhizhua](https://forum.cmascenter.org/u/Zhizhua)\ **Post date:** [October 24, 2020, 7:24am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/1 "2020-10-24T07:24:02Z")
Hi all, When activating process analysis, I find that the overall concentration change of particle nitrate (hourly difference of data from CONC file) does not equal to the sum of the contributions of the 12 processes (CHEM+EMIS+…). I check that the equality does hold for ozone concentration, where AERO (NPF + COND + COAG + GROW) contributes 0. I also check OH and HNO3 (g), the former AERO contributes 0 while the latter not, and the former holds equality while the latter not. Then I guess the problem may come from PA of AERO, and I find that CMAQ 5.3 updated PA for aerosol subprocesses. Is this problem a bug or due to my mistakes? Does anyone have any suggestion? ---
**Author:** ![cgnolte](https://yyz1.discourse-cdn.com/flex027/user_avatar/forum.cmascenter.org/cgnolte/32/53_2.png) [@cgnolte](https://forum.cmascenter.org/u/cgnolte)\ **Post date:** [October 25, 2020, 2:47am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/2 "2020-10-25T02:47:13Z")
Hi Zhizhua, It sounds like you have done the right sort of tests to demonstrate a problem. How large is the mass closure issue in a given time step, on an absolute and relative basis? Are the errors at the surface, or aloft? ---
**Author:** ![Zhizhua](https://avatars.discourse-cdn.com/v4/letter/z/90db22/32.png) [@Zhizhua](https://forum.cmascenter.org/u/Zhizhua)\ **Post date:** [October 25, 2020, 6:18am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/3 "2020-10-25T06:18:17Z")
Thank you for you reply, cgnolte. The concentration difference is very large for ANO3. I only calculate one day, and the absolute difference is from -19 to 20 ug/m3/hr, and relative difference is -1590% to 1786%. I also calculate “inert” species (Na and EC), the concentration differences are negligible since the contributions of AERO are negligible (not 0). The error is at the surface. I have not checked aloft yet. ---
**Author:** ![Ben\_Murphy](https://avatars.discourse-cdn.com/v4/letter/b/74df32/32.png) [@Ben\_Murphy](https://forum.cmascenter.org/u/Ben_Murphy)\ **Post date:** [October 26, 2020, 3:48pm UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/4 "2020-10-26T15:48:00Z")
Hi Zhizhua, Thanks for your description and sorry you are seeing these issues. May I ask you a few more clarifying questions to think more deeply about this problem? - What domain are you running your simulation for? - I’m wondering what you see for a pollutant like ASO4 which will have condensation but no evaporation? - Are you summing multiple aerosol modes for these values? That shouldn’t be a problem but useful to know when thinking these things through. The coagulation contribution should also be 0 for the sum of i and j mode variables of a particular species so that’s one additional thing to check. Best regards, Ben ---
**Author:** ![Zhizhua](https://avatars.discourse-cdn.com/v4/letter/z/90db22/32.png) [@Zhizhua](https://forum.cmascenter.org/u/Zhizhua)\ **Post date:** [October 27, 2020, 2:44am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/5 "2020-10-27T02:44:14Z")
Hi Ben, Thanks for your reply. The following are answers to your questions, and I hope these could make the problem clearer. 1. The domain is east and southeast Asia, but I only chose an area about 70,000 km2 in China (highly polluted area) to calculate. 2. The differences of ASO4 are also very large (data for another day: 0.014~1.014 ug/m3/hr, -623%~1493%; avg. 0.324 ug/m3/hr, 85%). By the way, only COND has non-negligible contribution for both ASO4 and ANO3. 3. Yes. I summed i and j modes. However, the contribution of COAG is not 0, both for ANO3 and ASO4 (though very small, of 10^-7 order of the net effect, while the contribution of NPF is exact 0 for ANO3). ---
**Author:** ![Ben\_Murphy](https://avatars.discourse-cdn.com/v4/letter/b/74df32/32.png) [@Ben\_Murphy](https://forum.cmascenter.org/u/Ben_Murphy)\ **Post date:** [October 28, 2020, 11:20am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/6 "2020-10-28T11:20:44Z")
Hi Zhizhua, Thank you for your answers. Unfortunately, this does not sound like a problem with an easy fix, and I’m uncertain why it is happening on your domain and not on ours. Can you please send me your process analysis input file and the output file as well? I will direct message you with my email address. Ben ---
**Author:** ![Ben\_Murphy](https://avatars.discourse-cdn.com/v4/letter/b/74df32/32.png) [@Ben\_Murphy](https://forum.cmascenter.org/u/Ben_Murphy)\ **Post date:** [November 21, 2020, 1:04pm UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/7 "2020-11-21T13:04:42Z")
We have successfully diagnosed this problem and are investigating a potential solution. Once it has cleared our internal code review process, we will post the update. The bottom line is that aerosol processes are being double-counted for the process analysis fields for all time steps beyond the first one of each output time step. This error does not impact the results of the core model – the predicted concentrations are not corrupted by this error. -Ben ---
**Author:** ![bmbmnb](https://avatars.discourse-cdn.com/v4/letter/b/94ad74/32.png) [@bmbmnb](https://forum.cmascenter.org/u/bmbmnb)\ **Post date:** [March 15, 2021, 3:40am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/8 "2021-03-15T03:40:47Z")
Hi Ben, I want to know if the problem has been solved for CMAQ 5.3.2. Best regards, Yunqing ---
**Author:** ![Ben\_Murphy](https://avatars.discourse-cdn.com/v4/letter/b/74df32/32.png) [@Ben\_Murphy](https://forum.cmascenter.org/u/Ben_Murphy)\ **Post date:** [March 15, 2021, 1:01pm UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/9 "2021-03-15T13:01:58Z")
Hi Yunqing, We have found a solution but it is still under review since it is lumped with a few other updates. I have tried to distill just the pertinent bug fix in the update to the file uploaded here. Please remove the .txt extension and add it to the CCTM/src/procan/pa folder. I have not tested it, so please let me know if it fails to compile or run. [pa\_update.F.txt](https://forum.cmascenter.org/uploads/short-url/gV95NUAC23aVKFA3VpjprhQ00De.txt) (30.6 KB) Please also modify the call to pa\_update\_aero in scriproc.F to the following: IF ( LIPR ) CALL PA\_UPDATE\_AERO ( CGRID, JDATE, JTIME ) Best wishes, Ben ---
**Author:** ![bmbmnb](https://avatars.discourse-cdn.com/v4/letter/b/94ad74/32.png) [@bmbmnb](https://forum.cmascenter.org/u/bmbmnb)\ **Post date:** [March 16, 2021, 1:45am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/10 "2021-03-16T01:45:32Z")
Hi Ben, Thanks for your help, I’ll give it a try. Best wishes, Yunqing ---
**Author:** ![bmbmnb](https://avatars.discourse-cdn.com/v4/letter/b/94ad74/32.png) [@bmbmnb](https://forum.cmascenter.org/u/bmbmnb)\ **Post date:** [March 16, 2021, 2:49am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/11 "2021-03-16T02:49:01Z")
Hi Ben, I tried to compile CMAQ 5.3.2 with the pa\_update.F file you gave to me, but it failed with some errors. (I can compile original CMAQ 5.3.2 successfully with the exact same configuration.) Attached is my build log.[build.txt](https://forum.cmascenter.org/uploads/short-url/42ghAM0fHebgP7EM2NGwvAKNTnV.txt) (338.3 KB) Yunqing ---
**Author:** ![Ben\_Murphy](https://avatars.discourse-cdn.com/v4/letter/b/74df32/32.png) [@Ben\_Murphy](https://forum.cmascenter.org/u/Ben_Murphy)\ **Post date:** [March 16, 2021, 11:56am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/12 "2021-03-16T11:56:05Z")
Hi Yunqing, My mistake, those files were not ready. I’ve compiled and tested these. Please see if they fix your issue. Cheers, Ben[pa\_update.F.txt](https://forum.cmascenter.org/uploads/short-url/cAKt79Asp02PaboLaQiggrIFnHw.txt) (31.3 KB) [sciproc.F.txt](https://forum.cmascenter.org/uploads/short-url/5FUWighJ67ytV8TP9Ac3ZKWC75N.txt) (14.0 KB) ---
**Author:** ![bmbmnb](https://avatars.discourse-cdn.com/v4/letter/b/94ad74/32.png) [@bmbmnb](https://forum.cmascenter.org/u/bmbmnb)\ **Post date:** [March 22, 2021, 3:48am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/13 "2021-03-22T03:48:53Z")
Hi Ben, Sorry for the late reply. I’ve been trying out the \*.F files you gave me. I can compile these files successfully, but there are still some problems with PA output. I checked nitrate and CO for the output, the same problem still exist(mass closure problem with AERO). Here are my pa.ctl, pa\_report and output.[mass\_closure.csv](https://forum.cmascenter.org/uploads/short-url/uMQk7F0W5lQuE4HX5t9AZErVKdK.csv) (51.6 KB) [pa\_cb6r3\_ae7\_aq.ctl.txt](https://forum.cmascenter.org/uploads/short-url/dyKuGukjijf6aSDOzbOPmNFp9ki.txt) (1.8 KB) [PA\_REPORT.20201210.txt](https://forum.cmascenter.org/uploads/short-url/fenRrPcSfF2eSOEShZ1WHTVL6na.txt) (9.2 KB) I ran the model for 96 hours and extract the output for two grids((1,1) & (7,14) for PA domain). Yunqing ---
**Author:** ![Zhizhua](https://avatars.discourse-cdn.com/v4/letter/z/90db22/32.png) [@Zhizhua](https://forum.cmascenter.org/u/Zhizhua)\ **Post date:** [March 25, 2021, 6:47am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/14 "2021-03-25T06:47:21Z")
Hi Yunqing, Ben gave me three files to replace the old ones, and I tested them successfully. Hope they are helpful to you. [sciproc.F.txt](https://forum.cmascenter.org/uploads/short-url/5bleAS7DN2X11FBqDpuFdx0x2mk.txt) (14.0 KB) [pa\_update.F.txt](https://forum.cmascenter.org/uploads/short-url/lbuxMk4yL5KCGwsGZEGLc0JzKFD.txt) (32.9 KB) [aero\_subs.F.txt](https://forum.cmascenter.org/uploads/short-url/cilLLz6uxnyHZxeyQDqwGrptWv4.txt) (136.0 KB) Zhizhua ---
**Author:** ![Ben\_Murphy](https://avatars.discourse-cdn.com/v4/letter/b/74df32/32.png) [@Ben\_Murphy](https://forum.cmascenter.org/u/Ben_Murphy)\ **Post date:** [March 25, 2021, 1:53pm UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/15 "2021-03-25T13:53:40Z")
Thanks for sending those Zhizhua. I had lost track of our exchange last fall. @Yunqing: please see if those files from Zhizhua help your case. I tried them myself and it didn’t change any of my results. It’s possible that the problems you are seeing are a result of the inherent flaw in using concentration-based process analysis (rather than a mass-based approach) to try to quantify mass closure. However, the extreme deviation that you see for ANO3J makes me think otherwise. For grid cell 7,14 you report pretty massive gains by PA that are unsupported by the concentration change in the early hours of the simulaiton. Can you determine which process(es) are dominating that increase? Best, Ben ---
**Author:** ![bmbmnb](https://avatars.discourse-cdn.com/v4/letter/b/94ad74/32.png) [@bmbmnb](https://forum.cmascenter.org/u/bmbmnb)\ **Post date:** [March 26, 2021, 1:30am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/16 "2021-03-26T01:30:46Z")
Thank you for sending me those files Zhizhua. I’ll give it a try. Yunqing ---
**Author:** ![bmbmnb](https://avatars.discourse-cdn.com/v4/letter/b/94ad74/32.png) [@bmbmnb](https://forum.cmascenter.org/u/bmbmnb)\ **Post date:** [March 26, 2021, 1:36am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/17 "2021-03-26T01:36:56Z")
Thanks for your reply Ben. The difference was mainly caused by process TRNM.[mass\_closure.csv](https://forum.cmascenter.org/uploads/short-url/2FiCDP4wYakMaZogAdCFqSdAKT1.csv) (98.2 KB) Best Yunqing ---
**Author:** ![Ben\_Murphy](https://avatars.discourse-cdn.com/v4/letter/b/74df32/32.png) [@Ben\_Murphy](https://forum.cmascenter.org/u/Ben_Murphy)\ **Post date:** [March 26, 2021, 4:09am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/18 "2021-03-26T04:09:48Z")
Can you tell me what advection option you are using? Is it ‘local’ or ‘wrf’? Also, if it’s not too much trouble, it might be useful to rerun with some of the sub-processes within TRNM separated out (e.g. VDIF, HDIF, HADV, and ZADV) to see which is contributing the apparent error. -Ben ---
**Author:** ![bmbmnb](https://avatars.discourse-cdn.com/v4/letter/b/94ad74/32.png) [@bmbmnb](https://forum.cmascenter.org/u/bmbmnb)\ **Post date:** [March 29, 2021, 2:18am UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/19 "2021-03-29T02:18:53Z")
Hi Ben, I’m using ‘wrf’ advection option and wrf meteorological output. I rerun the model (with files you gave me), it shows that ‘VDIF’ contribute the most.[mass\_closure.csv](https://forum.cmascenter.org/uploads/short-url/svwlBGfz5lOFHJfZoWQAJJRn70r.csv) (134.6 KB) Yunqing ---
**Author:** ![Ben\_Murphy](https://avatars.discourse-cdn.com/v4/letter/b/74df32/32.png) [@Ben\_Murphy](https://forum.cmascenter.org/u/Ben_Murphy)\ **Post date:** [April 8, 2021, 9:18pm UTC](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863/20 "2021-04-08T21:18:36Z")
Hi Yunqing, Just so you’re aware, I’m still working on this issue on my side in between other projects. I have some new updates I’ve made that will hopefully alleviate the issue you’re seeing. However, I’m still testing them. When they are ready, would you consider running some tests yourself to see if they resolve your issue? Best wishes, Ben [Next page](https://forum.cmascenter.org/t/cmaq-5-3-1-ipr-mass-closure-problem/1863.md?page=2)