Changeset 6491


Ignore:
Timestamp:
Mar 11, 2009, 2:17:50 PM (15 years ago)
Author:
myall
Message:
 
Location:
anuga_work/production
Files:
3 edited

Legend:

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

    r6486 r6491  
    8888
    8989# Number of boundary segments
    90 num_ocean_segments = len(event_sts) - 1
     90num_ocean_segments = len(event_sts) #- 1
    9191# Number of landward_boundary points
    9292num_land_points = 0 #file_length(project.landward_boundary)
     
    9494# Boundary tags refer to project.landward_boundary
    9595# 4 points equals 5 segments start at N
    96 boundary_tags={'back': range(num_ocean_segments+1,
    97                              num_ocean_segments+num_land_points),
    98                'side': [num_ocean_segments,
    99                         num_ocean_segments+num_land_points],
     96boundary_tags={#'back': range(num_ocean_segments+1,
     97               #              num_ocean_segments+num_land_points),
     98               #'side': [num_ocean_segments,
     99               #         num_ocean_segments+num_land_points],
    100100               'ocean': range(num_ocean_segments)}
    101101
     
    106106                                    interior_regions=project.interior_regions,
    107107                                    mesh_filename=project.meshes,
    108                                     use_cache=True,
     108                                    use_cache=False,
    109109                                    verbose=True)
    110110print domain.statistics()
     
    127127else:
    128128    IC = 0
    129 domain.set_quantity('stage', IC, use_cache=True, verbose=True)
     129domain.set_quantity('stage', IC, use_cache=False, verbose=True)
    130130domain.set_quantity('friction', project.friction)
    131131domain.set_quantity('elevation',
    132132                    filename=project.combined_elevation+'.pts',
    133                     use_cache=True,
     133                    use_cache=False,
    134134                    verbose=True,
    135135                    alpha=project.alpha)
     
    155155                        default_boundary=Bd,
    156156                        boundary_polygon=bounding_polygon_sts,                   
    157                         use_cache=True,
     157                        use_cache=False,
    158158                        verbose=True)
    159     domain.set_boundary({'back': Br,
    160                      'side': Bd,
    161                      'ocean': Bf})
     159    domain.set_boundary({'ocean': Bf})
    162160   
    163161elif project.wave == 'Tb':
     
    165163                       verbose=True)
    166164
    167     domain.set_boundary({'back': Br,
    168                          'side': Bd,
    169                          'ocean': Tb})
     165    domain.set_boundary({'ocean': Tb})
    170166else:
    171167    print 'No wave specified in project script (Bf or Tb)'
  • anuga_work/production/carnarvon/get_timeseries.py

    r6317 r6491  
    1919
    2020##time_dir1 = '20081202_084220_run_final_1_27283_alpha0.1_kvanputt'
    21 time_dir2 = '20081202_084202_run_final_1_27255_alpha0.1_kvanputt'
     21##time_dir2 = '20081202_084202_run_final_1_27255_alpha0.1_kvanputt'
    2222##time_dir3 = '20081202_084132_run_final_1_68693_alpha0.1_kvanputt'
    2323##time_dir4 = '20081202_084025_run_final_0_68693_alpha0.1_kvanputt'
    24 ##time_dir5 = '20081202_083932_run_final_0_27255_alpha0.1_kvanputt'
     24time_dir5 = '20081202_083932_run_final_0_27255_alpha0.1_kvanputt'
    2525##time_dir6 = '20081201_103449_run_final_0_27283_alpha0.1_kvanputt'
    2626
    27 time_dirs = [time_dir2] #1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
     27time_dirs = [time_dir5] #1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]
    2828
    2929for time_dir in time_dirs:
  • anuga_work/production/carnarvon/project.py

    r6317 r6491  
    9494
    9595# gauges - used in get_timeseries.py
    96 gauge_name = 'MH_gauges.csv' #'carnarvon.csv'
     96gauge_name = 'MH_gauges_deep.csv' #'carnarvon.csv'
    9797gauge_name2 = 'thinned_MGA50.csv'
    9898#buildings - used in run_building_inundation.py
Note: See TracChangeset for help on using the changeset viewer.