(vf, dim)
| 524 | |
| 525 | # Returns effective radius |
| 526 | def GetReff(vf, dim): |
| 527 | x1, y1, z1, hx, hy, hz = GetGeom(vf.shape) |
| 528 | # cell volume |
| 529 | vc = hx * hy * hz |
| 530 | # integral of vf |
| 531 | v = vf.sum() * vc |
| 532 | return GetReff2(v / hz) if dim == 2 else GetReff3(v) |
| 533 | |
| 534 | |
| 535 | # Returns path template from sample path |