what case this error?
WRF FILES DO NOT SEEM TO BE FROM SAME DOMAIN
- SUBROUTINE SETUP_WRFEM - READING WRF HEADER
WRF GRID DIMENSIONS (X,Y,Z) 85 76 33
*** SUBROUTINE: CHKWPSHDR
*** WRF FILES DO NOT SEEM TO BE FROM SAME DOMAIN
*** VARIABLE = DX
*** FIRST FILE = /sqwork/project/rc_cmaq/data/cq_rc_cmaq/wrf/wrfout_d02_2020-05-11_00:00:00
*** VALUE IN FIRST FILE = 1933.333
*** NEW FILE = /sqwork/project/rc_cmaq/data/cq_rc_cmaq/wrf/geo_em.d02.nc
*** VALUE IN NEW FILE = 1933.333
*** ERROR ABORT in subroutine CHKWPSHDR
ABNORMAL TERMINATION IN CHKWPSHDR
Date and time 0:00:00 May 11, 2020 (2020132:000000)
find in code:
rcode = nf90_get_att (cdfid, nf90_global, ‘DX’, rval)
IF ( rcode == nf90_noerr ) THEN
IF ( ABS( rval - met_resoln ) > rtol ) THEN
WRITE (*,f9000) TRIM(pname), 'DX'
WRITE (*,f9200) TRIM(fl1), met_resoln, TRIM(fl), rval
CALL graceful_stop (pname)
ENDIF
ELSE
WRITE (*,f9400) TRIM(pname), 'DX', TRIM(nf90_strerror(rcode))
CALL graceful_stop (pname)
ENDIF
Hi,
Do you have your met_em and geo_em files stored in separate directories? It looks like mcip is reading in a geo_em file where it expects to find a met_em.
Megan
@CN_CODER –
Try changing the value of RTOL in chkwpshdr.f90 from 1.0e-4 to 1.0e-3, and recompile. I think the precision is too strict on RTOL for your case.
–Tanya
Hy,
I’m facing following error while running MCIP using wrfout file:
==============================================================================
- SUBROUTINE SETUP_WRFEM - READING WRF HEADER
WRF GRID DIMENSIONS (X,Y,Z) 80 61 35
*** SUBROUTINE: LL2XY_LAM
*** TANGENT LAMBERT CONFORMAL PROJECTION DETECTED
*** TRUE LATITUDES = 37.430 37.430
*** MAY NOT WORK PROPERLY IN SPATIAL ALLOCATOR
*** ERROR ABORT in subroutine LL2XY_LAM
ABNORMAL TERMINATION IN LL2XY_LAM
Date and time 12:00:00 June 1, 2020 (2020153:120000)
Error running mcip
ERROR, file on 20200601 not produced
Also, please see my attached namelist.mcip file.
Please let me know how to fix this.namelist.mcip.txt (680 Bytes)
Thanks
Rasel
@munshimdrasel –
If you are not planning to use the Spatial Allocator tool, then you can find this section of the code in MCIP and comment out the call to graceful_stop. Recompile and proceed.
If you are planning to use the Spatial Allocator tool, then you would need a different map projection.
–Tanya
@tlspero
I’m sorry I’m a beginner in CMAQ. I didn’t understand what you’re saying. Can you explain more?
Also, I’ve used the same script to run for the whole US domain. It worked fine then. However, now I reduced my model domain in wrf and changed mcip script accordingly. Now I’m facing this error. I’m not sure what’s going on.
Thanks
Rasel
@munshimdrasel –
Sorry about that.
If you changed your WRF domain and reran WRF, then that might be the issue that is causing the problem. If you previously ran CMAQ and did not need to use the Spatial Allocator, then you may not need it for this application, either.
The error message prevents users from moving forward if there is a type of projection that cannot be handled by one of the downstream tools.
In your case and since you seem to have run CMAQ before, let’s assume that you will not need the Spatial Allocator (which is optional for some applications with CMAQ). CMAQ can handle your WRF domain. Edit MCIP subroutine ll2xy_lam.f90 by adding a Fortran comment mark ("!") to the first column of line 73, and delete the comment mark from line 74. Then recompile MCIP using its Makefile. Then rerun.
Please let me know if that solves your problem.
–Tanya
@tlspero
I’ve rebuilt mcip like the way you suggested. Now I’ve got this error:
==============================================================================
- SUBROUTINE SETUP_WRFEM - READING WRF HEADER
WRF GRID DIMENSIONS (X,Y,Z) 80 61 35
*** SUBROUTINE: LL2XY_LAM
*** TANGENT LAMBERT CONFORMAL PROJECTION DETECTED
*** TRUE LATITUDES = 37.430 37.430
*** MAY NOT WORK PROPERLY IN SPATIAL ALLOCATOR
+++ met_pcp_incr = 0
*** SUBROUTINE: SETGRIDDEFS
*** Requested MCIP X domain exceeds met domain
*** METCOL, METROW = 80 61
*** MCIP domain (col, row) = 83 64
*** ERROR ABORT in subroutine SETGRIDDEFS
ABNORMAL TERMINATION IN SETGRIDDEFS
Date and time 12:00:00 June 1, 2020 (2020153:120000)
Any suggestions?
Thanks
Rasel
@munshimdrasel –
This error means that the size of the MCIP output domain that you are requesting is larger than the WRF domain that you are supplying. In your script, you need to reduce the NCOLS and NROWS for the MCIP domain so that they are less than METCOL and METROW, respectively. Try NCOLS=77 and NROWS=58, or set BTRIM=0 and skip this section in the script.
Let me know if that works.
Thank you for the help. It worked.
1 Like