s_10to10.f90: Description of the input/output structure

Purpose: Program to generate synthetic synthetic spectra using the 10to10
         line list at specified temperatures. The spectrum is given either 
         as a stick diagram (cm/molecule), or in the form of cross-sections  
         with a Gaussian/Doppler profile on a grid of frequency points.
         (cm^2/molecule).  The partition function value can be also recomputed, 
         if it is not provided by the user.

Input structure:
 
    ! number of Transition files to read
    !
    read*,nfiles
 
    ! names of Transition files 
    !
    do i = 1,nfiles
      read*,intfilename(i)
    enddo

    ! name of the Energy file
    !
    read*,enrfilename

    ! temperature (K), partition function; 
    ! if partfunc <0 the partition function will be computed using the energies
     provided
    !
    read*,temp,partfunc

    ! frequency range (cm-1) and number of grid points; the latter is not
    referenced for proftype='stick'
    !
    read*,freql,freqr,npoints

    ! form of spectra (profiles): 'gauss', 'doppler', 'bin' or 'stick' 
    !
    read*,proftype
    !
    ! type of spectra:
    ! 'absorption': absorption coefficients (cm/mol) 
    ! 'emission': emission coefficients (Ergs/mol/Sr)
    ! 'partfunc': partition function for a series of temperatures together with
    derivatives
    !
    read*,specttype

    ! HWPar = HWHM (cm-1) if ProfType='gauss'
    ! HWPar = molecular mean molar mass (g) if ProfType='doppler'
    ! HWPar = a print our threshold for absorption coefficient (cm/mol) if 
    ProfType='stick'
    !
    read*, HWPar

    ! intensity cutoff (cm/mol or Ergs/mol/Sr) to spead up evaluation of 
    intensity, can be set negative for a zero-cutoff.
    ! not referenced for 'stick'
    !
    read*,thresh
    !

    !Output structure: 

    ! for 'gauss', 'doppler', 'bin'
    !
    print('(2(1x,es16.8))'), (Freq(i), Intens(i), i = 1, Npoints)
 
    ! for 'stick' the typical record  has the following format:

!     nu(cm-1)         I(cm/mol)   G_f   J_f  v1 v2 L2 v3 L3 M3 v4 L4 M4 Gvib  Kf tf Gr     E_f       <- G_i  J_i   v1 v2 L2 v3 L3 M3 v4 L4 M4 Gvib  Ki ti Gr     E_i
   3.50079450E+01  1.08612338E-24  F1    2    0  0  0  0  0  0  1  1  1 F2     1   0 F2     1351.3129 <- F2    1    0  0  0  0  0  0  1  1  1 F2     1   1 F1     1316.3049

    ! where i,f stand for initial and final states, respectively;
    ! v1,v2^L2,v3^{L3,M3},v4^{L4,M4} are the normal mode quantum numbers;
    ! G, Gv, and Gr are total, vibrational and rotational symmetries;
    ! E_i and E_f are the rovibrational term values;
    ! J,K,t (tau) are the rotational quantum numbers: angular momentum, projection of 
    ! the angular momentum onto the z-axis, and rotational parity.