# dissolution of nitrogen
a1 = -59.6274
a2 = +85.7661
a3 = +24.3696
b1 = +0.051580
b2 = +0.026320
b3 = +0.0037252
Sa = 65
Tw = seq(12,18,.5)
K1 = a1+a2*(100/Tw)+a3*log2(Tw/100)
K2 = a1+a2*(100/Tw)+a3*log2(Tw/100)+Sa*(b1+b2*(Tw/100)+b3*(Tw/100)^2)
plot(Tw,K1, col='blue', type='l')
lines(Tw, K2, col="red",lty=2)