Changeset 8273
- Timestamp:
- Dec 7, 2011, 5:02:30 PM (13 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_core/source/anuga/utilities/sww_merge.py
r8270 r8273 11 11 from anuga.file.sww import SWW_file, Write_sww 12 12 13 def sww_merge(domain_global_name, np, verbose=False): 13 14 14 def sww_merge(swwfiles, output, verbose = False): 15 output = domain_global_name+".sww" 16 swwfiles = [ domain_global_name+"_P"+str(v)+"_"+str(np)+".sww" for v in range(np)] 17 18 _sww_merge(swwfiles, output, verbose) 19 20 21 def _sww_merge(swwfiles, output, verbose): 15 22 """ 16 23 Merge a list of sww files into a single file. … … 28 35 verbose True to log output information 29 36 """ 30 37 38 if verbose: 39 print "MERGING SWW Files" 40 31 41 static_quantities = ['elevation'] 32 42 dynamic_quantities = ['stage', 'xmomentum', 'ymomentum'] … … 88 98 89 99 # Write out the SWW file 100 101 if verbose: 102 print 'Writing file ', output, ':' 90 103 fido = NetCDFFile(output, netcdf_mode_w) 91 104 sww = Write_sww(static_quantities, dynamic_quantities) … … 130 143 help='number of processors used to produce sww files') 131 144 parser.add_argument('-f', type=str, default="domain", 132 help=' base sww filename')145 help='domain global name') 133 146 parser.add_argument('-v', nargs='?', type=bool, const=True, default=False, 134 147 help='verbosity') … … 137 150 138 151 np = args.np 139 filebase = args.f152 domain_global_name = args.f 140 153 verbose = args.v 141 154 142 #print np143 #print filebase144 #print verbose145 146 147 148 output = filebase+".sww"149 swwfiles = [ filebase+"_P"+str(v)+"_"+str(np)+".sww" for v in range(np)]150 151 #print swwfiles152 155 153 156 try: 154 sww_merge( swwfiles, output, verbose)157 sww_merge(domain_global_name, np, verbose) 155 158 except: 156 msg = 'ERROR: When merging sww files '+" ".join(swwfiles)159 msg = 'ERROR: When merging sww files %s '% domain_global_name 157 160 print msg 158 161 raise
Note: See TracChangeset
for help on using the changeset viewer.