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
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. **
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).
Following up with the aforementioned question, I wonder if the 10 subroutines as mentioned above would be all the CCTM programs I should revise for the purpose of adding the USGS28 to CCTM in M3DRY option.
Following up with the aforementioned question, I wonder if the 10 subroutines as mentioned above are all the CCTM programs I should revise for the purpose of adding the USGS28 to CCTM in M3DRY option.
Also FYR, please find the attached 10 modified subroutines, in which I added USGS28 option.
Greetings everyone,
My apologies for the long-waiting time to this update. It’s been delayed due to work priority shifts.
I recently revised and modified the 10 subroutines, which were presented previously. The CCTM can eventually run successfully, with using meteorological data processed from WRF output data using USGS28 land use data, DDM being turned on, and setting M3DRY deposition. This is accomplished now by adding the USGS28 land use category to the 10 subroutines of the CCTM source codes.
Please find the following attachments of the ten modified CCTM subroutines FYR in case that you may need.
I have also rerun WRF with the MODIS land use data and regenerated the meteorological data, which is inherently compatible with the CMAQ5.5, as you know. I will compare the results between them, and QA/QC the data and modelling results. In case that any inconsistencies would be found, I will update the information to let you know.
Moreover, if you may find any issues in the codes I have modified, or have any questions or suggestions, please let me know so that they can get improved if possible.