Changeset 5369


Ignore:
Timestamp:
May 27, 2008, 4:07:41 PM (16 years ago)
Author:
duncan
Message:

comments

Location:
anuga_core/source/anuga/shallow_water
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/shallow_water/data_manager.py

    r5358 r5369  
    44384438    Given the info to replicate the URS grid and a polygon output
    44394439    a file that specifies the cloud of boundary points for URS.
     4440
     4441    This creates a .urs file.  This is in the format used by URS;
     4442    1st line is the number of points,
     4443    each line after represents a point,in lats and longs.
    44404444   
    44414445    Note: The polygon cannot cross zones or hemispheres.
     4446
     4447    A work-a-round for different zones or hemispheres is to run this twice,
     4448    once for each zone, and then combine the output.
    44424449   
    44434450    file_name - name of the urs file produced for David.
  • anuga_core/source/anuga/shallow_water/test_data_manager.py

    r5361 r5369  
    66406640        #print "seg_lat_long",  poly_lat_long
    66416641       
    6642         #geo=URS_points_needed_to_file('test_example_poly3', poly_lat_long,
     6642      #   geo=URS_points_needed_to_file('test_example_poly3', poly_lat_long,
     6643#                                   ZONE,
     6644#                                   LL_LAT, LL_LONG,
     6645#                                   GRID_SPACING,
     6646#                                   LAT_AMOUNT, LONG_AMOUNT,
     6647#                                   isSouthernHemisphere=False,
     6648#                                   export_csv=True,
     6649#                                   verbose=self.verbose)
     6650       
    66436651        geo=URS_points_needed(poly_lat_long,
    66446652                                  ZONE,
     
    66476655                                  LAT_AMOUNT, LONG_AMOUNT,
    66486656                                  isSouthHemisphere=False,
    6649                                   verbose=self.verbose)
     6657                                  verbose=self.verbose)
     6658       
    66506659        results = ImmutableSet(geo.get_data_points(as_lat_long=True,
    66516660                                  isSouthHemisphere=False))
Note: See TracChangeset for help on using the changeset viewer.