vc RFM Downloads

RFM Downloads

01DEC10 AD Add rfmrd2.pro
This web-page contains links to RFM input files. These are (mostly) text files so you should be able to view the files just by clicking on the filename, and your web-browser should also be able to download them.

Cross-Section Data (.xsc files)

These are the cross-sectional data released with HITRAN 1992, HITRAN 1996 and HITRAN 2000 with slight modifications to convert them to the RFM format (add extra records starting with '***')

For [52] (h2k), [54] (h2k), [5 6] (h2k), [70], [72]-[75], [78]-[79] and [81], the start/end wavenumbers for each range were lined up to within 0.01 wavenumber. (This resulted in a small amount of data, representing noise only, being removed from some of the ranges.)

Suffix '_h2k' refers to HITRAN 2000 release, suffix '_h96' refers to HITRAN 1996 release, '_h92' refers to HITRAN 1992 release.

rfm_xsc.tar.Z (4.6 Mb) Download all following XSC data as a single compressed tar file

TIPS Data

This are the HITRAN Partition Sum data for each molecule/isotope combined into a single file suitable for input to the RFM using the TPS option (default is to calculate partition sums using polynomial approximations).
rfm.tps.Z
(Compressed) TIPS file for all molecules/isotopes downloaded from the HITRAN ftp site on 13FEB03

IDL Procedures

Various useful IDL procedures (=subroutines) for reading RFM files
rfmrd.pro
Read RFM ASCII spectrum. Also rfmrd2.pro for version which also works with RFM binary spectra.
reaprf.pro
Read atmospheric profile

FORTRAN Programs

Some useful fortran programs for creating/examining HITRAN binary files
hitbin.f (v1.20)
Convert HITRAN ASCII file to binary format.
Note that this is a new version (30JUN09) which should be able to handle the change in the format of the self-broadening paramter in HITRAN 2008 (changed from F5.4 to F5.3) but I haven't yet had a chance to check it myself.
Previous version (v1.10), which definitely works with earlier versions of HITRAN data, is here: [hitbin_v1.10.f]
This is a simplified version of the GENLN2 program HITLIN with the same output format. Compile with "f77 hitbin.f -o hitbin", then dialogue is as follows:
R-HITBIN: Running HITBIN v1.01
Output identifying program and version being run
Input HITRAN ASCII file:
Insert name of file to be converted (plus directory paths), up to C*80
Wavenumber range (cm-1) [<CR>=all]:
Insert lower, upper wavenumber limits to be extracted from ASCII file. Default (pressing the carriage return key) is to use complete range in file.
HITRAN ID#s of gases to use (-ve=exclude) [<CR>=use all]:
Enter a list of indices (see Appendix E) of absorbers to be selected (if all positive) or excluded (if all negative) from ASCII file. Default is to include everything. NB: cannot mix positive and negative indices.
Record Length for binary file (22=DEC, 88=others) [<CR>=88]:
Binary files are created with 88-byte records. However in the FORTRAN OPEN statement, some machines/compilers define the keyword RECL to represent words of 4-bytes, whereas most other machines define RECL to be bytes.
If in doubt, I suggest using the default figure of 88 *but* check the size of the resulting output file. If it is significantly larger than the input ASCII file, then your set-up probably requires 22 instead.
See also RECLEN parameter in the RFM module reclen.inc.
New binary file:
Name of the new binary file to be created, up to C*80
Header for new file (up to 48 chars):
Comment to be inserted in Record#1 of new binary file, default=none.
Merge with existing binary file, filename [=none]:
Existing binary file in which to insert ASCII data (optional). Default is to perform simple ASCII to binary conversion.
Status of new ASCII data (.GE.10) [=10]: (no merge)
Status of new ASCII data (.GE.10) [=increment]: (merge)
Different versions of line transition data in the binary files are identified by a status number (LSTAT) 10 or more. If merging an ASCII file with existing binary data only the transition with the highest status number is preserved and the lower status line is removed from the output file (the GENLN2 program HITLIN acts slightly differently in retaining the lower status line but converting LSTAT to negative). This is to allow partial updates of a binary file. Default for simple ASCII to binary conversion is LSTAT=10, default for merging is to first check through the old binary file for the highest status and then assign new lines with a status of +1 higher.
Remaining outputs are diagnostic information. Successful completion is indicated by the line 
 R-HITBIN: Successful completion
Internally identified failure conditions are indicated by output beginning 
 F-xxxxxx
where xxxxxx is the subroutine in which the error occurred.

exabin.f
Examine binary line data file, eg HITBIN output. Prompts user for filename and either a range of record numbers or a wavenumber range. Output to terminal is then a simple listing of each record (converted back into HITRAN-like format for transition data records).
hitfmt.f
Reformat new HITRAN2004 as old format ASCII files, suitable for conversion to binary file by hitbin. This simply strips out much of the information in the HITRAN 2004 file (C*160 records) and writes the minimal set required for RFM calculations into the old-style C*100 records.

Latest version (30JUN09) fixes a potential bug when reading the HITRAN molecule ID (spotted by Mike Burton - thanks) and allows for the new HITRAN molecule# up to 42.