- Timestamp:
- May 6, 2009, 12:06:04 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/australia_ph2/sydney/export_results.py
r6926 r6982 29 29 ##cellsize = 5 30 30 31 #timestep = None # None means no timestep!32 timestep = 031 timestep = None # None means no timestep! 32 #timestep = 0 33 33 34 34 ###### … … 50 50 51 51 # one or more key strings from var_equations above 52 var = [' elevation']52 var = ['stage'] 53 53 54 54 ###### … … 79 79 for time_dir in time_dirs: 80 80 81 name1 = join(directory, time_dir, project.scenario_name) 82 ## name2 = join(directory, time_dir, project.scenario_name)+'_time_21060_0' 83 ## name3 = join(directory, time_dir, project.scenario_name)+'_time_42120_0' 84 ## name4 = join(directory, time_dir, project.scenario_name)+'_time_63180_0' 85 ## name5 = join(directory, time_dir, project.scenario_name)+'_time_53280_0' 86 ## name6 = join(directory, time_dir, project.scenario_name)+'_time_66600_0' 87 ## name7 = join(directory, time_dir, project.scenario_name)+'_time_79920_0' 81 name = join(directory, time_dir, project.scenario_name) 88 82 89 names = [name1]#, name2, name3,name4]#,name5,name6,name7] 90 91 asc_name = [] 92 93 for name in names: 94 95 outname = name + '_' + which_area + '_' + which_var 96 quantityname = var_equations[which_var] 83 outname = name + '_' + which_area + '_' + which_var 84 quantityname = var_equations[which_var] 85 86 print 'start sww2dem: time_dir=%s' % time_dir 87 88 sww2dem(name, basename_out = outname, 89 quantity = quantityname, 90 timestep = timestep, 91 cellsize = cellsize, 92 easting_min = easting_min, 93 easting_max = easting_max, 94 northing_min = northing_min, 95 northing_max = northing_max, 96 reduction = max, 97 verbose = True, 98 format = 'asc') 97 99 98 print 'start sww2dem: time_dir=%s' % time_dir99 100 sww2dem(name, basename_out = outname,101 quantity = quantityname,102 timestep = timestep,103 cellsize = cellsize,104 easting_min = easting_min,105 easting_max = easting_max,106 northing_min = northing_min,107 northing_max = northing_max,108 reduction = max,109 verbose = True,110 format = 'asc')111 112 ## asc_name.append(outname + '.asc')113 ##114 ## maxasc_outname = join(directory, time_dir, project.scenario_name)+'_'+which_area+'_'+which_var+'_max.asc'115 ##116 ## print 'max asc outname ', maxasc_outname117 ## print 'asc_name ', str(asc_name)118 ##119 ## MaxAsc(maxasc_outname, asc_name)
Note: See TracChangeset
for help on using the changeset viewer.