doc/test_doc/hemispherical_rind/hemispherical_grid.py

   1
   2 ####### Hemispherical Grid Test ##########
   3
   4 import hexablock
   5
   6
   7 # ======================================================= test_sphere
   8 def test_sphere () :
   9
  10     doc  = hexablock.addDocument ("default")
  11
  12     center  = doc.addVertex (0,0,0)
  13     radius  = 8
  14     orig    = doc.addVertex (0,0,0)
  15     vz      = doc.addVector (0,0,1)
  16     vx      = doc.addVector (1,0,0)
  17     radhole = 1
  18     nrad    = 3
  19     nang    = 16
  20     nhaut   = 8
  21
  22     doc.makeSphere (center, vx, vz, radius, radhole, orig, nrad, nang, nhaut)
  23     doc.saveVtk ("hemispherical.vtk")
  24     return doc
  25
  26 # ================================================================= Test
  27 doc = test_sphere  ()
  28