[Gmsh] cut high order postprocessing views

Christophe Geuzaine cgeuzaine at ulg.ac.be
Fri Oct 2 15:40:37 CEST 2009


Marcus Drosson wrote:
> Dear all,
> 
> is there any way to use the "adapt visualization grid" without the GUI? 
> I have defined a gmsh script, which cuts the solution field at a given 
> section and exports the result as a raw text file. The postprocessing 
> view is composed of high order polynomials.
> My problem: when I run the script, the cut section contains only a few 
> points, corresponding to the initial grid. When I use the GUI to first 
> adapt the visualization grid and then run "cut plane" from the GUI, a 
> much better representation is obtained. (I am particularly interested in 
> the gradient close to the wall. That's why one of the reasons why I need 
> to adapt the grid.)
> As I have to do the same operation many times, I really would like to 
> use a script.
> 

Hello Marcus - Here's an example:

Merge "adaptive_view.pos";

View[0].AdaptVisualizationGrid =1;

Plugin(CutPlane).RecurLevel = 4;
Plugin(CutPlane).TargetError = 0;
Plugin(CutPlane).D = -3;
Plugin(CutPlane).iView = 0;
Plugin(CutPlane).Run;


(Beware that all plugins are not yet ready for high-order maps. Extract, 
Evaluate & co will just take components "as is" : they won't perform any 
interpolation.)

Cheers,

Christophe



> Thanks in advance
> Marcus
> 


-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science
http://www.montefiore.ulg.ac.be/~geuzaine