Issues on sensinput file

CMAQ-DDM-3D
1

Hi,
I’m working on CMAQ-DDM-3D v5.4. Before that, I have run CCTM scripts without DDM-3D and ISAM greatly. The error log showed “forrtl: severe (24): end-of-file during read, unit -5, file Internal List-Directed Read”.
Is there anything wrong with my sensinput file?

The sensinput file is:
EVOCOT

  • EMIS*
  • OTHER*
  • SPECIES*
  • ACETONE, ACETYLENE, ACROLEIN, ALK1, ALK2, ALK3, ALK4, ALK5, ARO1, ARO2, BACL, BALD, BENZENE, BUTADIENE13, CCHO, CCOOH, CRES, ETHENE, ETOH, GLY, HCHO, HCOOH, IPRD, ISOPRENE, MACR, MEK, MEOH, MGLY, MVK, NAPHTHAL, OLE1, OLE2, OXYL, PRD2, RCOOH, RCHO, RNO3, TMBENZ124, TOLUENE, TERP*

EVOCCI

  • EMIS*
  • GRIDCI, POINTCI*
  • SPECIES*
  • ACETONE, ACETYLENE, ACROLEIN, ALK1, ALK2, ALK3, ALK4, ALK5, ARO1, ARO2, BACL, BALD, BENZENE, BUTADIENE13, CCHO, CCOOH, CRES, ETHENE, ETOH, GLY, HCHO, HCOOH, IPRD, ISOPRENE, MACR, MEK, MEOH, MGLY, MVK, NAPHTHAL, OLE1, OLE2, OXYL, PRD2, RCOOH, RCHO, RNO3, TMBENZ124, TOLUENE, TERP*

EVOCAL

  • EMIS*
  • GRIDCI, POINTCI, OTHER*
  • SPECIES*
  • ACETONE, ACETYLENE, ACROLEIN, ALK1, ALK2, ALK3, ALK4, ALK5, ARO1, ARO2, BACL, BALD, BENZENE, BUTADIENE13, CCHO, CCOOH, CRES, ETHENE, ETOH, GLY, HCHO, HCOOH, IPRD, ISOPRENE, MACR, MEK, MEOH, MGLY, MVK, NAPHTHAL, OLE1, OLE2, OXYL, PRD2, RCOOH, RCHO, RNO3, TMBENZ124, TOLUENE, TERP*

ENOX

  • EMIS*
  • OTHER*
  • SPECIES*
  • NO, NO2*

2EVOC

  • HIGH*
  • EVOCAL*
  • EVOCAL*

END

I hope someone would give me some suggentions. Thanks!

2

Its possible that the character string for the species is too long. I can show you how to allow for more text in that read statement.

But first, another problem could be the use of too many characters in the sensitivity parameter name. We are recommending only 3 characters to accommodate the much longer chemical species names in current mechanisms and still stay within the IOAPI species names character limitations.

So, where you have EVOCOT, EVOCCI, and EVOCAL, it could be reading all three as “EVO”

You could just try “COT” “CCI” and “CAL” and come back and let me know if that was enough to get through.

If not, please, post original sensinput.dat without the forum formatting and we can go from there.

Sergey

3

Hi,
Glad for your suggestions. I tried 3 characters names and the same errores existed. Then, I tried to use “ALL” as the species and it worked. So, the errors are possible caused by the length of the species string.
The sensinput file:
sensinput.2019ci.dat.txt (1023 Bytes)

4

Yes - your SPECIES string is around 280 characters and the file is coded to read a max of 256.

You can open the file sinput.F and search for the line:

  CHARACTER( 256 ) :: TXT, TXT2, EMISTYPE

You can change that 256 to 512, for example.

Additionally, find the lines:

    READ( JVUNIT, '(A256)', END = 999 ) TXT

Change that A256 to the same number (A512, for example).

After these 2 changes, you will need to recompile the code. Please, report back if this solves your problem.

Sergey

1 Like
5

It addressed the issues.
Thanks a lot!

1 Like