Ignore:
Timestamp:
Mar 13, 2009, 9:02:45 AM (15 years ago)
Author:
myall
Message:

darwin now working; LHI still has problems with boundary

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/australia_ph2/LHI/run_model.py

    r6491 r6502  
    8080# Reading the landward defined points, this incorporates the original clipping
    8181# polygon minus the 100m contour
    82 
    83 # NO LANDWARD BOUNDARY FOR ISLAND MODEL
    8482##landward_boundary = read_polygon(project.landward_boundary)
    8583
     
    8886
    8987# Number of boundary segments
    90 num_ocean_segments = len(event_sts) #- 1
     88num_ocean_segments = len(event_sts) - 1
    9189# Number of landward_boundary points
    9290num_land_points = 0 #file_length(project.landward_boundary)
     
    9492# Boundary tags refer to project.landward_boundary
    9593# 4 points equals 5 segments start at N
    96 boundary_tags={#'back': range(num_ocean_segments+1,
    97                #              num_ocean_segments+num_land_points),
    98                #'side': [num_ocean_segments,
    99                #         num_ocean_segments+num_land_points],
     94boundary_tags={'back': range(num_ocean_segments+1,
     95                             num_ocean_segments+num_land_points),
     96               'side': [num_ocean_segments,
     97                        num_ocean_segments+num_land_points],
    10098               'ocean': range(num_ocean_segments)}
    10199
     
    106104                                    interior_regions=project.interior_regions,
    107105                                    mesh_filename=project.meshes,
    108                                     use_cache=False,
     106                                    use_cache=True,
    109107                                    verbose=True)
    110108print domain.statistics()
     
    127125else:
    128126    IC = 0
    129 domain.set_quantity('stage', IC, use_cache=False, verbose=True)
     127domain.set_quantity('stage', IC, use_cache=True, verbose=True)
    130128domain.set_quantity('friction', project.friction)
    131129domain.set_quantity('elevation',
    132130                    filename=project.combined_elevation+'.pts',
    133                     use_cache=False,
     131                    use_cache=True,
    134132                    verbose=True,
    135133                    alpha=project.alpha)
     
    155153                        default_boundary=Bd,
    156154                        boundary_polygon=bounding_polygon_sts,                   
    157                         use_cache=False,
     155                        use_cache=True,
    158156                        verbose=True)
    159     domain.set_boundary({'ocean': Bf})
     157    domain.set_boundary({'back': Br,
     158                     'side': Bd,
     159                     'ocean': Bf})
    160160   
    161161elif project.wave == 'Tb':
     
    163163                       verbose=True)
    164164
    165     domain.set_boundary({'ocean': Tb})
     165    domain.set_boundary({'back': Br,
     166                         'side': Bd,
     167                         'ocean': Tb})
    166168else:
    167169    print 'No wave specified in project script (Bf or Tb)'
Note: See TracChangeset for help on using the changeset viewer.