Changeset 1920


Ignore:
Timestamp:
Oct 14, 2005, 11:33:57 AM (18 years ago)
Author:
ole
Message:

Karratha

Location:
production/karratha_2005
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • production/karratha_2005/export_results.py

    r1877 r1920  
    55
    66#name = project.outputname + '_0.0tide_notsunami'
    7 name = project.outputname + '_0.0tide'
     7#name = project.outputname + '_0.0tide'
    88
    9 sww2dem(name, basename_out = name+'_crop',
    10         quantity = 'depth',
    11         cellsize = 10,
    12         easting_min = 465646,
    13         easting_max = 480286,
    14         northing_min = 7712865,
    15         northing_max = 7724233,       
     9name = project.outputname
     10
     11sww2dem(name, basename_out = name + '_momentum',
     12        #quantity = 'stage-elevation',  #Depth
     13        quantity = '(xmomentum**2 + ymomentum**2)**0.5',  #Absolute momentum
     14        cellsize = 50,
     15        #cellsize = 10,       
     16        #easting_min = 465646,
     17        #easting_max = 480286,
     18        #northing_min = 7712865,
     19        #northing_max = 7724233,       
    1620        reduction = max,
    1721        verbose = True,
  • production/karratha_2005/project.py

    r1877 r1920  
    77import sys
    88
    9 #Making assumptions abot the location of scenario data
     9#Making assumptions about the location of scenario data
    1010scenario_dir_name = 'karratha_tsunami_scenario_2005'
    11 basename = 'karratha_100m'
     11#basename = 'karratha_100m'
     12basename = 'elevation50m'
    1213boundary_basename = 'SU-AU'
    1314
  • production/karratha_2005/run_karratha.py

    r1877 r1920  
    6262       'origin': project.mesh_origin,
    6363       'mean_stage': tide,
    64        'zscale': 2,                 #Enhance tsunami
     64       'zscale': 1,                 #Enhance tsunami
    6565       'fail_on_NaN': False,
    6666       'inverted_bathymetry': True},
     67      #evaluate = True,
    6768      verbose = True)
    6869#FIXME: Dependencies
     
    110111domain.store = True
    111112
    112 #domain.quantities_to_be_stored = ['stage', 'xmomentum', 'ymomentum']
    113 domain.quantities_to_be_stored = ['stage']
     113domain.quantities_to_be_stored = ['stage', 'xmomentum', 'ymomentum']
     114#domain.quantities_to_be_stored = ['stage']
    114115       
    115116print 'Number of triangles = ', len(domain)
     
    141142                   f=lambda t: [(60<t<660)*4, 0, 0])
    142143
    143 domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bf} )
     144
     145domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bf} ) #MOST tsunami
    144146
    145147#domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bd} )  #Nothing
     
    150152t0 = time.time()
    151153
    152 for t in domain.evolve(yieldstep = 60, finaltime = 15000):
    153     domain.write_time()
    154     domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
     154#for t in domain.evolve(yieldstep = 60, finaltime = 15000):
     155#    domain.write_time()
     156#    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
    155157
    156 for t in domain.evolve(yieldstep = 20, finaltime = 35000,
    157                        skip_initial_step = True):
    158     domain.write_time()
    159     domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
     158#for t in domain.evolve(yieldstep = 20, finaltime = 35000,
     159#                       skip_initial_step = True):
     160#    domain.write_time()
     161#    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')
    160162
    161 for t in domain.evolve(yieldstep = 60, finaltime = 40000,
    162                        skip_initial_step = True):
     163for t in domain.evolve(yieldstep = 60, finaltime = 40000):
     164                       #skip_initial_step = True):
    163165    domain.write_time()
    164166    domain.write_boundary_statistics(tags = 'ocean') #quantities = 'stage')       
Note: See TracChangeset for help on using the changeset viewer.