- Timestamp:
- Jun 21, 2012, 11:30:04 AM (12 years ago)
- Location:
- trunk/anuga_work/development/Busselton
- Files:
-
- 6 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_work/development/Busselton/Analysis/export_results_max.py
r8449 r8452 24 24 rootDir = r"/short/w85/Shane/Output" 25 25 scenario_name = "Busselton" 26 event = "B4" 27 #directory = join(rootDir,event) 26 event = "B2Full" 28 27 29 time_dir1 = "final_"+event+"_sfm547_np1" 30 time_dir2 = "final_"+event+"_sfm547_np1_From21600" 31 32 time_dirs = [time_dir1,time_dir2] 28 time_dirs = [#"final_"+event+"_sfm547_np1",\ 29 # "final_"+event+"_sfm547_np1_From28800",\ 30 # "final_"+event+"_sfm547_np1_From57600",\ 31 # "final_"+event+"_sfm547_np1_From86400",\ 32 # "final_"+event+"_sfm547_np1_From108000",\ 33 "final_"+event+"_sfm547_np1_From115200"] 33 34 34 35 #Modify the cellsize value to set the size of the raster you require … … 66 67 ###### 67 68 68 if myid == 0: 69 time_dir = time_dirs[0] 70 var = vars[0] 69 proc = 0 70 Needed = True 71 72 for dirID in range(len(time_dirs)): 71 73 72 elif myid == 1: 73 time_dir = time_dirs[1] 74 var = vars[1] 74 if myid == proc: 75 time_dir = time_dirs[dirID] 76 var = vars[0] 77 78 elif myid == proc + 1: 79 time_dir = time_dirs[dirID] 80 var = vars[1] 81 82 proc+=2 83 84 if myid >= len(time_dirs)*2: 85 print "I am not needed: Proc ", myid 86 Needed = False 87 88 if Needed: 89 90 if var not in var_equations: 91 print 'Unrecognized variable name: %s' % which_var 92 sys.exit() 75 93 76 elif myid == 2: 77 time_dir = time_dirs[0] 78 var = vars[1] 94 for which_area in area: 95 if which_area == 'All': 96 easting_min = None 97 easting_max = None 98 northing_min = None 99 northing_max = None 100 else: 101 try: 102 easting_min = eval('project.xmin%s' % which_area) 103 easting_max = eval('project.xmax%s' % which_area) 104 northing_min = eval('project.ymin%s' % which_area) 105 northing_max = eval('project.ymax%s' % which_area) 106 except AttributeError: 107 print 'Unrecognized area name: %s' % which_area 108 break 79 109 80 elif myid == 3: 81 time_dir = time_dirs[1]82 var = vars[0]110 names = [] 111 Name = glob.glob(join(rootDir,event,time_dir,scenario_name) +'*.sww') 112 names.append(Name) 83 113 84 if var not in var_equations: 85 print 'Unrecognized variable name: %s' % which_var 86 sys.exit() 87 88 for which_area in area: 89 if which_area == 'All': 90 easting_min = None 91 easting_max = None 92 northing_min = None 93 northing_max = None 94 else: 95 try: 96 easting_min = eval('project.xmin%s' % which_area) 97 easting_max = eval('project.xmax%s' % which_area) 98 northing_min = eval('project.ymin%s' % which_area) 99 northing_max = eval('project.ymax%s' % which_area) 100 except AttributeError: 101 print 'Unrecognized area name: %s' % which_area 102 break 103 104 names = [] 105 Name = glob.glob(join(rootDir,event,time_dir,scenario_name) +'*.sww') 106 names.append(Name) 107 108 for names in names: 109 110 asc_name = [] 111 print 'start sww2dem: time_dir = %s' % time_dir 112 for name in names: 113 114 test = name 115 name = test[:-4] 116 outname = name + '_' + which_area + '_' + var + str(cellsize) + "m" 117 quantityname = var_equations[var] 118 119 print 'start sww2dem: name = %s' % os.path.basename(name) 120 print outname 121 anuga.sww2dem(name+'.sww', outname+'.asc', 122 quantity = quantityname, 123 cellsize = cellsize, 124 easting_min = easting_min, 125 easting_max = easting_max, 126 northing_min = northing_min, 127 northing_max = northing_max, 128 reduction = timestep, 129 verbose = True) 130 131 asc_name.append(outname + '.asc') 114 for names in names: 132 115 116 asc_name = [] 117 print 'start sww2dem: time_dir = %s' % time_dir 118 for name in names: 119 120 test = name 121 name = test[:-4] 122 outname = name + '_' + which_area + '_' + var + str(cellsize) + "m" 123 quantityname = var_equations[var] 124 125 print 'start sww2dem: name = %s' % os.path.basename(name) 126 print outname 127 anuga.sww2dem(name+'.sww', outname+'.asc', 128 quantity = quantityname, 129 cellsize = cellsize, 130 easting_min = easting_min, 131 easting_max = easting_max, 132 northing_min = northing_min, 133 northing_max = northing_max, 134 reduction = timestep, 135 verbose = True) 136 137 asc_name.append(outname + '.asc') 138 133 139 barrier() 134 140 -
trunk/anuga_work/development/Busselton/Analysis/master.py
r8449 r8452 27 27 wall[32] = "07:00:00" 28 28 29 np = 429 np = 16 30 30 31 31 ScriptName = "/home/547/sfm547/Busselton/Analysis/export_results_max.py" -
trunk/anuga_work/development/Busselton/B1Full/anuga.log
r8449 r8452 1 2012-06-1 6 18:13:52,079INFO anuga.utilities.log:0 |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL2 2012-06-1 6 18:13:52,210INFO polygon:959 |--------------------------------------------------------------------------------3 2012-06-1 6 18:13:52,210INFO polygon:960 |Polygon Max triangle area (m^2) Total area (km^2) Estimated #triangles4 2012-06-1 6 18:13:52,211INFO polygon:962 |--------------------------------------------------------------------------------5 2012-06-1 6 18:13:52,211INFO polygon:976 |Interior 100 56.48 5647626 2012-06-1 6 18:13:52,212INFO polygon:984 |Bounding 10000 352.99 352997 2012-06-1 6 18:13:52,212INFO polygon:989 |Estimated total number of triangles: 8572308 2012-06-1 6 18:13:52,212INFO polygon:990 |Note: This is generally about 20% less than the final amount1 2012-06-18 16:19:21,783 INFO anuga.utilities.log:0 |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL 2 2012-06-18 16:19:21,816 INFO polygon:959 |-------------------------------------------------------------------------------- 3 2012-06-18 16:19:21,816 INFO polygon:960 |Polygon Max triangle area (m^2) Total area (km^2) Estimated #triangles 4 2012-06-18 16:19:21,816 INFO polygon:962 |-------------------------------------------------------------------------------- 5 2012-06-18 16:19:21,817 INFO polygon:976 |Interior 100 56.48 564762 6 2012-06-18 16:19:21,817 INFO polygon:984 |Bounding 10000 352.99 35299 7 2012-06-18 16:19:21,818 INFO polygon:989 |Estimated total number of triangles: 857230 8 2012-06-18 16:19:21,818 INFO polygon:990 |Note: This is generally about 20% less than the final amount -
trunk/anuga_work/development/Busselton/B2Full/anuga.log
r8449 r8452 1 2012-06-1 5 15:13:55,701INFO anuga.utilities.log:0 |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL2 2012-06-1 5 15:13:55,798INFO polygon:959 |--------------------------------------------------------------------------------3 2012-06-1 5 15:13:55,798INFO polygon:960 |Polygon Max triangle area (m^2) Total area (km^2) Estimated #triangles4 2012-06-1 5 15:13:55,798INFO polygon:962 |--------------------------------------------------------------------------------5 2012-06-1 5 15:13:55,802INFO polygon:976 |Interior 100 56.48 5647626 2012-06-1 5 15:13:55,802INFO polygon:984 |Bounding 10000 352.99 352997 2012-06-1 5 15:13:55,802INFO polygon:989 |Estimated total number of triangles: 8572308 2012-06-1 5 15:13:55,802INFO polygon:990 |Note: This is generally about 20% less than the final amount1 2012-06-18 16:31:17,239 INFO anuga.utilities.log:0 |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL 2 2012-06-18 16:31:17,265 INFO polygon:959 |-------------------------------------------------------------------------------- 3 2012-06-18 16:31:17,265 INFO polygon:960 |Polygon Max triangle area (m^2) Total area (km^2) Estimated #triangles 4 2012-06-18 16:31:17,265 INFO polygon:962 |-------------------------------------------------------------------------------- 5 2012-06-18 16:31:17,266 INFO polygon:976 |Interior 100 56.48 564762 6 2012-06-18 16:31:17,266 INFO polygon:984 |Bounding 10000 352.99 35299 7 2012-06-18 16:31:17,266 INFO polygon:989 |Estimated total number of triangles: 857230 8 2012-06-18 16:31:17,266 INFO polygon:990 |Note: This is generally about 20% less than the final amount -
trunk/anuga_work/development/Busselton/B3Full/anuga.log
r8449 r8452 1 2012-06-1 3 09:57:27,185INFO anuga.utilities.log:0 |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL2 2012-06-1 3 09:57:27,296 INFO polygon:959 |--------------------------------------------------------------------------------3 2012-06-1 3 09:57:27,297INFO polygon:960 |Polygon Max triangle area (m^2) Total area (km^2) Estimated #triangles4 2012-06-1 3 09:57:27,297INFO polygon:962 |--------------------------------------------------------------------------------5 2012-06-1 3 09:57:27,298INFO polygon:976 |Interior 100 56.48 5647626 2012-06-1 3 09:57:27,298INFO polygon:984 |Bounding 10000 352.99 352997 2012-06-1 3 09:57:27,298INFO polygon:989 |Estimated total number of triangles: 8572308 2012-06-1 3 09:57:27,299INFO polygon:990 |Note: This is generally about 20% less than the final amount1 2012-06-18 16:31:17,284 INFO anuga.utilities.log:0 |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL 2 2012-06-18 16:31:17,306 INFO polygon:959 |-------------------------------------------------------------------------------- 3 2012-06-18 16:31:17,306 INFO polygon:960 |Polygon Max triangle area (m^2) Total area (km^2) Estimated #triangles 4 2012-06-18 16:31:17,306 INFO polygon:962 |-------------------------------------------------------------------------------- 5 2012-06-18 16:31:17,307 INFO polygon:976 |Interior 100 56.48 564762 6 2012-06-18 16:31:17,307 INFO polygon:984 |Bounding 10000 352.99 35299 7 2012-06-18 16:31:17,307 INFO polygon:989 |Estimated total number of triangles: 857230 8 2012-06-18 16:31:17,307 INFO polygon:990 |Note: This is generally about 20% less than the final amount -
trunk/anuga_work/development/Busselton/B4Full/anuga.log
r8449 r8452 1 2012-06- 16 18:14:42,834INFO anuga.utilities.log:0 |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL2 2012-06- 16 18:14:42,868INFO polygon:959 |--------------------------------------------------------------------------------3 2012-06- 16 18:14:42,868INFO polygon:960 |Polygon Max triangle area (m^2) Total area (km^2) Estimated #triangles4 2012-06- 16 18:14:42,868INFO polygon:962 |--------------------------------------------------------------------------------5 2012-06- 16 18:14:42,869INFO polygon:976 |Interior 100 56.48 5647626 2012-06- 16 18:14:42,870INFO polygon:984 |Bounding 10000 352.99 352997 2012-06- 16 18:14:42,870INFO polygon:989 |Estimated total number of triangles: 8572308 2012-06- 16 18:14:42,870INFO polygon:990 |Note: This is generally about 20% less than the final amount1 2012-06-20 22:53:16,886 INFO anuga.utilities.log:0 |Logfile is './anuga.log' with logging level of INFO, console logging level is CRITICAL 2 2012-06-20 22:53:16,952 INFO polygon:959 |-------------------------------------------------------------------------------- 3 2012-06-20 22:53:16,953 INFO polygon:960 |Polygon Max triangle area (m^2) Total area (km^2) Estimated #triangles 4 2012-06-20 22:53:16,953 INFO polygon:962 |-------------------------------------------------------------------------------- 5 2012-06-20 22:53:16,953 INFO polygon:976 |Interior 100 56.48 564762 6 2012-06-20 22:53:16,953 INFO polygon:984 |Bounding 10000 352.99 35299 7 2012-06-20 22:53:16,954 INFO polygon:989 |Estimated total number of triangles: 857230 8 2012-06-20 22:53:16,954 INFO polygon:990 |Note: This is generally about 20% less than the final amount
Note: See TracChangeset
for help on using the changeset viewer.