(*************************************************************************
 * THIS SCRIPT PLOTTED 3D ARCTIC PH DATA ON TOP OF THE BATHYMETRIC SHELL * 
 * OF THE OCEAN FLOOR. IT PLOTS THE FIGURE AND EXPORTS A ROTATING        *
 * ANIMATION (IN AVI FORMAT) OF THE FIGURE. THE SAME SCRIPT WAS USED TO  *
 * PLOT PCO2 DATA ON A 3D MAP.                                           *
 *************************************************************************
 * SCRIPT WRITTEN BY MATT GOVE FOR USGS/JACOBS TECHNOLOGY IN 2010.       *
 *************************************************************************)

data = Import["Desktop/arctic_bathymetry_model/new_csv/all_bathy_mac.csv"];

spc = Import["Desktop/arctic_bathymetry_model/new_csv/3d_ph_in_situ.csv"];

s1pc = Import["Desktop/arctic_bathymetry_model/Real CTD Data/pH/s1ph.csv"];
s4pc = Import["Desktop/arctic_bathymetry_model/Real CTD Data/pH/s4ph.csv"];
s5pc = Import["Desktop/arctic_bathymetry_model/Real CTD Data/pH/s5ph.csv"];
s6pc = Import["Desktop/arctic_bathymetry_model/Real CTD Data/pH/s6ph.csv"];
s7pc = Import["Desktop/arctic_bathymetry_model/Real CTD Data/pH/s7ph.csv"];

p1 = ListPlot3D[data, 
  ColorFunction -> "BlueGreenYellow", 
  PlotRange -> {{-160, -120}, {70, 83}, {-5.059, 1.5}}];

pspc = ListPointPlot3D[spc, 
  PlotRange -> {{-160, -120}, {70, 83}, {.078, .083}}, 
  ColorFunction -> (Blend[{Yellow, Orange, Magenta, Cyan}, #3] &), 
  PlotStyle -> PointSize[0.0075]];

pcpl1 = ListPointPlot3D[s1pc, 
  ColorFunction -> (Blend[{Yellow, Orange, Magenta, Cyan}, #2] &), 
  PlotRange -> {{-160, -120}, {70, 83}, {-5.059, 1.5}}, 
  PlotStyle -> PointSize[0.0075]];
pcpl4 = ListPointPlot3D[s4pc, 
  ColorFunction -> (Blend[{Yellow, Orange, Magenta, Cyan}, #2] &), 
  PlotRange -> {{-160, -120}, {70, 83}, {-5.059, 1.5}}, 
  PlotStyle -> PointSize[0.0075]];
pcpl5 = ListPointPlot3D[s5pc, 
  ColorFunction -> (Blend[{Yellow, Orange, Magenta, Cyan}, #2] &), 
  PlotRange -> {{-160, -120}, {70, 83}, {-5.059, 1.5}}, 
  PlotStyle -> PointSize[0.0075]];
pcpl6 = ListPointPlot3D[s6pc, 
  ColorFunction -> (Blend[{Yellow, Orange, Magenta, Cyan}, #2] &), 
  PlotRange -> {{-160, -120}, {70, 83}, {-5.059, 1.5}}, 
  PlotStyle -> PointSize[0.0075]];
pcpl7 = ListPointPlot3D[s7pc, 
  ColorFunction -> (Blend[{Yellow, Orange, Magenta, Cyan}, #2] &), 
  PlotRange -> {{-160, -120}, {70, 83}, {-5.059, 1.5}}, 
  PlotStyle -> PointSize[0.0075]];

p3d = Show[{p1, pspc, pcpl1, pcpl4, pcpl5, pcpl6, pcpl7}];

Do[anim[i_] := 
  Graphics3D[Rotate[p3d[[1]], i*Degree, {0, 0, 1}], Boxed -> False, 
  ImageSize -> 1024, BoxRatios -> {2, 2, .5}, 
  ViewPoint -> {Pi, Pi/2, \[Pi]/2}], {i, 0, 1, 1}]

v = Table[anim[i], {i, 0, 360, 4}];

Export["Desktop/Arctic_ph_excl_final.avi", v, "FrameRate" -> 7]