Ignore:
Timestamp:
Nov 1, 2007, 11:12:32 AM (17 years ago)
Author:
ole
Message:

Cosmetics according to style guide.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/shallow_water/data_manager.py

    r4775 r4776  
    12901290    """
    12911291
    1292     #FIXME: Can this be written feasibly using write_pts?
     1292    # FIXME: Can this be written feasibly using write_pts?
    12931293
    12941294    import os
     
    12981298    root = basename_in
    12991299
    1300     #Get NetCDF
    1301     infile = NetCDFFile(root + '.dem', 'r')  #Open existing netcdf file for read
     1300    # Get NetCDF
     1301    infile = NetCDFFile(root + '.dem', 'r')  # Open existing netcdf file for read
    13021302
    13031303    if verbose: print 'Reading DEM from %s' %(root + '.dem')
     
    13051305    ncols = infile.ncols[0]
    13061306    nrows = infile.nrows[0]
    1307     xllcorner = infile.xllcorner[0]  #Easting of lower left corner
    1308     yllcorner = infile.yllcorner[0]  #Northing of lower left corner
     1307    xllcorner = infile.xllcorner[0]  # Easting of lower left corner
     1308    yllcorner = infile.yllcorner[0]  # Northing of lower left corner
    13091309    cellsize = infile.cellsize[0]
    13101310    NODATA_value = infile.NODATA_value[0]
     
    13151315    false_northing = infile.false_northing[0]
    13161316
    1317     #Text strings
     1317    # Text strings
    13181318    projection = infile.projection
    13191319    datum = infile.datum
     
    13211321
    13221322
    1323     #Get output file
     1323    # Get output file
    13241324    if basename_out == None:
    13251325        ptsname = root + '.pts'
     
    13311331    outfile = NetCDFFile(ptsname, 'w')
    13321332
    1333     #Create new file
     1333    # Create new file
    13341334    outfile.institution = 'Geoscience Australia'
    13351335    outfile.description = 'NetCDF pts format for compact and portable storage ' +\
    13361336                      'of spatial point data'
    1337     #assign default values
     1337    # Assign default values
    13381338    if easting_min is None: easting_min = xllcorner
    13391339    if easting_max is None: easting_max = xllcorner + ncols*cellsize
     
    13411341    if northing_max is None: northing_max = yllcorner + nrows*cellsize
    13421342
    1343     #compute offsets to update georeferencing
     1343    # Compute offsets to update georeferencing
    13441344    easting_offset = xllcorner - easting_min
    13451345    northing_offset = yllcorner - northing_min
    13461346
    1347     #Georeferencing
     1347    # Georeferencing
    13481348    outfile.zone = zone
    1349     outfile.xllcorner = easting_min #Easting of lower left corner
    1350     outfile.yllcorner = northing_min #Northing of lower left corner
     1349    outfile.xllcorner = easting_min # Easting of lower left corner
     1350    outfile.yllcorner = northing_min # Northing of lower left corner
    13511351    outfile.false_easting = false_easting
    13521352    outfile.false_northing = false_northing
     
    13571357
    13581358
    1359     #Grid info (FIXME: probably not going to be used, but heck)
     1359    # Grid info (FIXME: probably not going to be used, but heck)
    13601360    outfile.ncols = ncols
    13611361    outfile.nrows = nrows
     
    13641364    totalnopoints = nrows*ncols
    13651365
    1366     # calculating number of NODATA_values for each row in clipped region
    1367     #FIXME: use array operations to do faster
     1366    # Calculating number of NODATA_values for each row in clipped region
     1367    # FIXME: use array operations to do faster
    13681368    nn = 0
    13691369    k = 0
     
    13901390    index2 = thisj
    13911391
    1392     # dimension definitions
     1392    # Dimension definitions
    13931393    nrows_in_bounding_box = int(round((northing_max-northing_min)/cellsize))
    13941394    ncols_in_bounding_box = int(round((easting_max-easting_min)/cellsize))
     
    14071407    outfile.createDimension('number_of_dimensions', 2) #This is 2d data
    14081408
    1409     # variable definitions
     1409    # Variable definitions
    14101410    outfile.createVariable('points', Float, ('number_of_points',
    14111411                                             'number_of_dimensions'))
     
    14171417
    14181418    lenv = index2-index1+1
    1419     #Store data
     1419    # Store data
    14201420    global_index = 0
    1421     #for i in range(nrows):
     1421    # for i in range(nrows):
    14221422    for i in range(i1_0,thisi+1,1):
    14231423        if verbose and i%((nrows+10)/10)==0:
     
    14291429        no_NODATA = sum(v == NODATA_value)
    14301430        if no_NODATA > 0:
    1431             newcols = lenv - no_NODATA #ncols_in_bounding_box - no_NODATA
     1431            newcols = lenv - no_NODATA # ncols_in_bounding_box - no_NODATA
    14321432        else:
    1433             newcols = lenv #ncols_in_bounding_box
     1433            newcols = lenv # ncols_in_bounding_box
    14341434
    14351435        telev = zeros(newcols, Float)
Note: See TracChangeset for help on using the changeset viewer.