* LAGriT Example Input file * Carl Gable * gable@lanl.gov * * Wed Aug 21 22:15:45 UTC 2002 * * Triangulate a set of polygons, stitch the polygons together * and refine and smooth the final mesh. * *--------------------------------- * * Triangulate the first polygon * cmo / create / cmotri1 / / / tri read / avs / polygon1.inp / cmotri1 define / ICOLOR / 1 * * Call another LAGriT input file which has the triangulation commands * infile / lagrit.input_tri_macro dump / gmv / output_tri1.gmv /cmotri1 ********** * * Triangulate the next polygon * cmo / create / cmotri2 / / / tri read / avs / polygon2.inp / cmotri2 define / ICOLOR / 2 infile / lagrit.input_tri_macro dump / gmv / output_tri2.gmv / cmotri2 dump / avs / output_tri2.inp / cmotri2 **** * * Triangulate the next polygon * cmo / create / cmotri3 / / / tri read / avs / polygon3.inp / cmotri3 define / ICOLOR / 3 infile / lagrit.input_tri_macro dump / gmv / output_tri3.gmv / cmotri3 * * * Triangulate the next polygon * cmo / create / cmotri4 / / / tri read / avs / polygon4.inp / cmotri4 define / ICOLOR / 4 infile / lagrit.input_tri_macro dump / gmv / output_tri4.gmv / cmotri4 * * * Triangulate the next polygon * cmo / create / cmotri5 / / / tri read / avs / polygon5.inp / cmotri5 define / ICOLOR / 5 infile / lagrit.input_tri_macro dump / gmv / output_tri5.gmv / cmotri5 * * Merge all the triangulated polygons into a single MO * addmesh / merge / cmo_all / cmotri1 / cmotri2 addmesh / merge / cmo_all / cmo_all / cmotri3 addmesh / merge / cmo_all / cmo_all / cmotri4 addmesh / merge / cmo_all / cmo_all / cmotri5 * * Figure out where there are duplicate points and remove duplicates * filter / 1 0 0 * * Recompute the jtet element adjancy information. * geniee * * Recompute the itp array * resetpts / itp * * Create parent/child points at material interfaces * settets * * Reset the imt values of nodes to be the same as the element itetclr * resetpts / itetclr dump / gmv / output_all.gmv / cmo_all * * Refine all triangles using edge mid-point bisection * refine2d resetpts / itp * recon 0 recon 1 smooth / position / esug / 1 0 0 dump / gmv / output_smooth1.gmv / cmo_all * * Refine again * refine2d resetpts / itp smooth / position / esug / 1 0 0 * recon 1 recon 0 smooth / position / esug / 1 0 0 * recon 1 recon 0 dump / gmv / output_smooth2.gmv / cmo_all * * Refine again * refine2d resetpts / itp recon 0 smooth / position / esug / 1 0 0 recon 0 smooth / position / esug / 1 0 0 recon 0 smooth / position / esug / 1 0 0 recon 0 smooth / position / esug / 1 0 0 * recon 1 recon 0 dump / gmv / output_smooth3.gmv / cmo_all * * Refine again * refine2d resetpts / itp recon 0 smooth / position / esug / 1 0 0 recon 0 smooth / position / esug / 1 0 0 recon 0 smooth / position / esug / 1 0 0 recon 0 smooth / position / esug / 1 0 0 * recon 1 recon 0 dump / gmv / output_smooth4.gmv / cmo_all finish