library(plotrix) deg <- c(-180.52,-179.93,-179.34,-178.75,-178.16,-177.58,-176.98,-176.39,-175.81,-175.22,-174.63,-174.04,-173.45,-172.87,-172.28, -171.69,-171.10,-170.50,-169.91,-169.32,-168.74,-168.15,-167.56,-166.98,-166.39,-165.80,-165.21,-164.62,-164.02,-163.44, -162.86,-162.27,-161.67,-161.08,-160.48,-159.90,-159.31,-158.72,-158.14,-157.55,-156.96,-156.38,-155.78,-155.19,-154.61, -154.02,-153.43,-152.83,-152.24,-151.64,-151.06,-150.47,-149.89,-149.30,-148.71,-148.11,-147.53,-146.94,-146.35,-145.76, -145.17,-144.57,-143.99,-143.40,-142.81,-142.23,-141.64,-141.05,-140.46,-139.87,-139.28,-138.70,-138.11,-137.52,-136.93, -136.34,-135.75,-135.17,-134.57,-133.99,-133.39,-132.80,-132.22,-131.63,-131.04,-130.44,-129.85,-129.25,-128.67,-128.08, -127.50,-126.91,-126.32,-125.73,-125.15,-124.55,-123.97,-123.38,-122.80,-122.21,-121.61,-121.03,-120.44,-119.85,-119.26, -118.67,-118.08,-117.49,-116.90,-116.31,-115.72,-115.13,-114.55,-113.96,-113.38,-112.78,-112.20,-111.61,-111.02,-110.44, -109.85,-109.27,-108.67,-108.07,-107.49,-106.90,-106.31,-105.72,-105.13,-104.54,-103.95,-103.36,-102.77,-102.18,-101.60, -101.01,-100.41,-99.83,-99.24,-98.65,-98.06,-97.47,-96.89,-96.29,-95.71,-95.12,-94.53,-93.95,-93.36,-92.77, -92.19,-91.60,-91.01,-90.42,-89.84,-89.25,-88.67,-88.08,-87.49,-86.91,-86.32,-85.72,-85.14,-84.55,-83.96, -83.37,-82.78,-82.19,-81.60,-81.01,-80.42,-79.84,-79.25,-78.67,-78.08,-77.48,-76.89,-76.31,-75.72,-75.13, -74.55,-73.96,-73.38,-72.79,-72.20,-71.61,-71.03,-70.43,-69.85,-69.26,-68.67,-68.08,-67.49,-66.90,-66.32, -65.73,-65.14,-64.55,-63.96,-63.37,-62.79,-62.20,-61.61,-61.01,-60.43,-59.84,-59.24,-58.65,-58.06,-57.47, -56.89,-56.30,-55.72,-55.13,-54.53,-53.95,-53.35,-52.77,-52.18,-51.59,-51.01,-50.42,-49.82,-49.24,-48.65, -48.05,-47.47,-46.88,-46.29,-45.70,-45.11,-44.53,-43.94,-43.35,-42.77,-42.18,-41.59,-41.00,-40.42,-39.81, -39.22,-38.63,-38.05,-37.46,-36.87,-36.28,-35.69,-35.11,-34.52,-33.92,-33.34,-32.74,-32.15,-31.56,-30.97, -30.37,-29.79,-29.21,-28.62,-28.02,-27.43,-26.85,-26.26,-25.68,-25.09,-24.50,-23.92,-23.33,-22.74,-22.15, -21.56,-20.97,-20.38,-19.79,-19.20,-18.60,-18.01,-17.42,-16.83,-16.23,-15.64,-15.05,-14.45,-13.87,-13.28, -12.69,-12.09,-11.50,-10.91,-10.31,-9.73,-9.13,-8.55,-7.95,-7.37,-6.78,-6.19,-5.59,-5.00,-4.41, -3.82,-3.23,-2.65,-2.06,-1.47,-0.88,-0.29,0.29,0.88,1.47,2.06,2.64,3.23,3.82,4.42, 5.01,5.61,6.19,6.78,7.37,7.97,8.55,9.14,9.73,10.31,10.91,11.49,12.08,12.67,13.26, 13.85,14.43,15.02,15.61,16.20,16.79,17.38,17.97,18.55,19.15,19.73,20.32,20.91,21.49,22.08, 22.66,23.25,23.84,24.42,25.01,25.61,26.21,26.79,27.38,27.96,28.55,29.14,29.73,30.32,30.91, 31.49,32.08,32.67,33.26,33.85,34.44,35.04,35.62,36.21,36.80,37.39,37.99,38.58,39.17,39.75, 40.34,40.92,41.51,42.10,42.69,43.28,43.88,44.47,45.05,45.65,46.24,46.82,47.41,48.00,48.59, 49.17,49.77,50.36,50.94,51.53,52.12,52.71,53.29,53.88,54.47,55.05,55.64,56.23,56.82,57.41, 58.00,58.60,59.19,59.77,60.36,60.95,61.54,62.13,62.72,63.30,63.89,64.48,65.06,65.66,66.26, 66.85,67.44,68.02,68.61,69.20,69.80,70.39,70.97,71.56,72.15,72.74,73.33,73.91,74.50,75.10, 75.68,76.27,76.86,77.45,78.04,78.63,79.22,79.81,80.39,80.98,81.57,82.15,82.74,83.33,83.92, 84.51,85.10,85.69,86.28,86.86,87.46,88.05,88.64,89.23,89.83,90.42,91.00,91.59,92.18,92.78, 93.36,93.95,94.55,95.14,95.73,96.31,96.91,97.49,98.08,98.67,99.25,99.84,100.43,101.02,101.61, 102.20,102.79,103.38,103.97,104.55,105.14,105.73,106.32,106.91,107.50,108.08,108.68,109.27,109.86,110.44, 111.03,111.61,112.20,112.79,113.38,113.97,114.56,115.14,115.74,116.32,116.91,117.50,118.09,118.68,119.27, 119.85,120.44,121.03,121.61,122.20,122.79,123.38,123.97,124.55,125.14,125.74,126.32,126.91,127.50,128.10, 128.69,129.28,129.86,130.45,131.03,131.62,132.22,132.80,133.39,133.99,134.57,135.16,135.75,136.33,136.93, 137.53,138.12,138.70,139.29,139.88,140.47,141.07,141.65,142.24,142.83,143.41,144.00,144.58,145.18,145.77, 146.36,146.95,147.54,148.12,148.72,149.30,149.89,150.48,151.07,151.66,152.24,152.83,153.42,154.01,154.60, 155.20,155.78,156.38,156.96,157.56,158.15,158.73,159.32,159.91,160.50,161.09,161.67,162.26,162.86,163.45, 164.04,164.63,165.22,165.80,166.39,166.99,167.57,168.17,168.76,169.34,169.93,170.52,171.10,171.69,172.28, 172.86,173.46,174.04,174.63,175.22,175.80,176.39,176.99,177.57,178.16,178.75,179.34) mag <- c(-37.20,-37.30,-38.70,-43.40,-48.90,-50.40,-47.10,-43.80,-41.10,-38.20,-36.40,-35.50,-35.10,-35.20,-35.70, -36.80,-38.10,-39.00,-40.20,-40.60,-40.60,-40.50,-40.40,-40.40,-40.40,-40.40,-40.70,-41.10,-41.40,-42.00, -41.70,-40.50,-38.80,-37.20,-36.20,-35.60,-35.40,-35.80,-37.10,-38.50,-40.60,-43.80,-46.80,-50.30,-52.80, -52.70,-52.10,-50.60,-48.80,-48.10,-47.50,-47.60,-48.80,-51.00,-51.00,-49.30,-46.10,-44.80,-43.10,-42.70, -42.20,-42.70,-43.70,-45.40,-47.00,-50.50,-54.20,-52.80,-48.80,-46.00,-44.40,-43.60,-43.60,-45.00,-47.10, -49.30,-53.50,-63.50,-72.40,-61.00,-59.70,-55.70,-53.60,-49.80,-46.60,-44.60,-43.30,-42.60,-42.90,-43.80, -45.60,-49.20,-52.80,-53.90,-50.80,-48.20,-47.20,-46.40,-46.40,-46.50,-46.50,-46.80,-47.20,-47.80,-47.90, -45.30,-43.40,-42.50,-41.30,-40.50,-39.80,-39.90,-40.40,-41.30,-42.80,-43.90,-44.60,-44.10,-43.00,-41.90, -41.00,-40.50,-41.10,-41.90,-42.50,-42.90,-42.90,-42.60,-42.30,-42.50,-42.40,-42.00,-41.20,-40.00,-39.10, -38.80,-39.10,-40.00,-40.90,-42.40,-42.40,-41.60,-40.60,-40.00,-39.50,-39.40,-39.30,-38.90,-38.90,-38.50, -38.10,-37.70,-37.10,-36.80,-37.00,-37.60,-38.60,-38.90,-38.60,-37.80,-37.10,-36.80,-36.60,-36.60,-36.30, -35.70,-34.90,-34.10,-33.80,-33.70,-34.00,-34.70,-35.40,-35.90,-36.30,-36.40,-36.60,-37.10,-37.90,-38.40, -38.50,-38.00,-37.70,-37.60,-37.70,-37.70,-37.70,-37.50,-37.00,-36.90,-36.50,-36.30,-36.30,-36.30,-36.30, -36.40,-36.90,-37.00,-37.60,-38.20,-39.00,-39.50,-39.70,-39.70,-39.20,-38.30,-37.00,-35.90,-35.10,-34.60, -34.30,-34.00,-33.90,-33.80,-33.70,-33.50,-33.00,-32.20,-31.30,-30.50,-29.60,-28.90,-28.10,-27.50,-27.10, -26.90,-26.90,-27.00,-27.20,-27.40,-27.90,-28.30,-28.70,-29.10,-29.10,-29.10,-28.90,-28.70,-28.60,-28.70, -28.90,-29.00,-29.00,-28.60,-28.20,-27.30,-26.50,-25.60,-24.80,-24.20,-23.70,-23.40,-23.20,-23.20,-23.30, -23.60,-23.90,-24.50,-25.20,-26.20,-27.50,-28.80,-29.80,-29.90,-29.00,-27.70,-26.50,-25.40,-24.80,-24.80, -25.20,-26.20,-27.70,-29.50,-30.20,-28.50,-26.30,-24.70,-23.80,-23.40,-23.60,-24.30,-25.50,-26.20,-25.60, -24.10,-22.60,-21.50,-20.60,-20.30,-20.20,-20.60,-21.60,-23.50,-27.20,-28.70,-22.00,-16.20,-12.20,-8.80, -6.30,-4.40,-2.90,-1.80,-1.00,-0.30,0.00,0.00,-0.30,-0.80,-1.70,-2.90,-4.30,-6.40,-8.90, -12.00,-16.50,-21.20,-23.40,-21.70,-20.30,-19.70,-19.70,-20.00,-20.30,-20.60,-21.00,-21.70,-22.90,-24.30, -25.90,-26.60,-26.00,-24.80,-24.10,-23.80,-24.10,-25.20,-27.00,-29.70,-31.60,-30.30,-27.60,-25.70,-24.40, -23.80,-23.70,-24.20,-25.10,-26.40,-28.10,-29.70,-30.60,-30.10,-28.80,-27.30,-25.90,-24.70,-23.70,-22.90, -22.40,-22.10,-22.00,-22.30,-22.60,-23.40,-24.40,-25.80,-27.50,-29.30,-30.30,-30.10,-29.30,-28.20,-27.30, -26.90,-27.10,-27.40,-28.00,-28.70,-29.70,-30.70,-31.60,-31.60,-31.00,-30.00,-29.10,-28.50,-28.10,-27.90, -27.80,-27.50,-27.30,-27.10,-27.10,-27.30,-27.80,-28.60,-29.40,-30.30,-31.20,-31.60,-32.10,-32.20,-32.40, -33.00,-34.00,-35.40,-36.90,-38.70,-39.90,-40.60,-40.30,-39.50,-38.30,-37.40,-36.50,-35.80,-35.40,-35.40, -35.50,-36.20,-37.00,-38.30,-39.80,-41.40,-41.80,-40.80,-39.70,-38.90,-38.70,-38.40,-37.90,-37.40,-37.10, -37.00,-37.10,-37.00,-36.40,-36.00,-36.10,-36.60,-37.40,-38.20,-38.50,-38.30,-38.30,-38.10,-38.10,-37.60, -36.80,-35.90,-35.00,-34.40,-34.20,-34.20,-34.50,-35.00,-35.80,-37.20,-39.10,-41.80,-44.30,-44.80,-42.70, -40.90,-39.70,-38.70,-38.20,-37.80,-37.70,-38.00,-38.60,-39.50,-40.70,-41.70,-42.30,-42.40,-42.20,-41.90, -42.00,-43.30,-45.50,-46.90,-46.50,-43.30,-40.50,-39.20,-38.40,-38.00,-38.10,-38.40,-39.10,-40.80,-44.10, -49.60,-45.90,-41.60,-39.20,-38.20,-38.40,-39.50,-41.40,-43.30,-45.10,-47.10,-49.90,-53.00,-48.50,-43.60, -41.20,-40.20,-41.00,-43.10,-49.40,-53.20,-45.90,-42.00,-40.10,-39.30,-39.60,-41.20,-44.50,-49.00,-47.20, -43.60,-42.60,-43.40,-45.80,-49.80,-52.10,-52.30,-54.40,-51.30,-46.10,-43.20,-41.70,-40.80,-41.10,-42.10, -44.20,-47.90,-51.30,-53.80,-54.10,-57.80,-50.30,-45.00,-42.20,-41.40,-41.20,-41.50,-41.50,-41.00,-40.10, -38.80,-37.40,-36.60,-36.60,-37.60,-38.90,-40.70,-42.00,-41.70,-40.60,-39.00,-38.00,-37.30,-37.20,-37.40, -37.90,-37.40,-36.20,-35.30,-35.50,-37.30,-41.40,-44.40,-41.80,-38.50,-37.80,-38.90,-41.50,-44.50,-47.00, -46.90,-43.40,-39.50,-37.10,-36.10,-36.30,-37.60,-40.50,-46.00,-49.40,-46.10,-42.80,-41.10,-40.40,-40.90, -42.00,-46.60,-58.90,-44.60,-38.90,-36.20,-35.40,-35.90,-37.00,-37.60,-36.70,-35.70) polar.plot(lengths=mag, polar.pos=deg, start=90, clockwise=TRUE, rp.type="p", show.centroid=TRUE)