A couple of thoughts on Raymer's & Shevell's e estimation methods as I understand them...
Raymer uses a semi-emperical approach to estimate K from leading edge suction. The issues with this approach are: a) you use a look up chart for S based on Cl but this is based on knowing the design Cl of the wing, b) the technique was geared more toward thin wing, low aspect ratio jets, c) it doesn't account for the entire plane including the fuselage.
Shevell's approach is also a semi-emperical approach based on data obtained from data from Douglas aircraft and accounts for fuselage affects. However in the past when I've run the calculations using Shevell's estimating method it estimates e much higher than what the P-51 data posted by Pyro shows.
I believe both of their estimation methods are intended for estimating initial designs. They are great for that purpose. However to reverse engineer and accurately estimate lift dependent viscous drag for an airplane I think we're stuck with getting into nasty Navier-Stokes world of CFD...
...or go find some good drag polars for specific aircraft based on flight tests or wind tunnel tests

.