   mesh
    rectilinear
      nx = 4
      ny = 4
      nz = 4
      bx =  3
      by =  3
      bz =  3
      gmin = -1.0 -1.0 -1.0
      gmax =  1.0  1.0  1.0
    end
    user defined geometry transformation
     "
     double r = sqrt(inxcoord*inxcoord+inycoord*inycoord);
     double theta = atan2(inycoord,inxcoord);
     if(inzcoord > 0.0) 
      {
        theta = theta + (3.14159 / 4.0)*(inzcoord/1.0);
        r = r*(1.0-inzcoord/1.1);
        outxcoord = r*cos(theta);
	outycoord = r*sin(theta);
      }
     "
    end
  end

