Support interpolating to obtain the image values. This makes it

possible to use IOJpeg together with mesh refinement.

There is a new keyword parameter "gridpoints" with possible values
"hyperslab" and "interpolate".  Hyperslabbing works only on uniform
grids, whereas interpolation also works with AMR.  Additional
parameters select where in the domain the points should be
interpolated.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusIO/IOJpeg/trunk@130 eff87b29-5268-4891-90a3-a07138403961

1 files changed, 93 insertions, 5 deletions

diff --git a/param.ccl b/param.ccl
index 15d8e89..9fdad3f 100644
--- a/param.ccl
+++ b/param.ccl
@@ -37,6 +37,12 @@ KEYWORD mode "Output mode to use" STEERABLE = ALWAYS
   "standard" :: "Generate a file for each out_every timesteps"
 } "standard"
 
+KEYWORD gridpoints "How to access grid points" STEERABLE = RECOVER
+{
+  "hyperslab"   :: "use locations of grid points"
+  "interpolate" :: "interpolate to arbitrary points"
+} "hyperslab"
+
 
 ########################
 # Specific to jpegs
@@ -58,8 +64,9 @@ INT colormap_factor "How to scale float values to rgb color" STEERABLE = ALWAYS
 
 KEYWORD colormap "How to set the colormap" STEERABLE = ALWAYS
 {
-  "auto"   :: "Set automatically using min/max of grid variables"
-  "custom" :: "Set min/max manually"
+  "auto"     :: "Set automatically using min/max of grid variables"
+  "auto-old" :: "Set automatically using min/max of grid variables, using the old reduction interface which is still used by Carpet"
+  "custom"   :: "Set min/max manually"
 } "custom"
 
 REAL colormap_min "minimum value to be mapped to colors" STEERABLE = ALWAYS
@@ -77,9 +84,9 @@ INT refinement_factor "Refine each 2D slice by a certain factor (using interpola
 } 1
 
 
-################################
-# Choosing what planes to output
-################################
+###################################################
+# Choosing what planes to output when hyperslabbing
+###################################################
 REAL out2D_yzplane_x "x-coord for 2D planes in yz" STEERABLE = RECOVER
 {
   *:*      :: "A value between [xmin, xmax]"
@@ -113,6 +120,87 @@ INT out2D_xyplane_zi "z-index (from 0) for 2D planes in xy" STEERABLE = RECOVER
 } -1
 
 
+###################################################
+# Choosing what region to output when interpolating
+###################################################
+
+STRING interpolator_name "Name of the interpolator" STEERABLE=always
+{
+  ".*" :: "must be a registered interpolator"
+} "Lagrange polynomial interpolation"
+
+STRING interpolator_options "Options for the interpolator" STEERABLE=always
+{
+  ".*" :: "must be a valid option specification"
+} "order=2"
+
+STRING interpolator_coordinates "Coordinate system" STEERABLE=always
+{
+  ".*" :: "must be a registered coordinate system"
+} "cart3d"
+
+BOOLEAN multiply_by_radius "Multiply valus by r" STEERABLE=always
+{
+} "no"
+
+
+
+REAL array2d_x0 "Origin" STEERABLE=always
+{
+  *:* :: ""
+} 0.0
+
+REAL array2d_y0 "Origin" STEERABLE=always
+{
+  *:* :: ""
+} 0.0
+
+REAL array2d_z0 "Origin" STEERABLE=always
+{
+  *:* :: ""
+} 0.0
+
+INT array2d_npoints_i "Number of grid points for the 2D grid arrays in the i direction"
+{
+  0:* :: ""
+} 10
+
+REAL array2d_dx_i "Spacing" STEERABLE=always
+{
+  0.0:* :: ""
+} 0.0
+
+REAL array2d_dy_i "Spacing" STEERABLE=always
+{
+  0.0:* :: ""
+} 0.0
+
+REAL array2d_dz_i "Spacing" STEERABLE=always
+{
+  0.0:* :: ""
+} 0.0
+
+INT array2d_npoints_j "Number of grid points for the 2D grid arrays in the j direction"
+{
+  0:* :: ""
+} 10
+
+REAL array2d_dx_j "Spacing" STEERABLE=always
+{
+  0.0:* :: ""
+} 0.0
+
+REAL array2d_dy_j "Spacing" STEERABLE=always
+{
+  0.0:* :: ""
+} 0.0
+
+REAL array2d_dz_j "Spacing" STEERABLE=always
+{
+  0.0:* :: ""
+} 0.0
+
+
 #############################################################################
 ### import IOUtil parameters
 #############################################################################