User manual Sydney example crashes

[duncan]

When trying to run the Sydney example I got a crash. Does this work for anyone else? Is there a known reason why it doesn't work?

The crash;

Using attribute elevation from file m:\inundation\sydney_tsunami_examp
le_2006\topographies\sydneytopo.pts
Available attributes: ['elevation']
Caching: looking for cached files c:\.python_cache\fit_to_mesh[-132031
2002]_{Result,Args,Admin}.z
+----------------------------------------------------------
| Tue Aug 29 10:38:42 2006. Evaluating function fit_to_mesh
+----------------------------------------------------------
| Arguments:    (Array: (14961, 2), Array: (29311, 3), Array: (989056,
 2), Array: (989056,))
| Keyword Args: {'alpha': None, 'mesh_origin': (-1, 316000.0, 6238000.
0), 'verbose': True, 'data_origin': (56, 332090, 6246250)}
| Reason:       No cached result
+----------------------------------------------------------
Traceback (most recent call last):
  File "I:\documentation\user_manual\examples\runsydney.py", line 119,
 in ?
    verbose = True)
  File "i:\anuga_core\source\anuga\pyvolution\domain.py", line 301, in
 set_quantity
    self.quantities[name].set_values(*args, **kwargs)
  File "i:\anuga_core\source\anuga\pyvolution\quantity.py", line 352,
in set_values
    use_cache = use_cache)
  File "i:\anuga_core\source\anuga\pyvolution\quantity.py", line 745,
in set_values_from_file
    use_cache = use_cache)
  File "i:\anuga_core\source\anuga\pyvolution\quantity.py", line 606,
in set_values_from_geospatial_data
    use_cache = use_cache)
  File "i:\anuga_core\source\anuga\pyvolution\quantity.py", line 672,
in set_values_from_points
    compression=False)
  File "I:\anuga_core\source\anuga\caching\caching.py", line 349, in c
ache
    T = apply(func,args,kwargs)
  File "I:\anuga_core\source\anuga\fit_interpolate\fit.py", line 438,
in fit_to_mesh
    alpha=alpha)
  File "I:\anuga_core\source\anuga\fit_interpolate\fit.py", line 93, i
n __init__
    max_vertices_per_cell)
  File "i:\anuga_core\source\anuga\fit_interpolate\general_fit_interpo
late.py", line 83, in __init__
    self.mesh.check_integrity()
  File "I:\anuga_core\source\anuga\pyvolution\neighbour_mesh.py", line
 601, in check_integrity
    a0 = anglediff(v1, v0)
  File "i:\anuga_core\source\anuga\utilities\numerical_tools.py", line
 73, in anglediff
    a1 = angle(v1)
  File "i:\anuga_core\source\anuga\utilities\numerical_tools.py", line
 49, in angle
    theta = acos(p)
ValueError: math domain error

[ole]

Doesn't work for me either - but in my case a file can't be found:

X:\inundation\sandpits\onielsen\anuga_core\documentation\user_manual\examples>python runsydney.py refzone 56 Traceback (most recent call last):

File "runsydney.py", line 28, in ?

import project # Define file names and polygons

File "X:\inundation\sandpits\onielsen\anuga_core\documentation\user_manual\examples\project.py", line 59, in ?

northern_polygon = read_polygon(polygonptsfile0 + '.csv')

File "X:\inundation\sandpits\onielsen\anuga_core\source\anuga\utilities\polygon.py", line 624, in read_polygon

fid = open(filename)

IOError: [Errno 2] No such file or directory: 'X:
inundation
sydney_tsunami_example_2006
polygons
poly2.csv'

[ole]

With INUNDATIONHOME set, it worked fine for me.

[Duncan <duncan.gray@…>]

works for me as well now. Change set #6388 would'v fixed this.

[duncan]

Oops. dyslexic error Changeset [3688] by ole Fixed get_boundary polygon based on information from new unit test. It now works for all examples of meshes provided.