/out,g_area,out
 /com,g_area.mac
!  computes the area of the active elements
!  output:
!     t_area   =  total area of the active elements
!     x_xct    =  x centroid of the active elements
!
 *get,_elmx,elem,,num,max
  *set,_mskv
  *dim,_mskv,,_elmx
  *set,_area
  *dim,_area,,_elmx
  *set,_xrea
  *dim,_xrea,,_elmx
   csys,1
   *vget,_mskv(1),elem,,esel       !  12/15/94
   *voper,_mskv(1),_mskv(1),ge,0       !  12/15/94
   *vmask,_mskv(1)       !  12/15/94
   *vget,_area(1),elem,1,geom       !  12/15/94
   *vmask,_mskv(1)       !  12/15/94
   *vget,_xrea(1),elem,1,cent,x     !  12/15/94
   *vmask,_mskv(1)
   *voper,_xrea(1),_xrea(1),mult,_area(1)
   *vmask,_mskv(1)       !  12/15/94
   *vscfun,_tarea,sum,_area(1)       !  12/15/94
   *vmask,_mskv(1)
   *vscfun,_xtarea,sum,_xrea(1)
   _xct=_xtarea/_tarea

  /com, output variables
  t_area=_tarea
  x_xct=_xct


  _elmx=  $  *set,_mskv  $  *set,_area  $  *set,_xrea  $  _xct= $  _tarea=
  _xtrea=

/out