Changeset 3808


Ignore:
Timestamp:
Oct 17, 2006, 12:05:03 PM (17 years ago)
Author:
nick
Message:

testing urs2sww

Location:
anuga_work/production/dampier_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/dampier_2006/project.py

    r3802 r3808  
    1515# Location and naming of scenario data
    1616state = 'western_australia'
    17 #scenario_dir_name = 'dampier_tsunami_scenario_2006'
    18 scenario_dir_name = 'karratha_tsunami_scenario_2005' # Tmp location
     17scenario_dir_name = 'dampier_tsunami_scenario_2006'
     18#scenario_dir_name = 'karratha_tsunami_scenario_2005' # Tmp location
    1919
    2020basename = 'elevation50m'
    21 boundary_basename = 'SU-AU'
     21boundary_basename = 'dampier'
    2222
    2323
     
    4343gaugedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'gauges'+sep
    4444polygondir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'polygons'+sep
    45 boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep
     45boundarydir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'boundaries'+sep+'urs'+sep+'test'+sep
    4646outputdir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'output'+sep
    4747tidedir = home+sep+state+sep+scenario_dir_name+sep+'anuga'+sep+'tide_data'+sep
     
    5151meshname = meshdir + basename
    5252outputname = outputtimedir + basename  #Used by post processing
    53 #boundaryname = outputtimedir + basename  #Used by post processing
    54 
    55 
     53boundaryname = boundarydir + boundary_basename  #Used by post processing
    5654
    5755
     
    6058refzone = 50
    6159south = degminsec2decimal_degrees(-20,55,0)
    62 north = degminsec2decimal_degrees(-20,15,0)
     60north = degminsec2decimal_degrees(-19,15,0)
    6361west = degminsec2decimal_degrees(116,17,0)
    64 east = degminsec2decimal_degrees(117,10,0)
     62east = degminsec2decimal_degrees(118,10,0)
    6563
    6664p0 = [south, degminsec2decimal_degrees(116,32,0)]
  • anuga_work/production/dampier_2006/run_dampier.py

    r3802 r3808  
    4747meshname = project.meshname+'.msh'
    4848source_dir = project.boundarydir
     49boundary_file = project.boundaryname
     50
     51from anuga.shallow_water.data_manager import urs2sww
     52
     53urs2sww(boundary_file, verbose='true',
     54       minlat=project.south,
     55       maxlat=project.north,
     56       minlon=project.west,
     57       maxlon=project.east)
     58
     59import sys; sys.exit()
     60Bf = File_boundary(source_dir + project.boundary_basename + '.sww',
     61                   domain, verbose = True)
     62                   
     63Br = Reflective_boundary(domain)
     64Bd = Dirichlet_boundary([tide,0,0])
     65domain.set_boundary({'back': Br,
     66                     'side': Bd,
     67                     'ocean': Bf})
    4968
    5069
     
    7594
    7695
     96#---------------------------
     97# this is check that no interior polygon is outside the bounding poly
     98#------------------------------
     99count = 0
     100for i in range(len(interior_regions)):
     101    region = interior_regions[i]
     102    interior_polygon = region[0]
     103    if len(inside_polygon(interior_polygon, bounding_polygon,
     104                   closed = True, verbose = False)) <> len(interior_polygon):
     105        print 'WARNING: interior polygon %d is outside bounding polygon' %(i)
     106        count += 1
     107
     108if count == 0:
     109    print 'interior regions OK'
     110else:
     111    print 'check out your interior polygons'
     112    print 'check %s in production directory' %figname
     113    import sys; sys.exit()
     114
     115#-------------------------------------
    77116
    78117print 'start create mesh from regions'
     
    119158print 'Available boundary tags', domain.get_boundary_tags()
    120159
     160
     161
     162urs2sww(boundary_file,
     163       minlat=project.south,
     164       maxlat=project.north,
     165       minlon=project.west,
     166       maxlon=project.east)
     167
    121168Bf = File_boundary(source_dir + project.boundary_basename + '.sww',
    122169                   domain, verbose = True)
Note: See TracChangeset for help on using the changeset viewer.