Changeset 2293


Ignore:
Timestamp:
Jan 26, 2006, 2:25:09 PM (18 years ago)
Author:
sexton
Message:

Updated files for Sydney tsunami test case (better interior regions, change to slump parameters)

Location:
production/sydney_2006
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • production/sydney_2006/project.py

    r2185 r2293  
    11"""Common filenames and locations for topographic data, meshes and outputs.
     2Also includes origin for slump scenario.
    23"""
    34
     
    910#Making assumptions about the location of scenario data
    1011scenario_dir_name = 'sydney_tsunami_scenario_2006'
    11 basename = 'sydney_2006_100' # get from Neil (DEM or topo data)
     12# original 100m data
     13#basename = 'sydney_2006_100' # get from Neil/Ingo (DEM or topo data)
     14# revised 100m data
     15coarsename = 'bathyland100' # get from Neil/Ingo (DEM or topo data)
     16# revised 25m data
     17#finename = 'bathyland25' # get from Neil/Ingo (DEM or topo data)
     18#finename = 'bathyland100' # get from Neil/Ingo (DEM or topo data)
     19finename = 'bathy_dem25' # get from Neil/Ingo (DEM or topo data) Wed 25 Jan
     20
     21# creating easting and northing max and min for fine data - region of interest
     22eastingmin = 332090
     23eastingmax = 347500
     24northingmin = 6246250
     25northingmax = 6264100
     26#northingmax = 6280100
     27
     28# creating easting and northing max and min for export viz purposes
     29#eminviz = 332000
     30#emaxviz = 350000
     31#nminviz = 6231000
     32#nmaxviz = 6283000
     33
     34# version2
     35eminviz = 318000
     36emaxviz = 351000
     37nminviz = 6231000
     38nmaxviz = 6283000
     39   
     40basename = 'slump'
    1241
    1342if sys.platform == 'win32':
     
    2453
    2554meshname = meshdir + basename
    26 demname = datadir + basename
     55coarsedemname = datadir + coarsename
     56finedemname = datadir + finename
     57combineddemname = datadir + 'sydneytopo'
    2758#boundaryname = boundarydir + boundary_basename
    2859outputname = outputdir + basename  #Used by post processing
    29 
     60#newoutputname = outputdir+'Combined Data v1'+sep+'d400'+sep+basename
     61gauge_filename = outputdir + 'sydney_gauges.xya'
    3062
    3163#Georeferencing
     
    3365
    3466#Origin of existing dem (FIXME: Temporary measure)
    35 mesh_origin = (56, 314036, 6224951)  # input from Neil's data
     67#x0_origin = 314036 #original 100m data
     68#y0_origin = 6224951
     69refzone = 56
     70x0_origin = 314036.58727982 #revised 100m and 25m data
     71y0_origin = 6224951.2960092
     72mesh_origin = (refzone, x0_origin, y0_origin)  # input from Neil's data
    3673
    3774# define clipping polygon
    38 #south = degminsec2decimal_degrees(-20,45,0)
    39 #north = degminsec2decimal_degrees(-20,15,0)
    40 #west = degminsec2decimal_degrees(116,30,0)
    41 #east = degminsec2decimal_degrees(117,0,0)
    42 
    4375south = degminsec2decimal_degrees(-34,05,0)
    44 north = degminsec2decimal_degrees(-33,45,0)
    45 west = degminsec2decimal_degrees(151,0,0)
     76north = degminsec2decimal_degrees(-33,33,0)
     77west = degminsec2decimal_degrees(151,1,0)
    4678east = degminsec2decimal_degrees(151,30,0)
    47 
    4879p0 = [south, west]
    4980p1 = [south, east]
     
    5283   
    5384polygonall = [p0, p1, p2, p3]
    54 refzone = 56
    55 
     85
     86dsouth = degminsec2decimal_degrees(-34,05,0)
     87dnorth = degminsec2decimal_degrees(-33,33,0)
     88dnorth1 = degminsec2decimal_degrees(-33,40,0)
     89dnorth2 = degminsec2decimal_degrees(-33,58,30)
     90dnorth3 = degminsec2decimal_degrees(-33,46,0)
     91dwest = degminsec2decimal_degrees(151,2,20)
     92deast1 = degminsec2decimal_degrees(151,20,0)
     93deast2 = degminsec2decimal_degrees(151,48,0)
     94deast3 = degminsec2decimal_degrees(151,10,0)
     95deast4 = degminsec2decimal_degrees(151,9,0)
     96
     97dp0 = [dsouth, dwest]
     98dp1 = [dsouth, deast1]
     99dp2 = [dnorth2, deast2]
     100dp3 = [dnorth1, deast2]
     101dp4 = [dnorth, deast1]
     102dp5 = [dnorth, deast4]
     103dp6 = [dnorth3, deast3]
     104dp7 = [dnorth3, dwest]
     105dp8 = [dnorth, dwest]
     106   
     107#diffpolygonall = [dp0, dp1, dp2, dp3, dp4, dp5, dp6, dp7]
     108diffpolygonall = [dp0, dp1, dp2, dp3, dp4, dp5, dp6, dp7]
     109# to put chunk back in
     110#diffpolygonall = [dp0, dp1, dp2, dp3, dp4, dp8]
     111
     112#Interior regions - the Harbour - take 2
     113harbour_1x = degminsec2decimal_degrees(-33,51,0)
     114harbour_1y = degminsec2decimal_degrees(151,2,30)
     115harbour_12x = degminsec2decimal_degrees(-33,51,0)
     116harbour_12y = degminsec2decimal_degrees(151,5,0)
     117harbour_13x = degminsec2decimal_degrees(-33,52,15)
     118harbour_13y = degminsec2decimal_degrees(151,5,0)
     119harbour_2x = degminsec2decimal_degrees(-33,53,0)
     120harbour_2y = degminsec2decimal_degrees(151,17,20)
     121harbour_3x = degminsec2decimal_degrees(-33,47,0)
     122harbour_3y = degminsec2decimal_degrees(151,20,30)
     123#harbour_4x = degminsec2decimal_degrees(-33,48,0)
     124harbour_4x = degminsec2decimal_degrees(-33,47,50)
     125harbour_4y = degminsec2decimal_degrees(151,8,10)
     126#harbour_5x = degminsec2decimal_degrees(-33,49,0)
     127harbour_5x = degminsec2decimal_degrees(-33,48,10)
     128harbour_5y = degminsec2decimal_degrees(151,8,0)
     129harbour_6x = degminsec2decimal_degrees(-33,49,0)
     130harbour_6y = degminsec2decimal_degrees(151,2,30)
     131harbour_7x = degminsec2decimal_degrees(-33,34,30)
     132harbour_7y = degminsec2decimal_degrees(151,20,20)
     133harbour_8x = degminsec2decimal_degrees(-33,33,30)
     134harbour_8y = degminsec2decimal_degrees(151,15,0)
     135harbour_9x = degminsec2decimal_degrees(-33,45,30)
     136harbour_9y = degminsec2decimal_degrees(151,15,0)
     137harbour_10x = degminsec2decimal_degrees(-33,45,10)
     138harbour_10y = degminsec2decimal_degrees(151,11,40)
     139harbour_11x = degminsec2decimal_degrees(-33,45,10)
     140harbour_11y = degminsec2decimal_degrees(151,11,40)
     141harbour_14x = degminsec2decimal_degrees(-33,49,10)
     142harbour_14y = degminsec2decimal_degrees(151,11,40)
     143harbour_15x = degminsec2decimal_degrees(-33,48,55)
     144harbour_15y = degminsec2decimal_degrees(151,2,30)
     145
     146k02 = [harbour_1x, harbour_1y]
     147k12 = [harbour_2x, harbour_2y]
     148k22 = [harbour_3x, harbour_3y]
     149k32 = [harbour_4x, harbour_4y]
     150k42 = [harbour_5x, harbour_5y]
     151k52 = [harbour_6x, harbour_6y]
     152k62 = [harbour_7x, harbour_7y]
     153k72 = [harbour_8x, harbour_8y]
     154k82 = [harbour_9x, harbour_9y]
     155k92 = [harbour_10x, harbour_10y]
     156k102 = [harbour_11x, harbour_11y]
     157k112 = [harbour_12x, harbour_12y]
     158k122 = [harbour_13x, harbour_13y]
     159k132 = [harbour_14x, harbour_14y]
     160k142 = [harbour_15x, harbour_15y]
     161
     162#harbour_polygon_2 = [k02, k112, k122, k12, k22, k62, k72, k82, k92, k102, k32, k42, k52]
     163#harbour_polygon_2 = [k122, k12, k22, k62, k72, k102] #worked
     164harbour_polygon_2 = [k02, k112, k122, k12, k22, k62, k72, k82, k102, k42, k52] #worked
     165
     166#Interior region - Botany Bay - take 2
     167bb_1x = degminsec2decimal_degrees(-34,3,0)
     168bb_1y = degminsec2decimal_degrees(151,2,30)
     169bb_10x = degminsec2decimal_degrees(-34,3,0)
     170bb_10y = degminsec2decimal_degrees(151,8,0)
     171bb_2x = degminsec2decimal_degrees(-34,3,0)
     172bb_2y = degminsec2decimal_degrees(151,14,0)
     173bb_3x = degminsec2decimal_degrees(-33,53,30)
     174bb_3y = degminsec2decimal_degrees(151,17,20)
     175bb_4x = degminsec2decimal_degrees(-33,53,0)
     176bb_4y = degminsec2decimal_degrees(151,8,0)
     177bb_5x = degminsec2decimal_degrees(-33,57,30)
     178bb_5y = degminsec2decimal_degrees(151,8,0)
     179bb_6x = degminsec2decimal_degrees(-33,57,30)
     180bb_6y = degminsec2decimal_degrees(151,2,30) 
     181bb_7x = degminsec2decimal_degrees(-33,53,30)
     182bb_7y = degminsec2decimal_degrees(151,12,30) 
     183bb_8x = degminsec2decimal_degrees(-33,55,20)
     184bb_8y = degminsec2decimal_degrees(151,8,0) 
     185bb_9x = degminsec2decimal_degrees(-33,55,20)
     186bb_9y = degminsec2decimal_degrees(151,12,30) 
     187
     188j02 = [bb_1x, bb_1y]
     189j12 = [bb_2x, bb_2y]
     190j22 = [bb_3x, bb_3y]
     191j32 = [bb_4x, bb_4y]
     192j42 = [bb_5x, bb_5y]
     193j52 = [bb_6x, bb_6y]
     194j62 = [bb_7x, bb_7y]
     195j72 = [bb_8x, bb_8y]
     196j82 = [bb_9x, bb_9y]
     197j92 = [bb_10x, bb_10y]
     198
     199botanybay_polygon_2 = [j92, j12, j22, j62, j82, j72, j42, j02]
     200#botanybay_polygon_2 = [j02, j12, j22, j32, j42, j52] # worked
     201
     202# this section sets up source origin for slump scenario
     203# close to harbour opening (343000,6255000)
     204# x0 = 28964
     205# y0 = 30049
     206# around 10km harbour opening (353000,6255000)
     207#x0 = 38964
     208#y0 = 30049
     209# around 18km across from harbour(361000,6255000)
     210#x0 = 46964
     211#y0 = 30049
     212# around 42km across from harbour(385000,6255000)
     213x0 = 28964 + 42000
     214y0 = 30049
     215# close to botany bay opening (340000,6236000)
     216# x0 = 25964
     217# y0 = 11049
     218# around 10km from botany bay opening (350000,6236000)
     219# x0 = 35964
     220# y0 = 11049
     221# around 21km from botany bay opening (361000,6236000)
     222#x0 = 46964
     223#y0 = 11049
     224
     225# not used for sydney scenario, original interior regions listed though
    56226# setting up problem area for doing just around the harbour
    57227hsouth = degminsec2decimal_degrees(-33,54,0)
     
    66236polygon_h = [hp0, hp1, hp2, hp3]
    67237
    68 #Interior regions - the Harbour
     238#Interior regions - the Harbour - take 1
    69239harbour_south = degminsec2decimal_degrees(-33,53,0)
    70240harbour_north = degminsec2decimal_degrees(-33,47,0)
    71241harbour_west = degminsec2decimal_degrees(151,5,0)
    72242harbour_east = degminsec2decimal_degrees(151,19,0)
     243
     244#harbour_south1 = degminsec2decimal_degrees(-33,53,0)
     245#harbour_south2 = degminsec2decimal_degrees(-33,52,0)
     246#harbour_north1 = degminsec2decimal_degrees(-33,45,0)
     247#harbour_north2 = degminsec2decimal_degrees(-33,48,0)
     248#harbour_west = degminsec2decimal_degrees(151,5,0)
     249#harbour_east = degminsec2decimal_degrees(151,19,0)
    73250
    74251k0 = [harbour_south, harbour_west]
     
    91268polygon_bb = [bp0, bp1, bp2, bp3]
    92269
    93 #Interior region - Botany Bay
     270#Interior region - Botany Bay - take 1
    94271botanybay_south = degminsec2decimal_degrees(-33,58,0)
    95272botanybay_north = degminsec2decimal_degrees(-34,1,0)
     
    103280
    104281botanybay_polygon = [j0, j1, j2, j3]
    105 
  • production/sydney_2006/run_sydney_smf.py

    r2292 r2293  
    135135    # len = 4500.0
    136136    len = 15000.0
    137     thk = 670 #thk = 760.0
     137    #thk = 670 #thk = 760.0
     138    thk = 75.0
    138139    wid = 15000.0 # for slump scenario, wid=len
    139140    #dep = 400 #1000.0 #d = 400 doesn't make sense with thk = 670!
    140     dep = 1000.0
     141    dep = 400.0
    141142    rad = 3330
    142143    dp = 0.23
Note: See TracChangeset for help on using the changeset viewer.