- Timestamp:
- May 26, 2009, 3:52:17 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/gold_coast_2009/Arc_asc2raster_GDA94z56.py
r7004 r7089 25 25 output_dir = "anuga\\outputs\\" 26 26 27 time_dir1 = '20090327_135943_run_final_0_51469_lfountai' 28 ##time_dir2 = '20090327_153206_run_final_0_51348_lfountai' 29 ##time_dir3 = '20090327_153231_run_final_0_51253_lfountai' 30 ##time_dir4 = '20090327_153304_run_final_0_50863_lfountai' 31 ##time_dir5 = '20090327_153403_run_final_0_51423_lfountai' 32 ##time_dir6 = '20090327_153602_run_final_0_62833_lfountai' 27 time_dir1 = '20090508_150215_run_final_0_51469_lfountai' 28 time_dir2 = '20090511_165539_run_final_0_50863_lfountai' 29 time_dir3 = '20090511_161526_run_final_1.1_51469_kvanputt' 30 time_dir4 = '20090518_154710_run_final_0_50994_lfountai' 31 ##time_dir5 = '20090519_160510_run_final_1.1_50994_lfountai' 32 ##time_dir6 = '20090521_220101_run_final_1.1_50863_kvanputt' 33 ##time_dir7 = '20090522_164526_run_final_0_51392_lfountai' 34 ##time_dir8 = '20090522_164640_run_final_1.1_51392_lfountai' 35 ##time_dir9 = '20090522_164948_run_final_0_51423_lfountai' 36 ##time_dir10 = '20090522_165600_run_final_1.1_51423_lfountai' 33 37 34 time_dirs = [time_dir1 ] #[time_dir2, time_dir3, time_dir4, time_dir5, time_dir6]38 time_dirs = [time_dir1, time_dir2, time_dir3, time_dir4] #, time_dir5, time_dir6] 35 39 36 40 for time_dir in time_dirs: … … 40 44 raster_gbd = folder + 'raster.gdb' 41 45 #contour = raster_gbd + '\\contour_dep' 42 land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Coast_250"46 ## land = scenario_dir + "ArcGIS\\gold_coast.gdb\\Coast_250" 43 47 ## ocean = scenario_dir + "map_work\\Perth.gdb\\Outlines\\initial_condition_ocean" 44 48 … … 51 55 52 56 #replication dictionary 53 replicate = (('gold_coast', ''), ('_time_2 1960_0', 'b'), ('_time_43920_0', 'c'),54 (' time_65880_0', 'd'), ('_', ''), ('max','M_'), ('depth','_dep_'),57 replicate = (('gold_coast', ''), ('_time_29220_0', 'b'), ('_time_58440_0', 'c'), 58 ('_', ''), ('max','M_'), ('depth','_dep_'), 55 59 ('speed', '_spe_'), ('elevation', '_ele_'), ('stage','_stage')) 56 60 57 61 generate_filename = [] 58 input_ascii = glob.glob(folder + '* depth_max.asc')62 input_ascii = glob.glob(folder + '*_max.asc') 59 63 print time_dir 60 64 … … 80 84 81 85 82 output_extract = output_DEM + 'E'83 print 'Output Extract ',output_extract84 85 print 'Process: Extract by Mask'86 gp.ExtractByMask_sa(output_DEM, land, output_extract)87 86 ## output_extract = output_DEM + 'E' 87 ## print 'Output Extract ',output_extract 88 ## 89 ## print 'Process: Extract by Mask' 90 ## gp.ExtractByMask_sa(output_DEM, land, output_extract) 91 ## 88 92
Note: See TracChangeset
for help on using the changeset viewer.