Ignore:
Timestamp:
Sep 30, 2005, 4:07:56 PM (19 years ago)
Author:
ole
Message:

Comments

File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/pyvolution/util.py

    r1854 r1859  
    164164    #In fact, this is where origin should be converted to that of domain
    165165    #Also, check that file covers domain fully.
     166    #If we use the suggested Point_set class for interpolation points
     167    #here it would be easier
     168
     169    #Take into account:
     170    #- domain's georef
     171    #- sww file's georef
     172    #- interpolation points georef
     173
    166174
    167175    assert type(filename) == type(''),\
     
    285293        xllcorner = fid.xllcorner[0]
    286294        yllcorner = fid.yllcorner[0]
     295        zone = fid.zone[0]       
    287296
    288297        x = fid.variables['x'][:]
     
    315324            time = time - domain.starttime + starttime
    316325
    317 
     326        #FIXME Use method in geo to reconcile
     327        #if spatial:
     328        #assert domain.geo_reference.xllcorner == xllcorner
     329        #assert domain.geo_reference.yllcorner == yllcorner
     330        #assert domain.geo_reference.zone == zone       
     331       
    318332    if verbose:
    319333        print 'File_function data obtained from: %s' %filename
Note: See TracChangeset for help on using the changeset viewer.