static double dtExt(double y, double jsd) { double dy = (y - 1820) / 100; return -20 + jsd * dy * dy; }