Could not open ISAM_REGIONS

[bmbmnb]

Dear all,
I’m trying to run CMAQ_ISAM but I got an error message in my LOG file(see below), the model doesn’t seem to find ISAM region file. I make sure I provided the right path. Does anyone know why?

      Error reading IDs for  variable  CD CN CB CQ YX
 netCDF error number  -49  processing file "ISAM_REGIONS"
 NetCDF: Variable not found


 *** ERROR ABORT in subroutine SubhFile_Cell on PE 001
 Could not open ISAM_REGIONS

PM3EXIT: date&time specified as 0
Date&time specified as 0

Yunqing

run.txt (12.0 KB)

[hogrefe.christian]

The errors you posted point to a problem with reading variables from the ISAM_REGIONS file, not with not being able to find the file. Looking at that error message, the code is looking for a variable named “CD CN CB CQ YX” while you probably want it to be looking for five different variables named CD, CN, CB, CQ, and YX. To help us diagnose this further, please post your EmissCtrl file, your isam_control file, and the header of our ISAM_REGIONS file obtained with ncdump -h.

The final error message you encountered “Could not open ISAM_REGIONS” is misleading in this case, I will flag it for the developers to have a look at.

[bmbmnb]

Hello Hogrefe,
Sorry for the late reply. Attached are the files you asked for.
EmissCtrl_cb6r3_ae7_aq.txt (31.4 KB) isam_control.txt (4.4 KB) ISAM_REGIONS.txt (1.9 KB)
Yunqing

[hogrefe.christian]

Thanks for posting these files.

Looking at the header of the ISAM_REGIONS file you posted as ISAM_REGIONS.txt file, it seems the file does not fully comply with I/O API conventions, leading to an error when CMAQ tries to read from it.

Specifically, since variable names in I/O API have a length of 16 characters (including spaces), the VAR-LIST attribute cannot be " CD CN CB CQ YX YD OT " since it needs to contain 16 characters (including spaces) for each variable and have a length of NVARS x 16. You can look at the headers of MCIP files for examples of properly constructed VAR-LIST attributes. It also appears that a number of the attributes in your file that should be of type integer (FTYPE, GDTYP, NVARS, STIME, SDATE, etc.) are of type real which may also cause problems.

Generally, it is strongly advised to use Fortran code linking in the I/O API libraries or programs from I/O API’s M3tools to create CMAQ input files such as this one to ensure that they comply with I/O API standards. There are a number of other threads discussing this topic:

CMAQV5.3 ISAM Region not found

Yukui, A quick solution is to “borrow” the meta data. This solution assumes your new data is the same shape and consistent with the spatial domain in the borrowed meta data. For now, I assume you are using files with only one variable. This could be modified for multiple variables. Each step includes code that can be used from the command line. Make your mask with ncl (hereafter /path/to/nclmask.nc) Copy an existing IOAPI compliant mask file cp /path/to/existing/ioapimask.nc /path/to/ne…

Error creating GRIDMASK

Hi all~ I am trying to generate my own GRIDMASK since I use other emission inventory instead of EDGAR. I wrote a f90 file to generate a GRIDMASK file. But when SMOKE was reading the GRIDMASK, the error shows up: Value for PROMPTFLAG: N returning FALSE Error reading IDs for variable TZONES netCDF error number -49 processing file “GRIDMASK” NetCDF: Variable not found Could not open file “GRIDMASK”. I guess the error is caused by the global attribute I wrote into GRIDMASK, which is VAR…

[bmbmnb]

Thanks for your help, I will check my ISAM_REGION file for this problem.
Yunqing

[qinch18]

Hello，
Could you please share the specific method.