Add land use classification USGS28 to CMAQ 5.5

CMAQ
1

Hello everyone,
I was trying to add the USGS28 land classification to the CMAQ v5.5 CCTM codes. But encountered some issues need your insights and helps.
As we know, CMAQ is compatible with USGS24 or Modis land use scheme but incompatible with USGS28. However, attempting to utilize WRF output from a previous project which used USGS28 land use scheme, I need to use the MCIP output generated from the WRF output. To make things more complex, I need to use CMAQ DDM to conduct sensitivity analysis. Given currently CMAQ DDM is only compatible with M3DRY, I have no other choices than adding USGS28 to the CMAQ v5.5 codes.
The model/module versions used are:
WRF 4.5.1
GEOGRID 4.5
MCIP (5.4)
CMAQ 5.5
Compiler: GNU Fortran (GCC) 8.5.0
In the new version CMAQ v5.5, info associated with land use scheme scatters throughout 10 subroutines of the CCTM. If I read them correctly, they may include: aero_depv.F, centralized_io_module.F, CSQY_DATA.F, LSM_MOD.F, lus_data_module.F, LUS_DEFN.F, m3dry.F, PHOTOLYSIS_ALBEDO.F, HGSIM.F, PHOT_OPTICS.dat.
I added the USGS28 land use scheme to all the above codes. The codes were successfully compiled but the CCTM run resulted in error message as follows.
Thanks,
Wen Xu

  1. CCTM run error message:
    Program received signal SIGFPE: Floating-point exception - erroneous arithmetic operation.
    Backtrace for this error:
    #0 0x7f1003a4f171 in ???
    #1 0x7f1003a4e313 in ???
    #2 0x7f1002e0fd0f in ???
    #3 0xd77a4d in hlconst_
    at /data/CMAQ_5.5/CCTM/scripts/BLD_CCTM_v55_DDM3D_gcc/hlconst.F:523
    #4 0xb10411 in m3dry_
    at /data/CMAQ_5.5/CCTM/scripts/BLD_CCTM_v55_DDM3D_gcc/m3dry.F:583
    #5 0x9321ff in depv_defn_MOD_get_depv
    at /data/CMAQ_5.5/CCTM/scripts/BLD_CCTM_v55_DDM3D_gcc/DEPV_DEFN.F:544
    #6 0xaf32cf in vdiff
    at /data/CMAQ_5.5/CCTM/scripts/BLD_CCTM_v55_DDM3D_gcc/vdiffproc.F:480
    #7 0xa08fcb in sciproc
    at /data/CMAQ_5.5/CCTM/scripts/BLD_CCTM_v55_DDM3D_gcc/sciproc.F:231
    #8 0x9f5a59 in cmaq_driver_
    at /data/CMAQ_5.5/CCTM/scripts/BLD_CCTM_v55_DDM3D_gcc/driver.F:781
    #9 0x9ee7c6 in cmaq
    at /data/CMAQ_5.5/CCTM/scripts/BLD_CCTM_v55_DDM3D_gcc/cmaq_main.F:97
    #10 0x9eead5 in main
    at /data/CMAQ_5.5/CCTM/scripts/BLD_CCTM_v55_DDM3D_gcc/cmaq_main.F:32

Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.

mpirun noticed that process rank 94 with PID 0 on node azu-vnl0007 exited on signal 8 (Floating point exception).
2. Error from *.log:
Error: Invalid argument (22)
[azu-vnl0007:12308] [[3875,0],0] ORTE_ERROR_LOG: Data unpack would read past end of buffer in file util/show_help.c at line 501
|> Multiprocess control, output and error checking:
[azu-vnl0007:12308] Set MCA parameter “orte_base_help_aggregate” to 0 to see all help / error messages
0 /lib64/libucs.so.0(ucs_handle_error+0x2dc) [0x7f0fecf5413c]
Backtrace for this error:
** Runscript Detected an Error: CGRID file was not written. **

2

I guess you should seek comments from Tanya regarding this matter.

3

Hi Ryan,

Thank you very much for the info. The comments from Tanya is valuable and helpful. However, it’s for the CMAQ v5.3, in which only 3 subroutines are related to land use classifications. Since CMAQ v5.4, the associated info has been more decentralized (as in the 10 subroutines mentioned above).

Best regards,

Wen