Changeset 6797
- Timestamp:
- Apr 15, 2009, 8:51:36 AM (16 years ago)
- Location:
- anuga_work/production/australia_ph2/dampier
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/dampier/Arc_asc2raster_GDA94z50.py
r6605 r6797 25 25 output_dir="anuga\\outputs\\" 26 26 27 ##time_dir1 = '20081209_155330_run_final_0_27255_250m_none_lfountai' 28 ##time_dir2 = '20081209_161229_run_final_0_68693_250m_none_lfountai' 29 ##time_dirs = [time_dir1, time_dir2] 27 time_dir1 = '20090324_140112_run_final_0_70385_Bf_all_kvanputt' 30 28 31 #time_dir1 = '20090311_152848_run_final_0_70385_1203_Bf_jgriffin' 32 time_dir2 = '20090323_165436_run_final_0_70385_1203_Tb_jgriffin' 33 #time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin' 34 #time_dir4 = '20090311_153054_run_final_0_70339_1203_Bf_jgriffin' 35 ##time_dir5 = '20081211_095614_run_final_0.0_27283_alpha0.1_jgriffin' 36 ##time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt' 37 ## 38 time_dirs = [time_dir2]#, time_dir2, time_dir3, time_dir4] #, time_dir3, time_dir4, time_dir5, time_dir6] 29 time_dirs = [time_dir1] 39 30 40 31 for time_dir in time_dirs: … … 46 37 ## ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean" 47 38 48 ##print 'Process: Create File GDB'49 ##gp.CreateFileGDB_management(folder, "raster")39 print 'Process: Create File GDB' 40 gp.CreateFileGDB_management(folder, "raster") 50 41 51 42 gp.Workspace = raster_gbd -
anuga_work/production/australia_ph2/dampier/export_results_max.py
r6594 r6797 22 22 directory = project.output_folder 23 23 24 ##time_dir1 = '200903 11_152848_run_final_0_70385_1203_Bf_jgriffin'25 ##time_dir2 = '200903 11_152929_run_final_0_70385_1203_Tb_jgriffin'24 ##time_dir1 = '20090324_110442_run_final_0_70385_Bf_30km_kvanputt' 25 ##time_dir2 = '20090324_123534_run_final_0_70385_Bf_60km_kvanputt' 26 26 ##time_dir3 = '20090311_153009_run_final_0_70454_1203_Bf_jgriffin' 27 27 time_dir4 = '20090323_165436_run_final_0_70385_1203_Tb_jgriffin' … … 52 52 53 53 # one or more key strings from var_equations above 54 var = [' stage']54 var = ['elevation'] 55 55 56 56 ###### -
anuga_work/production/australia_ph2/dampier/project.py
r6605 r6797 29 29 tide = 0 # difference between MSL and HAT 30 30 zone = 50 # specify zone of model 31 event_number = 70385 # the event number or the mux file name 31 central_meridian = None 32 33 import sys 34 if len(sys.argv) > 1: 35 event_number = int(sys.argv[1]) 36 else: 37 event_number = 70844 # the event number or the mux file name 38 39 event_number_list = [70844, 27338, 70337] # To piggy back multiple events 40 32 41 alpha = 0.1 # smoothing parameter for mesh 33 42 friction=0.01 # manning's friction coefficient 34 43 starttime=0 # start time for simulation 35 finaltime=80000 # final time for simulation 44 finaltime=1000 # final time for simulation 45 46 # for sensetivity study 47 model = 'simple' #60km 90km 120km all 36 48 37 49 # index is only used when wave = Tb 38 index = 1 203# index from the PTHA39 wave = ' Bf' # Bf (sts wave) Tb (index wave)40 41 setup = ' final' # This can be one of three values50 index = 1307 # index from the PTHA 51 wave = 'Tb' # Bf (sts wave) Tb (index wave) 52 53 setup = 'trial' # This can be one of three values 42 54 # trial - coarsest mesh, fast 43 55 # basic - coarse mesh … … 55 67 #------------------------------------------------------------------------------- 56 68 57 output_comment = [setup, tide, event_number, index, wave]69 output_comment = [setup, tide, event_number, wave, model] 58 70 59 71 #------------------------------------------------------------------------------- … … 64 76 # Used in build_elevation.py 65 77 # Format for ascii grids, as produced in ArcGIS + a projection file 66 ascii_grid_filenames = [ ] # 250m grid 200578 ascii_grid_filenames = ['grid250m_simple'] # 250m grid 2005 67 79 68 80 # Format for point is x,y,elevation (with header) 69 point_filenames = [ 'grid250m.txt'] # 250m grid 200581 point_filenames = [] # 250m grid 2005 70 82 71 83 ### Add csv header list to all files in point_filenames … … 77 89 # Used in build_elevation.py 78 90 # Format for points easting,northing (no header) 79 bounding_polygon_filename = ' bounding_polygon.csv'91 bounding_polygon_filename = 'poly_all_' + model + '.csv' 80 92 bounding_polygon_maxarea = 100000 81 93 … … 106 118 # Thinned ordering file from Hazard Map (geographic) 107 119 # Format is index,latitude,longitude (with header) 108 urs_order_filename = 'urs_order .csv'120 urs_order_filename = 'urs_order_' + model + '.csv' 109 121 110 122 # Landward bounding points 111 123 # Format easting,northing (no header) 112 landward_boundary_filename = 'landward_boundary .csv'124 landward_boundary_filename = 'landward_boundary_' + model + '.csv' 113 125 114 126 # MUX input filename. … … 121 133 multi_mux = True 122 134 123 124 135 # Specify if share cache is to be used 136 # Whatever is specified here will be relative to INUNDATION_HOME/.cache 137 # If nothing is specified, local cache will be used. 138 cachedir = '.python_cache_phII' 125 139 126 140 ################################################################################ … … 164 178 home = join(os.getenv(ENV_INUNDATIONHOME), 'data') # Absolute path for data folder 165 179 muxhome = os.getenv(ENV_MUXHOME) 166 180 181 # Create absolute pathname for cache directory 182 # and change caching to use it 183 if 'cachedir' in dir(): 184 cachedir = join(os.getenv(ENV_INUNDATIONHOME), '.cache', cachedir) 185 from anuga.caching import caching 186 caching.set_option('cachedir', cachedir) 187 167 188 # check various directories/files that must exist 168 189 anuga_folder = join(home, state, scenario_folder, 'anuga') … … 201 222 202 223 # The absolute pathname for the .sts file, generated in build_boundary.py 203 event_sts = join(event_folder, scenario_name) 224 event_sts = join(event_folder, scenario_name)+ '_' +model 204 225 205 226 # The absolute pathname for the output folder names -
anuga_work/production/australia_ph2/dampier/run_model.py
r6594 r6797 45 45 from anuga.shallow_water.data_manager import copy_code_files 46 46 from anuga.shallow_water.data_manager import urs2sts 47 from anuga.geospatial_data.geospatial_data import Geospatial_data 47 48 from anuga.utilities.polygon import read_polygon, Polygon_function 48 49 … … 51 52 import build_urs_boundary as bub 52 53 import prepare_timeboundary as TB 54 55 ##Used to clip elevation pts file 56 ##elevation_geospatial = Geospatial_data(file_name=project.combined_elevation+'.pts', 57 ## verbose=True) 58 ##elevation_geospatial = elevation_geospatial.clip(project.bounding_polygon) 59 ##elevation = project.combined_elevation + '_' + project.model 60 ##elevation_geospatial.export_points_file(elevation + '.pts') 53 61 54 62 #------------------------------------------------------------------------------- … … 128 136 domain.set_quantity('friction', project.friction) 129 137 domain.set_quantity('elevation', 130 filename =project.combined_elevation+'.pts',138 filename = elevation + '.pts', 131 139 use_cache=True, 132 140 verbose=True, … … 156 164 verbose=True) 157 165 domain.set_boundary({'back': Br, 158 'side': B d,166 'side': Bt, 159 167 'ocean': Bf}) 160 168
Note: See TracChangeset
for help on using the changeset viewer.