- Timestamp:
- May 1, 2009, 1:37:13 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/perth/Arc_asc2raster_GDA94z50.py
r6842 r6942 25 25 output_dir="anuga\\outputs\\" 26 26 27 time_dir1 = '20090414_110305_run_final_0_26978_803_Tb__kvanputt' 28 time_dir2 = '20090415_090717_run_final_0_64479_803_Tb__kvanputt' 29 time_dir3 = '20090416_050507_run_final_0_70921_803_Tb__kvanputt' 27 time_dir1 = '20090415_090717_run_final_0_64479_803_Tb__kvanputt' 28 time_dir2 = '20090416_050507_run_final_0_70921_803_Tb__kvanputt' 30 29 31 events = [[time_dir1,0.658849],[time_dir2,0.543065],[time_dir3,0.52893]] 30 time_dir4 = '20090421_170526_run_final_0_64479_803_Tb_internal_mhingee' 32 31 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 32 ##time_dirs = [time_dir1, time_dir2]#, time_dir3] 33 time_dirs = [time_dir1] 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 + '* stage_max.asc')56 input_ascii = glob.glob(folder + '*elevation.asc') 60 57 61 58 for infile in input_ascii: … … 75 72 76 73 print 'Process: Define Projection' 77 # 74 #GDA_1994_MGA_Zone_54 78 75 gp.DefineProjection_management(output_DEM, "PROJCS['GDA_1994_MGA_Zone_50',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]" 79 76 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',117.0],PARAMETER['Scale_Factor',0.9996]" … … 83 80 ## print 'Process: Extract by Mask' 84 81 ## gp.ExtractByMask_sa(output_DEM, land, output_extract) 85 86 div_file = output_DEM.replace('stage','div_stage')87 if not div_file == output_DEM:88 print 'divide', output_DEM,' by', max_wave,' and call',div_file89 gp.Divide_sa(output_DEM,max_wave,div_file)90 82 83
Note: See TracChangeset
for help on using the changeset viewer.