[Gmsh] Could not orient surface normal ?

[Aleksejs Fomins]

Dear GMSH,

I have a model for a cavity simulation necessary for my PhD.

2D mesh works fine, 3D mesh throws a bunch of warnings like

  Warning : Converged for i=7 j=5 (err=1.81843e-08 iter=2) BUT xyz error = 12.5473 in point (-1.604927e+01,-5.104927e+01,4.445920e+01) on surface 123

and produces a completely rubbish mesh. I attach the geometry and the full error log.

Any advice would be appreciated.

Kind regards,
Aleksejs Fomins
-------------- next part --------------
/** \file
    \brief This is a mesh of a rectangular cavity in a gold sheet.
    
    This mesh is part of the series developed within the PhD thesis of Aleksejs Fomins.                       
 
    Copyright by LSPR AG, Aleksejs Fomins, 2002-2015. All rights reserved.
       
 */
 
 
 
/***************************************************/
/** given physical dimensions                     **/
/***************************************************/

cavitySizeX     = 50.0;			/** Length of the cavity  */
cavitySizeY     = 100.0;		/** Width of the cavity  */
cavityCurvTop   = 10.0;			/** Radius of curvature of the top edges */
cavitySizeCurvX = 10.0;			/** Radius of x-curvature of vertical edges */
cavitySizeCurvY = 10.0;			/** Radius of y-curvature of vertical edges */

layer1SizeZ = 50.0; 			/** Thickness of the metal layer (Au or Ag)  */
liquidSizeZ = 200.0; 			/** Thickness of the liquid filling  */

domainSizeX = 600.0;			/** Domain length */
domainSizeY = 600.0;			/** Domain width */
domainSizeZ = 600.0;			/** Domain height */

substrateSizeZ = 200.0;			/** Thickness of the glass substrate **/

/** Geometrical position of the center of the cavity within the domain **/
centerx = 0.0;
centery = 0.0;
centerz = 0.0;

/**************************************************/
/* Level Structure: z axis bottom up
 * 000  : Base surface of the geometry
 * 300  : Ceiling Glass Substrate
 * 400  : Ceiling Layer 1 (Curvature base)
 * 450  : Ceiling Layer 1 (Curvature top)
 * 1000 : Ceiling of the liquid filling
 * 2000 : Ceiling of the geometry
 * 
 * Layer Structure:
 * 50   : Cavity Boundary
 * 55   : Cavity Boundary Curvature (Top, Horizontal)
 * 56   : Cavity Boundary Curvature (Vertical)
 * 2000 : Domain Boundary
 * 2001 : Fictitious vertices on the boundary
 * 
/**************************************************/

/**************************************************/
/* characteristic length scales                   */
/**************************************************/

/** for the pillar mesh **/
clscavity	= 0.5 * layer1SizeZ;

clslv300ly50	= 0.5 * layer1SizeZ;
clslv400ly50	= 0.1 * layer1SizeZ;
clslv450ly50	= 0.1 * layer1SizeZ;

/** for the computational domain tessellation **/
clslv000ly2000      =   0.5 * substrateSizeZ;
clslv300ly2000      =   0.5 * substrateSizeZ;
clslv450ly2000      =   0.5 * layer1SizeZ;
clslv1000ly2000     =   0.5 * liquidSizeZ;
clslv2000ly2000     =   0.5 * (domainSizeZ - substrateSizeZ);

/***************************************************/
/** derived physical dimensions                   **/
/***************************************************/

cavityRadX = 0.5 * cavitySizeX;
cavityRadY = 0.5 * cavitySizeY;
cavityRadXext = cavityRadX + cavityCurvTop;
cavityRadYext = cavityRadY + cavityCurvTop;

/** Corners of the cavity base  */
cavityXMin1 = centerx - (cavityRadX - cavitySizeCurvX);  cavityXMin2 = centerx - cavityRadX;
cavityXMax1 = centerx + (cavityRadX - cavitySizeCurvX);  cavityXMax2 = centerx + cavityRadX;
cavityYMin1 = centery - (cavityRadY - cavitySizeCurvY);  cavityYMin2 = centery - cavityRadY;
cavityYMax1 = centery + (cavityRadY - cavitySizeCurvY);  cavityYMax2 = centery + cavityRadY;

/** Corners of the cavity curvature-extended  */
cavityXMin1ext = centerx - (cavityRadXext - cavitySizeCurvX);  cavityXMin2ext = centerx - cavityRadXext;
cavityXMax1ext = centerx + (cavityRadXext - cavitySizeCurvX);  cavityXMax2ext = centerx + cavityRadXext;
cavityYMin1ext = centery - (cavityRadYext - cavitySizeCurvY);  cavityYMin2ext = centery - cavityRadYext;
cavityYMax1ext = centery + (cavityRadYext - cavitySizeCurvY);  cavityYMax2ext = centery + cavityRadYext;

/** Corners of the domain  */
domainXmin = centerx - 0.5 * domainSizeX;
domainXmax = centerx + 0.5 * domainSizeX;
domainYmin = centery - 0.5 * domainSizeY;
domainYmax = centery + 0.5 * domainSizeY;

/** Additional fictitious points on the domain for the fictitious diagonal lines */
domainHalfX = domainXmin + domainSizeX * 0.5;
domainHalfY = domainYmin + domainSizeY * 0.5;

/** Heights of layers */
height000  = centerz - substrateSizeZ;
height300  = centerz;
height400  = centerz + layer1SizeZ - cavityCurvTop;
height450  = centerz + layer1SizeZ;
height1000 = centerz + layer1SizeZ + liquidSizeZ;
height2000 = centerz + domainSizeZ - substrateSizeZ;

/***************************************************/
/** physical domain tags - volume                 **/
/***************************************************/

physvollayer1		= 900;
physvolliquid		= 5000;
physvolsubstrate	= 6000;
physvolair		= 8000;
physsurfairboundary	= 9000;







/*****************************************************************************************************/
/** Points                                                            **/
/*****************************************************************************************************/

// level 300 - layer 50
vrtx100lv300ly50 = newp; Point(vrtx100lv300ly50) = { cavityXMin2, cavityYMin1, height300, clslv300ly50};
vrtx200lv300ly50 = newp; Point(vrtx200lv300ly50) = { cavityXMin1, cavityYMin2, height300, clslv300ly50};
vrtx300lv300ly50 = newp; Point(vrtx300lv300ly50) = { cavityXMax1, cavityYMin2, height300, clslv300ly50};
vrtx400lv300ly50 = newp; Point(vrtx400lv300ly50) = { cavityXMax2, cavityYMin1, height300, clslv300ly50};
vrtx500lv300ly50 = newp; Point(vrtx500lv300ly50) = { cavityXMax2, cavityYMax1, height300, clslv300ly50};
vrtx600lv300ly50 = newp; Point(vrtx600lv300ly50) = { cavityXMax1, cavityYMax2, height300, clslv300ly50};
vrtx700lv300ly50 = newp; Point(vrtx700lv300ly50) = { cavityXMin1, cavityYMax2, height300, clslv300ly50};
vrtx800lv300ly50 = newp; Point(vrtx800lv300ly50) = { cavityXMin2, cavityYMax1, height300, clslv300ly50};

// Level 300 - layer 56 - Interior curvature of vertical cavity edges
vrtx100lv300ly56 = newp; Point(vrtx100lv300ly56) = { cavityXMin1, cavityYMin1, height300, clslv300ly50};
vrtx300lv300ly56 = newp; Point(vrtx300lv300ly56) = { cavityXMax1, cavityYMin1, height300, clslv300ly50};
vrtx500lv300ly56 = newp; Point(vrtx500lv300ly56) = { cavityXMax1, cavityYMax1, height300, clslv300ly50};
vrtx700lv300ly56 = newp; Point(vrtx700lv300ly56) = { cavityXMin1, cavityYMax1, height300, clslv300ly50};

// level 400 - layer 50
vrtx100lv400ly50 = newp; Point(vrtx100lv400ly50) = { cavityXMin2, cavityYMin1, height400, clslv400ly50};
vrtx200lv400ly50 = newp; Point(vrtx200lv400ly50) = { cavityXMin1, cavityYMin2, height400, clslv400ly50};
vrtx300lv400ly50 = newp; Point(vrtx300lv400ly50) = { cavityXMax1, cavityYMin2, height400, clslv400ly50};
vrtx400lv400ly50 = newp; Point(vrtx400lv400ly50) = { cavityXMax2, cavityYMin1, height400, clslv400ly50};
vrtx500lv400ly50 = newp; Point(vrtx500lv400ly50) = { cavityXMax2, cavityYMax1, height400, clslv400ly50};
vrtx600lv400ly50 = newp; Point(vrtx600lv400ly50) = { cavityXMax1, cavityYMax2, height400, clslv400ly50};
vrtx700lv400ly50 = newp; Point(vrtx700lv400ly50) = { cavityXMin1, cavityYMax2, height400, clslv400ly50};
vrtx800lv400ly50 = newp; Point(vrtx800lv400ly50) = { cavityXMin2, cavityYMax1, height400, clslv400ly50};

// level 400 - layer 55 - Interior curvature of the cavity top edge
vrtx100lv400ly55 = newp; Point(vrtx100lv400ly55) = { cavityXMin2ext, cavityYMin1ext, height400, clslv400ly50};
vrtx200lv400ly55 = newp; Point(vrtx200lv400ly55) = { cavityXMin1ext, cavityYMin2ext, height400, clslv400ly50};
vrtx300lv400ly55 = newp; Point(vrtx300lv400ly55) = { cavityXMax1ext, cavityYMin2ext, height400, clslv400ly50};
vrtx400lv400ly55 = newp; Point(vrtx400lv400ly55) = { cavityXMax2ext, cavityYMin1ext, height400, clslv400ly50};
vrtx500lv400ly55 = newp; Point(vrtx500lv400ly55) = { cavityXMax2ext, cavityYMax1ext, height400, clslv400ly50};
vrtx600lv400ly55 = newp; Point(vrtx600lv400ly55) = { cavityXMax1ext, cavityYMax2ext, height400, clslv400ly50};
vrtx700lv400ly55 = newp; Point(vrtx700lv400ly55) = { cavityXMin1ext, cavityYMax2ext, height400, clslv400ly50};
vrtx800lv400ly55 = newp; Point(vrtx800lv400ly55) = { cavityXMin2ext, cavityYMax1ext, height400, clslv400ly50};

// Level 400 - layer 56 - Interior curvature of vertical cavity edges
vrtx100lv400ly56 = newp; Point(vrtx100lv400ly56) = { cavityXMin1, cavityYMin1, height400, clslv400ly50};
vrtx300lv400ly56 = newp; Point(vrtx300lv400ly56) = { cavityXMax1, cavityYMin1, height400, clslv400ly50};
vrtx500lv400ly56 = newp; Point(vrtx500lv400ly56) = { cavityXMax1, cavityYMax1, height400, clslv400ly50};
vrtx700lv400ly56 = newp; Point(vrtx700lv400ly56) = { cavityXMin1, cavityYMax1, height400, clslv400ly50};

// level 450 - layer 50
vrtx100lv450ly50 = newp; Point(vrtx100lv450ly50) = { cavityXMin2ext, cavityYMin1ext, height450, clslv450ly50};
vrtx200lv450ly50 = newp; Point(vrtx200lv450ly50) = { cavityXMin1ext, cavityYMin2ext, height450, clslv450ly50};
vrtx300lv450ly50 = newp; Point(vrtx300lv450ly50) = { cavityXMax1ext, cavityYMin2ext, height450, clslv450ly50};
vrtx400lv450ly50 = newp; Point(vrtx400lv450ly50) = { cavityXMax2ext, cavityYMin1ext, height450, clslv450ly50};
vrtx500lv450ly50 = newp; Point(vrtx500lv450ly50) = { cavityXMax2ext, cavityYMax1ext, height450, clslv450ly50};
vrtx600lv450ly50 = newp; Point(vrtx600lv450ly50) = { cavityXMax1ext, cavityYMax2ext, height450, clslv450ly50};
vrtx700lv450ly50 = newp; Point(vrtx700lv450ly50) = { cavityXMin1ext, cavityYMax2ext, height450, clslv450ly50};
vrtx800lv450ly50 = newp; Point(vrtx800lv450ly50) = { cavityXMin2ext, cavityYMax1ext, height450, clslv450ly50};

// Level 450 - layer 56 - Interior curvature of vertical cavity edges
vrtx100lv450ly56 = newp; Point(vrtx100lv450ly56) = { cavityXMin1ext, cavityYMin1ext, height450, clslv450ly50};
vrtx300lv450ly56 = newp; Point(vrtx300lv450ly56) = { cavityXMax1ext, cavityYMin1ext, height450, clslv450ly50};
vrtx500lv450ly56 = newp; Point(vrtx500lv450ly56) = { cavityXMax1ext, cavityYMax1ext, height450, clslv450ly50};
vrtx700lv450ly56 = newp; Point(vrtx700lv450ly56) = { cavityXMin1ext, cavityYMax1ext, height450, clslv450ly50};


// level 000 - layer 2000
vrtx100lv000ly2000 = newp; Point(vrtx100lv000ly2000) = { domainXmin, domainYmin, height000,  clslv000ly2000};
vrtx200lv000ly2000 = newp; Point(vrtx200lv000ly2000) = { domainXmax, domainYmin, height000,  clslv000ly2000};
vrtx300lv000ly2000 = newp; Point(vrtx300lv000ly2000) = { domainXmax, domainYmax, height000,  clslv000ly2000};
vrtx400lv000ly2000 = newp; Point(vrtx400lv000ly2000) = { domainXmin, domainYmax, height000,  clslv000ly2000};

// level 300 - layer 2000
vrtx100lv300ly2000 = newp; Point(vrtx100lv300ly2000) = { domainXmin, domainYmin, height300,  clslv300ly2000};
vrtx200lv300ly2000 = newp; Point(vrtx200lv300ly2000) = { domainXmax, domainYmin, height300,  clslv300ly2000};
vrtx300lv300ly2000 = newp; Point(vrtx300lv300ly2000) = { domainXmax, domainYmax, height300,  clslv300ly2000};
vrtx400lv300ly2000 = newp; Point(vrtx400lv300ly2000) = { domainXmin, domainYmax, height300,  clslv300ly2000};

// level 300 - layer 2001
vrtx100lv300ly2001 = newp; Point(vrtx100lv300ly2001) = { domainHalfX, domainYmin,  height300,  clslv300ly2000};
vrtx200lv300ly2001 = newp; Point(vrtx200lv300ly2001) = { domainXmax,  domainHalfY, height300,  clslv300ly2000};
vrtx300lv300ly2001 = newp; Point(vrtx300lv300ly2001) = { domainHalfX, domainYmax,  height300,  clslv300ly2000};
vrtx400lv300ly2001 = newp; Point(vrtx400lv300ly2001) = { domainXmin,  domainHalfY, height300,  clslv300ly2000};

// level 450 - layer 2000
vrtx100lv450ly2000 = newp; Point(vrtx100lv450ly2000) = { domainXmin, domainYmin, height450,  clslv450ly2000};
vrtx200lv450ly2000 = newp; Point(vrtx200lv450ly2000) = { domainXmax, domainYmin, height450,  clslv450ly2000};
vrtx300lv450ly2000 = newp; Point(vrtx300lv450ly2000) = { domainXmax, domainYmax, height450,  clslv450ly2000};
vrtx400lv450ly2000 = newp; Point(vrtx400lv450ly2000) = { domainXmin, domainYmax, height450,  clslv450ly2000};

// level 450 - layer 2001
vrtx100lv450ly2001 = newp; Point(vrtx100lv450ly2001) = { domainHalfX, domainYmin,  height450,  clslv450ly2000};
vrtx200lv450ly2001 = newp; Point(vrtx200lv450ly2001) = { domainXmax,  domainHalfY, height450,  clslv450ly2000};
vrtx300lv450ly2001 = newp; Point(vrtx300lv450ly2001) = { domainHalfX, domainYmax,  height450,  clslv450ly2000};
vrtx400lv450ly2001 = newp; Point(vrtx400lv450ly2001) = { domainXmin,  domainHalfY, height450,  clslv450ly2000};

// level 1000 - layer 2000
vrtx100lv1000ly2000 = newp; Point(vrtx100lv1000ly2000) = { domainXmin, domainYmin, height1000,  clslv1000ly2000};
vrtx200lv1000ly2000 = newp; Point(vrtx200lv1000ly2000) = { domainXmax, domainYmin, height1000,  clslv1000ly2000};
vrtx300lv1000ly2000 = newp; Point(vrtx300lv1000ly2000) = { domainXmax, domainYmax, height1000,  clslv1000ly2000};
vrtx400lv1000ly2000 = newp; Point(vrtx400lv1000ly2000) = { domainXmin, domainYmax, height1000,  clslv1000ly2000};

// level 2000 - layer 2000
vrtx100lv2000ly2000 = newp; Point(vrtx100lv2000ly2000) = { domainXmin, domainYmin, height2000,  clslv2000ly2000};
vrtx200lv2000ly2000 = newp; Point(vrtx200lv2000ly2000) = { domainXmax, domainYmin, height2000,  clslv2000ly2000};
vrtx300lv2000ly2000 = newp; Point(vrtx300lv2000ly2000) = { domainXmax, domainYmax, height2000,  clslv2000ly2000};
vrtx400lv2000ly2000 = newp; Point(vrtx400lv2000ly2000) = { domainXmin, domainYmax, height2000,  clslv2000ly2000};






/*************************************************************************/
/** Lines - cavity                                                      **/
/*************************************************************************/

// horizontal lines - level 300 - layer 50
hlin100lv300ly50 = newl; Ellipse ( hlin100lv300ly50 ) = {vrtx100lv300ly50, vrtx100lv300ly56, vrtx100lv300ly50, vrtx200lv300ly50};
hlin300lv300ly50 = newl; Ellipse ( hlin300lv300ly50 ) = {vrtx300lv300ly50, vrtx300lv300ly56, vrtx300lv300ly50, vrtx400lv300ly50};
hlin500lv300ly50 = newl; Ellipse ( hlin500lv300ly50 ) = {vrtx500lv300ly50, vrtx500lv300ly56, vrtx500lv300ly50, vrtx600lv300ly50};
hlin700lv300ly50 = newl; Ellipse ( hlin700lv300ly50 ) = {vrtx700lv300ly50, vrtx700lv300ly56, vrtx700lv300ly50, vrtx800lv300ly50};
hlin200lv300ly50 = newl; Line ( hlin200lv300ly50 ) = {vrtx200lv300ly50, vrtx300lv300ly50};
hlin400lv300ly50 = newl; Line ( hlin400lv300ly50 ) = {vrtx400lv300ly50, vrtx500lv300ly50};
hlin600lv300ly50 = newl; Line ( hlin600lv300ly50 ) = {vrtx600lv300ly50, vrtx700lv300ly50};
hlin800lv300ly50 = newl; Line ( hlin800lv300ly50 ) = {vrtx800lv300ly50, vrtx100lv300ly50};

// horizontal lines - level 400 - layer 50
hlin100lv400ly50 = newl; Ellipse ( hlin100lv400ly50 ) = {vrtx100lv400ly50, vrtx100lv400ly56, vrtx100lv400ly50, vrtx200lv400ly50};
hlin300lv400ly50 = newl; Ellipse ( hlin300lv400ly50 ) = {vrtx300lv400ly50, vrtx300lv400ly56, vrtx300lv400ly50, vrtx400lv400ly50};
hlin500lv400ly50 = newl; Ellipse ( hlin500lv400ly50 ) = {vrtx500lv400ly50, vrtx500lv400ly56, vrtx500lv400ly50, vrtx600lv400ly50};
hlin700lv400ly50 = newl; Ellipse ( hlin700lv400ly50 ) = {vrtx700lv400ly50, vrtx700lv400ly56, vrtx700lv400ly50, vrtx800lv400ly50};
hlin200lv400ly50 = newl; Line ( hlin200lv400ly50 ) = {vrtx200lv400ly50, vrtx300lv400ly50};
hlin400lv400ly50 = newl; Line ( hlin400lv400ly50 ) = {vrtx400lv400ly50, vrtx500lv400ly50};
hlin600lv400ly50 = newl; Line ( hlin600lv400ly50 ) = {vrtx600lv400ly50, vrtx700lv400ly50};
hlin800lv400ly50 = newl; Line ( hlin800lv400ly50 ) = {vrtx800lv400ly50, vrtx100lv400ly50};

// horizontal lines - level 450 - layer 50
hlin100lv450ly50 = newl; Ellipse ( hlin100lv450ly50 ) = {vrtx100lv450ly50, vrtx100lv450ly56, vrtx100lv450ly50, vrtx200lv450ly50};
hlin300lv450ly50 = newl; Ellipse ( hlin300lv450ly50 ) = {vrtx300lv450ly50, vrtx300lv450ly56, vrtx300lv450ly50, vrtx400lv450ly50};
hlin500lv450ly50 = newl; Ellipse ( hlin500lv450ly50 ) = {vrtx500lv450ly50, vrtx500lv450ly56, vrtx500lv450ly50, vrtx600lv450ly50};
hlin700lv450ly50 = newl; Ellipse ( hlin700lv450ly50 ) = {vrtx700lv450ly50, vrtx700lv450ly56, vrtx700lv450ly50, vrtx800lv450ly50};
hlin200lv450ly50 = newl; Line ( hlin200lv450ly50 ) = {vrtx200lv450ly50, vrtx300lv450ly50};
hlin400lv450ly50 = newl; Line ( hlin400lv450ly50 ) = {vrtx400lv450ly50, vrtx500lv450ly50};
hlin600lv450ly50 = newl; Line ( hlin600lv450ly50 ) = {vrtx600lv450ly50, vrtx700lv450ly50};
hlin800lv450ly50 = newl; Line ( hlin800lv450ly50 ) = {vrtx800lv450ly50, vrtx100lv450ly50};


// vertical lines - from level 300 -> 400 - layer 50
vlin100lv300ly50 = newl; Line ( vlin100lv300ly50 ) = {vrtx100lv300ly50, vrtx100lv400ly50};
vlin200lv300ly50 = newl; Line ( vlin200lv300ly50 ) = {vrtx200lv300ly50, vrtx200lv400ly50};
vlin300lv300ly50 = newl; Line ( vlin300lv300ly50 ) = {vrtx300lv300ly50, vrtx300lv400ly50};
vlin400lv300ly50 = newl; Line ( vlin400lv300ly50 ) = {vrtx400lv300ly50, vrtx400lv400ly50};
vlin500lv300ly50 = newl; Line ( vlin500lv300ly50 ) = {vrtx500lv300ly50, vrtx500lv400ly50};
vlin600lv300ly50 = newl; Line ( vlin600lv300ly50 ) = {vrtx600lv300ly50, vrtx600lv400ly50};
vlin700lv300ly50 = newl; Line ( vlin700lv300ly50 ) = {vrtx700lv300ly50, vrtx700lv400ly50};
vlin800lv300ly50 = newl; Line ( vlin800lv300ly50 ) = {vrtx800lv300ly50, vrtx800lv400ly50};

// vertical curved lines - from level 400 -> 450 - layer 50
vlincurv100lv400ly50 = newl; Ellipse ( vlincurv100lv400ly50 ) = {vrtx100lv400ly50, vrtx100lv400ly55, vrtx100lv400ly50, vrtx100lv450ly50};
vlincurv200lv400ly50 = newl; Ellipse ( vlincurv200lv400ly50 ) = {vrtx200lv400ly50, vrtx200lv400ly55, vrtx200lv400ly50, vrtx200lv450ly50};
vlincurv300lv400ly50 = newl; Ellipse ( vlincurv300lv400ly50 ) = {vrtx300lv400ly50, vrtx300lv400ly55, vrtx300lv400ly50, vrtx300lv450ly50};
vlincurv400lv400ly50 = newl; Ellipse ( vlincurv400lv400ly50 ) = {vrtx400lv400ly50, vrtx400lv400ly55, vrtx400lv400ly50, vrtx400lv450ly50};
vlincurv500lv400ly50 = newl; Ellipse ( vlincurv500lv400ly50 ) = {vrtx500lv400ly50, vrtx500lv400ly55, vrtx500lv400ly50, vrtx500lv450ly50};
vlincurv600lv400ly50 = newl; Ellipse ( vlincurv600lv400ly50 ) = {vrtx600lv400ly50, vrtx600lv400ly55, vrtx600lv400ly50, vrtx600lv450ly50};
vlincurv700lv400ly50 = newl; Ellipse ( vlincurv700lv400ly50 ) = {vrtx700lv400ly50, vrtx700lv400ly55, vrtx700lv400ly50, vrtx700lv450ly50};
vlincurv800lv400ly50 = newl; Ellipse ( vlincurv800lv400ly50 ) = {vrtx800lv400ly50, vrtx800lv400ly55, vrtx800lv400ly50, vrtx800lv450ly50};


/*************************************************************************/
/** Lines - substrate                                                   **/
/*************************************************************************/

// level 300
diaghlin100lv300ly2000 = newl; Line ( diaghlin100lv300ly2000 ) = {vrtx100lv300ly2000, vrtx200lv300ly50};
diaghlin200lv300ly2000 = newl; Line ( diaghlin200lv300ly2000 ) = {vrtx100lv300ly2001, vrtx300lv300ly50};
diaghlin300lv300ly2000 = newl; Line ( diaghlin300lv300ly2000 ) = {vrtx200lv300ly2000, vrtx400lv300ly50};
diaghlin400lv300ly2000 = newl; Line ( diaghlin400lv300ly2000 ) = {vrtx200lv300ly2001, vrtx500lv300ly50};
diaghlin500lv300ly2000 = newl; Line ( diaghlin500lv300ly2000 ) = {vrtx300lv300ly2000, vrtx600lv300ly50};
diaghlin600lv300ly2000 = newl; Line ( diaghlin600lv300ly2000 ) = {vrtx300lv300ly2001, vrtx700lv300ly50};
diaghlin700lv300ly2000 = newl; Line ( diaghlin700lv300ly2000 ) = {vrtx400lv300ly2000, vrtx800lv300ly50};
diaghlin800lv300ly2000 = newl; Line ( diaghlin800lv300ly2000 ) = {vrtx400lv300ly2001, vrtx100lv300ly50};

// level 450
diaghlin100lv450ly2000 = newl; Line ( diaghlin100lv450ly2000 ) = {vrtx100lv450ly2000, vrtx200lv450ly50};
diaghlin200lv450ly2000 = newl; Line ( diaghlin200lv450ly2000 ) = {vrtx100lv450ly2001, vrtx300lv450ly50};
diaghlin300lv450ly2000 = newl; Line ( diaghlin300lv450ly2000 ) = {vrtx200lv450ly2000, vrtx400lv450ly50};
diaghlin400lv450ly2000 = newl; Line ( diaghlin400lv450ly2000 ) = {vrtx200lv450ly2001, vrtx500lv450ly50};
diaghlin500lv450ly2000 = newl; Line ( diaghlin500lv450ly2000 ) = {vrtx300lv450ly2000, vrtx600lv450ly50};
diaghlin600lv450ly2000 = newl; Line ( diaghlin600lv450ly2000 ) = {vrtx300lv450ly2001, vrtx700lv450ly50};
diaghlin700lv450ly2000 = newl; Line ( diaghlin700lv450ly2000 ) = {vrtx400lv450ly2000, vrtx800lv450ly50};
diaghlin800lv450ly2000 = newl; Line ( diaghlin800lv450ly2000 ) = {vrtx400lv450ly2001, vrtx100lv450ly50};


/*************************************************************************/
/** Lines - computational domain                                        **/
/*************************************************************************/

// horizontal - level 000
hlin100lv000ly2000 = newl; Line ( hlin100lv000ly2000 ) = {vrtx100lv000ly2000,vrtx200lv000ly2000};
hlin200lv000ly2000 = newl; Line ( hlin200lv000ly2000 ) = {vrtx200lv000ly2000,vrtx300lv000ly2000};
hlin300lv000ly2000 = newl; Line ( hlin300lv000ly2000 ) = {vrtx300lv000ly2000,vrtx400lv000ly2000};
hlin400lv000ly2000 = newl; Line ( hlin400lv000ly2000 ) = {vrtx400lv000ly2000,vrtx100lv000ly2000};

// horizontal - level 300
hlin100lv300ly2000 = newl; Line ( hlin100lv300ly2000 ) = {vrtx100lv300ly2000,vrtx200lv300ly2000};
hlin200lv300ly2000 = newl; Line ( hlin200lv300ly2000 ) = {vrtx200lv300ly2000,vrtx300lv300ly2000};
hlin300lv300ly2000 = newl; Line ( hlin300lv300ly2000 ) = {vrtx300lv300ly2000,vrtx400lv300ly2000};
hlin400lv300ly2000 = newl; Line ( hlin400lv300ly2000 ) = {vrtx400lv300ly2000,vrtx100lv300ly2000};

// horizontal - level 300 - layer 2001
hlin100lv300ly2001 = newl; Line ( hlin100lv300ly2001 ) = {vrtx100lv300ly2000, vrtx100lv300ly2001};
hlin200lv300ly2001 = newl; Line ( hlin200lv300ly2001 ) = {vrtx100lv300ly2001, vrtx200lv300ly2000};
hlin300lv300ly2001 = newl; Line ( hlin300lv300ly2001 ) = {vrtx200lv300ly2000, vrtx200lv300ly2001};
hlin400lv300ly2001 = newl; Line ( hlin400lv300ly2001 ) = {vrtx200lv300ly2001, vrtx300lv300ly2000};
hlin500lv300ly2001 = newl; Line ( hlin500lv300ly2001 ) = {vrtx300lv300ly2000, vrtx300lv300ly2001};
hlin600lv300ly2001 = newl; Line ( hlin600lv300ly2001 ) = {vrtx300lv300ly2001, vrtx400lv300ly2000};
hlin700lv300ly2001 = newl; Line ( hlin700lv300ly2001 ) = {vrtx400lv300ly2000, vrtx400lv300ly2001};
hlin800lv300ly2001 = newl; Line ( hlin800lv300ly2001 ) = {vrtx400lv300ly2001, vrtx100lv300ly2000};

// horizontal - level 450
hlin100lv450ly2000 = newl; Line ( hlin100lv450ly2000 ) = {vrtx100lv450ly2000,vrtx200lv450ly2000};
hlin200lv450ly2000 = newl; Line ( hlin200lv450ly2000 ) = {vrtx200lv450ly2000,vrtx300lv450ly2000};
hlin300lv450ly2000 = newl; Line ( hlin300lv450ly2000 ) = {vrtx300lv450ly2000,vrtx400lv450ly2000};
hlin400lv450ly2000 = newl; Line ( hlin400lv450ly2000 ) = {vrtx400lv450ly2000,vrtx100lv450ly2000};

// horizontal - level 450 - layer 2001
hlin100lv450ly2001 = newl; Line ( hlin100lv450ly2001 ) = {vrtx100lv450ly2000, vrtx100lv450ly2001};
hlin200lv450ly2001 = newl; Line ( hlin200lv450ly2001 ) = {vrtx100lv450ly2001, vrtx200lv450ly2000};
hlin300lv450ly2001 = newl; Line ( hlin300lv450ly2001 ) = {vrtx200lv450ly2000, vrtx200lv450ly2001};
hlin400lv450ly2001 = newl; Line ( hlin400lv450ly2001 ) = {vrtx200lv450ly2001, vrtx300lv450ly2000};
hlin500lv450ly2001 = newl; Line ( hlin500lv450ly2001 ) = {vrtx300lv450ly2000, vrtx300lv450ly2001};
hlin600lv450ly2001 = newl; Line ( hlin600lv450ly2001 ) = {vrtx300lv450ly2001, vrtx400lv450ly2000};
hlin700lv450ly2001 = newl; Line ( hlin700lv450ly2001 ) = {vrtx400lv450ly2000, vrtx400lv450ly2001};
hlin800lv450ly2001 = newl; Line ( hlin800lv450ly2001 ) = {vrtx400lv450ly2001, vrtx100lv450ly2000};

// horizontal - level 1000
hlin100lv1000ly2000 = newl; Line ( hlin100lv1000ly2000 ) = {vrtx100lv1000ly2000,vrtx200lv1000ly2000};
hlin200lv1000ly2000 = newl; Line ( hlin200lv1000ly2000 ) = {vrtx200lv1000ly2000,vrtx300lv1000ly2000};
hlin300lv1000ly2000 = newl; Line ( hlin300lv1000ly2000 ) = {vrtx300lv1000ly2000,vrtx400lv1000ly2000};
hlin400lv1000ly2000 = newl; Line ( hlin400lv1000ly2000 ) = {vrtx400lv1000ly2000,vrtx100lv1000ly2000};

// horizontal - level 2000
hlin100lv2000ly2000 = newl; Line ( hlin100lv2000ly2000 ) = {vrtx100lv2000ly2000,vrtx200lv2000ly2000};
hlin200lv2000ly2000 = newl; Line ( hlin200lv2000ly2000 ) = {vrtx200lv2000ly2000,vrtx300lv2000ly2000};
hlin300lv2000ly2000 = newl; Line ( hlin300lv2000ly2000 ) = {vrtx300lv2000ly2000,vrtx400lv2000ly2000};
hlin400lv2000ly2000 = newl; Line ( hlin400lv2000ly2000 ) = {vrtx400lv2000ly2000,vrtx100lv2000ly2000};


// level 000 -> 300
vlin100lv000ly2000 = newl; Line ( vlin100lv000ly2000 ) = {vrtx100lv000ly2000,vrtx100lv300ly2000};
vlin200lv000ly2000 = newl; Line ( vlin200lv000ly2000 ) = {vrtx200lv000ly2000,vrtx200lv300ly2000};
vlin300lv000ly2000 = newl; Line ( vlin300lv000ly2000 ) = {vrtx300lv000ly2000,vrtx300lv300ly2000};
vlin400lv000ly2000 = newl; Line ( vlin400lv000ly2000 ) = {vrtx400lv000ly2000,vrtx400lv300ly2000};

// level 300 -> 450
vlin100lv300ly2000 = newl; Line ( vlin100lv300ly2000 ) = {vrtx100lv300ly2000,vrtx100lv450ly2000};
vlin200lv300ly2000 = newl; Line ( vlin200lv300ly2000 ) = {vrtx200lv300ly2000,vrtx200lv450ly2000};
vlin300lv300ly2000 = newl; Line ( vlin300lv300ly2000 ) = {vrtx300lv300ly2000,vrtx300lv450ly2000};
vlin400lv300ly2000 = newl; Line ( vlin400lv300ly2000 ) = {vrtx400lv300ly2000,vrtx400lv450ly2000};

// level 300 -> 450 - layer 2001
vlin100lv300ly2001 = newl; Line ( vlin100lv300ly2001 ) = {vrtx100lv300ly2001, vrtx100lv450ly2001};
vlin200lv300ly2001 = newl; Line ( vlin200lv300ly2001 ) = {vrtx200lv300ly2001, vrtx200lv450ly2001};
vlin300lv300ly2001 = newl; Line ( vlin300lv300ly2001 ) = {vrtx300lv300ly2001, vrtx300lv450ly2001};
vlin400lv300ly2001 = newl; Line ( vlin400lv300ly2001 ) = {vrtx400lv300ly2001, vrtx400lv450ly2001};

// level 450 -> 1000
vlin100lv450ly2000 = newl; Line ( vlin100lv450ly2000 ) = {vrtx100lv450ly2000,vrtx100lv1000ly2000};
vlin200lv450ly2000 = newl; Line ( vlin200lv450ly2000 ) = {vrtx200lv450ly2000,vrtx200lv1000ly2000};
vlin300lv450ly2000 = newl; Line ( vlin300lv450ly2000 ) = {vrtx300lv450ly2000,vrtx300lv1000ly2000};
vlin400lv450ly2000 = newl; Line ( vlin400lv450ly2000 ) = {vrtx400lv450ly2000,vrtx400lv1000ly2000};

// level 1000 -> 2000
vlin100lv1000ly2000 = newl; Line ( vlin100lv1000ly2000 ) = {vrtx100lv1000ly2000,vrtx100lv2000ly2000};
vlin200lv1000ly2000 = newl; Line ( vlin200lv1000ly2000 ) = {vrtx200lv1000ly2000,vrtx200lv2000ly2000};
vlin300lv1000ly2000 = newl; Line ( vlin300lv1000ly2000 ) = {vrtx300lv1000ly2000,vrtx300lv2000ly2000};
vlin400lv1000ly2000 = newl; Line ( vlin400lv1000ly2000 ) = {vrtx400lv1000ly2000,vrtx400lv2000ly2000};







/*************************************************************************/
/** Surfaces - cavity                                                   **/
/*************************************************************************/

// Horizontal
hlinloplv300ly50 = newreg; Line Loop(hlinloplv300ly50)  = {hlin100lv300ly50, hlin200lv300ly50, hlin300lv300ly50, hlin400lv300ly50, hlin500lv300ly50, hlin600lv300ly50, hlin700lv300ly50, hlin800lv300ly50};
hsurflv300ly50 = news; Plane Surface(hsurflv300ly50) = {hlinloplv300ly50};

// Vertical Level 300 -> 400
vlinlop100lv300ly50 = newreg; Line Loop(vlinlop100lv300ly50)  = { hlin100lv300ly50, vlin200lv300ly50, -hlin100lv400ly50, -vlin100lv300ly50 };
vlinlop200lv300ly50 = newreg; Line Loop(vlinlop200lv300ly50)  = { hlin200lv300ly50, vlin300lv300ly50, -hlin200lv400ly50, -vlin200lv300ly50 };
vlinlop300lv300ly50 = newreg; Line Loop(vlinlop300lv300ly50)  = { hlin300lv300ly50, vlin400lv300ly50, -hlin300lv400ly50, -vlin300lv300ly50 };
vlinlop400lv300ly50 = newreg; Line Loop(vlinlop400lv300ly50)  = { hlin400lv300ly50, vlin500lv300ly50, -hlin400lv400ly50, -vlin400lv300ly50 };
vlinlop500lv300ly50 = newreg; Line Loop(vlinlop500lv300ly50)  = { hlin500lv300ly50, vlin600lv300ly50, -hlin500lv400ly50, -vlin500lv300ly50 };
vlinlop600lv300ly50 = newreg; Line Loop(vlinlop600lv300ly50)  = { hlin600lv300ly50, vlin700lv300ly50, -hlin600lv400ly50, -vlin600lv300ly50 };
vlinlop700lv300ly50 = newreg; Line Loop(vlinlop700lv300ly50)  = { hlin700lv300ly50, vlin800lv300ly50, -hlin700lv400ly50, -vlin700lv300ly50 };
vlinlop800lv300ly50 = newreg; Line Loop(vlinlop800lv300ly50)  = { hlin800lv300ly50, vlin100lv300ly50, -hlin800lv400ly50, -vlin800lv300ly50 };
vsurf100lv300ly50 = news; Ruled Surface(vsurf100lv300ly50) = {vlinlop100lv300ly50};
vsurf300lv300ly50 = news; Ruled Surface(vsurf300lv300ly50) = {vlinlop300lv300ly50};
vsurf500lv300ly50 = news; Ruled Surface(vsurf500lv300ly50) = {vlinlop500lv300ly50};
vsurf700lv300ly50 = news; Ruled Surface(vsurf700lv300ly50) = {vlinlop700lv300ly50};
vsurf200lv300ly50 = news; Plane Surface(vsurf200lv300ly50) = {vlinlop200lv300ly50};
vsurf400lv300ly50 = news; Plane Surface(vsurf400lv300ly50) = {vlinlop400lv300ly50};
vsurf600lv300ly50 = news; Plane Surface(vsurf600lv300ly50) = {vlinlop600lv300ly50};
vsurf800lv300ly50 = news; Plane Surface(vsurf800lv300ly50) = {vlinlop800lv300ly50};

// Curved
vlinlop100lv400ly50 = newreg; Line Loop(vlinlop100lv400ly50)  = { hlin100lv400ly50, vlincurv200lv400ly50, -hlin100lv450ly50, -vlincurv100lv400ly50 };
vlinlop200lv400ly50 = newreg; Line Loop(vlinlop200lv400ly50)  = { hlin200lv400ly50, vlincurv300lv400ly50, -hlin200lv450ly50, -vlincurv200lv400ly50 };
vlinlop300lv400ly50 = newreg; Line Loop(vlinlop300lv400ly50)  = { hlin300lv400ly50, vlincurv400lv400ly50, -hlin300lv450ly50, -vlincurv300lv400ly50 };
vlinlop400lv400ly50 = newreg; Line Loop(vlinlop400lv400ly50)  = { hlin400lv400ly50, vlincurv500lv400ly50, -hlin400lv450ly50, -vlincurv400lv400ly50 };
vlinlop500lv400ly50 = newreg; Line Loop(vlinlop500lv400ly50)  = { hlin500lv400ly50, vlincurv600lv400ly50, -hlin500lv450ly50, -vlincurv500lv400ly50 };
vlinlop600lv400ly50 = newreg; Line Loop(vlinlop600lv400ly50)  = { hlin600lv400ly50, vlincurv700lv400ly50, -hlin600lv450ly50, -vlincurv600lv400ly50 };
vlinlop700lv400ly50 = newreg; Line Loop(vlinlop700lv400ly50)  = { hlin700lv400ly50, vlincurv800lv400ly50, -hlin700lv450ly50, -vlincurv700lv400ly50 };
vlinlop800lv400ly50 = newreg; Line Loop(vlinlop800lv400ly50)  = { hlin800lv400ly50, vlincurv100lv400ly50, -hlin800lv450ly50, -vlincurv800lv400ly50 };
vsurf100lv400ly50 = news; Ruled Surface(vsurf100lv400ly50) = {vlinlop100lv400ly50};
vsurf200lv400ly50 = news; Ruled Surface(vsurf200lv400ly50) = {vlinlop200lv400ly50};
vsurf300lv400ly50 = news; Ruled Surface(vsurf300lv400ly50) = {vlinlop300lv400ly50};
vsurf400lv400ly50 = news; Ruled Surface(vsurf400lv400ly50) = {vlinlop400lv400ly50};
vsurf500lv400ly50 = news; Ruled Surface(vsurf500lv400ly50) = {vlinlop500lv400ly50};
vsurf600lv400ly50 = news; Ruled Surface(vsurf600lv400ly50) = {vlinlop600lv400ly50};
vsurf700lv400ly50 = news; Ruled Surface(vsurf700lv400ly50) = {vlinlop700lv400ly50};
vsurf800lv400ly50 = news; Ruled Surface(vsurf800lv400ly50) = {vlinlop800lv400ly50};


/*************************************************************************/
/** Surfaces - substrate                                                **/
/*************************************************************************/

// Horizontal - Level 300
diaghlinlop100lv300ly2000 = newreg; Line Loop(diaghlinlop100lv300ly2000)  = { hlin100lv300ly2001, diaghlin200lv300ly2000, -hlin200lv300ly50, -diaghlin100lv300ly2000 };   // surf 100
diaghlinlop200lv300ly2000 = newreg; Line Loop(diaghlinlop200lv300ly2000)  = { hlin200lv300ly2001, diaghlin300lv300ly2000, -hlin300lv300ly50, -diaghlin200lv300ly2000 };   // surf 200
diaghlinlop300lv300ly2000 = newreg; Line Loop(diaghlinlop300lv300ly2000)  = { hlin300lv300ly2001, diaghlin400lv300ly2000, -hlin400lv300ly50, -diaghlin300lv300ly2000 };   // surf 300
diaghlinlop400lv300ly2000 = newreg; Line Loop(diaghlinlop400lv300ly2000)  = { hlin400lv300ly2001, diaghlin500lv300ly2000, -hlin500lv300ly50, -diaghlin400lv300ly2000 };   // surf 400
diaghlinlop500lv300ly2000 = newreg; Line Loop(diaghlinlop500lv300ly2000)  = { hlin500lv300ly2001, diaghlin600lv300ly2000, -hlin600lv300ly50, -diaghlin500lv300ly2000 };   // surf 500
diaghlinlop600lv300ly2000 = newreg; Line Loop(diaghlinlop600lv300ly2000)  = { hlin600lv300ly2001, diaghlin700lv300ly2000, -hlin700lv300ly50, -diaghlin600lv300ly2000 };   // surf 600
diaghlinlop700lv300ly2000 = newreg; Line Loop(diaghlinlop700lv300ly2000)  = { hlin700lv300ly2001, diaghlin800lv300ly2000, -hlin800lv300ly50, -diaghlin700lv300ly2000 };   // surf 700
diaghlinlop800lv300ly2000 = newreg; Line Loop(diaghlinlop800lv300ly2000)  = { hlin800lv300ly2001, diaghlin100lv300ly2000, -hlin100lv300ly50, -diaghlin800lv300ly2000 };   // surf 800
diaghsurf100lv300ly2000 = news; Plane Surface(diaghsurf100lv300ly2000) = {diaghlinlop100lv300ly2000};
diaghsurf300lv300ly2000 = news; Plane Surface(diaghsurf300lv300ly2000) = {diaghlinlop300lv300ly2000};
diaghsurf500lv300ly2000 = news; Plane Surface(diaghsurf500lv300ly2000) = {diaghlinlop500lv300ly2000};
diaghsurf700lv300ly2000 = news; Plane Surface(diaghsurf700lv300ly2000) = {diaghlinlop700lv300ly2000};
diaghsurf200lv300ly2000 = news; Ruled Surface(diaghsurf200lv300ly2000) = {diaghlinlop200lv300ly2000};
diaghsurf400lv300ly2000 = news; Ruled Surface(diaghsurf400lv300ly2000) = {diaghlinlop400lv300ly2000};
diaghsurf600lv300ly2000 = news; Ruled Surface(diaghsurf600lv300ly2000) = {diaghlinlop600lv300ly2000};
diaghsurf800lv300ly2000 = news; Ruled Surface(diaghsurf800lv300ly2000) = {diaghlinlop800lv300ly2000};

/*************************************************************************/
/** Surfaces - metal layer                                              **/
/*************************************************************************/

// Horizontal - Level 450
diaghlinlop100lv450ly2000 = newreg; Line Loop(diaghlinlop100lv450ly2000)  = { hlin100lv450ly2001, diaghlin200lv450ly2000, -hlin200lv450ly50, -diaghlin100lv450ly2000 };   // surf 100
diaghlinlop200lv450ly2000 = newreg; Line Loop(diaghlinlop200lv450ly2000)  = { hlin200lv450ly2001, diaghlin300lv450ly2000, -hlin300lv450ly50, -diaghlin200lv450ly2000 };   // surf 200
diaghlinlop300lv450ly2000 = newreg; Line Loop(diaghlinlop300lv450ly2000)  = { hlin300lv450ly2001, diaghlin400lv450ly2000, -hlin400lv450ly50, -diaghlin300lv450ly2000 };   // surf 300
diaghlinlop400lv450ly2000 = newreg; Line Loop(diaghlinlop400lv450ly2000)  = { hlin400lv450ly2001, diaghlin500lv450ly2000, -hlin500lv450ly50, -diaghlin400lv450ly2000 };   // surf 400
diaghlinlop500lv450ly2000 = newreg; Line Loop(diaghlinlop500lv450ly2000)  = { hlin500lv450ly2001, diaghlin600lv450ly2000, -hlin600lv450ly50, -diaghlin500lv450ly2000 };   // surf 500
diaghlinlop600lv450ly2000 = newreg; Line Loop(diaghlinlop600lv450ly2000)  = { hlin600lv450ly2001, diaghlin700lv450ly2000, -hlin700lv450ly50, -diaghlin600lv450ly2000 };   // surf 600
diaghlinlop700lv450ly2000 = newreg; Line Loop(diaghlinlop700lv450ly2000)  = { hlin700lv450ly2001, diaghlin800lv450ly2000, -hlin800lv450ly50, -diaghlin700lv450ly2000 };   // surf 700
diaghlinlop800lv450ly2000 = newreg; Line Loop(diaghlinlop800lv450ly2000)  = { hlin800lv450ly2001, diaghlin100lv450ly2000, -hlin100lv450ly50, -diaghlin800lv450ly2000 };   // surf 800
diaghsurf100lv450ly2000 = news; Plane Surface(diaghsurf100lv450ly2000) = {diaghlinlop100lv450ly2000};
diaghsurf300lv450ly2000 = news; Plane Surface(diaghsurf300lv450ly2000) = {diaghlinlop300lv450ly2000};
diaghsurf500lv450ly2000 = news; Plane Surface(diaghsurf500lv450ly2000) = {diaghlinlop500lv450ly2000};
diaghsurf700lv450ly2000 = news; Plane Surface(diaghsurf700lv450ly2000) = {diaghlinlop700lv450ly2000};
diaghsurf200lv450ly2000 = news; Ruled Surface(diaghsurf200lv450ly2000) = {diaghlinlop200lv450ly2000};
diaghsurf400lv450ly2000 = news; Ruled Surface(diaghsurf400lv450ly2000) = {diaghlinlop400lv450ly2000};
diaghsurf600lv450ly2000 = news; Ruled Surface(diaghsurf600lv450ly2000) = {diaghlinlop600lv450ly2000};
diaghsurf800lv450ly2000 = news; Ruled Surface(diaghsurf800lv450ly2000) = {diaghlinlop800lv450ly2000};


/*************************************************************************/
/** Surfaces - liquid                                                   **/
/*************************************************************************/

// Horizontal
hlinlop100lv1000ly2000  = newreg; Line Loop(hlinlop100lv1000ly2000)   = { hlin100lv1000ly2000,  hlin200lv1000ly2000,  hlin300lv1000ly2000,  hlin400lv1000ly2000 };   // level 000
hsurf100lv1000ly2000  = news; Plane Surface(hsurf100lv1000ly2000)  = {hlinlop100lv1000ly2000};

/*************************************************************************/
/** Surfaces - computational domain                                     **/
/*************************************************************************/

// Horizontal
hlinlop100lv000ly2000  = newreg; Line Loop(hlinlop100lv000ly2000)   = { hlin100lv000ly2000,  hlin200lv000ly2000,  hlin300lv000ly2000,  hlin400lv000ly2000 };   // level 000
hlinlop100lv2000ly2000 = newreg; Line Loop(hlinlop100lv2000ly2000)  = { hlin100lv2000ly2000, hlin200lv2000ly2000, hlin300lv2000ly2000, hlin400lv2000ly2000 };  // level 2000
hsurf100lv000ly2000  = news; Plane Surface(hsurf100lv000ly2000)  = {hlinlop100lv000ly2000};
hsurf100lv2000ly2000 = news; Plane Surface(hsurf100lv2000ly2000) = {hlinlop100lv2000ly2000};

// Vertical Level 000 -> 300
vlinlop100lv000ly2000 = newreg; Line Loop(vlinlop100lv000ly2000)  = {  hlin100lv000ly2000, vlin200lv000ly2000, -hlin100lv300ly2000, -vlin100lv000ly2000 };  // surf 100
vlinlop200lv000ly2000 = newreg; Line Loop(vlinlop200lv000ly2000)  = {  hlin200lv000ly2000, vlin300lv000ly2000, -hlin200lv300ly2000, -vlin200lv000ly2000 };  // surf 200
vlinlop300lv000ly2000 = newreg; Line Loop(vlinlop300lv000ly2000)  = {  hlin300lv000ly2000, vlin400lv000ly2000, -hlin300lv300ly2000, -vlin300lv000ly2000 };  // surf 300
vlinlop400lv000ly2000 = newreg; Line Loop(vlinlop400lv000ly2000)  = {  hlin400lv000ly2000, vlin100lv000ly2000, -hlin400lv300ly2000, -vlin400lv000ly2000 };  // surf 400
vsurf100lv000ly2000 = news; Plane Surface(vsurf100lv000ly2000) = {vlinlop100lv000ly2000};
vsurf200lv000ly2000 = news; Plane Surface(vsurf200lv000ly2000) = {vlinlop200lv000ly2000};
vsurf300lv000ly2000 = news; Plane Surface(vsurf300lv000ly2000) = {vlinlop300lv000ly2000};
vsurf400lv000ly2000 = news; Plane Surface(vsurf400lv000ly2000) = {vlinlop400lv000ly2000};

// Vertical Level 300 -> 450
vlinlop100lv300ly2000 = newreg; Line Loop(vlinlop100lv300ly2000)  = {  hlin100lv300ly2001, vlin100lv300ly2001, -hlin100lv450ly2001, -vlin100lv300ly2000 };  // surf 100
vlinlop200lv300ly2000 = newreg; Line Loop(vlinlop200lv300ly2000)  = {  hlin200lv300ly2001, vlin200lv300ly2000, -hlin200lv450ly2001, -vlin100lv300ly2001 };  // surf 100
vlinlop300lv300ly2000 = newreg; Line Loop(vlinlop300lv300ly2000)  = {  hlin300lv300ly2001, vlin200lv300ly2001, -hlin300lv450ly2001, -vlin200lv300ly2000 };  // surf 200
vlinlop400lv300ly2000 = newreg; Line Loop(vlinlop400lv300ly2000)  = {  hlin400lv300ly2001, vlin300lv300ly2000, -hlin400lv450ly2001, -vlin200lv300ly2001 };  // surf 200
vlinlop500lv300ly2000 = newreg; Line Loop(vlinlop500lv300ly2000)  = {  hlin500lv300ly2001, vlin300lv300ly2001, -hlin500lv450ly2001, -vlin300lv300ly2000 };  // surf 300
vlinlop600lv300ly2000 = newreg; Line Loop(vlinlop600lv300ly2000)  = {  hlin600lv300ly2001, vlin400lv300ly2000, -hlin600lv450ly2001, -vlin300lv300ly2001 };  // surf 300
vlinlop700lv300ly2000 = newreg; Line Loop(vlinlop700lv300ly2000)  = {  hlin700lv300ly2001, vlin400lv300ly2001, -hlin700lv450ly2001, -vlin400lv300ly2000 };  // surf 400
vlinlop800lv300ly2000 = newreg; Line Loop(vlinlop800lv300ly2000)  = {  hlin800lv300ly2001, vlin100lv300ly2000, -hlin800lv450ly2001, -vlin400lv300ly2001 };  // surf 400
vsurf100lv300ly2000 = news; Plane Surface(vsurf100lv300ly2000) = {vlinlop100lv300ly2000};
vsurf200lv300ly2000 = news; Plane Surface(vsurf200lv300ly2000) = {vlinlop200lv300ly2000};
vsurf300lv300ly2000 = news; Plane Surface(vsurf300lv300ly2000) = {vlinlop300lv300ly2000};
vsurf400lv300ly2000 = news; Plane Surface(vsurf400lv300ly2000) = {vlinlop400lv300ly2000};
vsurf500lv300ly2000 = news; Plane Surface(vsurf500lv300ly2000) = {vlinlop500lv300ly2000};
vsurf600lv300ly2000 = news; Plane Surface(vsurf600lv300ly2000) = {vlinlop600lv300ly2000};
vsurf700lv300ly2000 = news; Plane Surface(vsurf700lv300ly2000) = {vlinlop700lv300ly2000};
vsurf800lv300ly2000 = news; Plane Surface(vsurf800lv300ly2000) = {vlinlop800lv300ly2000};

// Vertical Level 450 -> 1000
vlinlop100lv450ly2000 = newreg; Line Loop(vlinlop100lv450ly2000)  = {  hlin100lv450ly2000, vlin200lv450ly2000, -hlin100lv1000ly2000, -vlin100lv450ly2000 };  // surf 100
vlinlop200lv450ly2000 = newreg; Line Loop(vlinlop200lv450ly2000)  = {  hlin200lv450ly2000, vlin300lv450ly2000, -hlin200lv1000ly2000, -vlin200lv450ly2000 };  // surf 200
vlinlop300lv450ly2000 = newreg; Line Loop(vlinlop300lv450ly2000)  = {  hlin300lv450ly2000, vlin400lv450ly2000, -hlin300lv1000ly2000, -vlin300lv450ly2000 };  // surf 300
vlinlop400lv450ly2000 = newreg; Line Loop(vlinlop400lv450ly2000)  = {  hlin400lv450ly2000, vlin100lv450ly2000, -hlin400lv1000ly2000, -vlin400lv450ly2000 };  // surf 400
vsurf100lv450ly2000 = news; Plane Surface(vsurf100lv450ly2000) = {vlinlop100lv450ly2000};
vsurf200lv450ly2000 = news; Plane Surface(vsurf200lv450ly2000) = {vlinlop200lv450ly2000};
vsurf300lv450ly2000 = news; Plane Surface(vsurf300lv450ly2000) = {vlinlop300lv450ly2000};
vsurf400lv450ly2000 = news; Plane Surface(vsurf400lv450ly2000) = {vlinlop400lv450ly2000};

// Vertical Level 1000 -> 2000
vlinlop100lv1000ly2000 = newreg; Line Loop(vlinlop100lv1000ly2000)  = {  hlin100lv1000ly2000, vlin200lv1000ly2000, -hlin100lv2000ly2000, -vlin100lv1000ly2000 };  // surf 100
vlinlop200lv1000ly2000 = newreg; Line Loop(vlinlop200lv1000ly2000)  = {  hlin200lv1000ly2000, vlin300lv1000ly2000, -hlin200lv2000ly2000, -vlin200lv1000ly2000 };  // surf 200
vlinlop300lv1000ly2000 = newreg; Line Loop(vlinlop300lv1000ly2000)  = {  hlin300lv1000ly2000, vlin400lv1000ly2000, -hlin300lv2000ly2000, -vlin300lv1000ly2000 };  // surf 300
vlinlop400lv1000ly2000 = newreg; Line Loop(vlinlop400lv1000ly2000)  = {  hlin400lv1000ly2000, vlin100lv1000ly2000, -hlin400lv2000ly2000, -vlin400lv1000ly2000 };  // surf 400
vsurf100lv1000ly2000 = news; Plane Surface(vsurf100lv1000ly2000) = {vlinlop100lv1000ly2000};
vsurf200lv1000ly2000 = news; Plane Surface(vsurf200lv1000ly2000) = {vlinlop200lv1000ly2000};
vsurf300lv1000ly2000 = news; Plane Surface(vsurf300lv1000ly2000) = {vlinlop300lv1000ly2000};
vsurf400lv1000ly2000 = news; Plane Surface(vsurf400lv1000ly2000) = {vlinlop400lv1000ly2000};





/*****************************************************************************/
/*****************************************************************************/
/** Volume                                                                  **/
/*****************************************************************************/
/*****************************************************************************/


/*****************************************************************************/
/** volume - substrate                                                      **/
/*****************************************************************************/

surflop000lv000ly2000 = newreg; Surface Loop(surflop000lv000ly2000) =    {hsurf100lv000ly2000,

                                                                          vsurf100lv000ly2000,
                                                                          vsurf200lv000ly2000,
                                                                          vsurf300lv000ly2000,
                                                                          vsurf400lv000ly2000,
                                                                          
                                                                          diaghsurf100lv300ly2000,
                                                                          diaghsurf200lv300ly2000,
                                                                          diaghsurf300lv300ly2000,
                                                                          diaghsurf400lv300ly2000,
									  diaghsurf500lv300ly2000,
									  diaghsurf600lv300ly2000,
									  diaghsurf700lv300ly2000,
									  diaghsurf800lv300ly2000,
                                                                          
                                                                          hsurflv300ly50 };         /** lower horizontal surface of the pillar **/
                                                                       
volume000lv000ly2000 = newv; Volume(volume000lv000ly2000) = {surflop000lv000ly2000};

Physical Volume (physvolsubstrate)      = {volume000lv000ly2000};



/*****************************************************************************/
/** volume - metal layer                                                    **/
/*****************************************************************************/

surflop000lv300ly2000 = newreg; Surface Loop(surflop000lv300ly2000) =    {diaghsurf100lv300ly2000,  /** Bottom surfaces of the metal layer  */
                                                                          diaghsurf200lv300ly2000,
                                                                          diaghsurf300lv300ly2000,
                                                                          diaghsurf400lv300ly2000,
									  diaghsurf500lv300ly2000,
									  diaghsurf600lv300ly2000,
									  diaghsurf700lv300ly2000,
									  diaghsurf800lv300ly2000,

									  diaghsurf100lv450ly2000,  /** Top surfaces of the metal layer  */
                                                                          diaghsurf200lv450ly2000,
                                                                          diaghsurf300lv450ly2000,
                                                                          diaghsurf400lv450ly2000,
									  diaghsurf500lv450ly2000,
									  diaghsurf600lv450ly2000,
									  diaghsurf700lv450ly2000,
									  diaghsurf800lv450ly2000,

                                                                          vsurf100lv300ly50,  /** Side walls of cavity  */
                                                                          vsurf200lv300ly50,
                                                                          vsurf300lv300ly50,
                                                                          vsurf400lv300ly50,
									  vsurf500lv300ly50,
									  vsurf600lv300ly50,
									  vsurf700lv300ly50,
									  vsurf800lv300ly50,

									  vsurf100lv400ly50,  /** Top edge curve of the cavity  */
                                                                          vsurf200lv400ly50,
                                                                          vsurf300lv400ly50,
                                                                          vsurf400lv400ly50,
									  vsurf500lv400ly50,
									  vsurf600lv400ly50,
									  vsurf700lv400ly50,
									  vsurf800lv400ly50,

                                                                          vsurf100lv300ly2000,  /** Side walls of domain boundary  */
                                                                          vsurf200lv300ly2000,
                                                                          vsurf300lv300ly2000,
                                                                          vsurf400lv300ly2000,
									  vsurf500lv300ly2000,
									  vsurf600lv300ly2000,
									  vsurf700lv300ly2000,
									  vsurf800lv300ly2000};
                                                                       
volume000lv300ly2000 = newv; Volume(volume000lv300ly2000) = {surflop000lv300ly2000};

Physical Volume (physvollayer1)      = {volume000lv300ly2000};




/*****************************************************************************/
/** volume - liquid fill                                                    **/
/*****************************************************************************/

surflop100lv300ly2000 = newreg; Surface Loop(surflop100lv300ly2000) =    {
									  /* diagonal surfaces of the top layer */
                                                                          diaghsurf100lv450ly2000,
                                                                          diaghsurf200lv450ly2000,
                                                                          diaghsurf300lv450ly2000,
                                                                          diaghsurf400lv450ly2000,
                                                                          diaghsurf500lv450ly2000,
									  diaghsurf600lv450ly2000,
									  diaghsurf700lv450ly2000,
									  diaghsurf800lv450ly2000,
									  
                                                                          /* vertical surfaces of the liquid **/                                                                    
                                                                          vsurf100lv450ly2000,
                                                                          vsurf200lv450ly2000,
                                                                          vsurf300lv450ly2000,
                                                                          vsurf400lv450ly2000,
    
                                                                          /** vertical surfaces of the cavity **/                                                                          
                                                                          vsurf100lv300ly50,
                                                                          vsurf200lv300ly50,
                                                                          vsurf300lv300ly50,
                                                                          vsurf400lv300ly50,
									  vsurf500lv300ly50,
									  vsurf600lv300ly50,
									  vsurf700lv300ly50,
									  vsurf800lv300ly50,
									  
									  vsurf100lv400ly50,  /** Top edge curve of the cavity  */
                                                                          vsurf200lv400ly50,
                                                                          vsurf300lv400ly50,
                                                                          vsurf400lv400ly50,
									  vsurf500lv400ly50,
									  vsurf600lv400ly50,
									  vsurf700lv400ly50,
									  vsurf800lv400ly50,
 
                                                                          /* bottom of the cavity */
                                                                          hsurflv300ly50,
 
                                                                          /* top of the liquid */
                                                                          hsurf100lv1000ly2000
                                                                       };


volume100lv300ly2000 = newv; Volume(volume100lv300ly2000) = {surflop100lv300ly2000};

Physical Volume (physvolliquid)      = {volume100lv300ly2000};





/*****************************************************************************/
/** volume - air everywhere, no substrate                                   **/
/*****************************************************************************/


surflop100lv1000ly2000 = newreg; Surface Loop(surflop100lv1000ly2000) =    {
									  /* top of the liquid */
                                                                          hsurf100lv1000ly2000,
                                                                          
                                                                          /* vertical surfaces of the computational domain **/                                                                    
                                                                          vsurf100lv1000ly2000,
                                                                          vsurf200lv1000ly2000,
                                                                          vsurf300lv1000ly2000,
                                                                          vsurf400lv1000ly2000,
 
                                                                          /* computational domain, upper horizontal */
                                                                          hsurf100lv2000ly2000
                                                                       };


volume100lv1000ly2000 = newv; Volume(volume100lv1000ly2000) = {surflop100lv1000ly2000};

Physical Volume (physvolair)      = {volume100lv1000ly2000};




/*****************************************************************************/
/** boundary - enclosing computational domain                               **/
/*****************************************************************************/


Physical Surface (physsurfairboundary)      = {hsurf100lv000ly2000,
                                       hsurf100lv2000ly2000,
                                       
                                       vsurf100lv000ly2000,
                                       vsurf200lv000ly2000,
                                       vsurf300lv000ly2000,
                                       vsurf400lv000ly2000,
                                       
                                       vsurf100lv300ly2000,
                                       vsurf200lv300ly2000,
                                       vsurf300lv300ly2000,
                                       vsurf400lv300ly2000,
				       vsurf500lv300ly2000,
				       vsurf600lv300ly2000,
				       vsurf700lv300ly2000,
				       vsurf800lv300ly2000,
                                       
                                       vsurf100lv450ly2000,
                                       vsurf200lv450ly2000,
                                       vsurf300lv450ly2000,
                                       vsurf400lv450ly2000,
				       
				       vsurf100lv1000ly2000,
                                       vsurf200lv1000ly2000,
                                       vsurf300lv1000ly2000,
                                       vsurf400lv1000ly2000

                                      };

-------------- next part --------------
Info    : -------------------------------------------------------
Info    : Gmsh version   : 2.10.1-svn
Info    : Build OS       : Linux
Info    : Build options  : 64Bit Ann Bamg Bfgs Blas(Generic) Blossom Chaco DIntegration Dlopen Fltk Gmm Jpeg(Fltk) Kbipack Lapack(Generic) LinuxJoystick MathEx Med Mesh Metis Mmg3d Mpeg Netgen ONELAB ONELABMetamodel OpenCascade OpenGL OptHom PETSc Parser Plugins Png(Fltk) Post SLEPc Salome Solver Tetgen(1.5) Voro3D Zlib
Info    : Build date     : 20150730
Info    : Build host     : debian
Info    : Packager       : geuzaine
Info    : Executable     : /opt/extlib/gmsh/2.10.1/gcc/4.8.4/bin/gmsh
Info    : Home directory : /home/fomins/
Info    : Launch date    : Wed Nov 11 10:44:13 2015
Info    : Command line   : gmsh cavity-rectangular-rev-1400.geo
Info    : -------------------------------------------------------
Info    : Reading 'cavity-rectangular-rev-1400.geo'...
Info    : Done reading 'cavity-rectangular-rev-1400.geo'
Info    : Meshing 1D...
Info    : Meshing curve 1 (Ellipse)
Info    : Meshing curve 2 (Ellipse)
Info    : Meshing curve 3 (Ellipse)
Info    : Meshing curve 4 (Ellipse)
Info    : Meshing curve 5 (Line)
Info    : Meshing curve 6 (Line)
Info    : Meshing curve 7 (Line)
Info    : Meshing curve 8 (Line)
Info    : Meshing curve 9 (Ellipse)
Info    : Meshing curve 10 (Ellipse)
Info    : Meshing curve 11 (Ellipse)
Info    : Meshing curve 12 (Ellipse)
Info    : Meshing curve 13 (Line)
Info    : Meshing curve 14 (Line)
Info    : Meshing curve 15 (Line)
Info    : Meshing curve 16 (Line)
Info    : Meshing curve 17 (Ellipse)
Info    : Meshing curve 18 (Ellipse)
Info    : Meshing curve 19 (Ellipse)
Info    : Meshing curve 20 (Ellipse)
Info    : Meshing curve 21 (Line)
Info    : Meshing curve 22 (Line)
Info    : Meshing curve 23 (Line)
Info    : Meshing curve 24 (Line)
Info    : Meshing curve 25 (Line)
Info    : Meshing curve 26 (Line)
Info    : Meshing curve 27 (Line)
Info    : Meshing curve 28 (Line)
Info    : Meshing curve 29 (Line)
Info    : Meshing curve 30 (Line)
Info    : Meshing curve 31 (Line)
Info    : Meshing curve 32 (Line)
Info    : Meshing curve 33 (Ellipse)
Info    : Meshing curve 34 (Ellipse)
Info    : Meshing curve 35 (Ellipse)
Info    : Meshing curve 36 (Ellipse)
Info    : Meshing curve 37 (Ellipse)
Info    : Meshing curve 38 (Ellipse)
Info    : Meshing curve 39 (Ellipse)
Info    : Meshing curve 40 (Ellipse)
Info    : Meshing curve 41 (Line)
Info    : Meshing curve 42 (Line)
Info    : Meshing curve 43 (Line)
Info    : Meshing curve 44 (Line)
Info    : Meshing curve 45 (Line)
Info    : Meshing curve 46 (Line)
Info    : Meshing curve 47 (Line)
Info    : Meshing curve 48 (Line)
Info    : Meshing curve 49 (Line)
Info    : Meshing curve 50 (Line)
Info    : Meshing curve 51 (Line)
Info    : Meshing curve 52 (Line)
Info    : Meshing curve 53 (Line)
Info    : Meshing curve 54 (Line)
Info    : Meshing curve 55 (Line)
Info    : Meshing curve 56 (Line)
Info    : Meshing curve 57 (Line)
Info    : Meshing curve 58 (Line)
Info    : Meshing curve 59 (Line)
Info    : Meshing curve 60 (Line)
Info    : Meshing curve 61 (Line)
Info    : Meshing curve 62 (Line)
Info    : Meshing curve 63 (Line)
Info    : Meshing curve 64 (Line)
Info    : Meshing curve 65 (Line)
Info    : Meshing curve 66 (Line)
Info    : Meshing curve 67 (Line)
Info    : Meshing curve 68 (Line)
Info    : Meshing curve 69 (Line)
Info    : Meshing curve 70 (Line)
Info    : Meshing curve 71 (Line)
Info    : Meshing curve 72 (Line)
Info    : Meshing curve 73 (Line)
Info    : Meshing curve 74 (Line)
Info    : Meshing curve 75 (Line)
Info    : Meshing curve 76 (Line)
Info    : Meshing curve 77 (Line)
Info    : Meshing curve 78 (Line)
Info    : Meshing curve 79 (Line)
Info    : Meshing curve 80 (Line)
Info    : Meshing curve 81 (Line)
Info    : Meshing curve 82 (Line)
Info    : Meshing curve 83 (Line)
Info    : Meshing curve 84 (Line)
Info    : Meshing curve 85 (Line)
Info    : Meshing curve 86 (Line)
Info    : Meshing curve 87 (Line)
Info    : Meshing curve 88 (Line)
Info    : Meshing curve 89 (Line)
Info    : Meshing curve 90 (Line)
Info    : Meshing curve 91 (Line)
Info    : Meshing curve 92 (Line)
Info    : Meshing curve 93 (Line)
Info    : Meshing curve 94 (Line)
Info    : Meshing curve 95 (Line)
Info    : Meshing curve 96 (Line)
Info    : Meshing curve 97 (Line)
Info    : Meshing curve 98 (Line)
Info    : Meshing curve 99 (Line)
Info    : Meshing curve 100 (Line)
Info    : Meshing curve 101 (Line)
Info    : Meshing curve 102 (Line)
Info    : Meshing curve 103 (Line)
Info    : Meshing curve 104 (Line)
Info    : Meshing curve 105 (Line)
Info    : Meshing curve 106 (Line)
Info    : Meshing curve 107 (Line)
Info    : Meshing curve 108 (Line)
Info    : Meshing curve 109 (Line)
Info    : Meshing curve 110 (Line)
Info    : Meshing curve 111 (Line)
Info    : Meshing curve 112 (Line)
Info    : Done meshing 1D (0.152998 s)
Info    : Meshing 2D...
Info    : Meshing surface 114 (Plane, Delaunay)
Info    : Meshing surface 123 (Ruled surface, MeshAdapt)
Info    : Meshing surface 124 (Ruled surface, MeshAdapt)
Info    : Meshing surface 125 (Ruled surface, MeshAdapt)
Info    : Meshing surface 126 (Ruled surface, MeshAdapt)
Info    : Meshing surface 127 (Plane, Delaunay)
Info    : Meshing surface 128 (Plane, Delaunay)
Info    : Meshing surface 129 (Plane, Delaunay)
Info    : Meshing surface 130 (Plane, Delaunay)
Info    : Meshing surface 139 (Ruled surface, MeshAdapt)
Info    : Meshing surface 140 (Ruled surface, MeshAdapt)
Info    : Meshing surface 141 (Ruled surface, MeshAdapt)
Info    : Meshing surface 142 (Ruled surface, MeshAdapt)
Info    : Meshing surface 143 (Ruled surface, MeshAdapt)
Info    : Meshing surface 144 (Ruled surface, MeshAdapt)
Info    : Meshing surface 145 (Ruled surface, MeshAdapt)
Info    : Meshing surface 146 (Ruled surface, MeshAdapt)
Info    : Meshing surface 155 (Plane, Delaunay)
Info    : Meshing surface 156 (Plane, Delaunay)
Info    : Meshing surface 157 (Plane, Delaunay)
Info    : Meshing surface 158 (Plane, Delaunay)
Info    : Meshing surface 159 (Ruled surface, MeshAdapt)
Info    : Meshing surface 160 (Ruled surface, MeshAdapt)
Info    : Meshing surface 161 (Ruled surface, MeshAdapt)
Info    : Meshing surface 162 (Ruled surface, MeshAdapt)
Info    : Meshing surface 171 (Plane, Delaunay)
Info    : Meshing surface 172 (Plane, Delaunay)
Info    : Meshing surface 173 (Plane, Delaunay)
Info    : Meshing surface 174 (Plane, Delaunay)
Info    : Meshing surface 175 (Ruled surface, MeshAdapt)
Info    : Meshing surface 176 (Ruled surface, MeshAdapt)
Info    : Meshing surface 177 (Ruled surface, MeshAdapt)
Info    : Meshing surface 178 (Ruled surface, MeshAdapt)
Info    : Meshing surface 180 (Plane, Delaunay)
Info    : Meshing surface 183 (Plane, Delaunay)
Info    : Meshing surface 184 (Plane, Delaunay)
Info    : Meshing surface 189 (Plane, Delaunay)
Info    : Meshing surface 190 (Plane, Delaunay)
Info    : Meshing surface 191 (Plane, Delaunay)
Info    : Meshing surface 192 (Plane, Delaunay)
Info    : Meshing surface 201 (Plane, Delaunay)
Info    : Meshing surface 202 (Plane, Delaunay)
Info    : Meshing surface 203 (Plane, Delaunay)
Info    : Meshing surface 204 (Plane, Delaunay)
Info    : Meshing surface 205 (Plane, Delaunay)
Info    : Meshing surface 206 (Plane, Delaunay)
Info    : Meshing surface 207 (Plane, Delaunay)
Info    : Meshing surface 208 (Plane, Delaunay)
Info    : Meshing surface 213 (Plane, Delaunay)
Info    : Meshing surface 214 (Plane, Delaunay)
Info    : Meshing surface 215 (Plane, Delaunay)
Info    : Meshing surface 216 (Plane, Delaunay)
Info    : Meshing surface 221 (Plane, Delaunay)
Info    : Meshing surface 222 (Plane, Delaunay)
Info    : Meshing surface 223 (Plane, Delaunay)
Info    : Meshing surface 224 (Plane, Delaunay)
Info    : Done meshing 2D (0.581919 s)
Info    : 3712 vertices 8210 elements
Info    : Meshing 3D...
Info    : Delaunay Meshing 4 volumes with 1 connected components
Info    : Meshing volume 226 (Delaunay)
Info    : Meshing volume 6002 (Delaunay)
Info    : Meshing volume 6004 (Delaunay)
Info    : Meshing volume 6006 (Delaunay)
Info    : 3580 points 21649 edges and 7316 faces in the initial mesh
Info    : Found region 6004
Info    : Found region 6002
Info    : Found region 226
Info    : Found region 6006
Info    : 0 points created - Worst tet radius is -nan (PTS removed 0 0)
Info    : 3D point insertion terminated (3585 points created):
Info    :  - 0 Delaunay cavities modified for star shapeness
Info    :  - 13 points could not be inserted
Info    :  - 17434 tetrahedra created in 0.012529 sec. (1391491 tets/sec.)
Info    : Done meshing 3D (0.41212 s)
Warning : Could not orient normal of surface 114
Info    : Reparameterizing edge 34 on face 123
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.75
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.5625
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.421875
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.316406
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.237305
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.177979
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.133484
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.100113
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.0750847
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.0563135
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.0422351
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.0316764
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.0237573
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.0178179
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.0133635
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.0100226
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.00751695
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.00563771
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.00422828
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.00317121
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.00237841
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.00178381
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.00133786
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.00100339
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.000752543
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.000564408
Info    : point -16.0493 -51.0493 44.4592 : Relaxation factor = 0.000423306
Warning : Converged for i=7 j=5 (err=1.81843e-08 iter=2) BUT xyz error = 12.5473 in point (-1.604927e+01,-5.104927e+01,4.445920e+01) on surface 123
Info    : Reparameterizing edge 84 on face 124
Info    : point -300 -125 50 : Relaxation factor = 0.75
Info    : point -300 -125 50 : Relaxation factor = 0.5625
Info    : point -300 -125 50 : Relaxation factor = 0.421875
Info    : point -300 -125 50 : Relaxation factor = 0.316406
Info    : point -300 -125 50 : Relaxation factor = 0.237305
Info    : point -300 -125 50 : Relaxation factor = 0.177979
Info    : point -300 -125 50 : Relaxation factor = 0.133484
Info    : point -300 -125 50 : Relaxation factor = 0.100113
Info    : point -300 -125 50 : Relaxation factor = 0.0750847
Info    : point -300 -125 50 : Relaxation factor = 0.0563135
Info    : point -300 -125 50 : Relaxation factor = 0.0422351
Info    : point -300 -125 50 : Relaxation factor = 0.0316764
Info    : point -300 -125 50 : Relaxation factor = 0.0237573
Info    : point -300 -125 50 : Relaxation factor = 0.0178179
Info    : point -300 -125 50 : Relaxation factor = 0.0133635
Info    : point -300 -125 50 : Relaxation factor = 0.0100226
Info    : point -300 -125 50 : Relaxation factor = 0.00751695
Info    : point -300 -125 50 : Relaxation factor = 0.00563771
Info    : point -300 -125 50 : Relaxation factor = 0.00422828
Info    : point -300 -125 50 : Relaxation factor = 0.00317121
Info    : point -300 -125 50 : Relaxation factor = 0.00237841
Info    : point -300 -125 50 : Relaxation factor = 0.00178381
Info    : point -300 -125 50 : Relaxation factor = 0.00133786
Info    : point -300 -125 50 : Relaxation factor = 0.00100339
Info    : point -300 -125 50 : Relaxation factor = 0.000752543
Info    : point -300 -125 50 : Relaxation factor = 0.000564408
Info    : point -300 -125 50 : Relaxation factor = 0.000423306
Info    : point -300 -125 50 : Relaxation factor = 0.000317479
Info    : point -300 -125 50 : Relaxation factor = 0.000238109
Info    : point -300 -125 50 : Relaxation factor = 0.000178582
Info    : point -300 -125 50 : Relaxation factor = 0.000133937
Info    : point -300 -125 50 : Relaxation factor = 0.000100452
Info    : point -300 -125 50 : Relaxation factor = 7.53393e-05
Info    : point -300 -125 50 : Relaxation factor = 5.65045e-05
Info    : point -300 -125 50 : Relaxation factor = 4.23784e-05
Info    : point -300 -125 50 : Relaxation factor = 3.17838e-05
Info    : point -300 -125 50 : Relaxation factor = 2.38378e-05
Warning : Converged for i=7 j=0 (err=1.69588e-08 iter=2) BUT xyz error = 337.268 in point (-3.000000e+02,-1.250000e+02,5.000000e+01) on surface 124
Info    : Reparameterizing edge 31 on face 124
Info    : point -15 50 34.2921 : Relaxation factor = 0.75
Info    : point -15 50 34.2921 : Relaxation factor = 0.5625
Info    : point -15 50 34.2921 : Relaxation factor = 0.421875
Info    : point -15 50 34.2921 : Relaxation factor = 0.316406
Info    : point -15 50 34.2921 : Relaxation factor = 0.237305
Info    : point -15 50 34.2921 : Relaxation factor = 0.177979
Info    : point -15 50 34.2921 : Relaxation factor = 0.133484
Info    : point -15 50 34.2921 : Relaxation factor = 0.100113
Info    : point -15 50 34.2921 : Relaxation factor = 0.0750847
Info    : point -15 50 34.2921 : Relaxation factor = 0.0563135
Info    : point -15 50 34.2921 : Relaxation factor = 0.0422351
Info    : point -15 50 34.2921 : Relaxation factor = 0.0316764
Info    : point -15 50 34.2921 : Relaxation factor = 0.0237573
Info    : point -15 50 34.2921 : Relaxation factor = 0.0178179
Info    : point -15 50 34.2921 : Relaxation factor = 0.0133635
Info    : point -15 50 34.2921 : Relaxation factor = 0.0100226
Info    : point -15 50 34.2921 : Relaxation factor = 0.00751695
Info    : point -15 50 34.2921 : Relaxation factor = 0.00563771
Info    : point -15 50 34.2921 : Relaxation factor = 0.00422828
Info    : point -15 50 34.2921 : Relaxation factor = 0.00317121
Info    : point -15 50 34.2921 : Relaxation factor = 0.00237841
Info    : point -15 50 34.2921 : Relaxation factor = 0.00178381
Warning : Converged for i=6 j=5 (err=1.71365e-08 iter=2) BUT xyz error = 104.893 in point (-1.500000e+01,5.000000e+01,3.429209e+01) on surface 124
Info    : Reparameterizing point 17 on face 124
Info    : point 25 40 40 : Relaxation factor = 0.75
Info    : point 25 40 40 : Relaxation factor = 0.5625
Info    : point 25 40 40 : Relaxation factor = 0.421875
Info    : point 25 40 40 : Relaxation factor = 0.316406
Info    : point 25 40 40 : Relaxation factor = 0.237305
Info    : point 25 40 40 : Relaxation factor = 0.177979
Info    : point 25 40 40 : Relaxation factor = 0.133484
Info    : point 25 40 40 : Relaxation factor = 0.100113
Info    : point 25 40 40 : Relaxation factor = 0.0750847
Info    : point 25 40 40 : Relaxation factor = 0.0563135
Info    : point 25 40 40 : Relaxation factor = 0.0422351
Info    : point 25 40 40 : Relaxation factor = 0.0316764
Info    : point 25 40 40 : Relaxation factor = 0.0237573
Info    : point 25 40 40 : Relaxation factor = 0.0178179
Info    : point 25 40 40 : Relaxation factor = 0.0133635
Info    : point 25 40 40 : Relaxation factor = 0.0100226
Info    : point 25 40 40 : Relaxation factor = 0.00751695
Info    : point 25 40 40 : Relaxation factor = 0.00563771
Info    : point 25 40 40 : Relaxation factor = 0.00422828
Info    : point 25 40 40 : Relaxation factor = 0.00317121
Info    : point 25 40 40 : Relaxation factor = 0.00237841
Info    : point 25 40 40 : Relaxation factor = 0.00178381
Info    : point 25 40 40 : Relaxation factor = 0.00133786
Info    : point 25 40 40 : Relaxation factor = 0.00100339
Info    : point 25 40 40 : Relaxation factor = 0.000752543
Info    : point 25 40 40 : Relaxation factor = 0.000564408
Info    : point 25 40 40 : Relaxation factor = 0.000423306
Info    : point 25 40 40 : Relaxation factor = 0.000317479
Info    : point 25 40 40 : Relaxation factor = 0.000238109
Info    : point 25 40 40 : Relaxation factor = 0.000178582
Warning : Converged for i=8 j=1 (err=1.80278e-08 iter=2) BUT xyz error = 91.9565 in point (2.500000e+01,4.000000e+01,4.000000e+01) on surface 124
Info    : Reparameterizing point 54 on face 126
Info    : point 300 0 0 : Relaxation factor = 0.75
Info    : point 300 0 0 : Relaxation factor = 0.5625
Info    : point 300 0 0 : Relaxation factor = 0.421875
Info    : point 300 0 0 : Relaxation factor = 0.316406
Info    : point 300 0 0 : Relaxation factor = 0.237305
Info    : point 300 0 0 : Relaxation factor = 0.177979
Info    : point 300 0 0 : Relaxation factor = 0.133484
Info    : point 300 0 0 : Relaxation factor = 0.100113
Info    : point 300 0 0 : Relaxation factor = 0.0750847
Info    : point 300 0 0 : Relaxation factor = 0.0563135
Info    : point 300 0 0 : Relaxation factor = 0.0422351
Info    : point 300 0 0 : Relaxation factor = 0.0316764
Info    : point 300 0 0 : Relaxation factor = 0.0237573
Info    : point 300 0 0 : Relaxation factor = 0.0178179
Info    : point 300 0 0 : Relaxation factor = 0.0133635
Info    : point 300 0 0 : Relaxation factor = 0.0100226
Info    : point 300 0 0 : Relaxation factor = 0.00751695
Info    : point 300 0 0 : Relaxation factor = 0.00563771
Info    : point 300 0 0 : Relaxation factor = 0.00422828
Info    : point 300 0 0 : Relaxation factor = 0.00317121
Info    : point 300 0 0 : Relaxation factor = 0.00237841
Info    : point 300 0 0 : Relaxation factor = 0.00178381
Info    : point 300 0 0 : Relaxation factor = 0.00133786
Info    : point 300 0 0 : Relaxation factor = 0.00100339
Info    : point 300 0 0 : Relaxation factor = 0.000752543
Info    : point 300 0 0 : Relaxation factor = 0.000564408
Info    : point 300 0 0 : Relaxation factor = 0.000423306
Info    : point 300 0 0 : Relaxation factor = 0.000317479
Info    : point 300 0 0 : Relaxation factor = 0.000238109
Info    : point 300 0 0 : Relaxation factor = 0.000178582
Info    : point 300 0 0 : Relaxation factor = 0.000133937
Info    : point 300 0 0 : Relaxation factor = 0.000100452
Info    : point 300 0 0 : Relaxation factor = 7.53393e-05
Info    : point 300 0 0 : Relaxation factor = 5.65045e-05
Info    : point 300 0 0 : Relaxation factor = 4.23784e-05
Info    : point 300 0 0 : Relaxation factor = 3.17838e-05
Warning : Converged for i=5 j=0 (err=1.46489e-08 iter=2) BUT xyz error = 327.969 in point (3.000000e+02,0.000000e+00,0.000000e+00) on surface 126
Info    : Reparameterizing point 6 on face 126
Info    : point 15 50 0 : Relaxation factor = 0.75
Info    : point 15 50 0 : Relaxation factor = 0.5625
Info    : point 15 50 0 : Relaxation factor = 0.421875
Info    : point 15 50 0 : Relaxation factor = 0.316406
Info    : point 15 50 0 : Relaxation factor = 0.237305
Info    : point 15 50 0 : Relaxation factor = 0.177979
Info    : point 15 50 0 : Relaxation factor = 0.133484
Info    : point 15 50 0 : Relaxation factor = 0.100113
Info    : point 15 50 0 : Relaxation factor = 0.0750847
Info    : point 15 50 0 : Relaxation factor = 0.0563135
Info    : point 15 50 0 : Relaxation factor = 0.0422351
Info    : point 15 50 0 : Relaxation factor = 0.0316764
Info    : point 15 50 0 : Relaxation factor = 0.0237573
Info    : point 15 50 0 : Relaxation factor = 0.0178179
Info    : point 15 50 0 : Relaxation factor = 0.0133635
Info    : point 15 50 0 : Relaxation factor = 0.0100226
Info    : point 15 50 0 : Relaxation factor = 0.00751695
Info    : point 15 50 0 : Relaxation factor = 0.00563771
Info    : point 15 50 0 : Relaxation factor = 0.00422828
Info    : point 15 50 0 : Relaxation factor = 0.00317121
Info    : point 15 50 0 : Relaxation factor = 0.00237841
Info    : point 15 50 0 : Relaxation factor = 0.00178381
Info    : point 15 50 0 : Relaxation factor = 0.00133786
Info    : point 15 50 0 : Relaxation factor = 0.00100339
Info    : point 15 50 0 : Relaxation factor = 0.000752543
Info    : point 15 50 0 : Relaxation factor = 0.000564408
Info    : point 15 50 0 : Relaxation factor = 0.000423306
Info    : point 15 50 0 : Relaxation factor = 0.000317479
Info    : point 15 50 0 : Relaxation factor = 0.000238109
Info    : point 15 50 0 : Relaxation factor = 0.000178582
Warning : Converged for i=7 j=2 (err=1.80278e-08 iter=2) BUT xyz error = 44.2267 in point (1.500000e+01,5.000000e+01,0.000000e+00) on surface 126
Warning : Could not orient normal of surface 127
Warning : Could not orient normal of surface 129
Warning : Could not orient normal of surface 130
Info    : Reparameterizing edge 54 on face 139
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.75
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.5625
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.421875
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.316406
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.237305
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.177979
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.133484
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.100113
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.0750847
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.0563135
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.0422351
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.0316764
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.0237573
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.0178179
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.0133635
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.0100226
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.00751695
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.00563771
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.00422828
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.00317121
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.00237841
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.00178381
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.00133786
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.00100339
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.000752543
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.000564408
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.000423306
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.000317479
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.000238109
Info    : point -21.1209 97.2394 50 : Relaxation factor = 0.000178582
Warning : Converged for i=7 j=8 (err=1.77676e-08 iter=2) BUT xyz error = 147.705 in point (-2.112090e+01,9.723940e+01,5.000000e+01) on surface 139
Info    : Reparameterizing edge 34 on face 141
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.75
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.5625
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.421875
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.316406
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.237305
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.177979
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.133484
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.100113
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.0750847
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.0563135
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.0422351
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.0316764
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.0237573
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.0178179
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.0133635
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.0100226
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.00751695
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.00563771
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.00422828
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.00317121
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.00237841
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.00178381
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.00133786
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.00100339
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.000752543
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.000564408
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.000423306
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.000317479
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.000238109
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.000178582
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.000133937
Info    : point -21.6565 -56.6565 49.4245 : Relaxation factor = 0.000100452
Warning : Converged for i=7 j=6 (err=1.74516e-08 iter=2) BUT xyz error = 50.2418 in point (-2.165651e+01,-5.665651e+01,4.942449e+01) on surface 141
Info    : Reparameterizing edge 81 on face 141
Info    : point 275 300 50 : Relaxation factor = 0.75
Info    : point 275 300 50 : Relaxation factor = 0.5625
Info    : point 275 300 50 : Relaxation factor = 0.421875
Info    : point 275 300 50 : Relaxation factor = 0.316406
Info    : point 275 300 50 : Relaxation factor = 0.237305
Info    : point 275 300 50 : Relaxation factor = 0.177979
Info    : point 275 300 50 : Relaxation factor = 0.133484
Info    : point 275 300 50 : Relaxation factor = 0.100113
Info    : point 275 300 50 : Relaxation factor = 0.0750847
Info    : point 275 300 50 : Relaxation factor = 0.0563135
Info    : point 275 300 50 : Relaxation factor = 0.0422351
Info    : point 275 300 50 : Relaxation factor = 0.0316764
Info    : point 275 300 50 : Relaxation factor = 0.0237573
Info    : point 275 300 50 : Relaxation factor = 0.0178179
Info    : point 275 300 50 : Relaxation factor = 0.0133635
Info    : point 275 300 50 : Relaxation factor = 0.0100226
Info    : point 275 300 50 : Relaxation factor = 0.00751695
Info    : point 275 300 50 : Relaxation factor = 0.00563771
Info    : point 275 300 50 : Relaxation factor = 0.00422828
Info    : point 275 300 50 : Relaxation factor = 0.00317121
Info    : point 275 300 50 : Relaxation factor = 0.00237841
Info    : point 275 300 50 : Relaxation factor = 0.00178381
Info    : point 275 300 50 : Relaxation factor = 0.00133786
Info    : point 275 300 50 : Relaxation factor = 0.00100339
Info    : point 275 300 50 : Relaxation factor = 0.000752543
Info    : point 275 300 50 : Relaxation factor = 0.000564408
Info    : point 275 300 50 : Relaxation factor = 0.000423306
Info    : point 275 300 50 : Relaxation factor = 0.000317479
Info    : point 275 300 50 : Relaxation factor = 0.000238109
Info    : point 275 300 50 : Relaxation factor = 0.000178582
Info    : point 275 300 50 : Relaxation factor = 0.000133937
Info    : point 275 300 50 : Relaxation factor = 0.000100452
Info    : point 275 300 50 : Relaxation factor = 7.53393e-05
Info    : point 275 300 50 : Relaxation factor = 5.65045e-05
Info    : point 275 300 50 : Relaxation factor = 4.23784e-05
Info    : point 275 300 50 : Relaxation factor = 3.17838e-05
Info    : point 275 300 50 : Relaxation factor = 2.38378e-05
Info    : point 275 300 50 : Relaxation factor = 1.78784e-05
Info    : point 275 300 50 : Relaxation factor = 1.34088e-05
Info    : point 275 300 50 : Relaxation factor = 1.00566e-05
Info    : point 275 300 50 : Relaxation factor = 7.54244e-06
Warning : Converged for i=8 j=8 (err=1.56089e-08 iter=2) BUT xyz error = 436.119 in point (2.750000e+02,3.000000e+02,5.000000e+01) on surface 141
Info    : Reparameterizing edge 80 on face 143
Info    : point 300 75 50 : Relaxation factor = 0.75
Info    : point 300 75 50 : Relaxation factor = 0.5625
Info    : point 300 75 50 : Relaxation factor = 0.421875
Info    : point 300 75 50 : Relaxation factor = 0.316406
Info    : point 300 75 50 : Relaxation factor = 0.237305
Info    : point 300 75 50 : Relaxation factor = 0.177979
Info    : point 300 75 50 : Relaxation factor = 0.133484
Info    : point 300 75 50 : Relaxation factor = 0.100113
Info    : point 300 75 50 : Relaxation factor = 0.0750847
Info    : point 300 75 50 : Relaxation factor = 0.0563135
Info    : point 300 75 50 : Relaxation factor = 0.0422351
Info    : point 300 75 50 : Relaxation factor = 0.0316764
Info    : point 300 75 50 : Relaxation factor = 0.0237573
Info    : point 300 75 50 : Relaxation factor = 0.0178179
Info    : point 300 75 50 : Relaxation factor = 0.0133635
Info    : point 300 75 50 : Relaxation factor = 0.0100226
Info    : point 300 75 50 : Relaxation factor = 0.00751695
Info    : point 300 75 50 : Relaxation factor = 0.00563771
Info    : point 300 75 50 : Relaxation factor = 0.00422828
Info    : point 300 75 50 : Relaxation factor = 0.00317121
Info    : point 300 75 50 : Relaxation factor = 0.00237841
Info    : point 300 75 50 : Relaxation factor = 0.00178381
Info    : point 300 75 50 : Relaxation factor = 0.00133786
Info    : point 300 75 50 : Relaxation factor = 0.00100339
Info    : point 300 75 50 : Relaxation factor = 0.000752543
Info    : point 300 75 50 : Relaxation factor = 0.000564408
Info    : point 300 75 50 : Relaxation factor = 0.000423306
Info    : point 300 75 50 : Relaxation factor = 0.000317479
Info    : point 300 75 50 : Relaxation factor = 0.000238109
Info    : point 300 75 50 : Relaxation factor = 0.000178582
Info    : point 300 75 50 : Relaxation factor = 0.000133937
Info    : point 300 75 50 : Relaxation factor = 0.000100452
Info    : point 300 75 50 : Relaxation factor = 7.53393e-05
Info    : point 300 75 50 : Relaxation factor = 5.65045e-05
Warning : Converged for i=8 j=8 (err=1.71452e-08 iter=2) BUT xyz error = 277.399 in point (3.000000e+02,7.500000e+01,5.000000e+01) on surface 143
Info    : Reparameterizing edge 4 on face 145
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.75
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.5625
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.421875
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.316406
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.237305
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.177979
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.133484
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.100113
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.0750847
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.0563135
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.0422351
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.0316764
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.0237573
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.0178179
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.0133635
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.0100226
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.00751695
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.00563771
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.00422828
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.00317121
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.00237841
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.00178381
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.00133786
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.00100339
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.000752543
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.000564408
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.000423306
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.000317479
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.000238109
Info    : point -22.0711 47.0711 0 : Relaxation factor = 0.000178582
Warning : Converged for i=0 j=7 (err=1.29252e-08 iter=2) BUT xyz error = 51.9615 in point (-2.207107e+01,4.707107e+01,0.000000e+00) on surface 145
Info    : Reparameterizing edge 14 on face 145
Info    : point 25 -35 40 : Relaxation factor = 0.75
Info    : point 25 -35 40 : Relaxation factor = 0.5625
Info    : point 25 -35 40 : Relaxation factor = 0.421875
Info    : point 25 -35 40 : Relaxation factor = 0.316406
Info    : point 25 -35 40 : Relaxation factor = 0.237305
Info    : point 25 -35 40 : Relaxation factor = 0.177979
Info    : point 25 -35 40 : Relaxation factor = 0.133484
Info    : point 25 -35 40 : Relaxation factor = 0.100113
Info    : point 25 -35 40 : Relaxation factor = 0.0750847
Info    : point 25 -35 40 : Relaxation factor = 0.0563135
Info    : point 25 -35 40 : Relaxation factor = 0.0422351
Info    : point 25 -35 40 : Relaxation factor = 0.0316764
Info    : point 25 -35 40 : Relaxation factor = 0.0237573
Info    : point 25 -35 40 : Relaxation factor = 0.0178179
Info    : point 25 -35 40 : Relaxation factor = 0.0133635
Info    : point 25 -35 40 : Relaxation factor = 0.0100226
Info    : point 25 -35 40 : Relaxation factor = 0.00751695
Info    : point 25 -35 40 : Relaxation factor = 0.00563771
Info    : point 25 -35 40 : Relaxation factor = 0.00422828
Info    : point 25 -35 40 : Relaxation factor = 0.00317121
Info    : point 25 -35 40 : Relaxation factor = 0.00237841
Info    : point 25 -35 40 : Relaxation factor = 0.00178381
Info    : point 25 -35 40 : Relaxation factor = 0.00133786
Info    : point 25 -35 40 : Relaxation factor = 0.00100339
Info    : point 25 -35 40 : Relaxation factor = 0.000752543
Info    : point 25 -35 40 : Relaxation factor = 0.000564408
Info    : point 25 -35 40 : Relaxation factor = 0.000423306
Info    : point 25 -35 40 : Relaxation factor = 0.000317479
Info    : point 25 -35 40 : Relaxation factor = 0.000238109
Info    : point 25 -35 40 : Relaxation factor = 0.000178582
Warning : Converged for i=2 j=8 (err=1.77635e-08 iter=2) BUT xyz error = 94.9145 in point (2.500000e+01,-3.500000e+01,4.000000e+01) on surface 145
Warning : Could not orient normal of surface 156
Warning : Could not orient normal of surface 158
Info    : Reparameterizing point 47 on face 159
Info    : point 300 300 -200 : Relaxation factor = 0.75
Info    : point 300 300 -200 : Relaxation factor = 0.5625
Info    : point 300 300 -200 : Relaxation factor = 0.421875
Info    : point 300 300 -200 : Relaxation factor = 0.316406
Info    : point 300 300 -200 : Relaxation factor = 0.237305
Info    : point 300 300 -200 : Relaxation factor = 0.177979
Info    : point 300 300 -200 : Relaxation factor = 0.133484
Info    : point 300 300 -200 : Relaxation factor = 0.100113
Info    : point 300 300 -200 : Relaxation factor = 0.0750847
Info    : point 300 300 -200 : Relaxation factor = 0.0563135
Info    : point 300 300 -200 : Relaxation factor = 0.0422351
Info    : point 300 300 -200 : Relaxation factor = 0.0316764
Info    : point 300 300 -200 : Relaxation factor = 0.0237573
Info    : point 300 300 -200 : Relaxation factor = 0.0178179
Info    : point 300 300 -200 : Relaxation factor = 0.0133635
Info    : point 300 300 -200 : Relaxation factor = 0.0100226
Info    : point 300 300 -200 : Relaxation factor = 0.00751695
Info    : point 300 300 -200 : Relaxation factor = 0.00563771
Info    : point 300 300 -200 : Relaxation factor = 0.00422828
Info    : point 300 300 -200 : Relaxation factor = 0.00317121
Info    : point 300 300 -200 : Relaxation factor = 0.00237841
Info    : point 300 300 -200 : Relaxation factor = 0.00178381
Info    : point 300 300 -200 : Relaxation factor = 0.00133786
Info    : point 300 300 -200 : Relaxation factor = 0.00100339
Info    : point 300 300 -200 : Relaxation factor = 0.000752543
Info    : point 300 300 -200 : Relaxation factor = 0.000564408
Info    : point 300 300 -200 : Relaxation factor = 0.000423306
Info    : point 300 300 -200 : Relaxation factor = 0.000317479
Info    : point 300 300 -200 : Relaxation factor = 0.000238109
Info    : point 300 300 -200 : Relaxation factor = 0.000178582
Info    : point 300 300 -200 : Relaxation factor = 0.000133937
Info    : point 300 300 -200 : Relaxation factor = 0.000100452
Info    : point 300 300 -200 : Relaxation factor = 7.53393e-05
Info    : point 300 300 -200 : Relaxation factor = 5.65045e-05
Warning : Converged for i=8 j=2 (err=1.91835e-08 iter=2) BUT xyz error = 613.544 in point (3.000000e+02,3.000000e+02,-2.000000e+02) on surface 159
Info    : Reparameterizing edge 43 on face 160
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.75
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.5625
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.421875
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.316406
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.237305
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.177979
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.133484
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.100113
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.0750847
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.0563135
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.0422351
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.0316764
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.0237573
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.0178179
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.0133635
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.0100226
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.00751695
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.00563771
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.00422828
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.00317121
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.00237841
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.00178381
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.00133786
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.00100339
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.000752543
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.000564408
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.000423306
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.000317479
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.000238109
Info    : point 45.0763 -58.9812 0 : Relaxation factor = 0.000178582
Warning : Converged for i=4 j=2 (err=1.94815e-08 iter=2) BUT xyz error = 194.896 in point (4.507625e+01,-5.898118e+01,0.000000e+00) on surface 160
Info    : Reparameterizing edge 81 on face 160
Info    : point 225 300 50 : Relaxation factor = 0.75
Info    : point 225 300 50 : Relaxation factor = 0.5625
Info    : point 225 300 50 : Relaxation factor = 0.421875
Info    : point 225 300 50 : Relaxation factor = 0.316406
Info    : point 225 300 50 : Relaxation factor = 0.237305
Info    : point 225 300 50 : Relaxation factor = 0.177979
Info    : point 225 300 50 : Relaxation factor = 0.133484
Info    : point 225 300 50 : Relaxation factor = 0.100113
Info    : point 225 300 50 : Relaxation factor = 0.0750847
Info    : point 225 300 50 : Relaxation factor = 0.0563135
Info    : point 225 300 50 : Relaxation factor = 0.0422351
Info    : point 225 300 50 : Relaxation factor = 0.0316764
Info    : point 225 300 50 : Relaxation factor = 0.0237573
Info    : point 225 300 50 : Relaxation factor = 0.0178179
Info    : point 225 300 50 : Relaxation factor = 0.0133635
Info    : point 225 300 50 : Relaxation factor = 0.0100226
Info    : point 225 300 50 : Relaxation factor = 0.00751695
Info    : point 225 300 50 : Relaxation factor = 0.00563771
Info    : point 225 300 50 : Relaxation factor = 0.00422828
Info    : point 225 300 50 : Relaxation factor = 0.00317121
Info    : point 225 300 50 : Relaxation factor = 0.00237841
Info    : point 225 300 50 : Relaxation factor = 0.00178381
Info    : point 225 300 50 : Relaxation factor = 0.00133786
Info    : point 225 300 50 : Relaxation factor = 0.00100339
Info    : point 225 300 50 : Relaxation factor = 0.000752543
Info    : point 225 300 50 : Relaxation factor = 0.000564408
Info    : point 225 300 50 : Relaxation factor = 0.000423306
Info    : point 225 300 50 : Relaxation factor = 0.000317479
Info    : point 225 300 50 : Relaxation factor = 0.000238109
Warning : Converged for i=1 j=8 (err=1.9315e-08 iter=2) BUT xyz error = 150.083 in point (2.250000e+02,3.000000e+02,5.000000e+01) on surface 160
Info    : Reparameterizing edge 69 on face 161
Info    : point 200 300 0 : Relaxation factor = 0.75
Info    : point 200 300 0 : Relaxation factor = 0.5625
Info    : point 200 300 0 : Relaxation factor = 0.421875
Info    : point 200 300 0 : Relaxation factor = 0.316406
Info    : point 200 300 0 : Relaxation factor = 0.237305
Info    : point 200 300 0 : Relaxation factor = 0.177979
Info    : point 200 300 0 : Relaxation factor = 0.133484
Info    : point 200 300 0 : Relaxation factor = 0.100113
Info    : point 200 300 0 : Relaxation factor = 0.0750847
Info    : point 200 300 0 : Relaxation factor = 0.0563135
Info    : point 200 300 0 : Relaxation factor = 0.0422351
Info    : point 200 300 0 : Relaxation factor = 0.0316764
Info    : point 200 300 0 : Relaxation factor = 0.0237573
Info    : point 200 300 0 : Relaxation factor = 0.0178179
Info    : point 200 300 0 : Relaxation factor = 0.0133635
Info    : point 200 300 0 : Relaxation factor = 0.0100226
Info    : point 200 300 0 : Relaxation factor = 0.00751695
Info    : point 200 300 0 : Relaxation factor = 0.00563771
Info    : point 200 300 0 : Relaxation factor = 0.00422828
Info    : point 200 300 0 : Relaxation factor = 0.00317121
Info    : point 200 300 0 : Relaxation factor = 0.00237841
Info    : point 200 300 0 : Relaxation factor = 0.00178381
Info    : point 200 300 0 : Relaxation factor = 0.00133786
Info    : point 200 300 0 : Relaxation factor = 0.00100339
Info    : point 200 300 0 : Relaxation factor = 0.000752543
Info    : point 200 300 0 : Relaxation factor = 0.000564408
Info    : point 200 300 0 : Relaxation factor = 0.000423306
Info    : point 200 300 0 : Relaxation factor = 0.000317479
Info    : point 200 300 0 : Relaxation factor = 0.000238109
Info    : point 200 300 0 : Relaxation factor = 0.000178582
Info    : point 200 300 0 : Relaxation factor = 0.000133937
Warning : Converged for i=4 j=8 (err=1.6763e-08 iter=2) BUT xyz error = 290 in point (2.000000e+02,3.000000e+02,0.000000e+00) on surface 161
Info    : Reparameterizing edge 61 on face 162
Info    : point -100 -300 0 : Relaxation factor = 0.75
Info    : point -100 -300 0 : Relaxation factor = 0.5625
Info    : point -100 -300 0 : Relaxation factor = 0.421875
Info    : point -100 -300 0 : Relaxation factor = 0.316406
Info    : point -100 -300 0 : Relaxation factor = 0.237305
Info    : point -100 -300 0 : Relaxation factor = 0.177979
Info    : point -100 -300 0 : Relaxation factor = 0.133484
Info    : point -100 -300 0 : Relaxation factor = 0.100113
Info    : point -100 -300 0 : Relaxation factor = 0.0750847
Info    : point -100 -300 0 : Relaxation factor = 0.0563135
Info    : point -100 -300 0 : Relaxation factor = 0.0422351
Info    : point -100 -300 0 : Relaxation factor = 0.0316764
Info    : point -100 -300 0 : Relaxation factor = 0.0237573
Info    : point -100 -300 0 : Relaxation factor = 0.0178179
Info    : point -100 -300 0 : Relaxation factor = 0.0133635
Info    : point -100 -300 0 : Relaxation factor = 0.0100226
Info    : point -100 -300 0 : Relaxation factor = 0.00751695
Info    : point -100 -300 0 : Relaxation factor = 0.00563771
Info    : point -100 -300 0 : Relaxation factor = 0.00422828
Info    : point -100 -300 0 : Relaxation factor = 0.00317121
Info    : point -100 -300 0 : Relaxation factor = 0.00237841
Info    : point -100 -300 0 : Relaxation factor = 0.00178381
Info    : point -100 -300 0 : Relaxation factor = 0.00133786
Info    : point -100 -300 0 : Relaxation factor = 0.00100339
Info    : point -100 -300 0 : Relaxation factor = 0.000752543
Info    : point -100 -300 0 : Relaxation factor = 0.000564408
Info    : point -100 -300 0 : Relaxation factor = 0.000423306
Info    : point -100 -300 0 : Relaxation factor = 0.000317479
Info    : point -100 -300 0 : Relaxation factor = 0.000238109
Info    : point -100 -300 0 : Relaxation factor = 0.000178582
Info    : point -100 -300 0 : Relaxation factor = 0.000133937
Warning : Converged for i=0 j=6 (err=1.99542e-08 iter=2) BUT xyz error = 223.506 in point (-1.000000e+02,-3.000000e+02,0.000000e+00) on surface 162
Info    : Reparameterizing edge 19 on face 162
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.75
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.5625
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.421875
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.316406
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.237305
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.177979
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.133484
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.100113
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.0750847
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.0563135
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.0422351
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.0316764
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.0237573
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.0178179
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.0133635
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.0100226
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.00751695
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.00563771
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.00422828
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.00317121
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.00237841
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.00178381
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.00133786
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.00100339
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.000752543
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.000564408
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.000423306
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.000317479
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.000238109
Info    : point 32.0711 57.0711 50 : Relaxation factor = 0.000178582
Warning : Converged for i=7 j=1 (err=1.91946e-08 iter=2) BUT xyz error = 267.062 in point (3.207107e+01,5.707107e+01,5.000000e+01) on surface 162
Warning : Could not orient normal of surface 189
Warning : Could not orient normal of surface 190
Warning : Could not orient normal of surface 192
Warning : Could not orient normal of surface 201
Warning : Could not orient normal of surface 204
Warning : Could not orient normal of surface 205
Warning : Could not orient normal of surface 206
Warning : Could not orient normal of surface 208
Warning : Could not orient normal of surface 221
Warning : Could not orient normal of surface 222
Warning : Could not orient normal of surface 223
Warning : Could not orient normal of surface 224
Info    : 3712 vertices 25644 elements
Warning : ------------------------------
Warning : Mesh generation error summary
Warning :    36 warnings
Warning :     0 errors
Warning : Check the full log for details
Warning : ------------------------------
Info    : Writing '/home/fomins/brr.txt'...