CubicSpline cs = new CubicSpline(concn85, lvisc);
for(int j=0; j<lengthC85; j++){
array[i][j] = cs.interpolate(conc85[j]);
}
}
BiCubicSpline bcs = new BiCubicSpline(temp85, conc85, array);
ret = Math.exp(bcs.interpolate(temperature, gPerL));
}
else{
if(temperature>100.0){
throw new IllegalArgumentException("Temperature, " + temperature + ", out of range");
}
else{
int lengthT100 = length60;
int lengthC100 = amC.indexOf(771.9);
double[] temp100 = amT.subarray_as_double(0, lengthT100-1);
double[] conc100 = amC.subarray_as_double(0, lengthC100);
lengthC100++;
double[][] array = new double[lengthT100][lengthC100];
for(int i=0; i<lengthT100; i++){
double[] lvisc = {Math.log(zero[i]), Math.log(forty[i]), Math.log(sixty[i])};
CubicSpline cs = new CubicSpline(concn100, lvisc);
for(int j=0; j<lengthC100; j++){
array[i][j] = cs.interpolate(conc100[j]);
}
}
BiCubicSpline bcs = new BiCubicSpline(temp100, conc100, array);
ret = Math.exp(bcs.interpolate(temperature, gPerL));
}
}
}
else{
if(temperature==20 && gPerL<=1208.2){