Changeset 6949
- Timestamp:
- May 5, 2009, 9:35:20 AM (15 years ago)
- Location:
- anuga_work/production/australia_ph2
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/carnarvon/Arc_asc2raster_GDA94z50.py
r6741 r6949 25 25 output_dir="anuga\\outputs\\" 26 26 27 time_dir1 = '20090 324_150124_run_final_0_27285_109_Bf_mhingee'28 time_dir2 = ' 20090324_161452_run_final_0_27285_109_Tb_mhingee'27 time_dir1 = '20090429_170742_run_trial_0_27285_84_Tb__kvanputt' 28 time_dir2 = '' 29 29 30 time_dirs = [time_dir 2] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6]30 time_dirs = [time_dir1] #, time_dir2] #4, time_dir5] #, time_dir3, time_dir4, time_dir5, time_dir6] 31 31 32 32 for time_dir in time_dirs: … … 52 52 53 53 generate_filename = [] 54 input_ascii = glob.glob(folder + '* .asc')54 input_ascii = glob.glob(folder + '*elevation.asc') 55 55 56 56 for infile in input_ascii: … … 70 70 71 71 print 'Process: Define Projection' 72 gp.DefineProjection_management(output_DEM, "PROJCS[' GDA_1994_MGA_Zone_49',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]"73 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',11 1.0],PARAMETER['Scale_Factor',0.9996]"72 gp.DefineProjection_management(output_DEM, "PROJCS['CM_113.5',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]" 73 ",PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',113.5],PARAMETER['Scale_Factor',0.9996]" 74 74 ",PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]") 75 75 ## output_extract = output_DEM + '_E' -
anuga_work/production/australia_ph2/carnarvon/project.py
r6856 r6949 40 40 friction=0.01 # manning's friction coefficient 41 41 starttime=0 # start time for simulation 42 finaltime= 1000 # final time for simulation43 setup = ' trial' # This can be one of three values42 finaltime=60000 # final time for simulation 43 setup = 'final' # This can be one of three values 44 44 # trial - coarsest mesh, fast 45 45 # basic - coarse mesh … … 128 128 # Thinned ordering file from Hazard Map (geographic) 129 129 # Format is index,latitude,longitude (with header) 130 urs_order_filename = 'urs_order _simple.csv'130 urs_order_filename = 'urs_order.csv' 131 131 132 132 # Landward bounding points -
anuga_work/production/australia_ph2/gold_coast/Arc_asc2raster_GDA94z50.py
r6842 r6949 29 29 time_dir3 = '20090409_175154_run_final_0_51453_3727_Tb__kvanputt' 30 30 31 time_dirs = [time_dir1, time_dir2, time_dir3] 31 ##time_dir1 = '20090421_150448_run_final_0_50883_3727_Tb_internal_mhingee' 32 33 time_dirs = [time_dir1]#, time_dir2, time_dir3] 32 34 33 35 for time_dir in time_dirs: … … 39 41 ## ocean = scenario_dir + "map_work\\\port_hedland.gdb\\outlines\\initial_conditions_ocean" 40 42 41 print 'Process: Create File GDB'42 gp.CreateFileGDB_management(folder, "raster")43 ## print 'Process: Create File GDB' 44 ## gp.CreateFileGDB_management(folder, "raster") 43 45 44 46 gp.Workspace = raster_gbd … … 53 55 54 56 generate_filename = [] 55 input_ascii = glob.glob(folder + '* max.asc')57 input_ascii = glob.glob(folder + '*elevation.asc') 56 58 57 59 for infile in input_ascii:
Note: See TracChangeset
for help on using the changeset viewer.