Figures from ftrue_higuchi.for program and also for Alsop's Figure 13
A. Figures for Higuchi Model

A.1 Exact solutions for R/T coefficients for HIGUCHI model when the "fundamental assumption" is applied.
a. (v_R_11)_Higuchi = (v_R_11)_latave Here, "EXACT" = "DIRECT" in the output.
b. (v_T_11)_Higuchi = (v_T_11)_latave Here, "EXACT" = "DIRECT" in the output.
The "scattering ratio" that you see besides the "Fractional energy" in the lowermost panel is actually the RATIOE(I) that you computed in ftrue_f13_enbal.for
RATIOE(I) = ( OMEGA(I) * K_LEFT(I) * LSHLFT(I) * RLFRHD(I)**2 +
1 OMEGA(I) * KRIGHT(I) * LSHRHT(I) * TLFRHD(I)**2 )
2 /
3 ( OMEGA(I) * K_LEFT(I) * LSHLFT(I) * 1.000D+00**2 )

A.2 Other parameters for HIGUCHI model when the "fundamental assumption" is applied. In the upper panel, the wavenumber k is always increasing as frequency increases. In the second and third panels, the subtraction of LEFT side and RIGHT will equal to parameter in the CENTER.

A.3 Ratio
B. Figures for Alsop's Figure 13

A.1 Exact solutions for R/T coefficients from laterally-averaged Alsop's Figure 13. The method where we used
( v_R_11 )_latave = (1 - a)/(1 + a ) and ( v_T_11 )_latave = 2/(1 + a)
produces the same results.
The "scattering ratio" that you see besides the "Fractional energy" in the lowermost panel is actually the RATIOE(I) that you computed in ftrue_f13_enbal.for
RATIOE(I) = ( OMEGA(I) * K_LEFT(I) * LSHLFT(I) * RLFRHD(I)**2 +
1 OMEGA(I) * KRIGHT(I) * LSHRHT(I) * TLFRHD(I)**2 )
2 /
3 ( OMEGA(I) * K_LEFT(I) * LSHLFT(I) * 1.000D+00**2 )

A.2 Other parameters. In the upper panel, the wavenumber k is always increasing as frequency increases. In the second and third panels, the subtraction of LEFT side and RIGHT will equal to parameter in the CENTER.

A.3

A.4 | Direct - Exact | .