Error in generating inline dust emissions while running CMAQ5.0.2

#1

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

#2

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?