MCIP and CCTM issue for CMAQ version 5.2.1


#1

Good afternoon CMAQ users:

I am trying to run CCTM, it shows these following error messages.

 XCENT_B:    -97.000000000000  XCENT3D (file):    -97.000000000000
 YCENT_B:     33.000000000000  YCENT3D (file):     40.000000000000
 XCELL_B:  12000.000000000000  XCELL3D (file):  12000.000000000000
 YCELL_B:  12000.000000000000  YCELL3D (file):  12000.000000000000


 *** ERROR ABORT in subroutine SubhFile_Cell on PE 000   
 File header inconsistent with GRID_CRO_2D
 Date&time specified as 0

However, on GRIDDESC, it shows
’ ’
‘LamCon_40N_97W’
2 33.000 45.000 -97.000 -97.000 33.000
’ ’
‘SE52BENCH’
‘LamCon_40N_97W’ -112000.000 406500.000 12000.000 12000.000 157 148 1
’ ’
which is the same as my emission file. And I check GRIDCRO2D, the global attributes are also the same as my emission file. In fact, all my mcip output (GRIDBDY2D, GRIDCRO2D, GRIDDOT2D, METBDY3D, METCRO3D, METCRO2D, METDOT3D) has these global attributes. I don’t know why these error messages were shown. Anyone has an idea?

            :P_ALP = 33. ;
            :P_BET = 45. ;
            :P_GAM = -97. ;
            :XCENT = -97. ;
            :YCENT = 33. ;

Best
Huan


#2

Hello Huan.
When you checked the GRIDCRO2D file, did you used m3edhdr iaopi tool?
This is a useful to edit values in ncf files in general. Maybe, if you edit the YCENT3d value in the GRIDCRO2d file, using that tool, from 40 to 45, you can run CCTM without error, but it is an suggestion. If you still have this error, maybe you are using different GRIDCRO2D files to get emission and process in CMAQ.

Tell us about it
Grettings


#3

Hi Grettings:

Thank you for your time!
I checked that this error message is from the default ocean file, which has different coordinate as my input. I used m3fake to create an ocean file that has the same coordinate as my input, but with 0, since my study area is inland. Thus, now, this error message is killed.
Currently, the error messages are INTERPX and PE 000. Actually, this PE 000 error always appears when I try to run CCTM. Do you have any idea how shall I solve this?

Best
Huan
>>—>> WARNING in subroutine INIT_MET:INTERPX
Variable “PURB” not in file GRID_CRO_2D
M3WARN: DTBUF 0:00:00 July 13, 2016 (2016195:000000)

 *** ERROR ABORT in subroutine INIT_MET on PE 000        
  Error interpolating variable PURB from GRID_CRO_2D
 Date and time 0:00:00   July 13, 2016  (2016195:000000)

application called MPI_Abort(MPI_COMM_WORLD, 37196448) - process 0
0.261u 0.244s 0:05.09 9.8% 0+0k 0+96io 0pf+0w
date
Wed Sep 19 15:55:03 EDT 2018
if ( ! -e /scratch/brown/fang63/CMAQ-5.2.1/data/output_CCTM_v521_intel_160713/LOGS ) then
mv: No match.
setenv NEW_START false
set TODAYG = date -ud "${TODAYG}+1days" +%Y-%m-%d
date -ud 2016-07-13+1days +%Y-%m-%d
set TODAYJ = date -ud "${TODAYG}" +%Y%j
date -ud 2016-07-14 +%Y%j
end
while ( 2016196 < = 2016195 )
exit


#4

Hi @fang63
I had to create an ocean file for my domain that was inland as well, however, what I did wasn’t correct apparently, would you mind sharing the script you used to create your OCEAN file?

This is the script I used for that:

Thank you!


#5

Hi Daniel(?):

Are you also planning to create a ocean file with zeros for the variables for an inland study area? Below is my script.

#!/bin/csh -f

#m3fake script to create a dummy ocean file
set APPL = 160713 #> Application Name (e.g. Gridname)
setenv GRIDDESC $CMAQ_HOME/data/met/mcip/$APPL/GRIDDESC
setenv GRID_NAME SE52BENCH
setenv OUTFILE $CMAQ_HOME/data/ocean/ocean_file.dummy.$GRID_NAME.ncf
/scratch/brown/fang63/CMAQ-5.2.1/lib/x86_64/intel/ioapi/Linux2_x86_64ifort/m3fake << EOF

2
$GRID_NAME
1
0
2
OPEN
UNKNOWN
OPEN
3
5

SURF
UNKNOWN
SURF
3
5

OUTFILE

EOF


#6

Thanks for sharing your script @fang63! Yes, I do also need to create an ocean file with zeros for the study area.

Quick question, what version of IOAPI are you using?


#7

I am using version 3.2


#8

Greetings,
We have corrected some errors in the sample run script we provide for creating a dummy ocean file. The updated tutorial can be found under DOCS/Tutorials on the CMAQ GitHub repository for CMAQv5.3b2:


Also CMAQv5.3 now allows a user to run CMAQ without an ocean input file. This is also explained in the tutorial. Hope this helps!