Hello, I have errors about “GRIDDESC”. I am sure the GRIDDESC is in the specified path and the GRIDNAME in it matches that I set in run_cctm.csh. Here are the logfile and GRIDDESC. Thank you so much!
Execution ID "mcip"
Grid name "US_d01_CROSS"
Dimensions: 299 rows, 459 cols, 23 lays, 16 vbles
NetCDF ID: 65536 opened as READONLY
Starting date and time 2023009:000000 (0:00:00 Jan. 9, 2023)
Timestep 010000 (1:00:00 hh:mm:ss)
Maximum current record number 25
MET data determined based on WRF ARW version 3.9.1
Closing file MET_CRO_3D
========================================
|>--- RETRIEVE HORIZONTAL GRID ---<|
========================================
File "GRIDDESC" opened for input on unit: 98
/share/atmoschem/nwu3/CMAQv5.3.2.1/data/mcip/GRIDDESC
>>--->> WARNING in subroutine HGRD_INIT on PE 000
Failure retrieving horizontal grid parameters
>>--->> WARNING in subroutine GRID_INIT
*** Failure defining horizontal domain
*** ERROR ABORT in subroutine DRIVER on PE 000
*** Failure defining domain configuration
PM3EXIT: date&time specified as 0
Date&time specified as 0
It would be interesting to see what happens when you use this file and m3tools program * latlonto try to generate a LATLON file – does that work, or else what are the error messages?
If the CCTM is properly using I/O API routine DSCGRID to read the GRIDDESC file, then you should be getting proper error-messages from it, describing just what the nature of the error is. But I don’t see such messages ;-(
C With GRID_NAME (only input) retrieve all horizontal grid parameters from
C the grid description file pointed to by the GRIDDESC env var:
WRITE( LOGDEV, * )
CALL LOG_HEADING( LOGDEV, "Retrieve Horizontal Grid" )
IF ( .NOT. DSCGRID ( GRID_NAME,
& COORD_SYS_NAME, GDTYP_GD,
& P_ALP_GD, P_BET_GD, P_GAM_GD,
& XCENT_GD, YCENT_GD,
& XORIG_GD, YORIG_GD, XCELL_GD, YCELL_GD,
& GL_NCOLS, GL_NROWS, NTHIK ) ) THEN
XMSG = 'Failure retrieving horizontal grid parameters'
CALL M3WARN ( 'HGRD_INIT', 0, 0, XMSG )
SUCCESS = .FALSE.; RETURN
END IF
If setenv GRID_NAME US_d01 and setenv GRIDDESC /share/atmoschem/nwu3/CMAQv5.3.2.1/data/mcip/GRIDDESC are both set in the run script, I am not quite sure why this error would occur if the GRIDDESC file is o.k. - the test suggested by @cjcoats would be an independent way of verifying that it is.
I am not sure if this is a problem or not, but the GRIDDESC file I’ve been using and the example GRIDDESC file provided with the I/O API documentation have a quoted space at the end of the coordinate and grid blocks, not at the beginning like in the one you posted above.
There are two parts to a GRIDDESC: the “map projections” section and the “grids” section. Both of these have lists of names and defining parameters, that are ended by a quoted blank in what would otherwise be the name-portion of the respective section. A blank leading the “map projections” section would say that there are zero map projections, and what follows should be the list of grid-names and grid-parameters. Which, of course, does not match up with the rest of the original poster’s file.
I’m sorry that I omit the last line in my GRIDDESC when I posted it. I will corrected it in the original post. It is generated by MCIP run, and have a quoted space at the end. The GRID_NAME and GRIDDESC are both specified in the run script. I’m quite confused about this error
