MCIP Error encontered

Dear @tlspero @cgnolte @lizadams @peng

I got the error below while running MCIPv4.4 with output from WRFv4.2:

2022-04-28 2022-04-29
1 35 28 1

***** SUBROUTINE: GETVERSION**
***** MM5v2-FORMATTED DATA ARE NO LONGER SUPPORTED**
***** - CONSIDER CONVERTING DATA SET**
***** TO MM5v3 FORMAT USING UTILITY CODE**

** ** *** ERROR ABORT in subroutine GETVERSION** ** ABNORMAL TERMINATION IN GETVERSION** ** Date and time 0:00:00 April 28, 2022 (2022118:000000)** **

Error running mcip

MCIP4.4 accompanied CMAQv5.2, the one I want to use. Kindly note that the MCIP4.4 runs well with WRF4.2 outputs on another server of ours. I just don’t know why encountering the error above on our new server. I need solution to it, please.
I will really appreciate your suggestions and help in fixing the above error.

Thanks and best regards

As suggested by @cjcoats, I used the recent version of MCIPV5.4 but still got the error below:

==============================================================================

              US EPA COMMUNITY MULTISCALE AIR QUALITY MODEL
                METEOROLOGY-CHEMISTRY INTERFACE PROCESSOR

                       MCIP V5.4 FROZEN 07/12/2022

==============================================================================


*** SUBROUTINE: SETUP
*** ERROR OPENING WRF NETCDF FILE
*** FILE = /p1g/dust-department/idsulaymon/MENA/may2022/wrf_output/wrfout_d02_2022-04-28_00:00:00
*** NCF: NetCDF: Unknown file format


 *** ERROR ABORT in subroutine SETUP
 ABNORMAL TERMINATION IN SETUP
 Date and time  0:00:00   April 28, 2022  (2022118:000000)

Error running mcip

Kindly help in proffering solution to it.

Thanks

The fact that are able to read the file on one system suggests to me that it is either something wrong with your NetCDF on the new system or that the file has become corrupted. NetCDF files ought to be platform-independent.
Are you able to ncdump -h this file to view its header?

Also, I will note that this version of WRF is from April 2020, several years after the versions of MCIP and CMAQ you are using. Beginning I think with WRF version 4.0, the default vertical coordinate changed. The older version of MCIP does not understand and properly handle this vertical coordinate. I do not know how much of an impact this would have on layer 1 pollutant concentrations. See the release notes in CMAQ v5.3.

Thanks for your response.

I could view the header.
I don’t think it has to do with the versions. Besides, I have used MCIPv5.4, the most recent version but still got errors.

Hi @Catalyst2682!

Did you do anything to post-process your WRF output files after you ran WRF? That is, is the wrfout file that you are providing to MCIP the exact same file that was generated by WRF?

If you are using the WRF output directly, then please post the list of the global attributes that you see in the WRF header from your “ncdump -h ” command.

–Tanya

Hello @tlspero,

I am sorry for responding lately.

No, I didn’t do anything to the WRF output and I am using it directly. Below is the list of the global attributes in the WRF header:

// global attributes:
:TITLE = " OUTPUT FROM WRF V4.2.1 MODEL" ;
:START_DATE = “2022-04-28_00:00:00” ;
:SIMULATION_START_DATE = “2022-04-28_00:00:00” ;
:WEST-EAST_GRID_DIMENSION = 241 ;
:SOUTH-NORTH_GRID_DIMENSION = 271 ;
:BOTTOM-TOP_GRID_DIMENSION = 44 ;
:DX = 12000.f ;
:DY = 12000.f ;
:AERCU_OPT = 0 ;
:AERCU_FCT = 1.f ;
:IDEAL_CASE = 0 ;
:DIFF_6TH_SLOPEOPT = 0 ;
:AUTO_LEVELS_OPT = 2 ;
:DIFF_6TH_THRESH = 0.1f ;
:DZBOT = 50.f ;
:DZSTRETCH_S = 1.3f ;
:DZSTRETCH_U = 1.1f ;
:SKEBS_ON = 0 ;
:SPEC_BDY_FINAL_MU = 1 ;
:USE_Q_DIABATIC = 0 ;
:GRIDTYPE = “C” ;
:DIFF_OPT = 1 ;
:KM_OPT = 4 ;
:DAMP_OPT = 0 ;
:DAMPCOEF = 0.2f ;
:KHDIF = 0.f ;
:KVDIF = 0.f ;
:MP_PHYSICS = 8 ;
:RA_LW_PHYSICS = 4 ;
:RA_SW_PHYSICS = 4 ;
:SF_SFCLAY_PHYSICS = 1 ;
:SF_SURFACE_PHYSICS = 2 ;
:BL_PBL_PHYSICS = 1 ;
:CU_PHYSICS = 3 ;
:SF_LAKE_PHYSICS = 0 ;
:SURFACE_INPUT_SOURCE = 1 ;
:SST_UPDATE = 0 ;
:GRID_FDDA = 1 ;
:GFDDA_INTERVAL_M = 360 ;
:GFDDA_END_H = 816 ;
:GRID_SFDDA = 0 ;
:SGFDDA_INTERVAL_M = 0 ;
:SGFDDA_END_H = 0 ;
:HYPSOMETRIC_OPT = 2 ;
:USE_THETA_M = 1 ;
:GWD_OPT = 0 ;
:SF_URBAN_PHYSICS = 0 ;
:SF_SURFACE_MOSAIC = 0 ;
:SF_OCEAN_PHYSICS = 0 ;
:SHCU_PHYSICS = 0 ;
:MFSHCONV = 0 ;
:FEEDBACK = 1 ;
:SMOOTH_OPTION = 0 ;
:SWRAD_SCAT = 1.f ;
:W_DAMPING = 0 ;
:DT = 20.f ;
:RADT = 30.f ;
:BLDT = 0.f ;
:CUDT = 5.f ;
:AER_OPT = 0 ;
:SWINT_OPT = 0 ;
:AER_TYPE = 1 ;
:AER_AOD550_OPT = 1 ;
:AER_ANGEXP_OPT = 1 ;
:AER_SSA_OPT = 1 ;
:AER_ASY_OPT = 1 ;
:AER_AOD550_VAL = 0.12f ;
:AER_ANGEXP_VAL = 1.3f ;
:AER_SSA_VAL = 0.85f ;
:AER_ASY_VAL = 0.9f ;
:MOIST_ADV_OPT = 1 ;
:SCALAR_ADV_OPT = 1 ;
:TKE_ADV_OPT = 1 ;
:DIFF_6TH_OPT = 0 ;
:DIFF_6TH_FACTOR = 0.12f ;
:FGDT = 0.f ;
:GUV = 0.0003f ;
:GT = 0.0003f ;
:GQ = 0.0001f ;
:IF_RAMPING = 0 ;
:DTRAMP_MIN = 60.f ;
:OBS_NUDGE_OPT = 0 ;
:BUCKET_MM = -1.f ;
:BUCKET_J = -1.f ;
:PREC_ACC_DT = 0.f ;
:ISFTCFLX = 0 ;
:ISHALLOW = 0 ;
:ISFFLX = 1 ;
:ICLOUD = 1 ;
:ICLOUD_CU = 0 ;
:TRACER_PBLMIX = 1 ;
:SCALAR_PBLMIX = 0 ;
:YSU_TOPDOWN_PBLMIX = 0 ;
:GRAV_SETTLING = 0 ;
:DFI_OPT = 0 ;
:NTASKS_X = 4 ;
:NTASKS_Y = 6 ;
:NTASKS_TOTAL = 24 ;
:SIMULATION_INITIALIZATION_TYPE = “REAL-DATA CASE” ;
:WEST-EAST_PATCH_START_UNSTAG = 1 ;
:WEST-EAST_PATCH_END_UNSTAG = 240 ;
:WEST-EAST_PATCH_START_STAG = 1 ;
:WEST-EAST_PATCH_END_STAG = 241 ;
:SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
:SOUTH-NORTH_PATCH_END_UNSTAG = 270 ;
:SOUTH-NORTH_PATCH_START_STAG = 1 ;
:SOUTH-NORTH_PATCH_END_STAG = 271 ;
:BOTTOM-TOP_PATCH_START_UNSTAG = 1 ;
:BOTTOM-TOP_PATCH_END_UNSTAG = 43 ;
:BOTTOM-TOP_PATCH_START_STAG = 1 ;
:BOTTOM-TOP_PATCH_END_STAG = 44 ;
:GRID_ID = 2 ;
:PARENT_ID = 1 ;
:I_PARENT_START = 39 ;
:J_PARENT_START = 55 ;
:PARENT_GRID_RATIO = 3 ;
:CEN_LAT = 23.96781f ;
:CEN_LON = 46.43536f ;
:TRUELAT1 = 0.f ;
:TRUELAT2 = 40.f ;
:MOAD_CEN_LAT = 25.00001f ;
:STAND_LON = 47.f ;
:POLE_LAT = 90.f ;
:POLE_LON = 0.f ;
:GMT = 0.f ;
:JULYR = 2022 ;
:JULDAY = 118 ;
:MAP_PROJ = 1 ;
:MAP_PROJ_CHAR = “Lambert Conformal” ;
:MMINLU = “MODIFIED_IGBP_MODIS_NOAH” ;
:NUM_LAND_CAT = 20 ;
:ISWATER = 17 ;
:ISLAKE = -1 ;
:ISICE = 15 ;
:ISURBAN = 13 ;
:ISOILWATER = 14 ;
:HYBRID_OPT = 2 ;
:ETAC = 0.2f ;

Thanks in anticipation as you help in fixing the error.

@cgnolte @tlspero

I recently got this error:

/g/home/models/CMAQv-5.2/PREP/mcip_new/src/mcip.exe: error while loading shared libraries: libnetcdff.so.5: cannot open shared object file: No such file or directory
Error running mcip

The libnetcdff.so.5 is in the environmental path of netcdf but I don’t know why getting such error again.

You need to set the LD_LIBRARY_PATH environment variable to point to the directory where libnetcdff.so.5 is located. Is it?

Yes, it is in the bashrc file.
Or do I need to also set it on my working environment before running run_mcip.csh script?

I have fixed this in MCIP4.4 and mcip.exe built successfully. However, when I attempted to run the code, I got GETVERSION error as shown below:


*** SUBROUTINE: GETVERSION
*** MM5v2-FORMATTED DATA ARE NO LONGER SUPPORTED
*** - CONSIDER CONVERTING DATA SET
*** TO MM5v3 FORMAT USING UTILITY CODE


 *** ERROR ABORT in subroutine GETVERSION
 ABNORMAL TERMINATION IN GETVERSION
 Date and time  0:00:00   April 28, 2022  (2022118:000000)

Error running mcip

In another instance, I built MCIP5.1 version but also got an error as shown below:

- SUBROUTINE SETUP_WRFEM - READING WRF HEADER
** WRF GRID DIMENSIONS (X,Y,Z) 241 271 43**


** *** SUBROUTINE: SETUP_WRFEM**
** *** ERROR RETRIEVING VARIABLE FROM WRF FILE**
** *** VARIABLE = ZNW**
** *** NCF: NetCDF: Index exceeds dimension bound**


** *** ERROR ABORT in subroutine SETUP_WRFEM**
** ABNORMAL TERMINATION IN SETUP_WRFEM**
** Date and time 0:00:00 April 28, 2022 (2022118:000000)**

Error running mcip

Thanks in anticipation in solving these problems.

I recommend creating a private module file and loading it prior to running mcip.

Instructions are here:
https://researchcomputing.princeton.edu/support/knowledge-base/custom-modules