added block example

examples/CMakeLists.txt

@@ -9,6 +9,11 @@ target_link_libraries(  EMSchur3D  "lemmacore" "fdem1d" "emschur3d")
 #    target_link_libraries( Interference ${VTK_LIBRARIES})
 #endif()
 
+install ( DIRECTORY mod  DESTINATION ${CMAKE_INSTALL_PREFIX}/share/EMSchur3D/ ) 
+install ( DIRECTORY grd  DESTINATION ${CMAKE_INSTALL_PREFIX}/share/EMSchur3D/ ) 
+install ( DIRECTORY src  DESTINATION ${CMAKE_INSTALL_PREFIX}/share/EMSchur3D/ ) 
+install ( FILES block.sh  DESTINATION ${CMAKE_INSTALL_PREFIX}/share/EMSchur3D/ ) 
+
 INSTALL_TARGETS( "/share/EMSchur3D/"
     EMSchur3D
 )

examples/block.sh

@@ -0,0 +1,2 @@
+./EMSchur3D grd/block50.grd  mod/cond.inp    mod/block_condhs.mod     src/1m.aem  test 
+#./EMSchur3D  grd/block50.grd   mod/res.inp  mod/block.mod   src/1m.aem  kest

examples/grd/block50.grd

@@ -0,0 +1,37 @@
+50	50	25		// nx ny nz
+-245	-245	-45		// ox oy oz
+ 
+10 10 10 10 10 10 10 10 10 10	// hx [0,1,... nz-1]
+10 10 10 10 10 10 10 10 10 10	//     |
+10 10 10 10 10 10 10 10 10 10	//     |
+10 10 10 10 10 10 10 10 10 10	//     |
+10 10 10 10 10 10 10 10 10 10	//     _
+
+10 10 10 10 10 10 10 10 10 10	// hy
+10 10 10 10 10 10 10 10 10 10
+10 10 10 10 10 10 10 10 10 10
+10 10 10 10 10 10 10 10 10 10
+10 10 10 10 10 10 10 10 10 10
+
+10 10 10 10 10 10 10 10 10 10	// hz
+10 10 10 10 10 10 10 10 10 10
+10 10 10 10 10 
+
+// Grid is complete, data fields 
+[cell data]
+$sigma  // 3D array
+$block
+0   .005			// air-earth interface,  sigma half space conductivity 
+nblocks=1
+10  60 				// block dimensions in x
+10  60				// y
+100 140				// z
+0.1 				// block conductivity
+[cell data]
+
+// need to be able to specify type of array scalar vector (# components), cell, edge face, name, type 
+// how to do this general but not too bloated?
+
+// Also, what if we want anisotropy, I guess just sigma_x etc.
+// This file should ONLY contain data that pertains to the grid, so that the background reference 
+// model probably doesn't fit? 

examples/mod/cond.inp

@@ -0,0 +1,10 @@
+1              // Number of NON-AIR layers
+20             // resistivity of bottom layer [Ohm-m]
+
+// 21 layer model
+// thick
+// 11, 5, 5, 4, 3.5, 4.5, 4.5, 5.5, 3, 6, 15, 3, 6, 3, 8, 9, 10,
+//        7, 8, 11;
+// rho 
+//  1, 16., 60, 25, 55, 20, 45, 5, 40, 20, 40, 20, 120, 40, 20, 45, 36, 30,
+//        48, 55, 40, 25;

examples/src/1m.aem

@@ -0,0 +1,108 @@
+// AEM Survey Specification file 
+
+     1    // N sources 
+
+100         // Source 1 frequency Hz
+MD           // Source 1 type
+0  0 0       // Source 1 relative position front/back, left/right, up/down 
+1            // Source 1 moment
+
+// Flight Plan (in same coordinates as grid)
+1 // points 
+// GPS location in x, y, z, polarisation [x,y,z]
+// line 0
+-0     -0      -3e1   0  0  1 
+-140   -100    -3e1   0  0  1 
+-130   -100    -3e1   0  0  1 
+-120   -100    -3e1   0  0  1 
+-110   -100    -3e1   0  0  1 
+-100   -100    -3e1   0  0  1 
+-90    -100    -3e1   0  0  1 
+-80    -100    -3e1   0  0  1 
+-70    -100    -3e1   0  0  1 
+-60    -100    -3e1   0  0  1 
+-50    -100    -3e1   0  0  1 
+-40    -100    -3e1   0  0  1 
+-30    -100    -3e1   0  0  1 
+-20    -100    -3e1   0  0  1 
+-10    -100    -3e1   0  0  1 
+0      -100    -3e1   0  0  1 
+10     -100    -3e1   0  0  1 
+20     -100    -3e1   0  0  1 
+30     -100    -3e1   0  0  1 
+40     -100    -3e1   0  0  1 
+50     -100    -3e1   0  0  1 
+60     -100    -3e1   0  0  1 
+70     -100    -3e1   0  0  1 
+80     -100    -3e1   0  0  1 
+90     -100    -3e1   0  0  1 
+100    -100    -3e1   0  0  1 
+110    -100    -3e1   0  0  1 
+120    -100    -3e1   0  0  1 
+130    -100    -3e1   0  0  1 
+140    -100    -3e1   0  0  1 
+150    -100    -3e1   0  0  1 
+// line 1
+-150    0    -3e1   0  0  1 
+-140    0    -3e1   0  0  1 
+-130    0    -3e1   0  0  1 
+-120    0    -3e1   0  0  1 
+-110    0    -3e1   0  0  1 
+-100    0    -3e1   0  0  1 
+-90     0    -3e1   0  0  1 
+-80     0    -3e1   0  0  1 
+-70     0    -3e1   0  0  1 
+-60     0    -3e1   0  0  1 
+-50     0    -3e1   0  0  1 
+-40     0    -3e1   0  0  1 
+-30     0    -3e1   0  0  1 
+-20     0    -3e1   0  0  1 
+-10     0    -3e1   0  0  1 
+0       0    -3e1   0  0  1 
+10      0    -3e1   0  0  1 
+20      0    -3e1   0  0  1 
+30      0    -3e1   0  0  1 
+40      0    -3e1   0  0  1 
+50      0    -3e1   0  0  1 
+60      0    -3e1   0  0  1 
+70      0    -3e1   0  0  1 
+80      0    -3e1   0  0  1 
+90      0    -3e1   0  0  1 
+100     0    -3e1   0  0  1 
+110     0    -3e1   0  0  1 
+120     0    -3e1   0  0  1 
+130     0    -3e1   0  0  1 
+140     0    -3e1   0  0  1 
+150     0    -3e1   0  0  1 
+// line 2
+-150   100    -3e1   0  0  1 
+-140   100    -3e1   0  0  1 
+-130   100    -3e1   0  0  1 
+-120   100    -3e1   0  0  1 
+-110   100    -3e1   0  0  1 
+-100   100    -3e1   0  0  1 
+-90    100    -3e1   0  0  1 
+-80    100    -3e1   0  0  1 
+-70    100    -3e1   0  0  1 
+-60    100    -3e1   0  0  1 
+-50    100    -3e1   0  0  1 
+-40    100    -3e1   0  0  1 
+-30    100    -3e1   0  0  1 
+-20    100    -3e1   0  0  1 
+-10    100    -3e1   0  0  1 
+0      100    -3e1   0  0  1 
+10     100    -3e1   0  0  1 
+20     100    -3e1   0  0  1 
+30     100    -3e1   0  0  1 
+40     100    -3e1   0  0  1 
+50     100    -3e1   0  0  1 
+60     100    -3e1   0  0  1 
+70     100    -3e1   0  0  1 
+80     100    -3e1   0  0  1 
+90     100    -3e1   0  0  1 
+100    100    -3e1   0  0  1 
+110    100    -3e1   0  0  1 
+120    100    -3e1   0  0  1 
+130    100    -3e1   0  0  1 
+140    100    -3e1   0  0  1 
+150    100    -3e1   0  0  1