Program can't load CGRID file properly


#1

Hello;

I have an issue with running CCTM. It can;’ tread or load the CGRID correctly. I tried different ICFILE or INIT_GASC_1 but they all did not work well. The error is shown below.

 >>--->> WARNING in subroutine RDTFLAG
 Time step error for file:  INIT_GASC_1
 M3WARN:  DTBUF 0:00:00   Jan. 1, 2014  (2014001:000000)
 
 >>--->> WARNING in subroutine XTRACT3
 Time step not available for file:  INIT_GASC_1
 M3WARN:  DTBUF 0:00:00   Jan. 1, 2014  (2014001:000000)


 *** ERROR ABORT in subroutine LOAD_CGRID on PE 000      
 Could not read NO2 from INIT_GASC_1

And the settings are here.

#> Initial conditions
if ($NEW_START == true || $NEW_START == TRUE ) then
setenv ICFILE ts20140101.bpch.ICON.nc
#setenv ICFILE ICON_v52_profile
setenv INITIAL_RUN Y #related to restart soil information file
rm -rf LOGDIR/CTM_LOG*{RUNID}* # Remove all Log Files Since this is a new start
mkdir -p $OUTDIR

Would you help me look at it? Appreciate your help!

Yijia


#2

Please try to use the following command to check that the ICFILE you are using is complete.
ncdump -h ts20140101.bpch.ICON.nc | grep NO2
You should see NO2 listed as follows:

float NO2(TSTEP, LAY, ROW, COL) ;
	NO2:long_name = "NO2             " ;
	NO2:units = "ppmV            " ;
	NO2:coordinates = "lon lat" ;
	NO2:var_desc = "Variable NO2                                                                    " ;
	:VAR-LIST = "AIRDEN          NO2             NO              O3   ...

`


#3

Hi Yijia,

I wonder how did you generate the IC file. According to the ncdump information you have provided, it does not look like a "normal" IOAPI file.

Cheers,
David


#4

Hi Yijia and David,

The *bpch.ICON.nc" file may need to be processed to put it in the format used for the initial condition for CMAQ. Until we determine how to do this, you could try using m3tshift (https://www.cmascenter.org/ioapi/documentation/all_versions/html/M3TSHIFT.html) to change the date on an ICON file for June 01, 2014 that is available on the google drive:
https://drive.google.com/drive/u/1/folders/1_6jVTUxAq9pBAZY-W4og6kxxEuqK4NNs

Liz


#5

Hi Yijia,

It looks like you need to use the geos2cmaq program available here:


This can be used to convert the geochem output file into a CMAQ ICON input file.

Liz