The CGRID file is the restart file that includes model species concentrations at the end of each simulation period.
I don’t think that this file would cause NaN values.
It looks like you are using a profile file for the initial conditions instead of a CGRID file.
setenv ICFILE ICON_v52_China_EV_profile
The log file does report no Initial conditions found for many species, but that should be ok as CMAQ sets the values to a low concentration.
No IC found for species XO2 in INIT_GASC_1; set to 1.00E-30
No IC found for species XO2N in INIT_GASC_1; set to 1.00E-30
No IC found for species NTROH in INIT_GASC_1; set to 1.00E-30
No IC found for species NTRALK in INIT_GASC_1; set to 1.00E-30
Please also check the mpi log files as you are running this case on 32 processors.
CTM_LOG_0[01-031].v52_gcc_China_EV_20121223
These per-processor log file likely contain more information
It looks like they are saved to the following directory according to your log file.
/projects/b1045/cmaq/cmaq_v5.2/data/output_CCTM_v52_gcc_China_EV_4xPROC/LOGS/
These files should all be similar/identical, so please attach one of them to this issue for us to review.
You are also getting this warning in your log file that may indicate a problem with the Emissions.
*WARNING:
One or more emissions surrogates assigned to model species
are not found in emissions sources but the CTM_EMISCHK
environment variable set to False so simulation will proceed.
You can also try to reduce the timestep by changing CTM_MINSYNC to 15
#> Sychronization Time Step and Tolerance Options
setenv CTM_MAXSYNC 300 #> max sync time step (sec) [ default: 720 ]
setenv CTM_MINSYNC 60 #> min sync time step (sec) [ default: 60 ]
recommend
- check WRF to see if wind speeds are too high.
- try compiling with the 5.2.1 version
- try different science modules? For Driver, switch between wrf and yamo
Just before the NaN values in the HDIFF module the ADV module is run.
Looking at your values and comparing them to the benchmark case:
after HADV G 1.0017700E+03 A 7.3521514E+12 N 4.5811871E-01
In the benchmark case, the values look like this
after HADV G 2.3950352E+03 A 6.2361959E+14 N 1.4660137E+01
The values being printed out are in the util/util/cksummer.F routine
GC_CKSUM / LCELLS, which I believe are values for the Gas Phase Species number concentration
AE_CKSUM / LCELLS, which I believe are values for the Aerosol Species number concentration
NR_CKSUM / LCELLS, which I believe are values for the Non-reactive Species number concentration
Please also attach your bldit_cctm.csh script, so that we can review your science process options and your compiler. You can also try recompiling the code with the -g option or do some sort of runtime trap for uninitialized values.
https://www.dursi.ca/post/stopping-your-program-at-the-first-nan.html
Another suggestion is to look at the PBL scheme used in the meteorology model and use the same scheme in the CMAQ model.