Changeset 6065
- Timestamp:
- Dec 11, 2008, 4:11:14 PM (15 years ago)
- Location:
- anuga_work/production/geraldton
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/geraldton/Arc_asc2raster_GDA94z50.py
r6041 r6065 22 22 gp.overwriteoutput = 1 23 23 24 output_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\"25 scenario_dir="geraldton_tsunami_scenario\\anuga\\outputs\\"24 scenario_dir = "\\\\nas2\\gemd\\georisk_models\\inundation\\data\\western_australia\\geraldton_tsunami_scenario\\" 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 30 28 time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'29 time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'31 ##time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt' 32 ##time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt' 30 33 time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt' 31 time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'34 ##time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt' 32 35 time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt' 33 time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'34 35 time_dirs = [time_dir 1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]36 ##time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt' 37 ## 38 time_dirs = [time_dir3, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6] 36 39 37 40 for time_dir in time_dirs: 38 41 39 42 # Local variables... 40 folder = output_dir + scenario_dir +time_dir +'\\'43 folder = scenario_dir + output_dir + time_dir +'\\' 41 44 raster_gbd = folder + 'raster.gdb' 42 land = scenario_dir + "map_work\\\Geraldton.gdb\\outlines\\Clip_Topography" 45 land = scenario_dir + "map_work\\\Geraldton.gdb\\outlines\\initial_condition" 46 ocean = scenario_dir + "map_work\\\Geraldton.gdb\\outlines\\initial_conditions_ocean" 43 47 44 # Process: Create File GDB... 45 #gp.CreateFileGDB_management(folder, "raster")48 ## print 'Process: Create File GDB' 49 ## gp.CreateFileGDB_management(folder, "raster") 46 50 47 51 gp.Workspace = raster_gbd … … 50 54 51 55 #replication dictionary 52 replicate = (('geraldton', ''),('_', ''),('max',' M_'),56 replicate = (('geraldton', ''),('_', ''),('max','_M'), 53 57 ('CBD', 'CDB'), 54 58 ('depth','_depth'),('speed', '_speed'), 55 ('elevation', '_ele_'), ('stage','_sta _'))59 ('elevation', '_ele_'), ('stage','_stage')) 56 60 57 61 generate_filename = [] 58 input_ascii = glob.glob(folder + '* max.asc')62 input_ascii = glob.glob(folder + '*stage_max.asc') 59 63 60 64 for infile in input_ascii: … … 69 73 70 74 print 'Output DEM ',output_DEM 71 output_extract = output_DEM + 'E'72 print 'Output Extract ',output_extract73 75 74 76 print 'Process: ASCII to Raster' … … 79 81 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]" 80 82 ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]") 81 82 print 'Process: Extract by Mask' 83 gp.ExtractByMask_sa(output_DEM, land, output_extract) 83 ## output_extract = output_DEM + 'E' 84 ## print 'Output Extract ',output_extract 85 ## print 'Process: Extract by Mask' 86 ## gp.ExtractByMask_sa(output_DEM, land, output_extract) 84 87 85 88 -
anuga_work/production/geraldton/export_results_max.py
r6040 r6065 20 20 21 21 directory = project.output_dir 22 ##time_dir1 = '20081209_155330_run_final_0_27255_250m_none_lfountai' 23 ##time_dir2 = '20081209_161229_run_final_0_68693_250m_none_lfountai' 24 ##time_dirs = [time_dir1, time_dir2] 22 25 23 time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt'24 time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt'26 ##time_dir1 = '20081117_141353_run_final_0_27283_alpha0.1_kvanputt' 27 ##time_dir2 = '20081117_141443_run_final_0.6_27283_alpha0.1_kvanputt' 25 28 time_dir3 = '20081117_141520_run_final_0_27255_alpha0.1_kvanputt' 26 time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt'29 ##time_dir4 = '20081117_141558_run_final_0.6_27255_alpha0.1_kvanputt' 27 30 time_dir5 = '20081117_141634_run_final_0_68693_alpha0.1_kvanputt' 28 time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt'29 30 time_dirs = [time_dir 1, time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]31 ##time_dir6 = '20081117_141717_run_final_0.6_68693_alpha0.1_kvanputt' 32 ## 33 time_dirs = [time_dir3, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6] 31 34 32 35 33 36 34 cellsize = 2035 #cellsize = 25037 #cellsize = 10 38 cellsize = 250 36 39 37 40 timestep = None # None means no timestep! … … 42 45 ###### 43 46 44 area = ['CBD'] # strings must match keys in var_equations below45 #area = ['All'] # 'All' means no special areas - the whole thing47 #area = ['CBD'] # strings must match keys in var_equations below 48 area = ['All'] # 'All' means no special areas - the whole thing 46 49 47 50 ###### … … 56 59 57 60 # one or more key strings from var_equations above 58 var = [' depth', 'speed']61 var = ['stage'] 59 62 60 63 ######
Note: See TracChangeset
for help on using the changeset viewer.