Merge pull request #263 from ljcarlin/feature-gmt-sphere-preprocessing

Feature gmt sphere preprocessing

example/gmt/Makefile.am

@@ -12,6 +12,11 @@ example_gmt_p4est_gmt_SOURCES = \
 
 LINT_CSOURCES += \
         $(example_gmt_p4est_gmt_SOURCES)
+
+bin_PROGRAMS += example/gmt/p4est_sphere_preprocessing
+
+example_gmt_p4est_sphere_preprocessing_SOURCES = example/gmt/sphere_preprocessing.c
+
 endif
 
 if P4EST_ENABLE_BUILD_3D

example/gmt/gmt_models.h

@@ -73,6 +73,18 @@ p4est_gmt_model_latlong_params_t;
 p4est_gmt_model_t  *p4est_gmt_model_latlong_new
   (p4est_gmt_model_latlong_params_t * params);
 
+/** Represents a segment of a geodesic in the sphere model.
+ * 
+ * Segments are restricted to lying on a single face of the cube-sphere.
+ * A segment is represented by its endpoints, given in tree-local
+ * reference coordinates.
+ */
+typedef struct p4est_gmt_sphere_geoseg
+{
+  p4est_topidx_t which_tree;
+  double p1x, p1y, p2x, p2y; /* Geodesic endpoints */
+} p4est_gmt_sphere_geoseg_t;
+
 /** Destroy model */
 void                p4est_gmt_model_destroy (p4est_gmt_model_t * model);
 

example/gmt/sphere_hello_world.csv

@@ -0,0 +1,15 @@
+0.0,50.0,0.0,130.0
+50.0,50.0,50.0,130.0
+0.0,90.0,50.0,90.0
+70.0,50.0,70.0,130.0
+70.0,50.0,110.0,50.0
+70.0,130.0,110.0,130.0
+70.0,90.0,110.0,90.0
+140.0,50.0,140.0,130.0
+140.0,130.0,180.0,130.0
+210.0,50.0,210.0,130.0
+210.0,130.0,250.0,130.0
+280.0,50.0,280.0,130.0
+320.0,50.0,320.0,130.0
+280.0,50.0,320.0,50.0
+280.0,130.0,320.0,130.0