PHOTON USE
  p
  p1
 
 
  msg  Drilling the sphere
  GR Z 1
  con mark x 1 fil;.01
  con mark x 20 fil;.01
  con mark z 20 fil;.01
  msg  Please wait .....
  surf mark x 0.99
  msg  Hit Enter to continue
  pause
  p
  p2
 
 
  msg  Drilling the cylinder
  gr z 1 y 1 19
  gr z m y 1 19
  msg  Please wait .....
  surf mark x 0.99
  msg  Hit Enter to continue
  pause
  p
  p3
 
 
  msg  Drilling the bullet
  gr z 5 y 1 19
  msg  Please wait .....
  surf mark x 0.99
  msg  Hit Enter to continue
  pause
  p
  p4
 
 
  msg  Drilling the annulus
  GR Z 1
  con mark Z 1 X 1 20 Y 1 M fil;.01
  con mark Z 10 X 1 20 Y 1 M fil;.01
  con mark Z m X 1 20 Y 1 M fil;.01
  GR OU Y M;GR OU Z M;GR OU Z M
  msg  Please wait .....
  surf mark x 0.99
  msg  Hit E to finish
  enduse
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
 
  DISPLAY
    This input demonstrates the "drilling-milling" technique
    available  in  PLANT  to generate marker distribution to
    fit the complex shapes on polar,  CARTES=F,  grids.  The
    method rests on PLANT ability to handle virtually "live"
    mathematical  expressions  as  the  arguments   of   its
    functions.
 
    All cases  of the file are wholly focused on calculation
    of marker disributions.  Therefore, no other actions are
    supported by input data.
 
    Four different   cases,   in   which   the  variants  of
    "drilling"  are  used  to  get  different  shapes,   are
    arranged  as  unsteady problem:  at each time moment the
    new marker distribution is PLANTed calculated and dumped
    into  specified  file  to  be  viewed by PHOTON commands
    supplied.
  ENDDIS
 
    PLANT information :
     * Data input groups used: 11, 19
     * Ground groups planted : 11, 19-6
     * Headings used  : INIT??, SC06??
     * Functions used : SPHERE, XYCIRC
     * Commands used  : IF, REGION
 
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
    GROUP 1. Run title and other preliminaries
TEXT( DRILling CYLINDER: CARTES=F
    GROUP 2. Transience; time-step specification
STEADY=F
GRDPWR(T,4,4.,1.0)
    GROUP 3. X-direction grid specification
CARTES=F;GRDPWR(X,40,2.*3.14,1.0)
    GROUP 4. Y-direction grid specification
GRDPWR(Y,40,10.,1.0)
    GROUP 5. Z-direction grid specification
GRDPWR(Z,40,20.,1.0)
    GROUP 7. Variables stored, solved & named
STORE(MARK)
    GROUP 11. Initialization of variable or porosity fields
iniadd=f;FIINIT(MARK)=0.0
 
NAMSAT=MOSG

   PLANTBEGIN
       1. Drilling the sphere
PATCH(INI1,INIVAL,1,NX,1,NY,1,NZ,1,1)
   VAL=SPHERE(1.,10.,10.,10.,5.)
INIT (INI1,MARK,0.,GRND)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    Ih the above settings the  sphere  is  "drilled-out"  in
    cylinder.   Crucial  faeture  of  this  example  is  the
    specification  of   center   coordinates:   unlike   for
    cartesian cases,  the line of X=0 and Y=0 in cylindrical
    polar coordinate system has got the coordinates  Xr=RV2D
    and  Yr=RV2D  in  reference  frame.  It is the latter in
    which the arguments of  geometrical  functions  must  be
    specified.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
       2. Drilling the cylinder
   MARK=0.
  IF(ISTEP.EQ.2.AND.LSWEEP.EQ.1)
   MARK=SPHERE(1.,10.,10.,ZGNZ,5.)
  IF(ISTEP.EQ.2.AND.LSWEEP.EQ.1)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    This example  seems  to  be  even simpler than foregoing
    one.  It deals  with  the  "drilling-out"  the  cylinder
    cavity  of 5 m diameter in coaxial cylindrical domain by
    spherical  "drill  bit".  The  "drilling  "  effect   is
    activated  by  using  the  distances  from  Z=0 plane as
    coordinates of sphere centre.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
       3. Drilling the bullet shape
   MARK=0.
  IF(ISTEP.EQ.3.AND.LSWEEP.EQ.1)
   MARK=SPHERE(1.,10.,10.,10.,5.)
  IF(ISTEP.EQ.3.AND.LSWEEP.EQ.1)
   MARK=XYCIRC(1.,10.,10.,5.)
  REGION(1,NX,1,NY,5,20,3,3) /LSWEEP.EQ.1
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    Another variant of  "drilling"  technique  is  used  for
    specification  of  bullet-shape  cavity.  The  shape  is
    combination  of  stationary   sphere   and   2D   circle
    "drill-bit".  The  "drilling" is activated by extents of
    REGION commands over which XYCIRC function  is  applied.
    The  command  PLACE  may  also be used with advantage to
    employ physical rather than cell numbers extents.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
       4. Drilling the annulus
   MARK=0.
  IF(ISTEP.EQ.4.AND.LSWEEP.EQ.1)
   MARK=SPHERE(1.,15.,10.,ZGNZ,4.)
  IF(ISTEP.EQ.4.AND.LSWEEP.EQ.1)
  >>>>>>>>>>>>>>>>>>>>>> Comment begins >>>>>>>>>>>>>>>>>>>>
    This examples differs from "drilling" the cylinder  case
    mainly  in  shifting  the axis of "drill bit" of 5 meter
    away from domain centre line.
  <<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
  PLANTEND

    GROUP 15. Termination of sweeps
LSWEEP=1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
    GROUP 23. Print-out & plot control
idispa=1; csg1=p; selref=t; resfac=1.e-2
dmpstk=t
DISTIL=T
EX(MARK)=1.812E-01
 LIBREF=603
STOP