[cig-commits] [commit] master: Add comment (c534a6e)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/a4d34fe709bd538894313ef1f1b5ab7a0e67791f...aad39993908bfd5cc020de4548464f6027d37173

>---------------------------------------------------------------

commit c534a6e5cd8ec095174a3e57c75660620b0eae66
Author: Siqi Zhang <siqi.zhang at mq.edu.au>
Date:   Tue Oct 28 13:26:07 2014 +1100

    Add comment


>---------------------------------------------------------------

c534a6e5cd8ec095174a3e57c75660620b0eae66
 source/velocity_boundary_conditions/gplates.cc | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/source/velocity_boundary_conditions/gplates.cc b/source/velocity_boundary_conditions/gplates.cc
index b0477bc..ba8cc7f 100644
--- a/source/velocity_boundary_conditions/gplates.cc
+++ b/source/velocity_boundary_conditions/gplates.cc
@@ -240,7 +240,10 @@ namespace aspect
             // Currently it would not be necessary to update the grid positions at every timestep
             // since they are not allowed to change. In case we allow this later, do it anyway.
             Tensor<1,3> spherical_position = get_grid_point_position(idx_theta,idx_phi,false);
-            if(idx_theta==0 || idx_theta==n_theta-1)spherical_position[1]=0.;
+            if(idx_theta==0 || idx_theta==n_theta-1)
+              // Gplate gpml files set the longitude of all points at poles to zero instead of vary from -180 to 180 degrees,
+              // we set this here as well to make sure that we get the correct cartesian velocity.
+              spherical_position[1]=0.;
             velocity_positions[idx_theta][idx_phi] = cartesian_surface_coordinates(spherical_position);
             (*velocity_values)[idx_theta][idx_phi] = sphere_to_cart_velocity(spherical_velocities,spherical_position)
                                                      / cmyr_si;