**************************************************************
** Runscript Detected an Error: CGRID file was not written. **
** This indicates that CMAQ was interrupted or an issue **
** exists with writing output. The runscript will now **
** abort rather than proceeding to subsequent days. **
**************************************************************
Processing Day/Time [YYYYDDD:HHMMSS]: 2020132:164320
Which is Equivalent to (UTC): 16:43:20 Monday, May 11, 2020
Time-Step Length (HHMMSS): 000140
After VDIFF : Gas 4.844E-01 | Aer 2.338E+10 | Non 4.978E-01
After COUPLE : Gas 4.585E+03 | Aer 2.563E+14 | Non 4.657E+03
After HADV : Gas 4.555E+03 | Aer 2.546E+14 | Non 4.604E+03
After ADV : Gas 4.558E+03 | Aer 2.546E+14 | Non 4.605E+03
After HDIFF : Gas 4.558E+03 | Aer 2.546E+14 | Non 4.605E+03
After DECOUPLE : Gas 4.817E-01 | Aer 2.315E+10 | Non 4.924E-01
After PHOT : Gas 4.817E-01 | Aer 2.315E+10 | Non 4.924E-01
*** ERROR ABORT in subroutine AQCHEM on PE 000
Maximum AQCHEM total iterations exceeded
PM3EXIT: DTBUF 16:43:20 May 11, 2020
Date and time 16:43:20 May 11, 2020 (2020132:164320)
AQCHEM errors can be difficult to debug. I do not see anything obvious.
Have you tried increasing the maximum allowable number of total iterations?
See this block beginning line 963 of aqchem.F:
IF ( ICNTAQ .GE. 60000 ) THEN
XMSG = 'Maximum AQCHEM total iterations exceeded'
CALL M3EXIT ( PNAME, JDATE, JTIME, XMSG, XSTAT2 )
END IF
Hi @CN_CODER, I was wondering if you were able to get rid of the AQCHEM error Maximum AQCHEM total iterations exceeded by trying @cgnolte’s suggestion to increase the maximum allowable number of total iterations from 60,000 to a higher number? If you did, what new number of iterations did you use? I am asking because I am encountering an error similar to yours on CMAQv5.3.2:
In my experience, problems in AQCHEM have always had bad data as their root cause.
Your log file contains these lines:
After VDIFF : Gas 1.061E-05 | Aer 8.467E-13 | Non 1.002E-29
After COUPLE : Gas 1.078E-01 | Aer 2.225E-08 | Non 1.006E-25
No BC's in file BNDY_CONC_1 for the following adv species: Set to 1.00E-30
APOCI
APOCJ
APNCOMI
APNCOMJ
ATTENTION: Applying fix to aerosol Boundary Conditions for aerosol modes.
Set verbose_rdbcon preprocessor flag to learn more.
After HADV : Gas 3.322E+00 | Aer 1.938E+11 | Non 1.447E-03
After ADV : Gas 3.322E+00 | Aer 1.937E+11 | Non 1.447E-03
Do you see the jump in aerosol species concentrations after the HADV step? That strongly indicates a problem with your BCON file, with either bad data or possibly incorrect units.
Thanks for your reply @cgnolte. I checked my BCON and ICON file. They used inconsistent profile files. Now I am using the avprofile_cb6r3m_ae7_kmtbr_hemi2016_v53beta2_m3dry_col051_row068 under the BLD directory, CCTM worked. Before ICON used my own profile file from my previous cmaq4.7.1. Actually, my previous profile file does not have so many layers (6 layers) in the vertical direction, which can be used to run cmaq4.7.1. I would like to ask if I can do some deletion of the vertical layers inside the avprofile(44 layers), because now CCTM is running very slowly, and I don’t know if there is any effect of this. Or maybe it’s because of my machine. I tried to delete some vertical layers but it reported an error that the variables could not be extracted. I uploaded my log file and part of my previous profile template. Thanks for your patience. buff_CMAQ_CCTMv531_dengt_20230324_122311_205505256.txt (27.4 KB)
Axueaha,
The number of vertical layers in your BCON file should match the number in your METCRO3D files.
You appear to be using CMAQv5.3.1, which is now over 3 years old. I recommend updating to CMAQv5.4.
Please refer to this post for suggestions on how to use this forum.
If you continue to have difficulties, please create a new thread to describe the specific problem you are experiencing.
