Changeset 2706 for production/onslow_2006/run_onslow.py
- Timestamp:
- Apr 12, 2006, 1:50:07 PM (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
production/onslow_2006/run_onslow.py
r2683 r2706 35 35 36 36 from geospatial_data import * 37 import sys 38 from pyvolution.util import Screen_Catcher 37 39 38 40 #------------------------------------------------------------------------------- … … 61 63 copy (project.codedir + 'run_onslow.py', project.outputtimedir + 'run_onslow.py') 62 64 print'output dir', project.outputtimedir 65 66 #normal screen output is stored in 67 screen_output_name = project.outputtimedir + "screen_output.txt" 68 69 #used to catch screen output to file 70 sys.stdout = Screen_Catcher(screen_output_name) 71 sys.stderr = Screen_Catcher(screen_output_name) 72 63 73 ''' 64 74 copied_files = False … … 75 85 #if access(project.MOST_dir + project.boundary_basename+'_ha.nc',F_OK) == 1 : 76 86 # print' most file', project.MOST_dir + project.boundary_basename 77 ''' 87 78 88 79 89 # fine data (clipping the points file to smaller area) … … 101 111 print'export G' 102 112 G.new_export_points_file(project.combined_dem_name + '.pts') 103 ''' 113 104 114 105 115 #------------------------------------------------------------------------------- … … 112 122 113 123 # original 114 interior_res = 1000115 interior_regions = [[project.poly_onslow, interior_res],116 [project.poly_thevenard, interior_res]]117 118 print 'number of interior regions' , len(interior_regions)124 #interior_res = 1000 125 #interior_regions = [[project.poly_onslow, interior_res], 126 # [project.poly_thevenard, interior_res]] 127 128 print 'number of interior regions'#, len(interior_regions) 119 129 120 130 from caching import cache … … 125 135 'bottom': [4], 'bottomright': [5], 126 136 'topright':[6]}, 127 'maximum_triangle_area': 500 00,128 'filename': meshname ,129 'interior_regions': interior_regions},130 verbose = True)137 'maximum_triangle_area': 500, 138 'filename': meshname}, 139 # 'interior_regions': interior_regions}, 140 verbose = True) 131 141 132 142 … … 181 191 use_cache = True, 182 192 verbose = True, 183 alpha = 0.0 5193 alpha = 0.0 184 194 ) 185 195 … … 245 255 t0 = time.time() 246 256 247 for t in domain.evolve(yieldstep = 500, finaltime = 3000):257 for t in domain.evolve(yieldstep = 10, finaltime = 10): 248 258 domain.write_time() 249 259 domain.write_boundary_statistics(tags = 'top') 250 260 ''' 251 261 for t in domain.evolve(yieldstep = 120, finaltime = 31000): 252 #,skip_initial_step = True):262 ,skip_initial_step = True): 253 263 domain.write_time() 254 264 domain.write_boundary_statistics(tags = 'top') 255 265 ''' 256 266 print 'That took %.2f seconds' %(time.time()-t0) 257 267
Note: See TracChangeset
for help on using the changeset viewer.