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?