*create a cylinder parallel to the y-axis inside a box

*use the coordsys command

cmo/create/3dmesh

surface/box/reflect /box/ 0.0,0.0, 0.0/ 1.0, 1.0, 1.0/

surface/h1/intrface /cylinder/ 0.2, -0.1, 0.8/ 0.2, 1.1, 0.8/.1/

region/H1/ le box and le h1 /

region/Fill/ le box and gt h1 /

mregion/Air/ le box and lt h1 /

mregion/Solid/Fill

rz/xyz/11,11,1/0.0,0.0,1.1/1.0,1.0,1.1/1,1,1/

pset/rays/seq/1,0,0/

regnpts/Fill/11/pset,get,rays/xyz/ 0.0,0.0,-0.1/1.0,1.0,-0.1/ &

0.0,1.0,-0.1/0,0/

coordsys/define/0.2,0.,0.8/1.,1.,0.8/0.2,2.0,0.8/

rz/rtz/1,7,11/5.,0.,0./5.,360.,1./0,1,1/0,0,0/

coordsys/normal

pset/ray1/seq/0,0,0/

regnpts/H1/3/pset,get,ray1/rtz/ 0.2,0.0, 0.8/0.2,1.0, 0.8/0,0/

filter/1,0,0/

zq/imt/1,0,0/1/

zq/itp/pset,get,ray1/dud/

zq/itp/pset,get,rays/dud/

setpts

search

settets

dump/gmv/gmvr.coordsys/3dmesh

finish

 

 

 

 

Return to LaGriT Home Page