Equation Eq. 3 of P.D.T Huibers, Applied Optics, Vol. 36, No. 16, pp. 3785-3787 (1997).

$$ n_{H2O} = 1.31279 + \frac{15.762}{\lambda} - \frac{4382.0}{\lambda^{2}} + \frac{1.1455x10^{6}}{\lambda^{3}}$$

Python Code

# Define the index of refraction function for water
def water_index(wavelength):
n_H2O = 1.31279 + 15.762/wavelength - 4382.0/(wavelength**2) + 1.1455e6/(wavelength**3)
return(n_H2O)

JavaScript Code

var n_H2O=[];
// Define the index of refraction function for water
function water_index(wavelength){
//creates the array for n_H20
for (var i =0;i < wavelength.length+1;i++){
n_H2O[i] = 1.31279 + 15.762/wavelength[i] - 4382.0/(wavelength[i]*wavelength[i]) + 1.1455e6/(wavelength[i]*wavelength[i]*wavelength[i]);
}
return(n_H2O)
}

Equation Eq. 3 of P.D.T Huibers, Applied Optics, Vol. 36, No. 16, pp. 3785-3787 (1997).

$$ n_{air}-1 = \frac{0.05792105}{238.0185- \frac{1x10^{6}}{\lambda^{2}}} + \frac{0.00167917}{57.362-\frac{1x10^6}{\lambda^{2}}} $$

Python Code

# Define the index of refraction function for air
# Note: wavelength is supposed to be in micrometers
def air_index_minus_one(wavelength):
term1 = 0.05792105/(238.0185 - 1.0e6/(wavelength**2))
term2 = 0.00167917/(57.362 - 1.0e6/(wavelength**2))
return(term1+term2)

JavaScript Code

var term1, term2;
var term=[];
// Define the index of refraction function for air
// Note: wavelength is supposed to be in micrometers
function air_index_minus_one(wavelength){
//create term1 and 2 arrays
for (var i=0;i < wavelength.length+1;i++){
term1 = 0.05792105/(238.0185 - 1.0e6/(wavelength[i]**2));
term2 = 0.00167917/(57.362 - 1.0e6/(wavelength[i]**2));
term[i] = term1+term2;
}
return(term)
}