[Kuljeet]

hi, i am trying to use this code to compute FSK distribution(this doesnt matter) from a file which has absorption coefficient and wavenumber. basically it is reading a file and in output it is creating some files. it is generating an error because of the format of input file. if someone can help i really appreciate it.

this is the reading file absco.dat

99999 000.0000-999.9900
0.00E+00 0.01
0.00E+00 0.02
0.00E+00 0.03
0.00E+00 0.04
0.00E+00 0.05
0.00E+00 0.06
0.00E+00 0.07
0.00E+00 0.08
0.00E+00 0.09
0.00E+00 0.1
0.00E+00 0.11
0.00E+00 0.12
... ...

there are two columns here.

and this is the format required

! absco.dat: input file for absorption coefficient absco
! 1st line: numk, Deta (in I5,F7.4 format)
! 2nd-(numk+1)th line: wvnm,absco (in e12.4 format)

thanks in advance

-K

[flotus1]

Which part do you want to change?
The format of the input file so it matches what the code expects?
Or the code so it can read the input file as-is?

[Kuljeet]

Thanks for the reply. I want to change the format of input file.

[flotus1]

No idea what to make of the first line. The code expects one integer length 5 and one floating point. The file contains one integer length 5 (this is fine) and a range? Change "000.0000-999.9900" to whatever single value "deta" is supposed to be. E.g. 99999 999.9900
Also note that 999.9900 is technically too large to fit into F7.4. That's 4 digits after the decimal point, leaving only 3 digits. One for the decimal point itself, one for a sign, and one for the digit before the decimal point. Compilers might read this without throwing an error, but some truncation WILL occur.

The other lines are more straightforward.
The code expects two E12.4 floating point numbers. In your example, it should look like this:

 0.0000E+00  0.1000E-01
 0.0000E+00  0.2000E-01
 0.0000E+00  0.3000E-01
 0.0000E+00  0.4000E-01
 0.0000E+00  0.5000E-01
 0.0000E+00  0.6000E-01
 0.0000E+00  0.7000E-01
 0.0000E+00  0.8000E-01
 0.0000E+00  0.9000E-01
 0.0000E+00  0.1000E+00
 0.0000E+00  0.1100E+00
 0.0000E+00  0.1200E+00

[Kuljeet]

got resolved. Thanks