------------------------- TEST INPUT FILE FOR green --------------------------

------- any line not starting with any of the recognized keywords is ignored
-------     keywords	PARAMETERS:

----------- window >>	EMIN EMAX (energy range where spectrum is computed)
----------- epsilon >>	EPSILON (imaginary part of energy for broadening)
----------- delta >>	DELTA (energy step: must be approx. EPSILON/4)
----------- tiers >>	N (number of tiers: N=0: only bright is included)
----------- tier_size >> NTIER SMAX (set to SMAX maximum size of tier NTIER 
-----------                          and of following tiers)
----------- bright >>	WEIGHT v1 v2 v3 .. (quantum numbers of bright state,
-----------                                implicitly padded with zeroes)
----------- coefficients >> FILENAME (where to find all potential parameters)
----------- frequencies >> FILENAME (where to find only the frequencies: this
-----------                         is in principle unnecessary, as the
-----------                         frequencies may be put in coefficients)
----------- turnoff $ KEYS where KEYS are zero or more of the following:
--------------------- zpe (set zero-point energy hbar omega/2 equal to zero)
--------------------- sum (max instead of sum in c.c.s. formula)
--------------------- square (do not perform square in c.c.s formula)
--------------------- spectrum (stop when hamiltonian matrix is generated)
--------------------- static (do dynamic Green's matrix allocation)
--------------------- alarm (skip beep at execution end)
----------- show $ KEYS where KEYS are zero or more of the following:
--------------------- tiers (show all states in the tiers basis)
--------------------- coefficients (show potential coefficients)
--------------------- weights (show cumulative coupling strengths of states)

------- the keywords may be called in almost any order


------------------------- Physical parameters --------------------------------
bright >> 1.0 2 0

tiers >>   5

tier_size >> 1 500
-- tier_size >> 3 20

------------------------- Parameters to set the calculation window -----------
window >> 1 9
epsilon >> 0.04
delta >>   0.01


------------------------- Options --------------------------------------------
turnoff $
-- zpe spectrum static square sum alarm

show $ tiers coefficients weights

------------------------- Data file names -----------------------------------

coefficients >>  test/test_phi2.txt test/test_phi3.txt
// all coefficients in form i1 i2 i3... param
// param must be all in the same units (e.g. all in cm^-1, or all in eV...)
// for example:
// 3 3 1234.5 means that the energy (h nu) of mode 3 is 1234.5
// 1 3 4 4 -233.1 means that the off diagonal coeff of q1 q3 q4 q4 is -233.1
