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?
Dear @lizadams,
I’ve solved this problem. For now, I prefer using CMAQv5.0.2 for some reasons.
Thanks so much and best regards.
Catalyst
Dear Catalyst2682,
I am getting the same error which you asked in the starting of this thread, can you please tell the solution of how you did it? It would be of great help.
Because as we are studying outside of the US, we do not have explicit DUST emissions input, and using GRIDCRO2D is generating the same error as yours
Regards,
Shivang
Hi Shivang,
Please start a new post, even though it is similar to an existing case and follow these best practices:
https://forum.cmascenter.org/t/please-read-before-posting/1321/3
You can reference an older issue in your new post by copying the link for that post into your issue, and use the @ command to tag someone such as @Catalyst2682 for assistance.