Original Gravity using Salinity Refractometer

#!/usr/bin/env python3

# ****************************************************************************
# Inputs: Original Salinity and Final Salinity (both in % solute by weight)
# Output: Final ABV
# Produced 2026-06-02
# Author: Sachin (Chino) Darji, adopted from other code.)
# Original author: unknown, from 2016-09-24
# ****************************************************************************

def salinity2brix(salinity):
    Brix = 1.187*salinity
    return Brix

OS = int(input('What was your original salinity reading?'))/10  # Original Salinity
FS = int(input('What is your final salinity reading?'))/10      # Final Salinity
OE = salinity2brix(OS)                                          # Original Extract (Brix)
AE = salinity2brix(FS)                                          # Apparent Extract (Brix)
OG = (OE / (258.6-((OE / 258.2)*227.1))) + 1
FG = 1.001843 - 0.002318474*OB - 0.000007775*OB*OB - 0.000000034*(OE**3) + 0.00574*AE + 0.00003344*(AE**2) + 0.000000086*(AE**3)
abv = (OG - FG) *131.25100
print("your ABV is ~", abv)
print("Your original extract (Brix) was: ", OE,"\nGiving an original gravity of ", OG, "\nYour apparent extract (Brix) was: ", AE, "\nGiving a final Gravity of ", FG)