# Make the ISMIP mesh B using awk
# Start from a square 1 x 1
# 
rm -rf square
rm -rf mesh_Bawk
ElmerGrid 1 2 square.grd
mv  square/mesh.nodes square/mesh0.nodes
awk '{print $1, $2, $3*20000.0, $4, $5, (-$3*20000.0*0.00872686779), (-1000.0+(500.0*sin(6.2831853*$3)))}' square/mesh0.nodes > square/mesh1.nodes
awk '{print $1, $2, $3, $6+((1.0-$4)*$7) , $5 }' square/mesh1.nodes > square/mesh.nodes
mv -f square mesh_Bawk
rm -f mesh_Bawk/mesh?.nodes
ElmerGrid 2 3 mesh_Bawk 

