:
h=$1

cat > square-${h}.cad << EOF
MeshVersionFormatted
  0
Dimension
  2
Vertices
  4
  0  0     1
  1  0     2
  1  1     3
  0  1     4
Edges
  4
  1  2     1
  2  3     2
  3  4     3
  4  1     4
hVertices
  ${h} ${h} ${h} ${h}
EOF

echo "! file \"square-${h}.cad\" created." 1>&2

echo "bamg -g square-${h}.cad -o square-${h}.bamg" 1>&2
bamg -g square-${h}.cad -o square-${h}.bamg

cat > square.dmn << EOF
EdgeDomainNames
4
bottom
right
top
left
EOF
echo "! file \"square.dmn\" created." 1>&2

echo "bamg2geo square-${h}.bamg square.dmn > square-${h}.geo" 1>&2
bamg2geo square-${h}.bamg square.dmn > square-${h}.geo


