USAGE:
             love-list [options] -f input-file
PURPOSE:

             Picks and lists Maxwell results, also mode summation.

Let # denote a number and ¤ denote a string.

OPTIONS:
 -b          - binary output
 -B0         - tidal (stress-free) Love-numbers
 -B1         - normal stress Love-numbers                                [-B1 is default]

 -C[C][#,#]  - copy (with comment) to output. -f and -o must precede.
               #,#: Degree range from to.                                      [0,999999]
                 The comment is the first line of the input mode file.

 -U          - unnormalize (normally listed = n l_n, n k_n)
 -l#[,#]     - Sph. harm. degree loop: from, to
 +[e][v][f]  - list elastic, viscous, fluid cases.

 -t[#[,#,#]] - Time loop limits for Impulse response (kyr)                            [0]
 -T[#[,#,#]] - Time loop limits for   Step  response (kyr)                            [0]
                 Specify begin,end,step.
               Elastic effect will be added if option +e is active.
 -Ti -Tir[p]   Special case -Ti requests Infinite time.
               Special case -Tir will also print mode sum residual.
               Special case -Tirp will print the partial sums term by term
 -Tir#         Add limiting parameter to print only big discrepancies.
 -logT       - Time loop species ₁₀log(time)
 -c#[,#]     - Minimum mode [H] strength to qualify for list                     [1.d-10]
               Minimum mode strength is also known as EPS
                 Optional: Maximum strength; only used to limit
                 the list of the modes (not used in mode sums).

 -cg{H|K|L|S}#   Select only the # greatest values of mode
                 strength or Laplace-s. Further selection may ensue.

 -H -K -L    - Select the component for modes scanned under -c#[#] option.
 -s#         - Maximum Laplace-s to qualify for list                                [1.d5]
 -H# -K# -L# - Print the largest # modes only
 -S#         - Print the slowest # modes only; -c# is in effect.
 -/s         - Print mode strength j divided by s_j
 -M          - Instead of mode strength, print mode magnitude minus log(EPS)  
               (cf. option -c)
 -M/s        - ... print mode magnitude j divided by s_j minus log(EPS)  
 -M[/s]q[o]  - ... print also/(o: only) magnitude range

 -g          - compute gravity-uplift ratio
 -f filename - input,  default file = prem.l120.um1.lm2.all
 -o filename - output, default is                                                 [stdout]
 +o filename - output, append. Blanks before filename are permitted

 +O¤ -O¤ t:¤ - -Oopenf-text with `_´ replacing ` ´, e.g.
                 -O22_B_Green/loadf.lov
               For reading elastic Love number files.
               Specify +O¤ to include the asymptotic values.
     t:¤ r:¤ - t:tag  for the data section, e.g.  t:L00EGBC
                 Specify r:¤ to list only the real part.

FORMAT of -Tirp output:
    n ____mark____    column  m    part.sum     p.s.-fluid     strength         s        strength/s
    2 Mode_partsum       1    1 -1.574503E+00  1.219400E+00 -1.528138E-01  2.618812E-01 -5.835236E-01
    2 Mode_partsum       1    2 -2.117496E+00  6.764063E-01 -1.518893E-02  2.797256E-02 -5.429939E-01

n: sph.harm. degree
m: mode index
column 1,2,3 designates Love number H,L,K

EXAMPLES:

 Example with Farrell type Love number input:
    love-list +e -l1,10000 -f prem.l120.um1.lm2.all -O22_B_Green/loadf.lov t:L00EGBC

 Example for mode discrepancy:
    love-list -c1.d-16 -Tir -l2,180 -f ../PAO/prem.l96_um0.5_lm10_n1-256


 Example for mode discrepancy convergence using a sorted mode list,
 sorted by strength/s of Love h:
    moded @ moded.ins :SRT; love-list -c1.d-16 -Tirp -l2,2 -f prem5.tst

.bye