Changeset 4063 for anuga_work/production


Ignore:
Timestamp:
Dec 6, 2006, 5:35:00 PM (18 years ago)
Author:
sexton
Message:

updates for nsw slide modelling and cairns demo

Location:
anuga_work/production
Files:
14 edited

Legend:

Unmodified
Added
Removed

  • anuga_work/production/broome_2006/export_results.py

    r3944 r4063  
    33
    44from anuga.shallow_water.data_manager import sww2dem
    5 from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
    65from os import sep
    76
    8 time_dir = '20061101_024119'
     7time_dir = '20061113_040953'
    98
    109directory = project.outputdir
    … …  
    3736sww2dem(name, basename_out = outname,
    3837            quantity = quantityname,
    39             cellsize = 100,      # would prefer this at 25
     38            cellsize = 25,      # would prefer this at 25
    4039            # define region for viz purposes
    4140            easting_min = project.eastingmin,

  • anuga_work/production/broome_2006/run_broome.py

    r3972 r4063  
    4343myid = 0
    4444numprocs = 1
    45 start_screen_catcher(project.outputtimedir, myid, numprocs)
     45#start_screen_catcher(project.outputtimedir, myid, numprocs)
    4646
    4747print 'USER:    ', project.user
    … …  
    104104G.export_points_file(project.combined_dem_name + '.pts')
    105105G.export_points_file(project.combined_dem_name + '.xya')
     106
    106107#----------------------------------------------------------------------------
    107108# Create the triangular mesh based on overall clipping polygon with a tagged

  • anuga_work/production/dampier_2006/project.py

    r4049 r4063  
    202202n_max_area = 7725000
    203203
     204e_min_area_mom = 442200
     205e_max_area_mom = 497028
     206n_min_area_mom = 7720095
     207n_max_area_mom = 7750500
     208
    204209poly_facility = read_polygon(polygons_dir+'facility.csv')
    205210

  • anuga_work/production/hobart_2006/run_hobart_clipdata_refine.py

    r3981 r4063  
    2626from anuga.abstract_2d_finite_volumes.combine_pts import combine_rectangular_points_files
    2727from anuga.geospatial_data.geospatial_data import *
    28 from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
    2928
    3029# Application specific imports
    … …  
    4241copy (project.codedir + 'run_hobart_clipdata_refine.py', project.outputtimedir + 'run_hobart_clipdata_refine.py')
    4342print'output dir', project.outputtimedir
    44 
    45 #normal screen output is stored in
    46 screen_output_name = project.outputtimedir + "screen_output.txt"
    47 screen_error_name = project.outputtimedir + "screen_error.txt"
    48 
    49 #used to catch screen output to file
    50 sys.stdout = Screen_Catcher(screen_output_name)
    51 #sys.stderr = Screen_Catcher(screen_output_name)
    52 sys.stderr = Screen_Catcher(screen_error_name)
    5343
    5444print 'USER:    ', project.user
    … …  
    123113# use 75 for onshore components (12.5m DEM)
    124114hobart_res = 7500
    125 bathy_res = 25000
    126 refine_res = 250
     115bathy_res = 50000
     116refine_res = 500
     117refine_res2 = 250
    127118interior_regions = [[project.poly_site13, refine_res],
    128119                    [project.poly_kingston, refine_res],
    … …  
    140131                              'e9': [9], 'e10': [10], 'e11': [11],
    141132                              'e12': [12], 'e13': [13], 'e14': [14]},
    142            'maximum_triangle_area': 700000,
     133           'maximum_triangle_area': 750000,
    143134           'filename': meshname,           
    144135           'interior_regions': interior_regions},

  • anuga_work/production/hobart_2006/run_hobart_test.py

    r3721 r4063  
    4040    mkdir (project.outputtimedir)
    4141copy (project.codedirname, project.outputtimedir + project.codename)
    42 copy (project.codedir + 'run_hobart.py', project.outputtimedir + 'run_hobart.py')
     42copy (project.codedir + 'run_hobart_test.py', project.outputtimedir + 'run_hobart_test.py')
    4343print'output dir', project.outputtimedir
    4444
    … …  
    157157
    158158# use 75 for onshore components (12.5m DEM)
    159 hobart_res = 50000
     159hobart_res = 2500
    160160interior_regions = [[project.poly_hobart1, hobart_res],
    161161                    [project.poly_hobart2, hobart_res],
    … …  
    172172                              'e9': [9], 'e10': [10], 'e11': [11],
    173173                              'e12': [12], 'e13': [13], 'e14': [14]},
    174            'maximum_triangle_area': 2500000,
     174           'maximum_triangle_area': 750000,
    175175           'filename': meshname,           
    176176           'interior_regions': interior_regions},
    … …  
    272272for t in domain.evolve(yieldstep = 240, finaltime = 7200):
    273273    domain.write_time()
    274     domain.write_boundary_statistics(tags = 'bottom')     
    275 
    276 for t in domain.evolve(yieldstep = 120, finaltime = 12600
     274    domain.write_boundary_statistics(tags = 'e13')     
     275
     276for t in domain.evolve(yieldstep = 30, finaltime = 16000
    277277                       ,skip_initial_step = True):
    278278    domain.write_time()
    279     domain.write_boundary_statistics(tags = 'bottom')     
     279    domain.write_boundary_statistics(tags = 'e13')     
     280
     281for t in domain.evolve(yieldstep = 240, finaltime = 25000
     282                       ,skip_initial_step = True):
     283    domain.write_time()
    280284   
    281285print 'That took %.2f seconds' %(time.time()-t0)

  • anuga_work/production/hobart_2006/run_hobart_usepts.py

    r3795 r4063  
    7878# clip 50m pts based on interior regions - want 50m data OUTSIDE of these polygons
    7979clip_regions = [project.poly_hobart1, project.poly_hobart2, \
    80                 project.poly_hobart3, project.poly_hobart4]
     80                project.poly_hobart3]
    8181
    8282# set up initial value as Geospatial data object
    … …  
    9696
    9797# creates pts file for onshore DEM - 25 and make a Geospatial data object
    98 dem2pts(onshore_dem_name_offshore_25,
     98dem2pts(onshore_offshore_dem_name_25,
    9999        easting_min=project.eastingmin25_3,
    100100        easting_max=project.eastingmax25_3,
    … …  
    147147# use 75 for onshore components (12.5m DEM)
    148148hobart_res = 7500
    149 bathy_res = 50000
     149bathy_res = 25000
    150150interior_regions = [[project.poly_hobart1, hobart_res],
    151151                    [project.poly_hobart2, hobart_res],
    … …  
    163163                              'e9': [9], 'e10': [10], 'e11': [11],
    164164                              'e12': [12], 'e13': [13], 'e14': [14],
    165                               'e15': [15]},
    166            'maximum_triangle_area': 500000,
     165                              },
     166           'maximum_triangle_area': 250000,
    167167           'filename': meshname,           
    168168           'interior_regions': interior_regions},
    … …  
    215215domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd,
    216216                      'e5': Bd,  'e6': Bd, 'e7': Bd, 'e8': Bd, 'e9': Bd,
    217                       'e10': Bd, 'e11': Bd, 'e12': Bf, 'e13': Bf, 'e14': Bf,
    218                       'e15': Bf} )
     217                      'e10': Bd, 'e11': Bf, 'e12': Bf, 'e13': Bf, 'e14': Bf,
     218                      } )
    219219
    220220#-------------------------------------------------------------------------------                                 
    … …  
    228228    domain.write_boundary_statistics(tags = 'e14')     
    229229
    230 for t in domain.evolve(yieldstep = 30, finaltime = 9000
     230for t in domain.evolve(yieldstep = 30, finaltime = 15000
    231231                       ,skip_initial_step = True):
    232232    domain.write_time()
    233233    domain.write_boundary_statistics(tags = 'e14')     
    234234
    235 for t in domain.evolve(yieldstep = 240, finaltime = 15000
     235for t in domain.evolve(yieldstep = 240, finaltime = 20000
    236236                       ,skip_initial_step = True):
    237237    domain.write_time()

  • anuga_work/production/newcastle_2006/project_slide.py

    r4058 r4063  
    6262
    6363# Necessary if using point datasets, rather than grid
    64 on_offshore_dem_name = datadir + on_offshore_name
     64nsw_dem_name = datadir + nsw100
    6565offshore_dem_name2 = datadir + offshore_name2
    6666offshore_dem_name3 = datadir + offshore_name3
    … …  
    7878
    7979# bounding polygon for study area
    80 polyAll = read_polygon(polygondir+'extentC.csv')
     80polyAll = read_polygon(polygondir+'extentC2.csv')
    8181
    8282print 'Area of bounding polygon', polygon_area(polyAll)/1000000.0
    … …  
    8787
    8888# newcastle digitized polygons
    89 poly_newcastle1 = read_polygon(polygondir+'.csv')
    90 poly_newcastle2 = read_polygon(polygondir+'.csv')
    91 poly_newcastle3 = read_polygon(polygondir+'.csv')
     89poly_local = read_polygon(polygondir+'local2.csv')
     90poly_newcastle = read_polygon(polygondir+'newcastle.csv')
    9291
    93 print 'Area of local polygon', polygon_area(poly_newcastle1)/1000000.0
    94 print 'Area of close polygon', polygon_area(poly_newcastle2)/1000000.0
    95 print 'Area of coastal polygon', polygon_area(poly_newcastle3)/1000000.0
     92print 'Area of local polygon', polygon_area(poly_local)/1000000.0
     93print 'Area of CBD polygon', polygon_area(poly_newcastle)/1000000.0
    9694
     95#plot_polygons([polyAll,poly_local,poly_newcastle],'fig',verbose=True)
    9796###################################################################
    9897# Clipping regions for export to asc and regions for clipping data
    … …  
    10099
    101100# clipping for nsw grid for conversion
    102 eastingmin_nsw = 546000
    103 eastingmax_nsw = 350390
     101eastingmin_nsw = 350390
     102eastingmax_nsw = 546000
    104103northingmin_nsw = 6253970
    105104northingmax_nsw = 6424600
    106105
     106poly_surveyclip = read_polygon(polygondir+'surveyclip.csv')
     107
    107108# exporting asc grid
    108 eastingmin =
    109 eastingmax =
    110 northingmin =
    111 northingmax =
     109#eastingmin =
     110#eastingmax =
     111#northingmin =
     112#northingmax =
    112113
    113114###################################################################

  • anuga_work/production/newcastle_2006/run_newcastle_slide.py

    r4058 r4063  
    22
    33Source data such as elevation and boundary data is assumed to be available in
    4 directories specified by project.py
    5 The output sww file is stored in project.outputtimedir
     4directories specified by project_slide.py
     5The output sww file is stored in project_slide.outputtimedir
    66
    7 The scenario is defined by a triangular mesh created from project.polygon,
     7The scenario is defined by a triangular mesh created from project_slide.polygon,
    88the elevation data and a tsunami wave generated by s submarine mass failure.
    99
    … …  
    3838
    3939# creates copy of code in output dir
    40 copy_code_files(project.outputtimedir,__file__,dirname(project.__file__)+sep+ project.__name__+'.py' )
     40copy_code_files(project_slide.outputtimedir,__file__,dirname(project_slide.__file__)+sep+ project_slide.__name__+'.py' )
    4141myid = 0
    4242numprocs = 1
    43 start_screen_catcher(project.outputtimedir, myid, numprocs)
     43start_screen_catcher(project_slide.outputtimedir, myid, numprocs)
    4444
    45 print 'USER:    ', project.user
     45print 'USER:    ', project_slide.user
    4646
    4747#-------------------------------------------------------------------------------
    … …  
    7979print 'add'
    8080G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\
    81     G12.clip(Geospatial_data(project_slide.polyAll))
    82     (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip))
     81    G12.clip(Geospatial_data(project_slide.polyAll)) +\
     82    (G4.clip(Geospatial_data(project_slide.polyAll)).clip_outside(Geospatial_data(project_slide.poly_surveyclip)))
    8383print 'export points'
    8484G.export_points_file(project_slide.combined_dem_name + '.pts')
    … …  
    8888#----------------------------------------------------------------------------
    8989# Create the triangular mesh based on overall clipping polygon with a tagged
    90 # boundary and interior regions defined in project.py along with
     90# boundary and interior regions defined in project_slide.py along with
    9191# resolutions (maximal area of per triangle) for each polygon
    9292#-------------------------------------------------------------------------------
    … …  
    9494from anuga.pmesh.mesh_interface import create_mesh_from_regions
    9595remainder_res = 500000
    96 local_res = 25000
    97 newcastle_res = 5000
    98 coast_res = 500
    99 interior_regions = [[project.poly_newcastle1, local_res],
    100                     [project.poly_newcastle2, newcastle_res],
    101                     [project.poly_newcastle3, coast_res]]
     96local_res = 50000
     97newcastle_res = 1000
     98interior_regions = [[project_slide.poly_local, local_res],
     99                    [project_slide.poly_newcastle, newcastle_res]]
    102100
    103101from caching import cache
    104102_ = cache(create_mesh_from_regions,
    105           project.polyAll,
     103          project_slide.polyAll,
    106104           {'boundary_tags': {'e0': [0], 'e1': [1], 'e2': [2],
    107                               'e3': [3], 'e4':[4], 'e5': [5],
    108                               'e6': [6]},
     105                              'e3': [3], 'e4':[4]},
    109106           'maximum_triangle_area': remainder_res,
    110107           'filename': meshname,
    … …  
    122119print domain.statistics()
    123120
    124 domain.set_name(project.basename)
    125 domain.set_datadir(project.outputtimedir)
     121domain.set_name(project_slide.basename)
     122domain.set_datadir(project_slide.outputtimedir)
    126123domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
    127124domain.set_minimum_storable_height(0.01)
    … …  
    135132domain.set_quantity('friction', 0.0)
    136133domain.set_quantity('elevation',
    137                     filename = project.combined_dem_name + '.pts',
     134                    filename = project_slide.combined_dem_name + '.pts',
    138135                    use_cache = True,
    139136                    verbose = True,
    … …  
    146143from smf import slide_tsunami
    147144
    148 tsunami_source = slide_tsunami(length=30000.0,
    149                                depth=400.0,
    150                                slope=6.0,
    151                                thickness=176.0,
    152                                radius=3330,
    153                                dphi=0.23,
    154                                x0=project.slump_origin[0],
    155                                y0=project.slump_origin[1],
    156                                alpha=0.0,
     145tsunami_source = slide_tsunami(length=project_slide.bulli_length,
     146                               width=project_slide.bulli_width,
     147                               depth=project_slide.bulli_depth,
     148                               slope=project_slide.bulli_slope,
     149                               thickness=project_slide.bulli_thickness,
     150                               x0=project_slide.slide_origin_c[0],
     151                               y0=project_slide.slide_origin_c[1],
     152                               alpha=project_slide.bulli_alpha,
    157153                               domain=domain)
    158154
    … …  
    165161Bd = Dirichlet_boundary([tide,0,0])
    166162
    167 domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd,
    168                       'e5': Bd,  'e6': Bd} )
     163domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd} )
    169164
    170165

  • anuga_work/production/onslow_2006/export_results.py

    r3650 r4063  
    5252
    5353if which_var == 4:  # Elevation
    54     outname = name + '_elevation'
     54    outname = name + '_elevation_25'
    5555    quantityname = 'elevation'  #Elevation
    5656
    … …  
    5959sww2dem(name, basename_out = outname,
    6060            quantity = quantityname,
    61             cellsize = 20,       # Trevor would like this at 25
     61            cellsize = 25,       
    6262            # define region for viz purposes
    63             easting_min = project.e_min_area,
    64             easting_max = project.e_max_area,
    65             northing_min = project.n_min_area,
    66             northing_max = project.n_max_area,       
     63            #easting_min = project.e_min_area,
     64            #easting_max = project.e_max_area,
     65            #northing_min = project.n_min_area,
     66            #northing_max = project.n_max_area,       
    6767            reduction = max, #this is because we want max quantityname
    6868            verbose = True,

  • anuga_work/production/onslow_2006/project.py

    r3788 r4063  
    166166polygons = [polyAll, export_region]
    167167figname = 'checking.png'
    168 from anuga.utilities.polygon import plot_polygons
    169 plot_polygons(polygons, figname, verbose = False)
    170 print figname
     168#from anuga.utilities.polygon import plot_polygons
     169#plot_polygons(polygons, figname, verbose = False)
     170#print figname
    171171#Interior region - Onslow town
    172172

  • anuga_work/production/pt_hedland_2006/export_results.py

    r3514 r4063  
    22import sys
    33
    4 from anuga.pyvolution.data_manager import sww2dem
    5 from anuga.pyvolution.ermapper_grids import read_ermapper_grid
    6 from anuga.pyvolution.util import Screen_Catcher
     4from anuga.shallow_water.data_manager import sww2dem
     5#from anuga.pyvolution.ermapper_grids import read_ermapper_grid
     6from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
    77from os import sep
    88
    9 time_dir = '20060706_235036' #test with coarse grid
    10 #time_dir = '20060707_001859' #MSL DLI data
     9#time_dir = '20060706_235036' #test with coarse grid
     10time_dir = '20060707_001859' #MSL DLI data
    1111#time_dir = '20060707_003301' #HAT DLI data
    1212#time_dir = '20060707_003424' #LAT DLI data
    … …  
    4949
    5050if which_var == 4:  # Elevation
    51     outname = name + '_elevation'
     51    outname = name + '_elevation_50'
    5252    quantityname = 'elevation'  #Elevation
    5353
    … …  
    5656sww2dem(name, basename_out = outname,
    5757            quantity = quantityname,
    58             cellsize = 20,       # Trevor would like this at 25
     58            cellsize = 50,       # Trevor would like this at 25
    5959            # define region for viz purposes
    60             easting_min = project.e_min_area,
    61             easting_max = project.e_max_area,
    62             northing_min = project.n_min_area,
    63             northing_max = project.n_max_area,       
     60            #easting_min = project.e_min_area,
     61            #easting_max = project.e_max_area,
     62            #northing_min = project.n_min_area,
     63            #northing_max = project.n_max_area,       
    6464            reduction = max, #this is because we want max quantityname
    6565            verbose = True,

  • anuga_work/production/sydney_2006/run_sydney_slide.py

    r4058 r4063  
    9898print 'add'
    9999G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\
    100     G12.clip(Geospatial_data(project_slide.polyAll))
     100    G12.clip(Geospatial_data(project_slide.polyAll)) +\
    101101    G2.clip(Geospatial_data(project_slide.poly_25mclip)) +\
    102102    G3.clip(Geospatial_data(project_slide.poly_origsyd)) +\
    … …  
    113113
    114114from anuga.pmesh.mesh_interface import create_mesh_from_regions
    115 remainder_res = 125000.
    116 local_res = 30000.
     115remainder_res = 150000.
     116local_res = 50000.
    117117coast_res = 500.
    118118interior_regions = [[project_slide.poly_syd1, local_res],
    … …  
    193193t0 = time.time()
    194194
    195 for t in domain.evolve(yieldstep = 30, finaltime = 480):
     195for t in domain.evolve(yieldstep = 30, finaltime = 5000):
    196196    domain.write_time()
    197197    domain.write_boundary_statistics(tags = 'e2')

  • anuga_work/production/wollongong_2006/project_slide.py

    r4058 r4063  
    6262
    6363# Necessary if using point datasets, rather than grid
    64 on_offshore_dem_name = datadir + on_offshore_name
     64on_offshore10_dem_name = datadir + on_offshore10_name
     65nsw_dem_name = datadir + nsw100
    6566offshore_dem_name1 = datadir + offshore_name1
    6667offshore_dem_name4 = datadir + offshore_name4
    … …  
    8788# areaA digitized polygons
    8889poly_local = read_polygon(polygondir+'local.csv')
    89 poly_gong = read_polygon(polygondir+'gong.csv')
     90poly_gong = read_polygon(polygondir+'gong2.csv')
    9091
    9192print 'Area of local polygon', polygon_area(poly_local)/1000000.0

  • anuga_work/production/wollongong_2006/run_gong_slide.py

    r4058 r4063  
    22
    33Source data such as elevation and boundary data is assumed to be available in
    4 directories specified by project.py
    5 The output sww file is stored in project.outputtimedir
     4directories specified by project_slide.py
     5The output sww file is stored in project_slide.outputtimedir
    66
    77The scenario is defined by a triangular mesh created from project.polygon,
    … …  
    3838
    3939# creates copy of code in output dir
    40 copy_code_files(project.outputtimedir,__file__,dirname(project.__file__)+sep+ project.__name__+'.py' )
     40copy_code_files(project_slide.outputtimedir,__file__,dirname(project_slide.__file__)+sep+ project_slide.__name__+'.py' )
    4141myid = 0
    4242numprocs = 1
    43 start_screen_catcher(project.outputtimedir, myid, numprocs)
     43start_screen_catcher(project_slide.outputtimedir, myid, numprocs)
    4444
    45 print 'USER:    ', project.user
     45print 'USER:    ', project_slide.user
    4646
    4747#-------------------------------------------------------------------------------
    … …  
    8181G4 = Geospatial_data(file_name = project_slide.nsw_dem_name + '.pts')
    8282print 'add'
    83 G = G11.clip(Geospatial(project_slide.poly_surveyclip)) +\
     83G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\
    8484    G12.clip(Geospatial_data(project_slide.polyAll)) +\
    8585    G2.clip(Geospatial_data(project_slide.poly_10mclip)) +\
    … …  
    9696
    9797from anuga.pmesh.mesh_interface import create_mesh_from_regions
    98 remainder_res = 500000
     98remainder_res = 100000
    9999local_res = 25000
    100 gong_res = 5000
    101 coast_res = 500
    102 interior_regions = [[project.poly_gong1, local_res],
    103                     [project.poly_gong2, gong_res],
    104                     [project.poly_gong3, coast_res]]
     100gong_res = 500
     101interior_regions = [[project_slide.poly_local, local_res],
     102                    [project_slide.poly_gong, gong_res]]
    105103
    106104from caching import cache
    107105_ = cache(create_mesh_from_regions,
    108           project.polyAll,
     106          project_slide.polyAll,
    109107           {'boundary_tags': {'e0': [0], 'e1': [1], 'e2': [2],
    110                               'e3': [3], 'e4':[4], 'e5': [5],
    111                               'e6': [6]},
     108                              'e3': [3], 'e4':[4]},
    112109           'maximum_triangle_area': remainder_res,
    113110           'filename': meshname,
    … …  
    125122print domain.statistics()
    126123
    127 domain.set_name(project.basename)
    128 domain.set_datadir(project.outputtimedir)
     124domain.set_name(project_slide.basename)
     125domain.set_datadir(project_slide.outputtimedir)
    129126domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
    130127domain.set_minimum_storable_height(0.01)
    … …  
    138135domain.set_quantity('friction', 0.0)
    139136domain.set_quantity('elevation',
    140                     filename = project.combined_dem_name + '.pts',
     137                    filename = project_slide.combined_dem_name + '.pts',
    141138                    use_cache = True,
    142139                    verbose = True,
    … …  
    149146from smf import slide_tsunami
    150147
    151 tsunami_source = slide_tsunami(length=30000.0,
    152                                depth=400.0,
    153                                slope=6.0,
    154                                thickness=176.0,
    155                                radius=3330,
    156                                dphi=0.23,
    157                                x0=project.slump_origin[0],
    158                                y0=project.slump_origin[1],
    159                                alpha=0.0,
     148tsunami_source = slide_tsunami(length=project_slide.bulli_length,
     149                               width=project_slide.bulli_width,
     150                               depth=project_slide.bulli_depth,
     151                               slope=project_slide.bulli_slope,
     152                               thickness=project_slide.bulli_thickness,
     153                               x0=project_slide.slide_origin_a[0],
     154                               y0=project_slide.slide_origin_a[1],
     155                               alpha=project_slide.bulli_alpha,
    160156                               domain=domain)
    161157
    … …  
    168164Bd = Dirichlet_boundary([tide,0,0])
    169165
    170 domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd,
    171                       'e5': Bd,  'e6': Bd} )
     166domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd} )
    172167
    173168
Note: See TracChangeset for help on using the changeset viewer.