Changeset 5373


Ignore:
Timestamp:
May 28, 2008, 3:26:37 PM (17 years ago)
Author:
jakeman
Message:

forgot file in last commit

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/abstract_2d_finite_volumes/generic_boundary_conditions.py

    r5221 r5373  
    190190
    191191    def __init__(self, filename, domain, time_thinning=1,
    192                  use_cache=False, verbose=False):
     192                 use_cache=False, verbose=False, boundary_polygon=None):
    193193        import time
    194194        from Numeric import array, zeros, Float
     
    210210        self.midpoint_coordinates = zeros( (len(domain.boundary), 2), Float)
    211211        boundary_keys = domain.boundary.keys()
    212 
    213212
    214213        xllcorner = domain.geo_reference.get_xllcorner()
     
    243242            self.boundary_indices[(vol_id, edge_id)] = i
    244243
    245 
    246244        if verbose: print 'Initialise file_function'
    247245        self.F = file_function(filename,
     
    251249                               time_thinning=time_thinning,
    252250                               use_cache=use_cache,
    253                                verbose=verbose)
     251                               verbose=verbose,
     252                               boundary_polygon=boundary_polygon)
     253
    254254        self.domain = domain
    255255
     
    277277        # Test that file function can be called
    278278        q = self.F(0, point_id=0)
    279 
    280279        d = len(domain.conserved_quantities)
    281280        msg = 'Values specified in file %s must be ' %filename
     
    294293
    295294        t = self.domain.time
    296 
    297295        if vol_id is not None and edge_id is not None:
    298296            i = self.boundary_indices[ vol_id, edge_id ]
    299297            res = self.F(t, point_id = i)
    300 
    301298            if res == NAN:
    302299                x,y=self.midpoint_coordinates[i,:]
Note: See TracChangeset for help on using the changeset viewer.