`L` <- c(1, 0, 0, 1, 1, 0, 2, 0, 0, 2, 0, 0, 0, 1, 0, 2, 2, 0, 1, 1, 0, 1, 0, 0, 0, 0, 2, 0, 1, 0, 0, 0, 1, 3, 0, 2, 1, 0, 1, 0, 1, 1, 1, 1, 2, 2, 0, 2, 1, 1, 1, 0, 0, 0, 3, 3, 0, 1, 0, 0, 1, 1, 0, 0, 3, 1, 0, 1, 0, 1, 1, 0, 0, 0, 1, 1, 0, 0, 0, 2, 1, 0, 0, 1, 0, 4, 0, 0, 0, 1, 0, 1, 1, 0, 0, 2, 1, 0, 2, 2, 0, 0, 1, 1, 4, 4, 0, 0) `nao` <- c(-0.9, -0.25, -1.355, -1.33, 0.44, -1.505, -1.715, -0.33, 1.375, -1.135, -1.285, 0.605, 0.36, 0.705, 1.38, -2.385, -1.875, -0.39, 0.77, 1.605, -0.43, -1.12, 1.575, 0.44, -1.32, -0.54, -1.49, -1.815, -2.395, 0.305, 0.735, -0.79, -1.07, -1.085, -0.54, -0.935, -0.79, 1.4, 0.31, -1.15, -0.725, -0.15, -0.64, 2.04, -1.18, -0.235, -0.07, -0.5, -0.75, -1.45, -0.235, -1.635, -0.46, -1.855, -0.925, 0.075, 2.9, -0.82, -0.17, -0.355, -0.17, 0.595, 0.655, 0.07, 0.33, 0.395, 1.165, 0.75, -0.275, -0.7, 0.88, -0.97, 1.155, 0.6, -0.075, -1.12, 1.48, -1.255, 0.255, 0.725, -1.23, -0.76, -0.38, -0.0150000000000000, -1.005, -1.605, 0.435, -0.695, -1.995, 0.315, -0.385, -0.175, -0.47, -1.215, 0.78, -1.86, -0.035, -2.7, -1.055, 1.21, 0.6, -0.71, 0.425, 0.155, -0.525, -0.565, -0.315, -1.36 ) `sst` <- c(-0.0533996666666676, -0.0793996666666672, -0.154733000000000, -0.523399666666667, -0.498066333333334, -0.347733000000000, -0.204066333333334, -0.331066333333333, -0.159066333333333, -0.205399666666666, -0.272066333333334, -0.183399666666668, -0.378399666666668, -0.494399666666668, -0.301399666666667, 0.0712669999999988, -0.103733000000001, -0.341066333333334, -0.403399666666668, -0.238399666666667, -0.342399666666668, -0.232399666666667, -0.375733000000001, -0.405066333333333, -0.241733, -0.159733000000001, 0.0459336666666654, 0.103267000000000, 0.0372669999999988, -0.127399666666667, 0.152600333333332, 0.213933666666667, 0.292600333333333, 0.259600333333334, 0.108600333333333, -0.0367330000000005, 0.146600333333333, 0.450933666666667, 0.308600333333333, 0.201600333333334, 0.0212669999999993, 0.329933666666667, 0.190267, 0.0606003333333322, 0.415600333333334, 0.107600333333333, -0.0527330000000001, -0.0640663333333341, -0.00373300000000043, 0.113933666666666, 0.00460033333333363, 0.298266999999998, 0.403267, 0.261267000000000, 0.0162669999999991, 0.257600333333334, -0.0353996666666667, 0.234267000000000, 0.220933666666666, 0.131600333333334, 0.331933666666665, 0.0766003333333328, 0.0589336666666671, -0.055733, -0.180066333333334, -0.162066333333334, 0.0572669999999983, -0.0503996666666673, -0.0860663333333336, -0.0553996666666676, -0.0797330000000009, -0.264733000000001, -0.239733000000001, -0.0727330000000006, -0.403399666666666, -0.198733000000001, -0.127399666666666, -0.101066333333334, -0.105732999999999, 0.00193366666666723, 0.0932669999999995, 0.00426699999999985, -0.204066333333334, -0.0683996666666668, -0.132733000000000, -0.127732999999999, -0.136732999999999, 0.308933666666667, 0.0479336666666665, 0.103600333333333, 0.244933666666666, 0.00693366666666613, -0.204066333333333, -0.063066333333334, 0.0129336666666664, 0.267266999999999, 0.102600333333333, 0.246933666666666, 0.603267, 0.344933666666666, 0.225267, 0.420933666666667, 0.280600333333332, 0.614600333333333, 0.448267, 0.551933666666667, 0.560267, 0.297600333333333) `soi` <- c(-0.996666666666667, -1.15666666666667, -1.20666666666667, 0.32, -0.0766666666666667, -0.83, 1.33, -0.516666666666667, 0.863333333333333, 0.383333333333333, 1.13, -1.19666666666667, -0.77, -0.963333333333333, -1.34666666666667, 1.19666666666667, 0.763333333333333, 2.61, -0.703333333333333, -0.913333333333333, 0.0733333333333333, 0.183333333333333, 0.243333333333333, -1.44, 0.75, -1.13333333333333, -0.186666666666667, -0.456666666666667, 0.77, 0.186666666666667, -0.216666666666667, -0.383333333333333, -0.816666666666667, 0.0366666666666667, -0.9, 0.396666666666667, -0.293333333333333, -0.05, 1.1, -0.923333333333333, -1.95333333333333, -1.77333333333333, 0.606666666666667, 0.67, -0.226666666666667, 0.69, -1.2, 0.47, -0.256666666666667, -0.03, 1.15, -1.12, -0.166666666666667, -1.14, 0.386666666666667, 1.43, 1.01, -0.726666666666667, -0.06, -0.116666666666667, 0.326666666666667, -0.26, 0.58, -0.83, 1.29, -1.33333333333333, -0.126666666666667, 0.236666666666667, -0.303333333333333, -1.05333333333333, 0.743333333333333, 1.56, -1.24333333333333, 1.1, 0.813333333333333, 1.95, -0.873333333333333, -1.25333333333333, -0.193333333333333, -0.316666666666667, -0.313333333333333, 0.07, -2.32666666666667, 0.293333333333333, -0.113333333333333, -0.0333333333333333, -0.28, -1.16333333333333, 1.54333333333333, 0.103333333333333, -0.553333333333333, -1.42, -0.66, -1.36, -1.77, -0.0566666666666667, 0.416666666666667, -1.85, 0.96, 0.283333333333333, 0.773333333333333, -0.426666666666667, -1.02333333333333, -0.263333333333333, -0.543333333333333, 0.163333333333333, -1.29, 0.19) `ssn` <- c(8.3, 0.6, 7.6, 11.1, 30.1, 55, 56.1, 85, 86.9, 38.8, 26.2, 4, 9.5, 1.2, 12.7, 49.5, 45.1, 129.4, 79.9, 54.7, 36.3, 17.8, 4.7, 13.2, 25.1, 60.2, 60.8, 68.4, 89.7, 34.4, 32.1, 19, 4, 5.1, 4, 42.1, 76, 100.7, 89.6, 112.6, 66.5, 65.9, 17.2, 10, 14.3, 34.9, 94.4, 169.4, 143.3, 145.3, 51.3, 83.1, 28.2, 19.3, 1.5, 42.7, 173.2, 235.8, 201.2, 145.2, 127.2, 63.6, 51.3, 38.8, 4.7, 16.8, 50.2, 76.8, 117.2, 91.3, 99.5, 50.2, 64, 59.3, 40.2, 13.9, 13.5, 44, 138.2, 188.4, 155, 167.3, 118.8, 50.3, 15.7, 3.9, 3.8, 33.9, 120.1, 176.7, 125.2, 125.3, 63.9, 22.4, 25.7, 11.8, 1.6, 51.3, 92.9, 71.5, 109.7, 150.7, 109.6, 48.7, 27.7, 21.9, 14.4, 2.4) `y` <- c(9.02563836163789, 9.03668418495465, 9.04259887197271, 9.0712255272393, 9.08454649181778, 9.08620066889547, 9.10356518572833, 9.11062839344156, 9.12203672681062, 9.12259671762722, 9.1338723654253, 9.138574701368, 9.15209858628223, 9.163552039194, 9.2224630017083, 9.23194793878745, 9.24147233666371, 9.25700822981192, 9.26400164483112, 9.2708452301019, 9.29268236663467, 9.34121517046173, 9.34830486304816, 9.38002006917346, 9.38405987061621, 9.41377564947338, 9.4155603368996, 9.442866440533, 9.45253033578189, 9.48669852938722, 9.49207834124204, 9.50538538802443, 9.51209262846602, 9.55280561463825, 9.55623410642682, 9.55870814967522, 9.57115571565614, 9.57633072264787, 9.57871653500318, 9.58083622827714, 9.60094804473258, 9.60714376726643, 9.61517600588662, 9.62400150482549, 9.66774860500008, 9.71654228699888, 9.74462785755593, 9.75733647762207, 9.75919674683235, 9.77849398387195, 9.79020942857908, 9.8020862282904, 9.81806678121376, 9.83168121558455, 9.83682204246685, 9.86016659307642, 9.87484489558616, 9.98448217601038, 9.99204005869188, 10, 10.0063790272329, 10.0621276657773, 10.0819820004496, 10.1292586441071, 10.1293177628760, 10.1642208732027, 10.1689801481721, 10.1785399555403, 10.190724112395, 10.2121089201684, 10.2426270645900, 10.2511744106878, 10.3138672203692, 10.3625703238466, 10.3663401552462, 10.3794481605134, 10.5030138625389, 10.5037720717920, 10.5512921071592, 10.5714023462808, 10.7350975551129, 10.7991305438384, 10.9084850188787, 10.9619676380509, 11.1453760041235 ) `Yr` <- c(2002, 1985, 1974, 1999, 1998, 1940, 1918, 1936, 1909, 1909, 1933, 1995, 1906, 1913, 1954, 1945, 1980, 1916, 1941, 1904, 1947, 1979, 2005, 1985, 1906, 1964, 1942, 1915, 1975, 1991, 1935, 1921, 1943, 1998, 1926, 1948, 1950, 1985, 1933, 1955, 2003, 1967, 1957, 1985, 1964, 1903, 1935, 1970, 1932, 1955, 1996, 1995, 1964, 1999, 1916, 1983, 2004, 2004, 1933, 2005, 1945, 1979, 1944, 1961, 1949, 1919, 1947, 1954, 2004, 2004, 1989, 1955, 2005, 1965, 1954, 1969, 1928, 1960, 1944, 1938, 1992, 1915, 2005, 1900, 1926) `u` <- 9 `Yr0` <- 1900 `M` <- 108 `N` <- 85 `C` <- 1e-04 `epsilon` <- 1e-12