% by Nirmal Jayaram, 4/27/07 % Stanford University % nirmalj@stanford.edu % % Chiou and Youngs attenuation equation, 2007 % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % % Input Variables % % T: Period (0.01 - 10s) % M: Magnitude % rrup: Closest distance to the rupture plane (km) % rjb: Joyner-Boore distance to the rupture plane (km) % delta: Dip % w: Rupture width (km) % ztor: Depth to top of rupture (km) % lambda: Rake angle % Vs30: Average shear velocity for top 30 m (m/s) % % Output variables % Sa: median spectral acceleration prediction % sigma: Intra-event standard error % tau: Inter-event standard error % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function [Sa sigma tau] = CY_2007_nga(T, M, rrup, rjb, delta, w, ztor, lambda, Vs30) ; % Coefficients period = [0.010 0.020 0.022 0.025 0.029 0.030 0.032 0.035 0.036 0.040 0.042 0.044 0.045 0.046 0.048 0.050 0.055 0.060 0.065 0.067 0.070 0.075 0.080 0.085 0.090 0.095 0.100 0.110 0.120 0.130 0.133 0.140 0.150 0.160 0.170 0.180 0.190 0.200 0.220 0.240 0.250 0.260 0.280 0.290 0.300 0.320 0.340 0.350 0.360 0.380 0.400 0.420 0.440 0.450 0.460 0.480 0.500 0.550 0.600 0.650 0.667 0.700 0.750 0.800 0.850 0.900 0.950 1.000 1.100 1.200 1.300 1.400 1.500 1.600 1.700 1.800 1.900 2.000 2.200 2.400 2.500 2.600 2.800 3.000 3.200 3.400 3.500 3.600 3.800 4.000 4.200 4.400 4.600 4.800 5.000 5.500 6.000 6.500 7.000 7.500 8.000 8.500 9.000 9.500 10.000]; c1 = [-1.2686 -1.2474 -1.2308 -1.2064 -1.1716 -1.1622 -1.1432 -1.1136 -1.1032 -1.0598 -1.0341 -1.0080 -0.9961 -0.9845 -0.9606 -0.9363 -0.8686 -0.8090 -0.7487 -0.7296 -0.7065 -0.6661 -0.6329 -0.5993 -0.5678 -0.5408 -0.5207 -0.4899 -0.4535 -0.4365 -0.4350 -0.4487 -0.4658 -0.4782 -0.4933 -0.5191 -0.5455 -0.5698 -0.6372 -0.7017 -0.7267 -0.7519 -0.8115 -0.8394 -0.8667 -0.9306 -0.9978 -1.0268 -1.0589 -1.1152 -1.1649 -1.2189 -1.2742 -1.2996 -1.3257 -1.3752 -1.4201 -1.5247 -1.6343 -1.7303 -1.7640 -1.8214 -1.8933 -1.9712 -2.0411 -2.1069 -2.1716 -2.2326 -2.3433 -2.4313 -2.5399 -2.6311 -2.7212 -2.8125 -2.8974 -2.9774 -3.0531 -3.1249 -3.2583 -3.3802 -3.4373 -3.4922 -3.5960 -3.6926 -3.7829 -3.8678 -3.9084 -3.9478 -4.0235 -4.0953 -4.1636 -4.2288 -4.2910 -4.3506 -4.4077 -4.5412 -4.6630 -4.7750 -4.8788 -4.9754 -5.0657 -5.1506 -5.2306 -5.3063 -5.3820]; c1a = [0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.1000 0.0994 0.0987 0.0978 0.0967 0.0960 0.0951 0.0931 0.0906 0.0828 0.0744 0.0652 0.0617 0.0547 0.0435 0.0323 0.0214 0.0112 0.0018 -0.0068 -0.0210 -0.0328 -0.0432 -0.0525 -0.0602 -0.0661 -0.0707 -0.0746 -0.0784 -0.0821 -0.0886 -0.0932 -0.0948 -0.0961 -0.0982 -0.0999 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000]; c1b = [-0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2554 -0.2543 -0.2527 -0.2511 -0.2495 -0.2478 -0.2442 -0.2405 -0.2385 -0.2366 -0.2327 -0.2308 -0.2288 -0.2250 -0.2211 -0.2192 -0.2173 -0.2135 -0.2098 -0.2061 -0.2025 -0.2007 -0.1990 -0.1955 -0.1920 -0.1839 -0.1768 -0.1705 -0.1686 -0.1651 -0.1603 -0.1562 -0.1525 -0.1491 -0.1461 -0.1433 -0.1384 -0.1342 -0.1306 -0.1274 -0.1247 -0.1221 -0.1199 -0.1177 -0.1158 -0.1140 -0.1107 -0.1080 -0.1067 -0.1056 -0.1036 -0.1019 -0.1006 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000 -0.1000]; c2 = [1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060 1.060]; c3 = [3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450 3.450]; cn = [2.996 3.292 3.352 3.429 3.501 3.514 3.533 3.551 3.555 3.563 3.563 3.561 3.559 3.557 3.553 3.547 3.531 3.513 3.493 3.484 3.471 3.448 3.423 3.397 3.369 3.341 3.312 3.255 3.199 3.145 3.129 3.093 3.044 2.997 2.952 2.910 2.870 2.831 2.760 2.692 2.658 2.626 2.564 2.533 2.505 2.449 2.397 2.372 2.348 2.303 2.261 2.222 2.185 2.167 2.150 2.118 2.087 2.017 1.957 1.904 1.887 1.855 1.812 1.773 1.737 1.704 1.675 1.648 1.605 1.572 1.546 1.526 1.511 1.498 1.489 1.481 1.474 1.470 1.463 1.458 1.456 1.456 1.455 1.456 1.457 1.458 1.459 1.461 1.463 1.465 1.468 1.470 1.473 1.475 1.478 1.483 1.488 1.492 1.496 1.498 1.499 1.500 1.501 1.501 1.502]; cm = [4.184 4.188 4.183 4.173 4.159 4.156 4.148 4.138 4.135 4.123 4.117 4.112 4.110 4.108 4.104 4.101 4.094 4.089 4.087 4.086 4.086 4.086 4.087 4.090 4.094 4.099 4.103 4.114 4.128 4.142 4.146 4.156 4.172 4.187 4.202 4.217 4.232 4.248 4.276 4.304 4.318 4.332 4.358 4.371 4.384 4.409 4.432 4.444 4.456 4.477 4.498 4.517 4.536 4.545 4.554 4.571 4.588 4.627 4.663 4.696 4.707 4.728 4.757 4.785 4.811 4.836 4.860 4.882 4.925 4.964 5.001 5.037 5.070 5.102 5.133 5.162 5.191 5.217 5.269 5.317 5.339 5.361 5.401 5.439 5.474 5.507 5.523 5.538 5.569 5.598 5.625 5.652 5.678 5.703 5.728 5.786 5.840 5.892 5.942 5.989 6.034 6.077 6.117 6.156 6.193]; c4 = [-2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10 -2.10]; c4a = [-0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50 -0.50]; crb = [50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50 50]; c5 = [6.1600 6.1580 6.1580 6.1590 6.1600 6.1600 6.1600 6.1600 6.1600 6.1600 6.1600 6.1600 6.1600 6.1600 6.1600 6.1600 6.1600 6.1570 6.1510 6.1490 6.1450 6.1370 6.1280 6.1180 6.1070 6.0940 6.0820 6.0550 6.0280 6.0000 5.9920 5.9720 5.9460 5.9200 5.8960 5.8730 5.8510 5.8300 5.7930 5.7580 5.7420 5.7260 5.6950 5.6800 5.6650 5.6360 5.6080 5.5950 5.5810 5.5550 5.5310 5.5070 5.4840 5.4730 5.4630 5.4440 5.4250 5.3820 5.3440 5.3110 5.3010 5.2820 5.2560 5.2330 5.2110 5.1920 5.1730 5.1550 5.1230 5.0960 5.0720 5.0520 5.0370 5.0240 5.0160 5.0100 5.0060 5.0030 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000 5.0000]; c6 = [0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4893 0.4892 0.4890 0.4886 0.4882 0.4878 0.4872 0.4867 0.4861 0.4849 0.4835 0.4822 0.4818 0.4810 0.4797 0.4784 0.4772 0.4760 0.4749 0.4739 0.4721 0.4707 0.4700 0.4694 0.4682 0.4676 0.4671 0.4660 0.4649 0.4644 0.4639 0.4629 0.4619 0.4609 0.4601 0.4597 0.4593 0.4586 0.4578 0.4564 0.4553 0.4544 0.4541 0.4536 0.4530 0.4525 0.4521 0.4517 0.4513 0.4511 0.4508 0.4506 0.4505 0.4505 0.4504 0.4504 0.4503 0.4503 0.4503 0.4503 0.4502 0.4502 0.4502 0.4502 0.4501 0.4501 0.4501 0.4501 0.4501 0.4501 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500 0.4500]; chm = [3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0]; c7 = [0.0512 0.0512 0.0512 0.0512 0.0512 0.0511 0.0511 0.0510 0.0509 0.0508 0.0507 0.0506 0.0506 0.0505 0.0505 0.0504 0.0502 0.0500 0.0499 0.0498 0.0497 0.0495 0.0494 0.0492 0.0491 0.0490 0.0488 0.0486 0.0484 0.0482 0.0482 0.0481 0.0479 0.0477 0.0476 0.0474 0.0473 0.0471 0.0468 0.0466 0.0464 0.0463 0.0460 0.0459 0.0458 0.0455 0.0453 0.0452 0.0450 0.0448 0.0445 0.0442 0.0439 0.0437 0.0436 0.0432 0.0429 0.0421 0.0412 0.0404 0.0401 0.0395 0.0387 0.0379 0.0372 0.0364 0.0357 0.0350 0.0336 0.0322 0.0308 0.0294 0.0280 0.0266 0.0253 0.0240 0.0226 0.0214 0.0188 0.0165 0.0154 0.0143 0.0124 0.0106 0.0090 0.0076 0.0069 0.0063 0.0052 0.0041 0.0033 0.0025 0.0019 0.0014 0.0010 0.0002 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000]; c9 = [1.0480 1.0830 1.0990 1.1250 1.1570 1.1650 1.1770 1.1940 1.1990 1.2150 1.2210 1.2270 1.2290 1.2320 1.2350 1.2380 1.2440 1.2480 1.2490 1.2500 1.2500 1.2490 1.2480 1.2450 1.2430 1.2390 1.2350 1.2260 1.2150 1.2040 1.2010 1.1930 1.1820 1.1700 1.1590 1.1470 1.1360 1.1240 1.1020 1.0810 1.0700 1.0590 1.0390 1.0290 1.0190 1.0000 0.9813 0.9722 0.9632 0.9458 0.9288 0.9123 0.8963 0.8885 0.8809 0.8658 0.8513 0.8167 0.7846 0.7542 0.7443 0.7256 0.6982 0.6721 0.6471 0.6229 0.5996 0.5771 0.5346 0.4951 0.4584 0.4244 0.3930 0.3639 0.3369 0.3120 0.2889 0.2675 0.2294 0.1965 0.1818 0.1680 0.1431 0.1210 0.1016 0.0846 0.0769 0.0698 0.0569 0.0459 0.0364 0.0283 0.0213 0.0154 0.0104 0.0012 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000]; cy1 = [-0.008040 -0.008113 -0.008155 -0.008230 -0.008351 -0.008387 -0.008454 -0.008564 -0.008599 -0.008754 -0.008833 -0.008906 -0.008945 -0.008984 -0.009050 -0.009121 -0.009288 -0.009433 -0.009561 -0.009603 -0.009656 -0.009733 -0.009782 -0.009808 -0.009809 -0.009788 -0.009753 -0.009627 -0.009460 -0.009269 -0.009205 -0.009054 -0.008832 -0.008610 -0.008398 -0.008183 -0.007976 -0.007776 -0.007397 -0.007044 -0.006877 -0.006715 -0.006408 -0.006263 -0.006123 -0.005859 -0.005614 -0.005497 -0.005386 -0.005175 -0.004980 -0.004799 -0.004632 -0.004552 -0.004477 -0.004332 -0.004199 -0.003903 -0.003652 -0.003435 -0.003368 -0.003246 -0.003078 -0.002929 -0.002795 -0.002674 -0.002564 -0.002464 -0.002287 -0.002138 -0.002010 -0.001898 -0.001802 -0.001717 -0.001643 -0.001578 -0.001519 -0.001467 -0.001381 -0.001311 -0.001281 -0.001255 -0.001209 -0.001172 -0.001142 -0.001117 -0.001106 -0.001096 -0.001079 -0.001065 -0.001052 -0.001042 -0.001032 -0.001024 -0.001016 -0.000999 -0.000987 -0.000977 -0.000968 -0.000961 -0.000955 -0.000950 -0.000945 -0.000941 -0.000937]; cy2 = [-0.007850 -0.007921 -0.007962 -0.008035 -0.008154 -0.008189 -0.008255 -0.008362 -0.008396 -0.008547 -0.008624 -0.008696 -0.008734 -0.008771 -0.008836 -0.008906 -0.009068 -0.009210 -0.009335 -0.009376 -0.009428 -0.009503 -0.009550 -0.009577 -0.009577 -0.009557 -0.009522 -0.009400 -0.009236 -0.009050 -0.008988 -0.008840 -0.008623 -0.008406 -0.008199 -0.007989 -0.007787 -0.007592 -0.007222 -0.006878 -0.006714 -0.006556 -0.006257 -0.006115 -0.005979 -0.005720 -0.005481 -0.005368 -0.005259 -0.005053 -0.004862 -0.004686 -0.004522 -0.004445 -0.004371 -0.004230 -0.004100 -0.003811 -0.003565 -0.003354 -0.003288 -0.003169 -0.003005 -0.002860 -0.002729 -0.002611 -0.002504 -0.002406 -0.002233 -0.002087 -0.001962 -0.001854 -0.001759 -0.001677 -0.001604 -0.001540 -0.001483 -0.001433 -0.001348 -0.001280 -0.001251 -0.001225 -0.001181 -0.001145 -0.001115 -0.001091 -0.001080 -0.001071 -0.001054 -0.001040 -0.001027 -0.001017 -0.001007 -0.001000 -0.000992 -0.000976 -0.000964 -0.000954 -0.000945 -0.000938 -0.000933 -0.000927 -0.000923 -0.000919 -0.000914]; cy3 = [4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00 4.00]; phi1 = [-0.482300 -0.472300 -0.469100 -0.463900 -0.456600 -0.454800 -0.451300 -0.446300 -0.444700 -0.438500 -0.435700 -0.433100 -0.431800 -0.430600 -0.428400 -0.426300 -0.422200 -0.419500 -0.418300 -0.418200 -0.418500 -0.419900 -0.422600 -0.426300 -0.431000 -0.436300 -0.442100 -0.454800 -0.467900 -0.481100 -0.485000 -0.493900 -0.506100 -0.517600 -0.528300 -0.538500 -0.547900 -0.556800 -0.573100 -0.587700 -0.594400 -0.600800 -0.612800 -0.618500 -0.623800 -0.633900 -0.643200 -0.647600 -0.651800 -0.659900 -0.667500 -0.674600 -0.681300 -0.684500 -0.687600 -0.693700 -0.699500 -0.712800 -0.724800 -0.735700 -0.739300 -0.745800 -0.755100 -0.763700 -0.771800 -0.779300 -0.786400 -0.793100 -0.805400 -0.816500 -0.826500 -0.835900 -0.844400 -0.852300 -0.859500 -0.866400 -0.872700 -0.878600 -0.889100 -0.898100 -0.902100 -0.905700 -0.912200 -0.917300 -0.921500 -0.924600 -0.925900 -0.926800 -0.928200 -0.928600 -0.928100 -0.926900 -0.924900 -0.922200 -0.918800 -0.907500 -0.892900 -0.875800 -0.856900 -0.837000 -0.816500 -0.795900 -0.775700 -0.756000 -0.737200]; phi2 = [-0.192800 -0.189500 -0.190900 -0.194100 -0.200600 -0.202600 -0.207100 -0.214400 -0.216900 -0.227500 -0.232900 -0.238300 -0.241000 -0.243700 -0.249100 -0.254400 -0.267100 -0.278900 -0.289500 -0.293400 -0.298700 -0.306600 -0.313200 -0.318700 -0.323100 -0.326600 -0.329300 -0.332800 -0.334300 -0.334300 -0.334100 -0.333300 -0.331400 -0.328800 -0.325700 -0.322400 -0.318600 -0.314700 -0.306000 -0.296700 -0.291900 -0.286800 -0.276500 -0.271300 -0.266000 -0.255400 -0.244900 -0.239600 -0.234400 -0.224000 -0.213900 -0.203900 -0.194300 -0.189600 -0.185000 -0.176100 -0.167600 -0.148000 -0.131100 -0.116400 -0.112000 -0.103900 -0.093100 -0.083770 -0.075640 -0.068540 -0.062260 -0.056690 -0.047240 -0.039500 -0.033080 -0.027680 -0.023140 -0.019290 -0.016030 -0.013270 -0.010920 -0.008925 -0.005808 -0.003580 -0.002728 -0.002017 -0.000930 -0.000190 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000]; phi3 = [-0.005911 -0.005864 -0.005817 -0.005730 -0.005578 -0.005533 -0.005451 -0.005316 -0.005274 -0.005087 -0.005002 -0.004917 -0.004873 -0.004834 -0.004758 -0.004682 -0.004517 -0.004374 -0.004266 -0.004228 -0.004182 -0.004120 -0.004079 -0.004058 -0.004050 -0.004058 -0.004075 -0.004136 -0.004228 -0.004348 -0.004383 -0.004476 -0.004617 -0.004762 -0.004903 -0.005047 -0.005185 -0.005316 -0.005572 -0.005799 -0.005905 -0.006006 -0.006195 -0.006289 -0.006371 -0.006532 -0.006678 -0.006745 -0.006812 -0.006943 -0.007055 -0.007162 -0.007263 -0.007314 -0.007365 -0.007447 -0.007529 -0.007712 -0.007868 -0.008003 -0.008043 -0.008115 -0.008213 -0.008304 -0.008379 -0.008438 -0.008497 -0.008548 -0.008626 -0.008695 -0.008747 -0.008791 -0.008818 -0.008835 -0.008862 -0.008871 -0.008880 -0.008880 -0.008888 -0.008888 -0.008888 -0.008888 -0.008888 -0.008888 -0.008888 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880 -0.008880]; phi4 = [0.100359 0.103066 0.104269 0.106074 0.109583 0.110415 0.112089 0.114716 0.115598 0.120311 0.122617 0.125023 0.126226 0.127429 0.129935 0.132442 0.138858 0.145375 0.151892 0.153998 0.158309 0.164525 0.170540 0.176355 0.181769 0.186983 0.191795 0.200518 0.207736 0.213752 0.215456 0.218464 0.222073 0.224580 0.226284 0.227086 0.227187 0.226685 0.224279 0.220469 0.218013 0.215557 0.209942 0.207035 0.204027 0.197811 0.191494 0.188386 0.185278 0.179163 0.173147 0.167332 0.161718 0.159011 0.156304 0.151090 0.146177 0.134547 0.124221 0.114997 0.112190 0.106776 0.099377 0.092739 0.086744 0.081330 0.076427 0.071956 0.064146 0.057579 0.051984 0.047192 0.043041 0.039412 0.036234 0.033426 0.030940 0.028714 0.024924 0.021826 0.020553 0.019270 0.017114 0.015300 0.013745 0.012402 0.011821 0.011239 0.010226 0.009334 0.008552 0.007850 0.007229 0.006677 0.006176 0.005133 0.004311 0.003659 0.003128 0.002697 0.002346 0.002045 0.001795 0.001584 0.001404]; tau1 = [0.3024 0.3034 0.3089 0.3151 0.3235 0.3257 0.3303 0.3339 0.3346 0.3384 0.3432 0.3468 0.3476 0.3493 0.3551 0.3634 0.3856 0.3913 0.4036 0.4092 0.4065 0.4101 0.4120 0.4131 0.4200 0.4243 0.4230 0.4126 0.4243 0.4088 0.4055 0.4012 0.3929 0.3778 0.3652 0.3526 0.3380 0.3234 0.3189 0.3087 0.3075 0.3018 0.2871 0.2876 0.2902 0.2836 0.2879 0.2931 0.2959 0.3009 0.3088 0.3175 0.3232 0.3215 0.3173 0.3063 0.3038 0.3095 0.3170 0.3205 0.3224 0.3203 0.3029 0.3086 0.2934 0.2999 0.3106 0.3165 0.3421 0.3510 0.3494 0.3469 0.3513 0.3591 0.3655 0.3735 0.3810 0.3818 0.3764 0.3889 0.3896 0.4008 0.4099 0.4045 0.4127 0.4195 0.4162 0.4158 0.4215 0.4287 0.4648 0.4520 0.4443 0.4466 0.4405 0.4379 0.4245 0.4414 0.4585 0.4888 0.5105 0.3588 0.3597 0.3528 0.2934]; sigma1 = [0.4720 0.4736 0.4747 0.4769 0.4815 0.4822 0.4835 0.4865 0.4879 0.4925 0.4945 0.4957 0.4962 0.4968 0.4985 0.4989 0.4998 0.5029 0.5027 0.5012 0.5036 0.5043 0.5068 0.5091 0.5092 0.5084 0.5104 0.5125 0.5148 0.5150 0.5141 0.5109 0.5096 0.5111 0.5117 0.5168 0.5182 0.5154 0.5110 0.5123 0.5129 0.5159 0.5256 0.5303 0.5343 0.5426 0.5429 0.5428 0.5408 0.5370 0.5373 0.5390 0.5395 0.5397 0.5398 0.5419 0.5461 0.5511 0.5614 0.5638 0.5652 0.5661 0.5641 0.5621 0.5605 0.5578 0.5540 0.5477 0.5442 0.5381 0.5383 0.5374 0.5379 0.5440 0.5447 0.5440 0.5413 0.5439 0.5460 0.5442 0.5431 0.5381 0.5321 0.5397 0.5480 0.5524 0.5532 0.5549 0.5547 0.5563 0.5521 0.5609 0.5709 0.5782 0.5837 0.5894 0.5997 0.6110 0.6038 0.6048 0.6156 0.6110 0.6264 0.6357 0.6445]; % interpolate between periods if neccesary if (length(find(period == T)) == 0) T_low = max(period(find(periodT))); [sa_low, sigma_low, tau_low] = youngs_atten(T_low, M, rrup, rjb, delta, w, ztor, lambda, Vs30); [sa_hi, sigma_hi, tau_hi] = youngs_atten(T_hi, M, rrup, rjb, delta, w, ztor, lambda, Vs30); x = [log(T_low) log(T_hi)]; Y_sa = [log(sa_low) log(sa_hi)]; Y_sigma = [sigma_low sigma_hi]; Y_tau = [tau_low tau_hi]; Sa = exp(interp1(x,Y_sa,log(T))); sigma = interp1(x,Y_sigma,log(T)); tau = interp1(x,Y_tau,log(T)); else i = find(period == T); % Identify the period frv = lambda >= 30 & lambda <= 150; % frv: 1 for lambda between 30 and 150, 0 otherwise fnm = lambda >= -120 & lambda <= -60; % fnm: 1 for lambda between -120 and -60, 0 otherwise term1 = c1(i); term2 = c1a(i) * frv; term3 = c1b(i) * fnm; term4 = c7(i) * (ztor - 4); term5 = c2(i) * (M - 6); term6 = ((c2(i) - c3(i))/cn(i)) * log (1 + exp(cn(i) * (cm(i) - M))); term7 = c4(i) * log(rrup + c5(i) * cosh(c6(i) * max(M - chm(i),0))); term8 = (c4a(i) - c4(i)) * log (sqrt(rrup^2 + (crb(i))^2)); term9 = (cy1(i) + cy2(i)/cosh(max(M-cy3(i),0)))*rrup; term10 = c9(i) * (cos(delta*pi()/180))^2 * tanh(rrup/2) * atan((w*cos(delta*pi()/180))/(2*(ztor+1))) * (2/pi()) * (1 - (rjb/(rrup + 0.001))); Sa1130 = exp(term1 + term2 + term3 + term4 + term5 + term6 + term7 + term8 + term9 + term10); term11 = phi1(i) * min(log(Vs30/1130),0); term12 = phi2(i) * (exp(phi3(i) * (min(Vs30,1130) - 360)) - exp(phi3(i) * (1130 - 360))) * log((Sa1130 + phi4(i))/phi4(i)); % Compute median Sa = exp(log(Sa1130) + term11 + term12); % Compute standard deviation sigma = sigma1(i); tau = tau1(i); end