Changeset 7022


Ignore:
Timestamp:
May 12, 2009, 3:17:41 PM (11 years ago)
Author:
jgriffin
Message:
 
File:
1 edited

Legend:

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

    r7002 r7022  
    8787# Number of boundary segments
    8888num_ocean_segments = len(event_sts) - 1
     89print num_ocean_segments
    8990# Number of landward_boundary points
    9091num_land_points = file_length(project.landward_boundary)
    91 
     92print num_land_points
    9293# Boundary tags refer to project.landward_boundary
    9394# 4 points equals 5 segments start at N
    94 boundary_tags={'back': [num_ocean_segments+1],
     95boundary_tags={'back': range(num_ocean_segments+2,
     96                             num_ocean_segments+5),
    9597               'side': [num_ocean_segments,
    96                         num_ocean_segments+2,
    97                         num_ocean_segments+3,
    98                         num_ocean_segments+num_land_points],
     98                        num_ocean_segments+1,
     99                        num_ocean_segments+5,
     100                        num_ocean_segments+6,
     101                        num_ocean_segments+7,
     102                        num_ocean_segments+8],
     103                        #num_ocean_segments+9],
    99104               'ocean': range(num_ocean_segments)}
     105##boundary_tags = {'back':[3],
     106##                 'side':[2,4],
     107##                 'ocean':[0,1,5]}
     108
     109print 'boundary tags', boundary_tags
    100110
    101111# Build mesh and domain
     
    105115                                    interior_regions=project.interior_regions,
    106116                                    mesh_filename=project.meshes,
    107                                     use_cache=True,
     117                                    use_cache=False,
    108118                                    verbose=True)
    109119print domain.statistics()
Note: See TracChangeset for help on using the changeset viewer.