CCTM running errors

CMAQ
1

Dear All,

I’m running CCTM and I’m getting the “Inconsistent Base/File Xorig” error; the related part of log file is as follows; I’ve tried several times, but I haven’t been able to solve it; as a beginner in CMAQ, can you tell me what I should do?

Thanks in advance

 ========================================
 |>---   OPENING CMAQ INPUT FILES   ---<|
 ========================================
 INIT_MEDC_1     :notused
 
 >>--->> WARNING in subroutine OPEN3
 File not available.
 
  
 "GRID_CRO_2D" opened as OLD:READ-ONLY   
 File name "/depot/gmichals/apps/cmaq/cmaq-5.4/data/test_2024/met/mcipv5.4/GRIDCRO2D_240311.nc"
 File type GRDDED3 
 Execution ID "mcip"
 Grid name "2024_test_CROSS"
 Dimensions: 132 rows, 142 cols, 1 lays, 7 vbles
 NetCDF ID:    131072  opened as READONLY            
 Time-independent data.
  
 "GRID_DOT_2D" opened as OLD:READ-ONLY   
 File name "/depot/gmichals/apps/cmaq/cmaq-5.4/data/test_2024/met/mcipv5.4/GRIDDOT2D_240311.nc"
 File type GRDDED3 
 Execution ID "mcip"
 Grid name "2024_test_DOT"
 Dimensions: 133 rows, 143 cols, 1 lays, 9 vbles
 NetCDF ID:    196608  opened as READONLY            
 Time-independent data.
  
 "LUFRAC_CRO" opened as OLD:READ-ONLY   
 File name "/depot/gmichals/apps/cmaq/cmaq-5.4/data/test_2024/met/mcipv5.4/LUFRAC_CRO_240311.nc"
 File type GRDDED3 
 Execution ID "mcip"
 Grid name "2024_test_CROSS"
 Dimensions: 132 rows, 142 cols, 21 lays, 1 vbles
 NetCDF ID:    262144  opened as READONLY            
 Time-independent data.
  
 "MET_CRO_2D" opened as OLD:READ-ONLY   
 File name "/depot/gmichals/apps/cmaq/cmaq-5.4/data/test_2024/met/mcipv5.4/METCRO2D_240311.nc"
 File type GRDDED3 
 Execution ID "mcip"
 Grid name "2024_test_CROSS"
 Dimensions: 132 rows, 142 cols, 1 lays, 35 vbles
 NetCDF ID:    327680  opened as READONLY            
 Starting date and time  2008002:000000 (0:00:00   Jan. 2, 2008)
 Timestep                          010000 (1:00:00 hh:mm:ss)
 Maximum current record number        25
  
 "MET_DOT_3D" opened as OLD:READ-ONLY   
 File name "/depot/gmichals/apps/cmaq/cmaq-5.4/data/test_2024/met/mcipv5.4/METDOT3D_240311.nc"
 File type GRDDED3 
 Execution ID "mcip"
 Grid name "2024_test_DOT"
 Dimensions: 133 rows, 143 cols, 44 lays, 6 vbles
 NetCDF ID:    393216  opened as READONLY            
 Starting date and time  2008002:000000 (0:00:00   Jan. 2, 2008)
 Timestep                          010000 (1:00:00 hh:mm:ss)
 Maximum current record number        25
 

 |> Reading Emission Control Namelist:
 +======================================
 
 Performing Basic Error Checks for Emission Scaling Rules

GR_EM_SYM_DATE_0 | T

 "GR_EMIS_001" opened as OLD:READ-ONLY   
 File name "/depot/gmichals/apps/cmaq/cmaq-5.4/data/test_2024/emis/emis_mole_all_2008002.ncf"
 File type GRDDED3 
 Execution ID "????????????????"
 Grid name "SE52BENCH"
 Dimensions: 132 rows, 142 cols, 6 lays, 46 vbles
 NetCDF ID:    458752  opened as READONLY            
 Starting date and time  2008002:000000 (0:00:00   Jan. 2, 2008)
 Timestep                          010000 (1:00:00 hh:mm:ss)
 Maximum current record number        25
 RELOFFX:  -0.458000000000000     
 XORIG_C, XORIG_F:   -71.0000000000000       -73.5420000000000
2

This message is saying that your emissions file is for a rather different region than your met and grid files: the SW corner for your emissions is about 250KM west of the SW corner for everything else.

So you need to get the right input files for your modeling domain.

3

Thanks for your reply so much! I am trying to feed CMAQ with WRF’s outcome, so I guess the domain setup differences between these two models cause my errors, according to your guidance. The emission data comes from a third party that cannot be modified; the related header is below.
##############################################
:NCOLS = 142 ;
:NROWS = 132 ;
:NLAYS = 6 ;
:NVARS = 46 ;
:GDTYP = 2 ;
:P_ALP = 0. ;
:P_BET = 60. ;
:P_GAM = 104.32 ;
:XCENT = 104.32 ;
:YCENT = 38.814 ;
:XORIG = -2647500. ;
:YORIG = -2375500. ;
:XCELL = 36000. ;
:YCELL = 36000. ;
###################################################
So my question is how to modify the domain for WRF; my previous namelist.wps is below. Looking forward to your further guidance.

&geogrid
parent_id = 1
parent_grid_ratio = 1
i_parent_start = 1
j_parent_start = 1
e_we = 145
e_sn = 135
geog_data_res = ‘default’
dx = 36000
dy = 36000
map_proj = ‘lambert’
ref_lat = 38.814
ref_lon = 104.32
truelat1 = 36.00
truelat2 = 37.00
pole_lat = 90
pole_lon = 0
stand_lon = 105.117
geog_data_path = ‘geog’
opt_geogrid_tbl_path = ‘./geogrid/’
/

4

For your emissions data, the map-projection parameters

:P_ALP = 0. ;
:P_BET = 60. ;
:P_GAM = 104.32 ;

correspond to the respective WRF map projection parameters:

truelat1 = 36.00
truelat2 = 37.00
stand_lon = 105.117

inidicating that the map projections for the two are quite different. You need a WRF run for which truelat1 = 0.0. truelat2 = 60.0, stand_lon = 104.32.

See https://www2.mmm.ucar.edu/wrf/users/docs/user_guide_v4/v4.1/users_guide_chap3.html

1 Like
5

Thank you very much for your help; I modified namelist.wps as you requested, as follows; and then used the generated WRF results to input into CMAQ, but still got exactly the same error as before; I am sure that I fed the modified WRF outcome into CMAQ again, and here I am posting the GRIDDESC file generated after the modification; You can see that -2556000/36000 = -71, still unchanged. Whereas we need -73.542.

I would like to know how I should further adjustments?

Thank you very much!

###########################################################
&share
wrf_core = ‘ARW’,
max_dom = 1,
start_date = ‘2008-01-01_00:00:00’,‘2008-01-01_00:00:00’,
end_date = ‘2008-01-03_00:00:00’,‘2008-01-03_00:00:00’,
interval_seconds = 3600
/

&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 12,
j_parent_start = 1, 11,

s_we = 1,
s_sn = 1,

e_we = 145, 142,
e_sn = 135, 112,
geog_data_res = ‘default’,‘default’,
dx = 36000,
dy = 36000,
map_proj = ‘lambert’,
ref_lat = 38.814,
ref_lon = 104.32,
truelat1 = 0.0,
truelat2 = 60.0,
stand_lon = 104.320,
################################################################
GRIDDESC

‘LamCon_40N_97W’
2 0.000 60.000 104.320 104.320 38.814
’ ’
‘2024_test’
‘LamCon_40N_97W’ -2556000.000 -2376000.000 36000.000 36000.000 142 132 1
################################################################

6

You probably know this already, but it is a difficult problem to reverse-engineer what WRF domain was used to generate a particular grid. It would be helpful if you could find the MCIP log file (or the wrfout files) that were used to prepare the meteorological inputs for CMAQ.

7

Thanks for your reply so much! I’ve been told that the center point in the header info of emission file seems to be out of alignment with the actual center point of that data; is this possible? It sounds crazy. Isn’t the use of latitude and longitude to indicate the center point intended to avoid a similar situation?

8

I’m not sure what you mean.
For a Lambert conformal projection, XCENT and YCENT are the longitude and latitude where X and Y are 0. XORIG and YORIG are the x- and y-coordinates (in meters, relative to (XCENT, YCENT)) corresponding to the lower left (southwest) corner of the domain. If you window (subset) the domain, then XORIG and YORIG change, but the values of X and Y (and LON and LAT) corresponding to a given geographic location do not change.

10

It seems that I need more time to learn as a beginner and thanks for your help!