* Create this plane at 0.5 km between points to be compatible with
* the mesh point distribution of 0.5 km.
* hextotet / 4 / will and points at cell centers
*

define/NX/33
define/NY/33
define / CMONAME / cmo
*
* Create the quad surfaces
*
* surface/f12/intrface/plane/0.0,0.0,-4.0/0.0,-4.0,-4.0/-4.0,-4.0,-8.0/
define / PI / 3.14159265358979

cmo/create/CMONAME///quad                                                         
quadxy NX NY/4.0,-4.0,0.0/4.0,4.0,0.0/-4.0,4.0,-8.0/-4.0,-4.0,-8.0
rzbrick/xyz/NX,NY,1/1,0,0/connect
hextotet / 4 / cmo_flt_surf / cmo
cmo / delete / cmo
define / CMONAME / cmo_flt_surf

* Use the math command to set the z coordinates to:
*
* cmo_1   z(i)=sin(y(i))
* cmo_2   z(i)=sin(x(i))
*
cmo / addatt / CMONAME / tempvar1 / vdouble / scalar / nnodes
cmo / addatt / CMONAME / tempvar2 / vdouble / scalar / nnodes
cmo / addatt / CMONAME / tempvar3 / vdouble / scalar / nnodes
cmo / copyatt / CMONAME / CMONAME / tempvar1 / xic /
* Scale the variable to go -1 to  1
math / multiply / CMONAME / tempvar1 / 1 0 0 / CMONAME / tempvar1 / .25
cmo / printatt / CMONAME / tempvar1 / minmax
* Scale the variable to go -PI to PI
math / multiply / CMONAME / tempvar1 / 1 0 0 / CMONAME / tempvar1 / PI
cmo / printatt / CMONAME / tempvar1 / minmax
* Now scale to determine how many cycles of -PI to PI
math / multiply / CMONAME / tempvar1 / 1 0 0 / CMONAME / tempvar1 / 2.0
cmo / printatt / CMONAME / tempvar1 / minmax
math/sin/ CMONAME / tempvar1 /1,0,0/ CMONAME / tempvar1

cmo / copyatt / CMONAME / CMONAME / tempvar2 / yic /
* Scale the variable to go -1 to  1
math / multiply / CMONAME / tempvar2 / 1 0 0 / CMONAME / tempvar2 / .25
* Scale the variable to go -PI to PI
math / multiply / CMONAME / tempvar2 / 1 0 0 / CMONAME / tempvar2 / PI
* Now scale to determine how many cycles of -PI to PI
math / multiply / CMONAME / tempvar2 / 1 0 0 / CMONAME / tempvar2 / 2.0

math/sin/ CMONAME / tempvar2 /1,0,0/ CMONAME / tempvar2
math / multiply / CMONAME / tempvar3 / 1 0 0 / &
                  CMONAME / tempvar1 / CMONAME / tempvar2 /
*
* The magnitude of tempvar3 is -1 to 1. For this problem the
* units are km. 
* Scale to variation is 0.01 which is 10m.
* Scale to variation is 0.02 which is 20m.
* Scale to variation is 0.05 which is 50m.
*
math / multiply / CMONAME / tempvar3 / 1 0 0 / CMONAME / tempvar3 / 0.05
math / add / CMONAME / zic / 1 0 0 / CMONAME / zic / CMONAME / tempvar3
dump / gmv / tmp.gmv / CMONAME

cmo / select / cube1
surface/f12/intrface/sheet/CMONAME
cmo / select / cube1

finish