SPC SHIP equation is incorrect about sharppy HOT 3 CLOSED

For the purpose of records...

Greg used the SHIP values from the SARS database and then used the SPC equation to calculate SHIP based on the parameters in the SARS hail file. This graph he made demonstrates the spread between SHIP in the database and SHIP calculated using it's own variables.

It obviously does not stay along the line as it should, leading the conclusion to be that the equation documented on the SPC webpage is in error, or the database is in error.

from sharppy.

August 23, 2014...wrote Ryan Jewell regarding these issues and sent him the figure attached to this thread by Kelton.

from sharppy.

Ryan Jewell wrote back a few hours later with some corrections to the SHIP equation. He mentioned that the denominator was incorrect. It is not 44,000,000 but rather 42,000,000:
"""
in addition to 42,000,000 in denominator, apply these checks before calculating and replot:

if 6km shear > 27 m/s, set to 27.
if 6km shear < 7 m/s, set to 7.
if mumixr > 13.6, set to 13.6
if mumixr < 11, set to 11
if 500 mb temp > -5.5, set to -5.5

Then calculate. THEN...

if mucape < 1300, ship = ship_(mucape/1300)
if 7-5 lr < 5.8, ship = ship_(lr75/5.8)
if freezing level < 2400m, ship = ship*(fzl/2400)

Then you'll have the final number nsharp has.
"""

Upon changes to the code, I loaded up a couple of different soundings (including the test sounding seen below) and the SHIP values compared perfectly. There may be some differences at some point, as I'm uncertain if the virtual temperature correction needs to be found in order to find the freezing level, but the consistency between the numbers I saw were spot on.

from sharppy.