subroutine int_Bz_LargeFld(val,eval) implicit none integer i,n parameter(n=1101) double precision temp(n), norm(n) double precision val,eval,slop,y_int c write (*,*) 'Give z Position (mm)' c read (*,*) val data temp/-50.0, !z postion in mm -KG & -49.0, & -48.0, & -47.0, & -46.0, & -45.0, & -44.0, & -43.0, & -42.0, & -41.0, & -40.0, & -39.0, & -38.0, & -37.0, & -36.0, & -35.0, & -34.0, & -33.0, & -32.0, & -31.0, & -30.0, & -29.0, & -28.0, & -27.0, & -26.0, & -25.0, & -24.0, & -23.0, & -22.0, & -21.0, & -20.0, & -19.0, & -18.0, & -17.0, & -16.0, & -15.0, & -14.0, & -13.0, & -12.0, & -11.0, & -10.0, & -9.0, & -8.0, & -7.0, & -6.0, & -5.0, & -4.0, & -3.0, & -2.0, & -1.0, & 0.0, & 1.0, & 2.0, & 3.0, & 4.0, & 5.0, & 6.0, & 7.0, & 8.0, & 9.0, & 10.0, & 11.0, & 12.0, & 13.0, & 14.0, & 15.0, & 16.0, & 17.0, & 18.0, & 19.0, & 20.0, & 21.0, & 22.0, & 23.0, & 24.0, & 25.0, & 26.0, & 27.0, & 28.0, & 29.0, & 30.0, & 31.0, & 32.0, & 33.0, & 34.0, & 35.0, & 36.0, & 37.0, & 38.0, & 39.0, & 40.0, & 41.0, & 42.0, & 43.0, & 44.0, & 45.0, & 46.0, & 47.0, & 48.0, & 49.0, & 50.0, & 51.0, & 52.0, & 53.0, & 54.0, & 55.0, & 56.0, & 57.0, & 58.0, & 59.0, & 60.0, & 61.0, & 62.0, & 63.0, & 64.0, & 65.0, & 66.0, & 67.0, & 68.0, & 69.0, & 70.0, & 71.0, & 72.0, & 73.0, & 74.0, & 75.0, & 76.0, & 77.0, & 78.0, & 79.0, & 80.0, & 81.0, & 82.0, & 83.0, & 84.0, & 85.0, & 86.0, & 87.0, & 88.0, & 89.0, & 90.0, & 91.0, & 92.0, & 93.0, & 94.0, & 95.0, & 96.0, & 97.0, & 98.0, & 99.0, & 100.0, & 101.0, & 102.0, & 103.0, & 104.0, & 105.0, & 106.0, & 107.0, & 108.0, & 109.0, & 110.0, & 111.0, & 112.0, & 113.0, & 114.0, & 115.0, & 116.0, & 117.0, & 118.0, & 119.0, & 120.0, & 121.0, & 122.0, & 123.0, & 124.0, & 125.0, & 126.0, & 127.0, & 128.0, & 129.0, & 130.0, & 131.0, & 132.0, & 133.0, & 134.0, & 135.0, & 136.0, & 137.0, & 138.0, & 139.0, & 140.0, & 141.0, & 142.0, & 143.0, & 144.0, & 145.0, & 146.0, & 147.0, & 148.0, & 149.0, & 150.0, & 151.0, & 152.0, & 153.0, & 154.0, & 155.0, & 156.0, & 157.0, & 158.0, & 159.0, & 160.0, & 161.0, & 162.0, & 163.0, & 164.0, & 165.0, & 166.0, & 167.0, & 168.0, & 169.0, & 170.0, & 171.0, & 172.0, & 173.0, & 174.0, & 175.0, & 176.0, & 177.0, & 178.0, & 179.0, & 180.0, & 181.0, & 182.0, & 183.0, & 184.0, & 185.0, & 186.0, & 187.0, & 188.0, & 189.0, & 190.0, & 191.0, & 192.0, & 193.0, & 194.0, & 195.0, & 196.0, & 197.0, & 198.0, & 199.0, & 200.0, & 201.0, & 202.0, & 203.0, & 204.0, & 205.0, & 206.0, & 207.0, & 208.0, & 209.0, & 210.0, & 211.0, & 212.0, & 213.0, & 214.0, & 215.0, & 216.0, & 217.0, & 218.0, & 219.0, & 220.0, & 221.0, & 222.0, & 223.0, & 224.0, & 225.0, & 226.0, & 227.0, & 228.0, & 229.0, & 230.0, & 231.0, & 232.0, & 233.0, & 234.0, & 235.0, & 236.0, & 237.0, & 238.0, & 239.0, & 240.0, & 241.0, & 242.0, & 243.0, & 244.0, & 245.0, & 246.0, & 247.0, & 248.0, & 249.0, & 250.0, & 251.0, & 252.0, & 253.0, & 254.0, & 255.0, & 256.0, & 257.0, & 258.0, & 259.0, & 260.0, & 261.0, & 262.0, & 263.0, & 264.0, & 265.0, & 266.0, & 267.0, & 268.0, & 269.0, & 270.0, & 271.0, & 272.0, & 273.0, & 274.0, & 275.0, & 276.0, & 277.0, & 278.0, & 279.0, & 280.0, & 281.0, & 282.0, & 283.0, & 284.0, & 285.0, & 286.0, & 287.0, & 288.0, & 289.0, & 290.0, & 291.0, & 292.0, & 293.0, & 294.0, & 295.0, & 296.0, & 297.0, & 298.0, & 299.0, & 300.0, & 301.0, & 302.0, & 303.0, & 304.0, & 305.0, & 306.0, & 307.0, & 308.0, & 309.0, & 310.0, & 311.0, & 312.0, & 313.0, & 314.0, & 315.0, & 316.0, & 317.0, & 318.0, & 319.0, & 320.0, & 321.0, & 322.0, & 323.0, & 324.0, & 325.0, & 326.0, & 327.0, & 328.0, & 329.0, & 330.0, & 331.0, & 332.0, & 333.0, & 334.0, & 335.0, & 336.0, & 337.0, & 338.0, & 339.0, & 340.0, & 341.0, & 342.0, & 343.0, & 344.0, & 345.0, & 346.0, & 347.0, & 348.0, & 349.0, & 350.0, & 351.0, & 352.0, & 353.0, & 354.0, & 355.0, & 356.0, & 357.0, & 358.0, & 359.0, & 360.0, & 361.0, & 362.0, & 363.0, & 364.0, & 365.0, & 366.0, & 367.0, & 368.0, & 369.0, & 370.0, & 371.0, & 372.0, & 373.0, & 374.0, & 375.0, & 376.0, & 377.0, & 378.0, & 379.0, & 380.0, & 381.0, & 382.0, & 383.0, & 384.0, & 385.0, & 386.0, & 387.0, & 388.0, & 389.0, & 390.0, & 391.0, & 392.0, & 393.0, & 394.0, & 395.0, & 396.0, & 397.0, & 398.0, & 399.0, & 400.0, & 401.0, & 402.0, & 403.0, & 404.0, & 405.0, & 406.0, & 407.0, & 408.0, & 409.0, & 410.0, & 411.0, & 412.0, & 413.0, & 414.0, & 415.0, & 416.0, & 417.0, & 418.0, & 419.0, & 420.0, & 421.0, & 422.0, & 423.0, & 424.0, & 425.0, & 426.0, & 427.0, & 428.0, & 429.0, & 430.0, & 431.0, & 432.0, & 433.0, & 434.0, & 435.0, & 436.0, & 437.0, & 438.0, & 439.0, & 440.0, & 441.0, & 442.0, & 443.0, & 444.0, & 445.0, & 446.0, & 447.0, & 448.0, & 449.0, & 450.0, & 451.0, & 452.0, & 453.0, & 454.0, & 455.0, & 456.0, & 457.0, & 458.0, & 459.0, & 460.0, & 461.0, & 462.0, & 463.0, & 464.0, & 465.0, & 466.0, & 467.0, & 468.0, & 469.0, & 470.0, & 471.0, & 472.0, & 473.0, & 474.0, & 475.0, & 476.0, & 477.0, & 478.0, & 479.0, & 480.0, & 481.0, & 482.0, & 483.0, & 484.0, & 485.0, & 486.0, & 487.0, & 488.0, & 489.0, & 490.0, & 491.0, & 492.0, & 493.0, & 494.0, & 495.0, & 496.0, & 497.0, & 498.0, & 499.0, & 500.0, & 501.0, & 502.0, & 503.0, & 504.0, & 505.0, & 506.0, & 507.0, & 508.0, & 509.0, & 510.0, & 511.0, & 512.0, & 513.0, & 514.0, & 515.0, & 516.0, & 517.0, & 518.0, & 519.0, & 520.0, & 521.0, & 522.0, & 523.0, & 524.0, & 525.0, & 526.0, & 527.0, & 528.0, & 529.0, & 530.0, & 531.0, & 532.0, & 533.0, & 534.0, & 535.0, & 536.0, & 537.0, & 538.0, & 539.0, & 540.0, & 541.0, & 542.0, & 543.0, & 544.0, & 545.0, & 546.0, & 547.0, & 548.0, & 549.0, & 550.0, & 551.0, & 552.0, & 553.0, & 554.0, & 555.0, & 556.0, & 557.0, & 558.0, & 559.0, & 560.0, & 561.0, & 562.0, & 563.0, & 564.0, & 565.0, & 566.0, & 567.0, & 568.0, & 569.0, & 570.0, & 571.0, & 572.0, & 573.0, & 574.0, & 575.0, & 576.0, & 577.0, & 578.0, & 579.0, & 580.0, & 581.0, & 582.0, & 583.0, & 584.0, & 585.0, & 586.0, & 587.0, & 588.0, & 589.0, & 590.0, & 591.0, & 592.0, & 593.0, & 594.0, & 595.0, & 596.0, & 597.0, & 598.0, & 599.0, & 600.0, & 601.0, & 602.0, & 603.0, & 604.0, & 605.0, & 606.0, & 607.0, & 608.0, & 609.0, & 610.0, & 611.0, & 612.0, & 613.0, & 614.0, & 615.0, & 616.0, & 617.0, & 618.0, & 619.0, & 620.0, & 621.0, & 622.0, & 623.0, & 624.0, & 625.0, & 626.0, & 627.0, & 628.0, & 629.0, & 630.0, & 631.0, & 632.0, & 633.0, & 634.0, & 635.0, & 636.0, & 637.0, & 638.0, & 639.0, & 640.0, & 641.0, & 642.0, & 643.0, & 644.0, & 645.0, & 646.0, & 647.0, & 648.0, & 649.0, & 650.0, & 651.0, & 652.0, & 653.0, & 654.0, & 655.0, & 656.0, & 657.0, & 658.0, & 659.0, & 660.0, & 661.0, & 662.0, & 663.0, & 664.0, & 665.0, & 666.0, & 667.0, & 668.0, & 669.0, & 670.0, & 671.0, & 672.0, & 673.0, & 674.0, & 675.0, & 676.0, & 677.0, & 678.0, & 679.0, & 680.0, & 681.0, & 682.0, & 683.0, & 684.0, & 685.0, & 686.0, & 687.0, & 688.0, & 689.0, & 690.0, & 691.0, & 692.0, & 693.0, & 694.0, & 695.0, & 696.0, & 697.0, & 698.0, & 699.0, & 700.0, & 701.0, & 702.0, & 703.0, & 704.0, & 705.0, & 706.0, & 707.0, & 708.0, & 709.0, & 710.0, & 711.0, & 712.0, & 713.0, & 714.0, & 715.0, & 716.0, & 717.0, & 718.0, & 719.0, & 720.0, & 721.0, & 722.0, & 723.0, & 724.0, & 725.0, & 726.0, & 727.0, & 728.0, & 729.0, & 730.0, & 731.0, & 732.0, & 733.0, & 734.0, & 735.0, & 736.0, & 737.0, & 738.0, & 739.0, & 740.0, & 741.0, & 742.0, & 743.0, & 744.0, & 745.0, & 746.0, & 747.0, & 748.0, & 749.0, & 750.0, & 751.0, & 752.0, & 753.0, & 754.0, & 755.0, & 756.0, & 757.0, & 758.0, & 759.0, & 760.0, & 761.0, & 762.0, & 763.0, & 764.0, & 765.0, & 766.0, & 767.0, & 768.0, & 769.0, & 770.0, & 771.0, & 772.0, & 773.0, & 774.0, & 775.0, & 776.0, & 777.0, & 778.0, & 779.0, & 780.0, & 781.0, & 782.0, & 783.0, & 784.0, & 785.0, & 786.0, & 787.0, & 788.0, & 789.0, & 790.0, & 791.0, & 792.0, & 793.0, & 794.0, & 795.0, & 796.0, & 797.0, & 798.0, & 799.0, & 800.0, & 801.0, & 802.0, & 803.0, & 804.0, & 805.0, & 806.0, & 807.0, & 808.0, & 809.0, & 810.0, & 811.0, & 812.0, & 813.0, & 814.0, & 815.0, & 816.0, & 817.0, & 818.0, & 819.0, & 820.0, & 821.0, & 822.0, & 823.0, & 824.0, & 825.0, & 826.0, & 827.0, & 828.0, & 829.0, & 830.0, & 831.0, & 832.0, & 833.0, & 834.0, & 835.0, & 836.0, & 837.0, & 838.0, & 839.0, & 840.0, & 841.0, & 842.0, & 843.0, & 844.0, & 845.0, & 846.0, & 847.0, & 848.0, & 849.0, & 850.0, & 851.0, & 852.0, & 853.0, & 854.0, & 855.0, & 856.0, & 857.0, & 858.0, & 859.0, & 860.0, & 861.0, & 862.0, & 863.0, & 864.0, & 865.0, & 866.0, & 867.0, & 868.0, & 869.0, & 870.0, & 871.0, & 872.0, & 873.0, & 874.0, & 875.0, & 876.0, & 877.0, & 878.0, & 879.0, & 880.0, & 881.0, & 882.0, & 883.0, & 884.0, & 885.0, & 886.0, & 887.0, & 888.0, & 889.0, & 890.0, & 891.0, & 892.0, & 893.0, & 894.0, & 895.0, & 896.0, & 897.0, & 898.0, & 899.0, & 900.0, & 901.0, & 902.0, & 903.0, & 904.0, & 905.0, & 906.0, & 907.0, & 908.0, & 909.0, & 910.0, & 911.0, & 912.0, & 913.0, & 914.0, & 915.0, & 916.0, & 917.0, & 918.0, & 919.0, & 920.0, & 921.0, & 922.0, & 923.0, & 924.0, & 925.0, & 926.0, & 927.0, & 928.0, & 929.0, & 930.0, & 931.0, & 932.0, & 933.0, & 934.0, & 935.0, & 936.0, & 937.0, & 938.0, & 939.0, & 940.0, & 941.0, & 942.0, & 943.0, & 944.0, & 945.0, & 946.0, & 947.0, & 948.0, & 949.0, & 950.0, & 951.0, & 952.0, & 953.0, & 954.0, & 955.0, & 956.0, & 957.0, & 958.0, & 959.0, & 960.0, & 961.0, & 962.0, & 963.0, & 964.0, & 965.0, & 966.0, & 967.0, & 968.0, & 969.0, & 970.0, & 971.0, & 972.0, & 973.0, & 974.0, & 975.0, & 976.0, & 977.0, & 978.0, & 979.0, & 980.0, & 981.0, & 982.0, & 983.0, & 984.0, & 985.0, & 986.0, & 987.0, & 988.0, & 989.0, & 990.0, & 991.0, & 992.0, & 993.0, & 994.0, & 995.0, & 996.0, & 997.0, & 998.0, & 999.0, &1000.0, &1001.0, &1002.0, &1003.0, &1004.0, &1005.0, &1006.0, &1007.0, &1008.0, &1009.0, &1010.0, &1011.0, &1012.0, &1013.0, &1014.0, &1015.0, &1016.0, &1017.0, &1018.0, &1019.0, &1020.0, &1021.0, &1022.0, &1023.0, &1024.0, &1025.0, &1026.0, &1027.0, &1028.0, &1029.0, &1030.0, &1031.0, &1032.0, &1033.0, &1034.0, &1035.0, &1036.0, &1037.0, &1038.0, &1039.0, &1040.0, &1041.0, &1042.0, &1043.0, &1044.0, &1045.0, &1046.0, &1047.0, &1048.0, &1049.0, &1050.0/ ! (mm) data norm/!integral Bz*dz large field systematic runs -KG & 0.0000, & 1.8140, & 3.6440, & 5.4910, & 7.3550, & 9.2360, & 11.1300, & 13.0490, & 14.9820, & 16.9320, & 18.8990, & 20.8840, & 22.8880, & 24.9110, & 26.9520, & 29.0120, & 31.0900, & 33.1880, & 35.3050, & 37.4410, & 39.5970, & 41.7730, & 43.9700, & 46.1860, & 48.4230, & 50.6810, & 52.9600, & 55.2610, & 57.5830, & 59.9260, & 62.2910, & 64.6790, & 67.0880, & 69.5200, & 71.9760, & 74.4540, & 76.9550, & 79.4800, & 82.0290, & 84.6010, & 87.1980, & 89.8190, & 92.4650, & 95.1360, & 97.8330, & 100.5540, & 103.3010, & 106.0740, & 108.8730, & 111.6980, & 114.5500, & 117.4300, & 120.3360, & 123.2700, & 126.2310, & 129.2210, & 132.2380, & 135.2830, & 138.3580, & 141.4610, & 144.5940, & 147.7560, & 150.9480, & 154.1700, & 157.4220, & 160.7050, & 164.0180, & 167.3620, & 170.7380, & 174.1440, & 177.5830, & 181.0530, & 184.5570, & 188.0920, & 191.6600, & 195.2620, & 198.8960, & 202.5640, & 206.2660, & 210.0020, & 213.7710, & 217.5750, & 221.4140, & 225.2880, & 229.1970, & 233.1410, & 237.1210, & 241.1370, & 245.1900, & 249.2780, & 253.4020, & 257.5630, & 261.7610, & 265.9950, & 270.2680, & 274.5780, & 278.9250, & 283.3110, & 287.7340, & 292.1960, & 296.6960, & 301.2350, & 305.8110, & 310.4280, & 315.0820, & 319.7760, & 324.5090, & 329.2810, & 334.0940, & 338.9450, & 343.8370, & 348.7680, & 353.7390, & 358.7500, & 363.8010, & 368.8920, & 374.0230, & 379.1950, & 384.4060, & 389.6590, & 394.9510, & 400.2840, & 405.6570, & 411.0710, & 416.5240, & 422.0190, & 427.5530, & 433.1280, & 438.7430, & 444.3980, & 450.0940, & 455.8300, & 461.6070, & 467.4230, & 473.2790, & 479.1750, & 485.1100, & 491.0860, & 497.1000, & 503.1540, & 509.2470, & 515.3800, & 521.5510, & 527.7620, & 534.0100, & 540.2970, & 546.6230, & 552.9860, & 559.3870, & 565.8270, & 572.3030, & 578.8170, & 585.3680, & 591.9560, & 598.5790, & 605.2400, & 611.9360, & 618.6690, & 625.4370, & 632.2390, & 639.0770, & 645.9490, & 652.8560, & 659.7970, & 666.7730, & 673.7810, & 680.8230, & 687.8980, & 695.0060, & 702.1460, & 709.3190, & 716.5230, & 723.7590, & 731.0260, & 738.3230, & 745.6520, & 753.0100, & 760.3990, & 767.8180, & 775.2670, & 782.7440, & 790.2510, & 797.7860, & 805.3490, & 812.9410, & 820.5590, & 828.2050, & 835.8780, & 843.5780, & 851.3050, & 859.0580, & 866.8380, & 874.6430, & 882.4730, & 890.3290, & 898.2100, & 906.1150, & 914.0450, & 921.9990, & 929.9780, & 937.9800, & 946.0060, & 954.0550, & 962.1270, & 970.2220, & 978.3400, & 986.4800, & 994.6420, & 1002.8300, & 1011.0300, & 1019.2600, & 1027.5100, & 1035.7800, & 1044.0700, & 1052.3800, & 1060.7200, & 1069.0700, & 1077.4500, & 1085.8400, & 1094.2500, & 1102.6900, & 1111.1400, & 1119.6200, & 1128.1100, & 1136.6200, & 1145.1500, & 1153.7000, & 1162.2700, & 1170.8500, & 1179.4600, & 1188.0800, & 1196.7200, & 1205.3800, & 1214.0600, & 1222.7600, & 1231.4700, & 1240.2000, & 1248.9500, & 1257.7100, & 1266.4900, & 1275.2900, & 1284.1000, & 1292.9300, & 1301.7800, & 1310.6400, & 1319.5200, & 1328.4100, & 1337.3200, & 1346.2500, & 1355.1900, & 1364.1400, & 1373.1100, & 1382.1000, & 1391.1000, & 1400.1100, & 1409.1400, & 1418.1800, & 1427.2400, & 1436.3100, & 1445.4000, & 1454.4900, & 1463.6000, & 1472.7300, & 1481.8700, & 1491.0100, & 1500.1800, & 1509.3500, & 1518.5400, & 1527.7400, & 1536.9500, & 1546.1700, & 1555.4000, & 1564.6500, & 1573.9100, & 1583.1700, & 1592.4500, & 1601.7400, & 1611.0400, & 1620.3500, & 1629.6700, & 1639.0000, & 1648.3400, & 1657.6900, & 1667.0500, & 1676.4200, & 1685.8000, & 1695.1900, & 1704.5900, & 1713.9900, & 1723.4100, & 1732.8300, & 1742.2600, & 1751.7000, & 1761.1500, & 1770.6100, & 1780.0800, & 1789.5500, & 1799.0300, & 1808.5200, & 1818.0200, & 1827.5200, & 1837.0300, & 1846.5500, & 1856.0800, & 1865.6100, & 1875.1500, & 1884.7000, & 1894.2600, & 1903.8200, & 1913.3900, & 1922.9700, & 1932.5500, & 1942.1400, & 1951.7400, & 1961.3500, & 1970.9600, & 1980.5700, & 1990.2000, & 1999.8300, & 2009.4700, & 2019.1100, & 2028.7600, & 2038.4200, & 2048.0800, & 2057.7600, & 2067.4300, & 2077.1200, & 2086.8000, & 2096.5000, & 2106.2000, & 2115.9100, & 2125.6300, & 2135.3500, & 2145.0700, & 2154.8100, & 2164.5500, & 2174.2900, & 2184.0400, & 2193.8000, & 2203.5700, & 2213.3400, & 2223.1100, & 2232.8900, & 2242.6800, & 2252.4700, & 2262.2700, & 2272.0800, & 2281.8900, & 2291.7100, & 2301.5300, & 2311.3600, & 2321.1900, & 2331.0300, & 2340.8700, & 2350.7200, & 2360.5800, & 2370.4400, & 2380.3000, & 2390.1800, & 2400.0500, & 2409.9300, & 2419.8200, & 2429.7100, & 2439.6000, & 2449.5000, & 2459.4000, & 2469.3100, & 2479.2300, & 2489.1400, & 2499.0600, & 2508.9900, & 2518.9200, & 2528.8500, & 2538.7900, & 2548.7300, & 2558.6800, & 2568.6300, & 2578.5800, & 2588.5300, & 2598.4900, & 2608.4500, & 2618.4200, & 2628.3900, & 2638.3600, & 2648.3300, & 2658.3100, & 2668.2900, & 2678.2800, & 2688.2600, & 2698.2500, & 2708.2400, & 2718.2400, & 2728.2300, & 2738.2300, & 2748.2300, & 2758.2400, & 2768.2400, & 2778.2500, & 2788.2600, & 2798.2700, & 2808.2900, & 2818.3100, & 2828.3200, & 2838.3400, & 2848.3700, & 2858.3900, & 2868.4200, & 2878.4500, & 2888.4800, & 2898.5100, & 2908.5400, & 2918.5800, & 2928.6200, & 2938.6600, & 2948.7000, & 2958.7500, & 2968.7900, & 2978.8400, & 2988.8900, & 2998.9400, & 3009.0000, & 3019.0500, & 3029.1100, & 3039.1700, & 3049.2300, & 3059.2900, & 3069.3600, & 3079.4300, & 3089.5000, & 3099.5700, & 3109.6500, & 3119.7200, & 3129.8000, & 3139.8900, & 3149.9700, & 3160.0500, & 3170.1400, & 3180.2300, & 3190.3200, & 3200.4200, & 3210.5200, & 3220.6100, & 3230.7200, & 3240.8200, & 3250.9200, & 3261.0300, & 3271.1400, & 3281.2600, & 3291.3700, & 3301.4900, & 3311.6100, & 3321.7300, & 3331.8500, & 3341.9800, & 3352.1100, & 3362.2400, & 3372.3700, & 3382.5000, & 3392.6400, & 3402.7800, & 3412.9200, & 3423.0600, & 3433.2000, & 3443.3400, & 3453.4900, & 3463.6400, & 3473.7900, & 3483.9400, & 3494.0900, & 3504.2500, & 3514.4000, & 3524.5600, & 3534.7200, & 3544.8800, & 3555.0400, & 3565.2000, & 3575.3600, & 3585.5200, & 3595.6800, & 3605.8500, & 3616.0100, & 3626.1800, & 3636.3500, & 3646.5100, & 3656.6800, & 3666.8500, & 3677.0200, & 3687.1900, & 3697.3500, & 3707.5200, & 3717.6900, & 3727.8600, & 3738.0300, & 3748.2000, & 3758.3700, & 3768.5400, & 3778.7100, & 3788.8700, & 3799.0400, & 3809.2100, & 3819.3800, & 3829.5500, & 3839.7200, & 3849.8800, & 3860.0500, & 3870.2200, & 3880.3800, & 3890.5500, & 3900.7200, & 3910.8800, & 3921.0500, & 3931.2100, & 3941.3800, & 3951.5400, & 3961.7100, & 3971.8700, & 3982.0300, & 3992.2000, & 4002.3600, & 4012.5200, & 4022.6900, & 4032.8500, & 4043.0100, & 4053.1700, & 4063.3400, & 4073.5000, & 4083.6600, & 4093.8200, & 4103.9800, & 4114.1500, & 4124.3100, & 4134.4700, & 4144.6300, & 4154.8000, & 4164.9600, & 4175.1200, & 4185.2900, & 4195.4500, & 4205.6100, & 4215.7800, & 4225.9400, & 4236.1100, & 4246.2700, & 4256.4400, & 4266.6000, & 4276.7700, & 4286.9300, & 4297.1000, & 4307.2700, & 4317.4300, & 4327.6000, & 4337.7700, & 4347.9400, & 4358.1000, & 4368.2700, & 4378.4400, & 4388.6100, & 4398.7800, & 4408.9500, & 4419.1200, & 4429.2900, & 4439.4500, & 4449.6200, & 4459.7900, & 4469.9600, & 4480.1300, & 4490.3000, & 4500.4700, & 4510.6300, & 4520.8000, & 4530.9700, & 4541.1400, & 4551.3000, & 4561.4700, & 4571.6300, & 4581.7900, & 4591.9600, & 4602.1200, & 4612.2800, & 4622.4400, & 4632.6000, & 4642.7600, & 4652.9200, & 4663.0700, & 4673.2300, & 4683.3800, & 4693.5300, & 4703.6800, & 4713.8300, & 4723.9800, & 4734.1200, & 4744.2700, & 4754.4100, & 4764.5500, & 4774.6900, & 4784.8200, & 4794.9600, & 4805.0900, & 4815.2200, & 4825.3500, & 4835.4800, & 4845.6000, & 4855.7200, & 4865.8400, & 4875.9600, & 4886.0700, & 4896.1900, & 4906.3000, & 4916.4100, & 4926.5100, & 4936.6200, & 4946.7200, & 4956.8200, & 4966.9200, & 4977.0100, & 4987.1100, & 4997.2000, & 5007.2900, & 5017.3700, & 5027.4600, & 5037.5400, & 5047.6200, & 5057.7000, & 5067.7800, & 5077.8500, & 5087.9200, & 5097.9900, & 5108.0600, & 5118.1200, & 5128.1900, & 5138.2500, & 5148.3000, & 5158.3600, & 5168.4200, & 5178.4700, & 5188.5200, & 5198.5700, & 5208.6100, & 5218.6600, & 5228.7000, & 5238.7400, & 5248.7800, & 5258.8200, & 5268.8500, & 5278.8900, & 5288.9200, & 5298.9500, & 5308.9800, & 5319.0000, & 5329.0300, & 5339.0500, & 5349.0700, & 5359.0800, & 5369.1000, & 5379.1100, & 5389.1200, & 5399.1300, & 5409.1400, & 5419.1400, & 5429.1400, & 5439.1400, & 5449.1400, & 5459.1300, & 5469.1300, & 5479.1200, & 5489.1000, & 5499.0900, & 5509.0700, & 5519.0500, & 5529.0300, & 5539.0000, & 5548.9700, & 5558.9300, & 5568.9000, & 5578.8600, & 5588.8200, & 5598.7700, & 5608.7200, & 5618.6700, & 5628.6100, & 5638.5500, & 5648.4800, & 5658.4200, & 5668.3400, & 5678.2700, & 5688.1900, & 5698.1000, & 5708.0100, & 5717.9200, & 5727.8200, & 5737.7100, & 5747.6100, & 5757.4900, & 5767.3800, & 5777.2500, & 5787.1300, & 5796.9900, & 5806.8600, & 5816.7100, & 5826.5700, & 5836.4100, & 5846.2500, & 5856.0900, & 5865.9200, & 5875.7400, & 5885.5600, & 5895.3800, & 5905.1800, & 5914.9800, & 5924.7800, & 5934.5700, & 5944.3500, & 5954.1300, & 5963.9100, & 5973.6700, & 5983.4300, & 5993.1900, & 6002.9300, & 6012.6800, & 6022.4100, & 6032.1400, & 6041.8700, & 6051.5800, & 6061.2900, & 6071.0000, & 6080.7000, & 6090.3900, & 6100.0700, & 6109.7500, & 6119.4300, & 6129.0900, & 6138.7500, & 6148.4000, & 6158.0500, & 6167.6900, & 6177.3200, & 6186.9500, & 6196.5700, & 6206.1900, & 6215.7900, & 6225.3900, & 6234.9900, & 6244.5700, & 6254.1500, & 6263.7300, & 6273.2900, & 6282.8500, & 6292.4000, & 6301.9500, & 6311.4800, & 6321.0100, & 6330.5400, & 6340.0500, & 6349.5600, & 6359.0600, & 6368.5500, & 6378.0400, & 6387.5100, & 6396.9800, & 6406.4400, & 6415.8900, & 6425.3400, & 6434.7700, & 6444.2000, & 6453.6200, & 6463.0300, & 6472.4300, & 6481.8200, & 6491.2000, & 6500.5800, & 6509.9400, & 6519.2900, & 6528.6400, & 6537.9700, & 6547.3000, & 6556.6100, & 6565.9200, & 6575.2100, & 6584.4900, & 6593.7600, & 6603.0200, & 6612.2700, & 6621.5100, & 6630.7400, & 6639.9600, & 6649.1600, & 6658.3500, & 6667.5300, & 6676.7000, & 6685.8500, & 6694.9900, & 6704.1200, & 6713.2400, & 6722.3400, & 6731.4300, & 6740.5100, & 6749.5700, & 6758.6200, & 6767.6500, & 6776.6700, & 6785.6800, & 6794.6700, & 6803.6500, & 6812.6100, & 6821.5500, & 6830.4900, & 6839.4000, & 6848.3000, & 6857.1900, & 6866.0500, & 6874.9100, & 6883.7500, & 6892.5700, & 6901.3700, & 6910.1600, & 6918.9300, & 6927.6800, & 6936.4200, & 6945.1400, & 6953.8400, & 6962.5200, & 6971.1900, & 6979.8400, & 6988.4700, & 6997.0800, & 7005.6700, & 7014.2500, & 7022.8100, & 7031.3500, & 7039.8700, & 7048.3700, & 7056.8500, & 7065.3100, & 7073.7500, & 7082.1700, & 7090.5700, & 7098.9600, & 7107.3200, & 7115.6600, & 7123.9800, & 7132.2800, & 7140.5600, & 7148.8200, & 7157.0500, & 7165.2700, & 7173.4600, & 7181.6300, & 7189.7800, & 7197.9100, & 7206.0100, & 7214.0900, & 7222.1500, & 7230.1800, & 7238.2000, & 7246.1800, & 7254.1500, & 7262.0900, & 7270.0000, & 7277.8900, & 7285.7600, & 7293.6000, & 7301.4200, & 7309.2100, & 7316.9700, & 7324.7100, & 7332.4200, & 7340.1000, & 7347.7600, & 7355.3900, & 7362.9900, & 7370.5600, & 7378.1100, & 7385.6300, & 7393.1200, & 7400.5800, & 7408.0100, & 7415.4100, & 7422.7800, & 7430.1200, & 7437.4300, & 7444.7100, & 7451.9600, & 7459.1800, & 7466.3600, & 7473.5100, & 7480.6300, & 7487.7200, & 7494.7800, & 7501.8000, & 7508.7900, & 7515.7400, & 7522.6600, & 7529.5500, & 7536.4000, & 7543.2200, & 7550.0000, & 7556.7500, & 7563.4600, & 7570.1300, & 7576.7700, & 7583.3700, & 7589.9400, & 7596.4700, & 7602.9600, & 7609.4100, & 7615.8300, & 7622.2100, & 7628.5500, & 7634.8500, & 7641.1100, & 7647.3400, & 7653.5300, & 7659.6800, & 7665.7800, & 7671.8500, & 7677.8900, & 7683.8800, & 7689.8300, & 7695.7400, & 7701.6100, & 7707.4500, & 7713.2400, & 7718.9900, & 7724.7000, & 7730.3700, & 7736.0100, & 7741.6000, & 7747.1500, & 7752.6600, & 7758.1300, & 7763.5600, & 7768.9500, & 7774.3000, & 7779.6100, & 7784.8700, & 7790.1000, & 7795.2900, & 7800.4400, & 7805.5400, & 7810.6100, & 7815.6400, & 7820.6200, & 7825.5700, & 7830.4800, & 7835.3500, & 7840.1700, & 7844.9600, & 7849.7100, & 7854.4200, & 7859.0900, & 7863.7200, & 7868.3200, & 7872.8700, & 7877.3900, & 7881.8600, & 7886.3000, & 7890.7000, & 7895.0700, & 7899.3900, & 7903.6800, & 7907.9300, & 7912.1400, & 7916.3200, & 7920.4600, & 7924.5600, & 7928.6200, & 7932.6500, & 7936.6500, & 7940.6100, & 7944.5300, & 7948.4200, & 7952.2700, & 7956.0900, & 7959.8700, & 7963.6200, & 7967.3400, & 7971.0200, & 7974.6700, & 7978.2800, & 7981.8600, & 7985.4100, & 7988.9300, & 7992.4100, & 7995.8600, & 7999.2800, & 8002.6700, & 8006.0300, & 8009.3500, & 8012.6500, & 8015.9100, & 8019.1500, & 8022.3500, & 8025.5200, & 8028.6700, & 8031.7800, & 8034.8700, & 8037.9300, & 8040.9600, & 8043.9600, & 8046.9300, & 8049.8700, & 8052.7900, & 8055.6800, & 8058.5400, & 8061.3800, & 8064.1900, & 8066.9800, & 8069.7300, & 8072.4600, & 8075.1700, & 8077.8500, & 8080.5100, & 8083.1400, & 8085.7400, & 8088.3300, & 8090.8900, & 8093.4200, & 8095.9300, & 8098.4200, & 8100.8800, & 8103.3200, & 8105.7400, & 8108.1400, & 8110.5100, & 8112.8600, & 8115.1900, & 8117.5000, & 8119.7900, & 8122.0600, & 8124.3000, & 8126.5300, & 8128.7300, & 8130.9200, & 8133.0800, & 8135.2300, & 8137.3500, & 8139.4600, & 8141.5400, & 8143.6100, & 8145.6600, & 8147.6900, & 8149.7000, & 8151.6900, & 8153.6700, & 8155.6200, & 8157.5600, & 8159.4800, & 8161.3900, & 8163.2800, & 8165.1500, & 8167.0000, & 8168.8400, & 8170.6600, & 8172.4600, & 8174.2500, & 8176.0300, & 8177.7800, & 8179.5200, & 8181.2500, & 8182.9600, & 8184.6600, & 8186.3400, & 8188.0100, & 8189.6600, & 8191.2900, & 8192.9200, & 8194.5300, & 8196.1200, & 8197.7000, & 8199.2700, & 8200.8200, & 8202.3600, & 8203.8900, & 8205.4000, & 8206.9100, & 8208.3900, & 8209.8700, & 8211.3300, & 8212.7800, & 8214.2200, & 8215.6500, & 8217.0600, & 8218.4600, & 8219.8500, & 8221.2300, & 8222.6000, & 8223.9600, & 8225.3000, & 8226.6300, & 8227.9600, & 8229.2700, & 8230.5700, & 8231.8600, & 8233.1400, & 8234.4100, & 8235.6600, & 8236.9100, & 8238.1500, & 8239.3800, & 8240.6000, & 8241.8000/ * do i=1,n-1 if(val.lt.temp(1).or.val.gt.temp(n)) then write(*,*) val,temp(1),temp(n) write(*,*) 'z value is out of range ... intBz code' stop elseif(val.eq.temp(i)) then eval = norm(i) elseif(val.eq.temp(i+1)) then eval = norm(i+1) elseif(val.gt.temp(i).and.val.lt.temp(i+1)) then slop = (norm(i+1)-norm(i))/(temp(i+1)-temp(i)) y_int = norm(i) - slop*temp(i) eval = slop*val+y_int goto 100 !-KG else continue endif enddo c 100 write (*,*) 'Bz_INT = ', eval 100 return end