Second pass at finding points on the level set surface. It now uses

grid arrays and works in parallel. Still needs to be generalized to work
in other modes than full mode and an algorithm has to be implemented
to figure out how many surfaces are actually present in the level set
function.


git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@78 2a26948c-0e4f-0410-aee8-f1d3e353619c

1 files changed, 21 insertions, 7 deletions

diff --git a/param.ccl b/param.ccl
index 16da9a2..60f1b2f 100644
--- a/param.ccl
+++ b/param.ccl
@@ -30,22 +30,22 @@ KEYWORD initial_f "Initial surface choice"
 
 REAL initial_rad "Initial radius of surface"
 {
-(0.0:            :: "Positive please"
+0.0:            :: "Positive please"
 } 1.0
 
 REAL initial_a "Initial a coefficient of ellipsoid"
 {
-(0.0:            :: "Positive please"
+0.0:            :: "Positive please"
 } 1.0
 
 REAL initial_b "Initial b coefficient of ellipsoid"
 {
-(0.0:            :: "Positive please"
+0.0:            :: "Positive please"
 } 1.0
 
 REAL initial_c "Initial c coefficient of ellipsoid"
 {
-(0.0:            :: "Positive please"
+0.0:            :: "Positive please"
 } 1.0
 
 REAL rotation_alpha "Rotation angle around z-axis of ellipsoid"
@@ -90,7 +90,7 @@ REAL cas_b "Initial b coefficient of ovaloid of cassini"
 
 REAL shell_width "Width of the evolution region in units of the grid spacing"
 {
-(0.0:            :: "Positive please"
+0.0:            :: "Positive please"
 } 7.0
 
 BOOLEAN normalize "normalize the derivatives of the level set function"
@@ -133,7 +133,7 @@ KEYWORD re_param_int_method "Integration method in pde re-parametrization"
 
 INT re_param_max_iter "maximum number of iteration in the re-parametrization"
 {
-(0:        :: "Positive please"
+0:        :: "Positive please"
 } 800
 
 KEYWORD pde_differences "Type of finite diffencing used in pde re-parametrization"
@@ -155,7 +155,7 @@ INT reparametrize_every_approx "Re-parametrize every using approx method"
 
 INT last_iteration_number "The starting iteration number for the EH_Finder (last iteration number of the simulation)"
 {
-(0:              :: "Positive please"
+0:              :: "Positive please"
 } 0
 
 INT saved_iteration_every "How often was the numerical data saved"
@@ -163,6 +163,20 @@ INT saved_iteration_every "How often was the numerical data saved"
 1:               :: "Positive please"
 } 1
 
+INT ntheta "Number of points in the theta direction when finding points on the surface"
+{
+1:               :: "Positive and odd please"
+} 51
+
+INT nphi "Number of points in the phi direction when finding points on the surface"
+{
+1:               :: "Positive and odd please"
+} 51
+INT n_array_ghosts "Number of ghost points in the surface grid array"
+{
+1:               :: "Positive please"
+} 1
+
 shares: grid
 
 USES KEYWORD domain