Error in generating inline dust emissions while running CMAQ5.0.2

Dear @tlspero , @lizadams, @thayna.gomes, @cjcoats,

Please, I’m having this problem.

I ran cctm using CMAQ5.0.2 with WRF4.1.2 outputs processed with MCIP4.5. Also, I tried WRF3.8.1 outputs processed with MCIP4.2 but still got same error. In the GRICRO2D, there are open shrublands (in LUFRAC_07) and closed shrublands (in LUFRAC_06) variables but got this error while running cctm:

Aerosol Emissions Processing in Vertical diffusion …
— Aero Species Mapped —
— PM Emis Species Mapped —
— Precursor Species Mapped —
Value for CTM_SSEMDIAG not defined;returning default: FALSE
Flag for writing the sea-salt-emission diagnostic file
Environment variable not set … Using default: 0
Value for CTM_WB_DUST: Y returning TRUE
Value for CTM_DUSTEM_DIAG not defined;returning default: FALSE
Flag for writing the windblown dust emission diagnostic file
Environment variable not set … Using default: 0

—>> WARNING in subroutine XTRACT3
Requested variable USGS_shrubland not available in file DUST_LU_1which contains variables:

 LAT
 LON
 MSFX2
 HT
 DLUSE
 LWMASK
 PURB
 LUFRAC_01
 LUFRAC_02
 LUFRAC_03
 LUFRAC_04
 LUFRAC_05
 LUFRAC_06
 LUFRAC_07
 LUFRAC_08
 LUFRAC_09
 LUFRAC_10
 LUFRAC_11
 LUFRAC_12
 LUFRAC_13
 LUFRAC_14
 LUFRAC_15
 LUFRAC_16
 LUFRAC_17
 LUFRAC_18
 LUFRAC_19
 LUFRAC_20
 LUFRAC_21
 
 *** ERROR ABORT in subroutine DUST_EMIS_INIT
 Could not read USGS_shrubland from DUST_LU_1
 Date and time  0:00:00   Jan. 31, 2016   (2016031:000000)

Command exited with non-zero status 1
0.23user 0.27system 0:09.60elapsed 5%CPU (0avgtext+0avgdata 163372maxresident)k
8800inputs+212960outputs (0major+39977minor)pagefaults 0swaps
mkdir -p /home/catalyst/Desktop/Build_WRF/CMAQv5.0.2/scripts/cctm/cctm_36km/logs/2016031
mv: No match.
@ i = 30 + 1
end
while ( 31 < 31 )
date
Thu Dec 03 17:25:32 CST 2019
exit

In my cctm script, I used:

setenv DUST_LU_1 MetPath/GRIDCRO2D{EXTN}
setenv DUST_LU_2 MetPath/GRIDCRO2D{EXTN}

Is there any setting to be done or modification to be inserted to make it working?

Your prompt help would be really appreciated.

Thanks in anticipation.

Catalyst

Hi Catalyst,
In the benchmark run script for CMAQv5.3, the environment variables are set as follows:

#> Windblown dust emissions configuration
if ( $CTM_WB_DUST == ‘Y’ ) then
# Input variables for BELD3 Landuse option
setenv DUST_LU_1 $LUpath/beld3_12US1_459X299_output_a_bench.nc
setenv DUST_LU_2 $LUpath/beld4_12US1_459X299_output_tot_bench.nc
endif

I examined the file using ncdump and see that the USGS_shrubland variable is in that file.
ncdump beld3_12US1_459X299_output_a_bench.nc | more
netcdf beld3_12US1_459X299_output_a_bench {
dimensions:
TSTEP = 1 ;
DATE-TIME = 2 ;
LAY = 1 ;
VAR = 120 ;
ROW = 80 ;
COL = 100 ;
variables:
int TFLAG(TSTEP, VAR, DATE-TIME) ;
TFLAG:units = “<YYYYDDD,HHMMSS>” ;
TFLAG:long_name = "TFLAG " ;
TFLAG:var_desc = "Timestep-valid flags: (1) YYYYDDD or (2) HHMMSS " ;
float USGS_urban(TSTEP, LAY, ROW, COL) ;
USGS_urban:long_name = "USGS_urban " ;
USGS_urban:units = "percent " ;
USGS_urban:var_desc = "BELD3 land use " ;
float USGS_drycrop(TSTEP, LAY, ROW, COL) ;
USGS_drycrop:long_name = "USGS_drycrop " ;
USGS_drycrop:units = "percent " ;
USGS_drycrop:var_desc = "BELD3 land use " ;
float USGS_irrcrop(TSTEP, LAY, ROW, COL) ;
USGS_irrcrop:long_name = "USGS_irrcrop " ;
USGS_irrcrop:units = "percent " ;
USGS_irrcrop:var_desc = "BELD3 land use " ;
float USGS_cropgrass(TSTEP, LAY, ROW, COL) ;
USGS_cropgrass:long_name = "USGS_cropgrass " ;
USGS_cropgrass:units = "percent " ;
USGS_cropgrass:var_desc = "BELD3 land use " ;
float USGS_cropwdlnd(TSTEP, LAY, ROW, COL) ;
USGS_cropwdlnd:long_name = "USGS_cropwdlnd " ;
USGS_cropwdlnd:units = "percent " ;
USGS_cropwdlnd:var_desc = "BELD3 land use " ;
float USGS_grassland(TSTEP, LAY, ROW, COL) ;
USGS_grassland:long_name = "USGS_grassland " ;
USGS_grassland:units = "percent " ;
USGS_grassland:var_desc = "BELD3 land use " ;
float USGS_shrubland(TSTEP, LAY, ROW, COL) ;
USGS_shrubland:long_name = "USGS_shrubland " ;
USGS_shrubland:units = "percent " ;
USGS_shrubland:var_desc = "BELD3 land use " ;
float USGS_shrubgrass(TSTEP, LAY, ROW, COL) ;
USGS_shrubgrass:long_name = "USGS_shrubgrass " ;
USGS_shrubgrass:units = "percent " ;
USGS_shrubgrass:var_desc = "BELD3 land use " ;

I am not sure if this option is available with CMAQ5.0.2. Is there a reason that you are not using CMAQv5.3 for your runs?