Changeset 8277


Ignore:
Timestamp:
Dec 8, 2011, 6:40:28 PM (8 years ago)
Author:
steve
Message:

Trying to fix sww_merge.py

Location:
trunk/anuga_core/source
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/source/anuga/utilities/sww_merge.py

    r8273 r8277  
    77
    88from Scientific.IO.NetCDF import NetCDFFile
    9 from anuga.config import netcdf_mode_r, netcdf_mode_w, \
    10                             netcdf_mode_a, netcdf_float
     9from anuga.config import netcdf_mode_r, netcdf_mode_w, netcdf_mode_a
     10from anuga.config import netcdf_float, netcdf_float32, netcdf_int
    1111from anuga.file.sww import SWW_file, Write_sww
    1212
     
    4949   
    5050        fid = NetCDFFile(filename, netcdf_mode_r)
     51
     52       
     53       
    5154        tris = fid.variables['volumes'][:]       
    5255         
     
    5861            out_s_quantities = {}
    5962            out_d_quantities = {}
     63
     64
     65            xllcorner = fid.xllcorner
     66            yllcorner = fid.yllcorner
     67
     68            print 'zone',fid.zone
     69
     70            print """
     71            fid.order
     72            fid.xllcorner;
     73            fid.yllcorner ;
     74            fid.zone;
     75            fid.false_easting;
     76            fid.false_northing;
     77            fid.datum;
     78            fid.projection;
     79            """
     80
     81
     82            print fid.order
     83            print fid.xllcorner;
     84            print fid.yllcorner ;
     85            print fid.zone;
     86            print fid.false_easting;
     87            print fid.false_northing;
     88            print fid.datum;
     89            print fid.projection;
     90
     91            order      = fid.order
     92            xllcorner  = fid.xllcorner;
     93            yllcorner  = fid.yllcorner ;
     94            zone       = fid.zone;
     95            false_easting  = fid.false_easting;
     96            false_northing = fid.false_northing;
     97            datum      = fid.datum;
     98            projection = fid.projection;
     99
    60100           
    61101            for quantity in static_quantities:
     
    95135    # Mash all points into a single big list   
    96136    points = [[xx, yy] for xx, yy in zip(x, y)]
     137
     138    points = num.asarray(points, dtype='f')
     139    print points
    97140    fid.close()
    98    
     141
     142    #---------------------------
    99143    # Write out the SWW file
     144    #---------------------------
    100145
    101146    if verbose:
     
    107152                             len(points),
    108153                             description=description,
    109                              sww_precision=netcdf_float)
    110 
    111 
     154                             sww_precision=netcdf_float32)
     155
     156
     157   
     158   
    112159    sww.store_triangulation(fido, points, out_tris)
     160
     161    fido.order      = order
     162    fido.xllcorner  = xllcorner;
     163    fido.yllcorner  = yllcorner ;
     164    fido.zone       = zone;
     165    fido.false_easting  = false_easting;
     166    fido.false_northing = false_northing;
     167    fido.datum      = datum;
     168    fido.projection = projection;
    113169       
    114170    sww.store_static_quantities(fido, verbose=verbose, **out_s_quantities)
     
    118174        q_values = out_d_quantities[q]
    119175        for i, time_slice in enumerate(q_values):
    120             fido.variables[q][i] = num.array(time_slice, netcdf_float)
     176            fido.variables[q][i] = num.array(time_slice, netcdf_float32)
    121177       
    122178        # This updates the _range values
  • trunk/anuga_core/source/anuga_parallel/run_parallel_sw_merimbula.py

    r8272 r8277  
    7878
    7979if myid == 0 and verbose: print 'DISTRIBUTING DOMAIN'
    80 domain = distribute(domain,verbose=True)
     80domain = distribute(domain)
    8181
    8282#domain.smooth = False
     
    119119# Merge the individual sww files into one file
    120120#--------------------------------------------------
    121 domain.sww_merge(verbose=True)
     121domain.sww_merge()
    122122
    123123finalize()
Note: See TracChangeset for help on using the changeset viewer.