Changeset 4564
- Timestamp:
- Jun 27, 2007, 10:33:33 AM (17 years ago)
- Location:
- anuga_work/production/broome_2006
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/broome_2006/build_broome.py
r4559 r4564 68 68 offshore_dir_name1 = project_urs.offshore_dir_name1 69 69 offshore_dir_name2 = project_urs.offshore_dir_name2 70 70 ''' 71 71 # creates DEM from asc data 72 72 print "creates DEMs from asc data" … … 102 102 103 103 print'clip combined geospatial object by bounding polygon' 104 G_clipped = G.clip(project_urs.poly_all)104 #G_clipped = G.clip(project_urs.poly_all) 105 105 #FIXME: add a clip function to pts 106 106 #print'shape of clipped data', G_clipped.get_data_points().shape … … 113 113 114 114 print'split' 115 G_small, G_other = G _clipped.split(.10,verbose=True)115 G_small, G_other = G.split(.10,verbose=True) 116 116 117 117 118 118 print 'export',project_urs.combined_dir_name + '.txt' 119 G _clipped.export_points_file(project_urs.combined_dir_name + '.txt')119 G.export_points_file(project_urs.combined_dir_name + '.txt') 120 120 print 'export', project_urs.combined_small_dir_name + '.txt' 121 121 G_small.export_points_file(project_urs.combined_small_dir_name + '.txt') 122 122 #G_clipped.export_points_file(project_urs.combined_dir_name + '.xya') 123 ''' 123 124 124 125 125 ''' 126 print'project_urs.combined_dir_name + 1.xya',project_urs.combined_dir_name + '1.xya' 127 G_all=Geospatial_data(file_name = project_urs.combined_dir_name + '1.xya') 126 print'export',project_urs.combined_dir_name+ '.txt' 127 G_all=Geospatial_data(file_name = project_urs.combined_dir_name+ '.txt') 128 128 print'split' 129 129 G_all_1, G_all_2 = G_all.split(.10) 130 130 print'export 1' 131 G_all_1.export_points_file(project_urs.combined_dir_name+'_small 1' + '.xya')131 G_all_1.export_points_file(project_urs.combined_dir_name+'_small' + '.xya') 132 132 print'export 2' 133 133 G_all_2.export_points_file(project_urs.combined_dir_name+'_other1' + '.xya') 134 134 ''' 135 135 136 136 #------------------------------------------------------------------------- -
anuga_work/production/broome_2006/project_urs.py
r4559 r4564 32 32 33 33 #tide = -5.3 34 tide = 035 #tide = 4.934 #tide = 0 35 tide = 4.9 36 36 37 37 #Maybe will try to make project a class to allow these parameters to be passed in. … … 41 41 finaltime=25000 42 42 setup='final' 43 source=' onslow'43 source='exmouth' 44 44 45 45 if setup =='trial': … … 70 70 71 71 #final topo name 72 combined_name ='broome_combined_elevation_z51 '72 combined_name ='broome_combined_elevation_z51_original' 73 73 combined_name1 ='broome_combined_elevation1' 74 74 combined_name_unclipped1 ='broome_combined_elevation_unclipped1' … … 124 124 output_dir = anuga_dir+'outputs'+sep 125 125 output_build_time_dir = output_dir+build_time+'_'+source+sep 126 #output_run_time_dir = output_dir +run_time+dir_comment+sep127 output_run_time_dir = output_dir +'20070622_004008_run_final_0_onslow_nbartzis'+sep126 output_run_time_dir = output_dir +run_time+dir_comment+sep 127 #output_run_time_dir = output_dir +'20070622_004008_run_final_0_onslow_nbartzis'+sep 128 128 output_run_time_dir_name = output_run_time_dir + scenario_name #Used by post processing 129 129 -
anuga_work/production/broome_2006/run_broome_urs.py
r4559 r4564 63 63 64 64 kwargs['output_dir']=project_urs.output_run_time_dir 65 kwargs['bathy_file']=project_urs.combined_dir_name _unclipped1+'.txt'65 kwargs['bathy_file']=project_urs.combined_dir_name+'.txt' 66 66 # kwargs['bathy_file']=project_urs.combined_small_dir_name + '.pts' 67 67 kwargs['boundary_file']=project_urs.boundaries_in_dir_name + '.sww' 68 68 69 start_screen_catcher(kwargs['output_dir'], myid, numprocs,extra_info='export',print_to_screen=True)70 '''71 69 print 'output_dir',kwargs['output_dir'] 72 70 if myid == 0: … … 143 141 domain.set_quantity('friction', kwargs['friction']) 144 142 145 print 'Start Set quantity' 143 print 'Start Set quantity',kwargs['bathy_file'] 146 144 147 145 domain.set_quantity('elevation', … … 230 228 kwargs['completed']=str(time.time()-t0) 231 229 232 ''' 233 if myid == 0: 234 # store_parameters(**kwargs) 230 if myid == 0: 231 store_parameters(**kwargs) 235 232 236 233 print 'memory usage before del domain',mem_usage() 237 #del domain234 del domain 238 235 print 'memory usage after del domain',mem_usage() 239 236 … … 254 251 verbose = False, 255 252 origin = None, 256 datum = ' WGS84',253 datum = 'GDA94', 257 254 format = 'asc') 258 255
Note: See TracChangeset
for help on using the changeset viewer.