Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
mesh:mshdem [2013/05/30 22:51]
gag
mesh:mshdem [2014/03/26 10:37] (current)
gag
Line 3: Line 3:
  
 * **Source Fortran File:** ''MshGlacierDEM.f90'' * **Source Fortran File:** ''MshGlacierDEM.f90''
 +
 +**NEW:** since version 6626, works for parallel meshes (Warning, the ''mesh_input.dat'' has to be modified, also for serial mesh)
  
 This tool deforms vertically a 1m thick 3D mesh to conform with the bedrock and surface elevations given as DEMs. Informations (name and grid definition) of the bedrock and surface DEMs are given is the file ''mesh_input.dat''. This file reads:\\ This tool deforms vertically a 1m thick 3D mesh to conform with the bedrock and surface elevations given as DEMs. Informations (name and grid definition) of the bedrock and surface DEMs are given is the file ''mesh_input.dat''. This file reads:\\
Line 25: Line 27:
 600.0 350.0 600.0 350.0
 ! Minimum ice thickness ! Minimum ice thickness
 +1.0
 +! Number of Partitions
 1.0 1.0
 </code>   </code>  
-The last parameter (''Minimum ice thickness''sets a minimum ice thickness to ensure strictly positive volume of the elements. +The parameter ''Minimum ice thickness'' sets a minimum ice thickness to ensure strictly positive volume of the elements. The parameter ''Number of Partitions'' is 1 for a serial mesh and equals to the number of partitions for a partitioned mesh.
  
 Nodal value of elevation are estimated from the bedrock and surface DEMs using bilinear interpolation. The DEMs are given as 2D regular grid in an ASCII file (x y z).  The data in the ASCII file have to be ordered such that    Nodal value of elevation are estimated from the bedrock and surface DEMs using bilinear interpolation. The DEMs are given as 2D regular grid in an ASCII file (x y z).  The data in the ASCII file have to be ordered such that   
 <code> <code>
-  x1 y1 val11 +  x1 y1 z11 
-  x2 y1 val21+  x2 y1 z21
   ...   ...
-  xn y1 valn1 +  xn y1 zn1 
-  x1 y2 val12+  x1 y2 z12
   ...   ...
-  xn yn valnn +  xn yn znn 
 </code> </code>
  
Line 45: Line 49:
   * (Nx, Ny) the number of cells in x and y directions    * (Nx, Ny) the number of cells in x and y directions 
   * No data are given by -9999   * No data are given by -9999
 +  
  
- +An example using ''MshGlacierDEM'' to construct Tête Rousse glacier mesh can be found in ''[ELMER_TRUNK]/elmerice/examples/Test_MshGlacierDEM'' for both a Serial and a Partitioned meshes.
-Partitioning, if needed, has to be performed after ''MshGlacierDEM''  +
- +
-An example using ''MshGlacierDEM'' to construct Tête Rousse glacier mesh can be downloaded {{:mesh:test_mshglacierdem.tar.gz|here}}.+
  
  
mesh/mshdem.1369954304.txt.gz · Last modified: 2013/05/30 22:51 by gag
CC Attribution-Share Alike 4.0 International
www.chimeric.de Valid CSS Driven by DokuWiki do yourself a favour and use a real browser - get firefox!! Recent changes RSS feed Valid XHTML 1.0