- Timestamp:
- May 6, 2009, 4:06:54 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/portland/Arc_asc2raster_GDA94z50.py
r6842 r6985 29 29 time_dir3 = '20090409_001857_run_final_0_64322_1969_Tb__mhingee' 30 30 31 events = [[time_dir1,0.430643],[time_dir2,0.462943],[time_dir3,0.418183]] 31 time_dir4 = '20090421_153138_run_final_0_64322_1969_Tb_internal_mhingee' 32 32 33 for event in events: 34 ##for time_dir in time_dirs: 35 time_dir = event[0] 36 max_wave = event[1] 37 print time_dir 38 print max_wave 33 time_dirs = [time_dir3] 34 35 for time_dir in time_dirs: 39 36 # Local variables... 40 37 folder = scenario_dir + output_dir + time_dir +'\\' … … 43 40 ## ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean" 44 41 45 print 'Process: Create File GDB'46 gp.CreateFileGDB_management(folder, "raster")42 ## print 'Process: Create File GDB' 43 ## gp.CreateFileGDB_management(folder, "raster") 47 44 48 45 gp.Workspace = raster_gbd … … 57 54 58 55 generate_filename = [] 59 input_ascii = glob.glob(folder + '* .asc')56 input_ascii = glob.glob(folder + '*elevation.asc') 60 57 61 58 for infile in input_ascii: … … 82 79 ## print 'Process: Extract by Mask' 83 80 ## gp.ExtractByMask_sa(output_DEM, land, output_extract) 84 # do this bit only if there are only stage asc files85 div_file = output_DEM.replace('stage','div_stage')86 print 'divide', output_DEM,' by', max_wave,' and call',div_file87 gp.Divide_sa(output_DEM,max_wave,div_file)88
Note: See TracChangeset
for help on using the changeset viewer.