Changeset 5781
- Timestamp:
- Sep 24, 2008, 11:56:50 AM (15 years ago)
- Location:
- anuga_work/production/perth
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/perth/build_boundary_27255.py
r5761 r5781 4 4 5 5 Input: order_filename from project.py 6 event_number needs to be reflected in the project file 6 7 Output: creates a sts file and csv files stored in project.boundaries_dir_event. 7 8 The run_perth.py is reliant on the output of this script. -
anuga_work/production/perth/build_boundary_27283.py
r5761 r5781 4 4 5 5 Input: order_filename from project.py 6 event_number needs to be reflected in the project file 6 7 Output: creates a sts file and csv files stored in project.boundaries_dir_event. 7 8 The run_perth.py is reliant on the output of this script. -
anuga_work/production/perth/build_perth.py
r5761 r5781 130 130 if access(project.topographies_dir,F_OK) == 0: 131 131 mkdir (project.topographies_dir) 132 G_clipped.export_points_file(project.combined_dir_name + '.pts')133 #G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC132 #G_clipped.export_points_file(project.combined_dir_name + '.pts') 133 G_clipped.export_points_file(project.combined_dir_name + '.txt') #Use for comparision in ARC 134 134 135 135 import sys -
anuga_work/production/perth/export_results_all.py
r5761 r5781 22 22 #time_dir = '20080526_104946_run_final_0.6_test_kvanputt' 23 23 #time_dir = '20080530_170833_run_final_0.6_exmouth_kvanputt' 24 time_dir1 = '20080815_103442_run_final_0.0_polyline_alpha0.1_kvanputt'25 time_dir 2 = '20080909_151438_run_final_0.0_polyline_alpha0.1_kvanputt'24 #time_dir1 = '20080815_103442_run_final_0.0_polyline_alpha0.1_kvanputt' 25 time_dir1 = '20080912_154439_run_final_0.6_27255_alpha0.1_kvanputt' 26 26 27 27 #cellsize = 20 … … 29 29 #timestep = 0 30 30 #area = ['Geordie', 'Sorrento', 'Fremantle', 'Rockingham'] 31 area = ['All'] 32 #var = [1,2] # Absolute momentum and depth 31 area = 'All' 32 which_area = area 33 var = [1,2] # Absolute momentum and depth 33 34 #var = [2] # depth 34 var = [0,4] #stage and elevation35 #var = [0,4] #stage and elevation 35 36 36 time_dir = [time_dir1, time_dir2]37 time_dirs = [time_dir1] 37 38 for time_dir in time_dirs: 38 39 39 40 name1 = directory+time_dir+sep+project.scenario_name 40 name2 = directory+time_dir+sep+'sww2'+sep+project.scenario_name+'_time_39900_0' #need to get assistance on how to make this into anything41 #name2 = directory+time_dir+sep+'sww2'+sep+project.scenario_name+'_time_39900_0' #need to get assistance on how to make this into anything 41 42 42 names = [name1 , name2]43 names = [name1] 43 44 for name in names: 44 45 45 for which_area in area:46 if which_area == 'All':47 48 if which_area == 'Geordie':49 easting_min = project.xminGeordie50 easting_max = project.xmaxGeordie51 northing_min = project.yminGeordie52 northing_max = project.ymaxGeordie53 54 if which_area == 'Sorrento':55 easting_min = project.xminSorrento56 easting_max = project.xmaxSorrento57 northing_min = project.yminSorrento58 northing_max = project.ymaxSorrento59 60 if which_area == 'Fremantle':61 easting_min = project.xminFremantle62 easting_max = project.xmaxFremantle63 northing_min = project.yminFremantle64 northing_max = project.ymaxFremantle65 66 if which_area == 'Rockingham':67 easting_min = project.xminRockingham68 easting_max = project.xmaxRockingham69 northing_min = project.yminRockingham70 northing_max = project.ymaxRockingham46 ## for which_area in area: 47 ## if which_area == 'All': 48 ## 49 ## if which_area == 'Geordie': 50 ## easting_min = project.xminGeordie 51 ## easting_max = project.xmaxGeordie 52 ## northing_min = project.yminGeordie 53 ## northing_max = project.ymaxGeordie 54 ## 55 ## if which_area == 'Sorrento': 56 ## easting_min = project.xminSorrento 57 ## easting_max = project.xmaxSorrento 58 ## northing_min = project.yminSorrento 59 ## northing_max = project.ymaxSorrento 60 ## 61 ## if which_area == 'Fremantle': 62 ## easting_min = project.xminFremantle 63 ## easting_max = project.xmaxFremantle 64 ## northing_min = project.yminFremantle 65 ## northing_max = project.ymaxFremantle 66 ## 67 ## if which_area == 'Rockingham': 68 ## easting_min = project.xminRockingham 69 ## easting_max = project.xmaxRockingham 70 ## northing_min = project.yminRockingham 71 ## northing_max = project.ymaxRockingham 71 72 72 73 … … 93 94 quantityname = 'elevation' #Elevation 94 95 95 if which_area == 'All'96 print 'start sww2dem', which_area96 else: 97 print 'start sww2dem',area 97 98 sww2dem(name, basename_out = outname, 98 99 quantity = quantityname, … … 102 103 verbose = True, 103 104 format = 'asc') 104 105 else:106 print 'start sww2dem',which_area, easting_min107 sww2dem(name, basename_out = outname,108 quantity = quantityname,109 #timestep = timestep,110 cellsize = cellsize,111 easting_min = easting_min,112 easting_max = easting_max,113 northing_min = northing_min,114 northing_max = northing_max,115 reduction = max,116 verbose = True,117 format = 'asc')118 105 ## 106 ## else: 107 ## print 'start sww2dem',which_area, easting_min 108 ## sww2dem(name, basename_out = outname, 109 ## quantity = quantityname, 110 ## #timestep = timestep, 111 ## cellsize = cellsize, 112 ## easting_min = easting_min, 113 ## easting_max = easting_max, 114 ## northing_min = northing_min, 115 ## northing_max = northing_max, 116 ## reduction = max, 117 ## verbose = True, 118 ## format = 'asc') 119 ## -
anuga_work/production/perth/project.py
r5761 r5781 41 41 scenario = 'perth_tsunami_scenario' 42 42 ##One or all can be changed each time the run_scenario script is excuted 43 tide = 0 .644 event_number = 2725545 #event_number = 2728343 tide = 0 #0.6 44 #event_number = 27255 45 event_number = 27283 46 46 alpha = 0.1 47 47 friction=0.01 -
anuga_work/production/perth/run_perth.py
r5761 r5781 130 130 # Setup initial conditions 131 131 #------------------------------------------------------------------------- 132 print 'Setup initial conditions'132 print 'Setup initial conditions' 133 133 134 134 #following sets the stage/water to be offcoast only
Note: See TracChangeset
for help on using the changeset viewer.