anuga 1.3.5__zip

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/__config__.py

@@ -0,0 +1,21 @@
+# This file is generated by C:\Users\steve\anuga\anuga_core\setup.py
+# It contains system_info results at the time of building this package.
+__all__ = ["get_info","show"]
+
+
+def get_info(name):
+    g = globals()
+    return g.get(name, g.get(name + "_info", {}))
+
+def show():
+    for name,info_dict in globals().items():
+        if name[0] == "_" or type(info_dict) is not type({}): continue
+        print(name + ":")
+        if not info_dict:
+            print("  NOT AVAILABLE")
+        for k,v in info_dict.items():
+            v = str(v)
+            if k == "sources" and len(v) > 200:
+                v = v[:60] + " ...\n... " + v[-60:]
+            print("    %s = %s" % (k,v))
+    

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/__init__.py

@@ -0,0 +1,363 @@
+""" ANUGA models the effect of tsunamis and flooding upon a terrain mesh.
+    In typical usage, a Domain class is created for a particular piece of
+    terrain. Boundary conditions are specified for the domain, such as inflow
+    and outflow, and then the simulation is run.
+
+    This is the public API to ANUGA. It provides a toolkit of often-used
+    modules, which can be used directly by including the following line in
+    the user's code:
+
+    import anuga
+        
+    This usage pattern abstracts away the internal heirarchy of the ANUGA
+    system, allowing the user to concentrate on writing simulations without
+    searching through the ANUGA source tree for the functions that they need.
+    
+    Also, it isolates the user from "under-the-hood" refactorings.
+"""
+
+#-----------------------------------------------------
+# Make selected classes available directly
+#-----------------------------------------------------
+
+from numpy.testing import Tester
+test = Tester().test
+
+
+# PEP0440 compatible formatted version, see:
+# https://www.python.org/dev/peps/pep-0440/
+#
+# Generic release markers:
+# X.Y
+# X.Y.Z # For bugfix releases
+#
+# Admissible pre-release markers:
+# X.YaN # Alpha release
+# X.YbN # Beta release
+# X.YrcN # Release Candidate
+# X.Y # Final release
+#
+# Dev branch marker is: 'X.Y.dev' or 'X.Y.devN' where N is an integer.
+#
+__version__ = '1.3.5'
+
+__svn_revision__ = filter(str.isdigit, "$Revision: 9641 $")
+
+#from anuga.__metadata__ import  __date__, __author__
+
+#from .version import git_revision as __git_revision__
+#from .version import svn_revision as __svn_revision__
+
+
+
+
+
+# We first need to detect if we're being called as part of the anuga setup
+# procedure itself in a reliable manner.
+try:
+    __ANUGA_SETUP__
+except NameError:
+    __ANUGA_SETUP__ = False
+
+    
+if __ANUGA_SETUP__:
+    import sys as _sys
+    _sys.stderr.write('Running from anuga source directory.\n')
+    del _sys
+else:
+
+    try:
+        from anuga.__config__ import show as show_config
+    except ImportError:
+        msg = """Error importing anuga: you should not try to import anuga from
+        its source directory; please exit the anuga source tree, and relaunch
+        your python interpreter from there."""
+        raise ImportError(msg)
+       
+    #--------------------------------
+    # Important basic classes
+    #--------------------------------
+    from anuga.shallow_water.shallow_water_domain import Domain
+    from anuga.abstract_2d_finite_volumes.quantity import Quantity
+    from anuga.abstract_2d_finite_volumes.region import Region
+    from anuga.geospatial_data.geospatial_data import Geospatial_data
+    from anuga.operators.base_operator import Operator
+    from anuga.structures.structure_operator import Structure_operator
+
+
+    from anuga.abstract_2d_finite_volumes.generic_domain import Generic_Domain
+    from anuga.abstract_2d_finite_volumes.neighbour_mesh import Mesh
+    #------------------------------------------------------------------------------ 
+    # Miscellaneous
+    #------------------------------------------------------------------------------ 
+    from anuga.abstract_2d_finite_volumes.util import file_function, \
+                                            sww2timeseries, sww2csv_gauges, \
+                                            csv2timeseries_graphs
+
+    from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross, \
+                                                        rectangular
+
+    from anuga.file.csv_file import load_csv_as_building_polygons,  \
+                                    load_csv_as_polygons
+
+    from anuga.file.sts import create_sts_boundary
+
+    from anuga.file.ungenerate import load_ungenerate
+
+    from anuga.geometry.polygon import read_polygon
+    from anuga.geometry.polygon import plot_polygons
+    from anuga.geometry.polygon import inside_polygon
+    from anuga.geometry.polygon import polygon_area
+    from anuga.geometry.polygon_function import Polygon_function
+
+    from anuga.abstract_2d_finite_volumes.pmesh2domain import \
+                                                pmesh_to_domain_instance
+
+    from anuga.utilities.system_tools import file_length
+    from anuga.utilities.sww_merge import sww_merge_parallel as sww_merge
+    from anuga.utilities.file_utils import copy_code_files
+    from anuga.utilities.numerical_tools import safe_acos as acos
+    import anuga.utilities.plot_utils as plot_utils
+
+
+    from anuga.caching import cache
+    from os.path import join
+    from anuga.config import indent
+
+
+
+    #----------------------------
+    # Parallel api 
+    #----------------------------
+    ## from anuga_parallel.parallel_api import distribute
+    ## from anuga_parallel.parallel_api import myid, numprocs, get_processor_name
+    ## from anuga_parallel.parallel_api import send, receive
+    ## from anuga_parallel.parallel_api import pypar_available, barrier, finalize
+
+    ## if pypar_available:
+    ##     from anuga_parallel.parallel_api import sequential_distribute_dump
+    ##     from anuga_parallel.parallel_api import sequential_distribute_load
+
+    from anuga.parallel.parallel_api import distribute
+    from anuga.parallel.parallel_api import myid, numprocs, get_processor_name
+    from anuga.parallel.parallel_api import send, receive
+    from anuga.parallel.parallel_api import pypar_available, barrier, finalize
+
+    if pypar_available:
+        from anuga.parallel.parallel_api import sequential_distribute_dump
+        from anuga.parallel.parallel_api import sequential_distribute_load
+
+
+    #-----------------------------
+    # Checkpointing
+    #-----------------------------
+    from anuga.shallow_water.checkpoint import load_checkpoint_file
+
+
+    #-----------------------------
+    # SwW Standard Boundaries
+    #-----------------------------
+    from anuga.shallow_water.boundaries import File_boundary
+    from anuga.shallow_water.boundaries import Reflective_boundary
+    from anuga.shallow_water.boundaries import Field_boundary
+    from anuga.shallow_water.boundaries import \
+                        Time_stage_zero_momentum_boundary
+    from anuga.shallow_water.boundaries import \
+                        Transmissive_stage_zero_momentum_boundary
+    from anuga.shallow_water.boundaries import \
+                        Transmissive_momentum_set_stage_boundary
+    from anuga.shallow_water.boundaries import \
+                        Transmissive_n_momentum_zero_t_momentum_set_stage_boundary
+    from anuga.shallow_water.boundaries import \
+                        Flather_external_stage_zero_velocity_boundary
+    from anuga.abstract_2d_finite_volumes.generic_boundary_conditions import \
+                        Compute_fluxes_boundary
+
+
+    #-----------------------------
+    # General Boundaries
+    #-----------------------------
+    from anuga.abstract_2d_finite_volumes.generic_boundary_conditions \
+                                import Dirichlet_boundary
+    from anuga.abstract_2d_finite_volumes.generic_boundary_conditions \
+                                import Time_boundary
+    from anuga.abstract_2d_finite_volumes.generic_boundary_conditions \
+                                import Time_space_boundary
+    from anuga.abstract_2d_finite_volumes.generic_boundary_conditions \
+                                import Transmissive_boundary
+
+
+
+    #-----------------------------
+    # Shallow Water Tsunamis
+    #-----------------------------
+
+    from anuga.tsunami_source.smf import slide_tsunami, slump_tsunami
+
+
+
+    #-----------------------------
+    # Forcing
+    # These are old, should use operators
+    #-----------------------------
+    from anuga.shallow_water.forcing import Inflow, Rainfall, Wind_stress
+
+
+
+    #-----------------------------
+    # File conversion utilities
+    #-----------------------------
+    from anuga.file_conversion.file_conversion import sww2obj, \
+                        timefile2netcdf, tsh2sww
+    from anuga.file_conversion.urs2nc import urs2nc
+    from anuga.file_conversion.urs2sww import urs2sww  
+    from anuga.file_conversion.urs2sts import urs2sts
+    from anuga.file_conversion.dem2pts import dem2pts                    
+    from anuga.file_conversion.esri2sww import esri2sww   
+    from anuga.file_conversion.sww2dem import sww2dem, sww2dem_batch 
+    from anuga.file_conversion.asc2dem import asc2dem
+    from anuga.file_conversion.xya2pts import xya2pts     
+    from anuga.file_conversion.ferret2sww import ferret2sww     
+    from anuga.file_conversion.dem2dem import dem2dem
+    from anuga.file_conversion.sww2array import sww2array
+
+    #-----------------------------
+    # Parsing arguments
+    #-----------------------------
+    from anuga.utilities.argparsing import create_standard_parser
+    from anuga.utilities.argparsing import parse_standard_args
+
+
+    def get_args():
+        """ Explicitly parse the argument list using standard anuga arguments
+
+        Don't use this if you want to setup your own parser
+        """
+        parser = create_standard_parser()
+        return parser.parse_args()
+
+
+    #-----------------------------
+    # Running Script
+    #-----------------------------
+    from anuga.utilities.run_anuga_script import run_script as run_anuga_script
+
+
+    #-----------------------------
+    # Mesh API
+    #-----------------------------
+    from anuga.pmesh.mesh_interface import create_mesh_from_regions
+
+    #-----------------------------
+    # SWW file access
+    #-----------------------------
+    from anuga.shallow_water.sww_interrogate import get_flow_through_cross_section
+
+    #---------------------------
+    # Operators
+    #---------------------------
+    from anuga.operators.kinematic_viscosity_operator import Kinematic_viscosity_operator
+
+    from anuga.operators.rate_operators import Rate_operator
+    from anuga.operators.set_friction_operators import Depth_friction_operator 
+
+    from anuga.operators.set_elevation_operator import Set_elevation_operator
+    from anuga.operators.set_quantity_operator import Set_quantity_operator
+    from anuga.operators.set_stage_operator import Set_stage_operator
+
+    from anuga.operators.set_elevation import Set_elevation
+    from anuga.operators.set_quantity import Set_quantity
+
+    from anuga.operators.erosion_operators import Bed_shear_erosion_operator
+    from anuga.operators.erosion_operators import Flat_slice_erosion_operator
+    from anuga.operators.erosion_operators import Flat_fill_slice_erosion_operator
+
+    #---------------------------
+    # Structure Operators
+    #---------------------------
+
+
+    if pypar_available:
+        from anuga.parallel.parallel_operator_factory import Inlet_operator
+        from anuga.parallel.parallel_operator_factory import Boyd_box_operator
+        from anuga.parallel.parallel_operator_factory import Boyd_pipe_operator
+        from anuga.parallel.parallel_operator_factory import Weir_orifice_trapezoid_operator
+        from anuga.parallel.parallel_operator_factory import Internal_boundary_operator
+    else:
+        from anuga.structures.inlet_operator import Inlet_operator
+        from anuga.structures.boyd_box_operator import Boyd_box_operator
+        from anuga.structures.boyd_pipe_operator import Boyd_pipe_operator
+        from anuga.structures.weir_orifice_trapezoid_operator import Weir_orifice_trapezoid_operator
+        from anuga.structures.internal_boundary_operator import Internal_boundary_operator
+
+
+    #----------------------------
+    # Parallel distribute
+    #----------------------------
+
+
+    #----------------------------
+    # 
+    #Added by Petar Milevski 10/09/2013
+    #import time, os
+
+    from anuga.utilities.model_tools import get_polygon_from_single_file
+    from anuga.utilities.model_tools import get_polygons_from_Mid_Mif
+    from anuga.utilities.model_tools import get_polygon_list_from_files
+    from anuga.utilities.model_tools import get_polygon_dictionary
+    from anuga.utilities.model_tools import get_polygon_value_list
+    from anuga.utilities.model_tools import read_polygon_dir
+    from anuga.utilities.model_tools import read_hole_dir_multi_files_with_single_poly
+    from anuga.utilities.model_tools import read_multi_poly_file
+    from anuga.utilities.model_tools import read_hole_dir_single_file_with_multi_poly
+    from anuga.utilities.model_tools import read_multi_poly_file_value
+    from anuga.utilities.model_tools import Create_culvert_bridge_Operator
+
+
+    #---------------------------
+    # User Access Functions
+    #---------------------------
+
+    from anuga.utilities.system_tools import get_user_name, get_host_name, \
+        get_revision_number
+    from anuga.utilities.mem_time_equation import estimate_time_mem
+
+
+    #-------------------------
+    # create domain functions
+    #-------------------------
+    from anuga.extras import create_domain_from_regions
+    from anuga.extras import create_domain_from_file
+    from anuga.extras import rectangular_cross_domain
+
+    
+    #import logging as log
+    from anuga.utilities import log
+
+    from anuga.config import g
+    from anuga.config import velocity_protection
+    
+## if use_psyco:
+##     # try using psyco if available
+##     try:
+##         import psyco
+##     except:
+##         import os
+##         import sys
+##         if os.name == 'posix' and os.uname()[4] in ['x86_64', 'ia64']:
+##             pass
+##             # Psyco isn't supported on 64 bit systems, but it doesn't matter
+##         elif sys.version[:3] == '2.7' :
+##             pass
+##             # Psyco isn't available for python 2.7 (16/05/2011)
+##         else:
+##             log.critical('WARNING: psyco (speedup) could not be imported, '
+##                          'you may want to consider installing it')
+##     else:
+##         psyco.full() # aggressively compile everything
+##         #psyco.background() # attempt to profile code - only compile most used
+        
+
+
+
+

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/__metadata__.py

@@ -0,0 +1,6 @@
+# Meta data for anuga
+__version__ = '1.3'
+__date__ = 'October 2011'
+__author__ =  'Ole Nielsen, Stephen Roberts, Duncan Gray'
+
+

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/abstract_2d_finite_volumes/__init__.py

@@ -0,0 +1,17 @@
+"""
+    A set of modules which is used to build a geometrical representation
+    of the test scenario domain.
+    
+    It contains functions for fitting data across an arbitrary terrain height
+    mesh. 
+"""
+
+from numpy.testing import Tester
+test = Tester().test
+
+
+
+
+
+
+

Users/steve/Miniconda/envs/anuga/Lib/site-packages/anuga/abstract_2d_finite_volumes/ermapper_grids.py

@@ -0,0 +1,251 @@
+#!/usr/bin/env python
+
+import numpy as num
+
+celltype_map = {'IEEE4ByteReal': num.float32, 'IEEE8ByteReal': num.float64}
+
+
+def write_ermapper_grid(ofile, data, header = {}):
+    """
+    write_ermapper_grid(ofile, data, header = {}):
+    
+    Function to write a 2D numeric array to an ERMapper grid.  There are a series of conventions adopted within
+    this code, specifically:
+    1)  The registration coordinate for the data is the SW (or lower-left) corner of the data
+    2)  The registration coordinates refer to cell centres
+    3)  The data is a 2D numeric array with the NW-most data in element (0,0) and the SE-most data in element (N,M)
+        where N is the last line and M is the last column
+    4)  There has been no testng of the use of a rotated grid.  Best to keep data in an NS orientation
+    
+    Input Parameters:
+    ofile:      string - filename for output (note the output will consist of two files
+                ofile and ofile.ers.  Either of these can be entered into this function
+    data:       array - 2D array containing the data to be output to the grid
+    header:     dictionary - contains spatial information about the grid, in particular:
+                    header['datum'] datum for the data ('"GDA94"')
+                    header['projection'] - either '"GEOGRAPHIC"' or '"PROJECTED"' 
+                    header['coordinatetype'] - either 'EN' (for eastings/northings) or
+                                                      'LL' (for lat/long)
+                    header['rotation'] - rotation of grid ('0:0:0.0')
+                    header['celltype'] - data type for writing data ('IEEE4ByteReal')
+                    header['nullcellvalue'] - value for null cells ('-99999')
+                    header['xdimension'] - cell size in x-dir in units dictated by 'coordinatetype' ('100')
+                    header['registrationcellx'] == '0'
+                    header['ydimension'] - cell size in y-dir in units dictated by 'coordinatetype' ('100')
+                    header['longitude'] - co-ordinate of registration cell ('0:0:0')
+                    header['latitude'] - co-ordinate of registration line ('0:0:0')
+                    header['nrofbands'] - number of bands ('1')
+                    header['value'] - name of grid ('"Default_Band"')
+
+                    Some entries are determined automatically from the data
+                    header['nroflines'] - number of lines in data
+                    header['nrofcellsperline'] - number of columns in data
+                    header['registrationcelly'] == last line of data
+
+    Written by Trevor Dhu, Geoscience Australia 2005 
+    """
+    # extract filenames for header and data files from ofile
+    ers_index = ofile.find('.ers')
+    if ers_index > 0:
+        data_file = ofile[0:ers_index]
+        header_file = ofile
+    else:
+        data_file = ofile
+        header_file = ofile + '.ers'
+
+
+    # Check that the data is a 2 dimensional array
+    data_size = num.shape(data)
+    assert len(data_size) == 2
+    
+    header['nroflines'] = str(data_size[0])
+    header['nrofcellsperline'] = str(data_size[1])
+
+
+    header = create_default_header(header)
+    write_ermapper_header(header_file, header)
+    write_ermapper_data(data, data_file, data_format = header['celltype'])
+
+    
+def read_ermapper_grid(ifile):
+    ers_index = ifile.find('.ers')
+    if ers_index > 0:
+        data_file = ifile[0:ers_index]
+        header_file = ifile
+    else:
+        data_file = ifile
+        header_file = ifile + '.ers'
+        
+    header = read_ermapper_header(header_file)
+
+    nroflines = int(header['nroflines'])
+    nrofcellsperlines = int(header['nrofcellsperline'])
+    data = read_ermapper_data(data_file)
+    data = num.reshape(data,(nroflines,nrofcellsperlines))
+    return data
+    
+
+def write_ermapper_header(ofile, header = {}):
+    
+    header = create_default_header(header)
+    # Determine if the dataset is in lats/longs or eastings/northings and set header parameters
+    # accordingly
+    if header['coordinatetype'] == 'LL':
+        X_Class = 'Longitude'
+        Y_Class = 'Latitude'
+    elif header['coordinatetype'] == 'EN':
+        X_Class = 'Eastings'
+        Y_Class = 'Northings'
+
+    # open the header file for writing to
+    fid = open(ofile,'wt')
+
+    # Begin writing the header
+    fid.write('DatasetHeader Begin\n')
+    fid.write('\tVersion\t\t= "6.4"\n')
+    fid.write('\tDatasetType\t= ERStorage\n')
+    fid.write('\tDataType\t= Raster\n')
+    fid.write('\tByteOrder\t= LSBFirst\n')
+
+    # Write the coordinate space information
+    fid.write('\tCoordinateSpace Begin\n')
+    fid.write('\t\tDatum\t\t\t = ' + header['datum'] + '\n')
+    fid.write('\t\tProjection\t\t = ' + header['projection'] + '\n')
+    fid.write('\t\tCoordinateType\t = ' + header['coordinatetype'] + '\n')
+    fid.write('\t\tRotation\t\t = ' + header['rotation'] + '\n')
+    fid.write('\t\tUnits\t\t = ' + header['units'] + '\n')
+    fid.write('\tCoordinateSpace End\n')
+
+    # Write the raster information
+    fid.write('\tRasterInfo Begin\n')
+    fid.write('\t\tCellType\t\t\t = ' + header['celltype'] + '\n')
+    fid.write('\t\tNullCellValue\t\t = ' + header['nullcellvalue'] + '\n')
+    fid.write('\t\tRegistrationCellX\t\t = ' + header['registrationcellx'] +'\n')
+    fid.write('\t\tRegistrationCellY\t\t = ' + header['registrationcelly'] +'\n')
+    # Write the cellsize information
+    fid.write('\t\tCellInfo Begin\n')
+    fid.write('\t\t\tXDimension\t\t\t = ' + header['xdimension'] + '\n')
+    fid.write('\t\t\tYDimension\t\t\t = ' + header['ydimension'] + '\n')
+    fid.write('\t\tCellInfo End\n')
+    # Continue with wrting the raster information
+    fid.write('\t\tNrOfLines\t\t\t = ' + header['nroflines'] + '\n')
+    fid.write('\t\tNrOfCellsPerLine\t = ' + header['nrofcellsperline'] + '\n')
+    # Write the registration coordinate information
+    fid.write('\t\tRegistrationCoord Begin\n')
+    fid.write('\t\t\t' + X_Class + '\t\t\t = ' + header[X_Class.lower()] + '\n')
+    fid.write('\t\t\t' + Y_Class + '\t\t\t = ' + header[Y_Class.lower()] + '\n')
+    fid.write('\t\tRegistrationCoord End\n')
+    # Continue with wrting the raster information
+    fid.write('\t\tNrOfBands\t\t\t = ' + header['nrofbands'] + '\n')
+    fid.write('\t\tBandID Begin\n')
+    fid.write('\t\t\tValue\t\t\t\t = ' + header['value'] + '\n')
+    fid.write('\t\tBandID End\n')
+    fid.write('\tRasterInfo End\n')
+    fid.write('DatasetHeader End\n')
+    
+    fid.close
+
+def read_ermapper_header(ifile):
+    # function for reading an ERMapper header from file
+    header = {}
+
+    fid = open(ifile,'rt')
+    header_string = fid.readlines()
+    fid.close()
+
+    for line in header_string:
+        if line.find('=') > 0:
+            tmp_string = line.strip().split('=')
+            header[tmp_string[0].strip().lower()]= tmp_string[1].strip()
+
+    return header                      
+
+def write_ermapper_data(grid, ofile, data_format=num.float32):
+
+
+    try:
+        data_format = celltype_map[data_format]
+    except:
+        pass
+    
+
+    #if isinstance(data_format, basestring):
+    #    #celltype_map is defined at top of code
+    #    if celltype_map.has_key(data_format):
+    #        data_format = celltype_map[data_format]
+    #    else:
+    #        msg = 'Format %s is not yet defined by celltype_map' %data_format
+    #        raise Exception(msg)
+        
+    
+    # Convert the array to data_format (default format is Float32)
+    grid_as_float = grid.astype(data_format)
+
+    # Convert array to a string for writing to output file
+    output_string = grid_as_float.tostring()
+
+    # open output file in a binary format and write the output string
+    fid = open(ofile,'wb')
+    fid.write(output_string)
+    fid.close()
+
+
+def read_ermapper_data(ifile, data_format = num.float32):
+    # open input file in a binary format and read the input string
+    fid = open(ifile,'rb')
+    input_string = fid.read()
+    fid.close()
+
+    # convert input string to required format (Note default format is num.float32)
+    grid_as_float = num.fromstring(input_string,data_format)
+    return grid_as_float
+
+def create_default_header(header = {}):
+    # fill any blanks in a header dictionary with default values
+    # input parameters:
+    # header:   a dictionary containing fields that are not meant
+    #           to be filled with default values
+    
+ 
+    if not header.has_key('datum'):
+        header['datum'] = '"GDA94"'
+    if not header.has_key('projection'):
+        header['projection'] = '"GEOGRAPHIC"'        
+    if not header.has_key('coordinatetype'):
+        header['coordinatetype'] = 'LL'
+    if not header.has_key('rotation'):
+        header['rotation'] = '0:0:0.0'
+    if not header.has_key('units'):
+        header['units'] = '"METERS"'
+    if not header.has_key('celltype'):
+        header['celltype'] = 'IEEE4ByteReal'
+    if not header.has_key('nullcellvalue'):
+        header['nullcellvalue'] = '-99999'
+    if not header.has_key('xdimension'):
+        header['xdimension'] = '100'
+    if not header.has_key('latitude'):
+        header['latitude'] = '0:0:0'
+    if not header.has_key('longitude'):
+        header['longitude'] = '0:0:0'        
+    if not header.has_key('ydimension'):
+        header['ydimension'] = '100'
+    if not header.has_key('nroflines'):
+        header['nroflines'] = '3'    
+    if not header.has_key('nrofcellsperline'):
+        header['nrofcellsperline'] = '4'
+    if not header.has_key('registrationcellx'):
+        header['registrationcellx'] = '0'
+    if not header.has_key('registrationcelly'):
+        header['registrationcelly'] = str(int(header['nroflines'])-1)
+    if not header.has_key('nrofbands'):
+        header['nrofbands'] = '1'
+    if not header.has_key('value'):
+        header['value'] = '"Default_Band"'
+
+
+    return header    
+                          
+    
+
+
+