
<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;"><font style="font-family: Arial;" size="2"><span style="font-size: 10pt;">Dear Gmsh users.<br>

<br>

I am part of a research team in the Applied Modeling and Computational <br>

Group at Imperial College London investigating the meshing of realistic <br>

domains. We are developing Python plugins for Quantum GIS that allow us <br>

to create Geo files from Shapefile polygon layers and use NetCDF files <br>

as a mesh-size metric followed by the creation of a mesh.<br>

<br>

For example, using GSHHS Shapefile data it is possible to create a <br>

domain of the sea surrounding the British Isles and using data stored in <br>

a NetCDF file to produce a finer mesh size in areas of shallow water and <br>

coarser resolution where the depth is greater. This has been achieved on <br>

a planar surface by creating a PostView (.pos) file from the NetCDF <br>

file, merging this to the Geo file, running the 'Triangulate' plugin and <br>

then using this as the background field.<br>

<br>

We would like to be able to perform a similar operation on the sphere in <br>

Cartesian space. The current code allows the projection of the domain <br>

onto a sphere and creates a functional Geo file. However, it has not <br>

been possible (in our experience) to implement the same process of using <br>

a PostView file containing the data to define the mesh size at a <br>

specific location on the sphere, as we want mesh size data on the <br>

surface of the sphere, but the PostView and triangulate use the 3D space <br>

(hence we get a non-usable PostView). Instead, we have used a Structured <br>

Field and used the LonLat operator to put that on the sphere surface.<br>

<br>

We have written a Python script which converts a NetCDF file into a <br>

*.dat file for use as a 'Structured' field inside Gmsh. This is imported <br>

and a 'LonLat' field is added. In a global domain, this works fine. <br>

However, when using a small area (e.g. around the UK) when visualised, <br>

the LonLat field is heavily distorted, appearing to have been stretched <br>

from the north pole, and indeed appears centered over the North pole.<br>

<br>

Does anybody have experience of getting limited (i.e. not global) fields <br>

into GMSH on the sphere (rather than in stereographic space)? Are there <br>

any other suitable methods of getting spatial data into GMSH for use as <br>

a mesh background field, other than Structured field or as a PostView?<br>

<br>

Thank you in advance,<br>

Shaun Lee</span></font></div>

</body>

</html>