Changeset 2273


Ignore:
Timestamp:
Jan 24, 2006, 5:40:50 PM (18 years ago)
Author:
ole
Message:

Updated generation of test files

File:
1 edited

Legend:

Unmodified
Added
Removed
  • inundation/pyvolution/test_data_manager.py

    r2263 r2273  
    7272
    7373
    74         #Write A testfile
    75 
    76         self.test_MOST_file = 'xxxsmall_ha.nc'
     74        #Write A testfile (not realistic. Values aren't realistic)
     75
     76        self.test_MOST_file = 'most_small'
    7777
    7878        longitudes = [150.66667, 150.83334, 151., 151.16667]
     
    8787
    8888
    89         fid = NetCDFFile(self.test_MOST_file, 'w')
    90        
    91         fid.createDimension(long_name,nx)
    92         fid.createVariable(long_name,'d',(long_name,))
    93         fid.variables[long_name].point_spacing='uneven'
    94         fid.variables[long_name].units='degrees_east'
    95         fid.variables[long_name].assignValue(longitudes)
    96 
    97         fid.createDimension(lat_name,ny)
    98         fid.createVariable(lat_name,'d',(lat_name,))
    99         fid.variables[lat_name].point_spacing='uneven'
    100         fid.variables[lat_name].units='degrees_north'
    101         fid.variables[lat_name].assignValue(latitudes)
    102 
    103         fid.createDimension('TIME',six)
    104         fid.createVariable('TIME','d',('TIME',))
    105         fid.variables['TIME'].point_spacing='uneven'
    106         fid.variables['TIME'].units='seconds'
    107         fid.variables['TIME'].assignValue([0.0, 0.1, 0.6, 1.1, 1.6, 2.1])
    108 
    109         fid.createVariable('HA','d',('TIME', lat_name, long_name))
    110         #fid.createVariable('HA','d',(lat_name, long_name))
    111         fid.variables['HA'].assignValue([[[0.3400644, 0, -46.63519, -6.50198],
    112                                           [-0.1214216, 0, 0, 0],
    113                                           [0, 0, 0, 0],
    114                                           [0, 0, 0, 0]],
    115                                          [[0.3400644, 2.291054e-005, -23.33335, -6.50198],
    116                                           [-0.1213987, 4.581959e-005, -1.594838e-007, 1.421085e-012],
    117                                           [2.291054e-005, 4.582107e-005, 4.581715e-005, 1.854517e-009],
    118                                           [0, 2.291054e-005, 2.291054e-005, 0]],
    119                                          [[0.3400644, 0.0001374632, -23.31503, -6.50198],
    120                                           [-0.1212842, 0.0002756907, 0.006325484, 1.380492e-006],
    121                                           [0.0001374632, 0.0002749264, 0.0002742863, 6.665601e-008],
    122                                           [0, 0.0001374632, 0.0001374632, 0]],
    123                                          [[0.3400644, 0.0002520159, -23.29672, -6.50198],
    124                                           [-0.1211696, 0.0005075303, 0.01264618, 6.208276e-006],
    125                                           [0.0002520159, 0.0005040318, 0.0005027961, 2.23865e-007],
    126                                           [0, 0.0002520159, 0.0002520159, 0]],
    127                                          [[0.3400644, 0.0003665686, -23.27842, -6.50198],
    128                                           [-0.1210551, 0.0007413362, 0.01896192, 1.447638e-005],
    129                                           [0.0003665686, 0.0007331371, 0.0007313463, 4.734126e-007],
    130                                           [0, 0.0003665686, 0.0003665686, 0]],
    131                                          [[0.3400644, 0.0004811212, -23.26012, -6.50198],
    132                                           [-0.1209405, 0.0009771062, 0.02527271, 2.617787e-005],
    133                                           [0.0004811212, 0.0009622425, 0.0009599366, 8.152277e-007],
    134                                           [0, 0.0004811212, 0.0004811212, 0]]])
     89        for ext in ['_ha.nc', '_ua.nc', '_va.nc', '_e.nc']:
     90            fid = NetCDFFile(self.test_MOST_file + ext, 'w')
     91       
     92            fid.createDimension(long_name,nx)
     93            fid.createVariable(long_name,'d',(long_name,))
     94            fid.variables[long_name].point_spacing='uneven'
     95            fid.variables[long_name].units='degrees_east'
     96            fid.variables[long_name].assignValue(longitudes)
     97
     98            fid.createDimension(lat_name,ny)
     99            fid.createVariable(lat_name,'d',(lat_name,))
     100            fid.variables[lat_name].point_spacing='uneven'
     101            fid.variables[lat_name].units='degrees_north'
     102            fid.variables[lat_name].assignValue(latitudes)
     103
     104            fid.createDimension('TIME',six)
     105            fid.createVariable('TIME','d',('TIME',))
     106            fid.variables['TIME'].point_spacing='uneven'
     107            fid.variables['TIME'].units='seconds'
     108            fid.variables['TIME'].assignValue([0.0, 0.1, 0.6, 1.1, 1.6, 2.1])
     109           
     110
     111            name = ext[1:3].upper()
     112            if name == 'E.': name = 'ELEVATION'
     113            fid.createVariable(name,'d',('TIME', lat_name, long_name))
     114            fid.variables[name].units='CENTIMETERS'
     115            fid.variables[name].missing_value=-1.e+034       
     116           
     117            fid.variables[name].assignValue([[[0.3400644, 0, -46.63519, -6.50198],
     118                                              [-0.1214216, 0, 0, 0],
     119                                              [0, 0, 0, 0],
     120                                              [0, 0, 0, 0]],
     121                                             [[0.3400644, 2.291054e-005, -23.33335, -6.50198],
     122                                              [-0.1213987, 4.581959e-005, -1.594838e-007, 1.421085e-012],
     123                                              [2.291054e-005, 4.582107e-005, 4.581715e-005, 1.854517e-009],
     124                                              [0, 2.291054e-005, 2.291054e-005, 0]],
     125                                             [[0.3400644, 0.0001374632, -23.31503, -6.50198],
     126                                              [-0.1212842, 0.0002756907, 0.006325484, 1.380492e-006],
     127                                              [0.0001374632, 0.0002749264, 0.0002742863, 6.665601e-008],
     128                                              [0, 0.0001374632, 0.0001374632, 0]],
     129                                             [[0.3400644, 0.0002520159, -23.29672, -6.50198],
     130                                              [-0.1211696, 0.0005075303, 0.01264618, 6.208276e-006],
     131                                              [0.0002520159, 0.0005040318, 0.0005027961, 2.23865e-007],
     132                                              [0, 0.0002520159, 0.0002520159, 0]],
     133                                             [[0.3400644, 0.0003665686, -23.27842, -6.50198],
     134                                              [-0.1210551, 0.0007413362, 0.01896192, 1.447638e-005],
     135                                              [0.0003665686, 0.0007331371, 0.0007313463, 4.734126e-007],
     136                                              [0, 0.0003665686, 0.0003665686, 0]],
     137                                             [[0.3400644, 0.0004811212, -23.26012, -6.50198],
     138                                              [-0.1209405, 0.0009771062, 0.02527271, 2.617787e-005],
     139                                              [0.0004811212, 0.0009622425, 0.0009599366, 8.152277e-007],
     140                                              [0, 0.0004811212, 0.0004811212, 0]]])
    135141           
    136142       
     
    18621868
    18631869
    1864     def test_ferret2sww(self):
     1870    def test_ferret2sww1(self):
    18651871        """Test that georeferencing etc works when converting from
    18661872        ferret format (lat/lon) to sww format (UTM)
     
    18811887        #Read
    18821888        from coordinate_transforms.redfearn import redfearn
    1883         fid = NetCDFFile(self.test_MOST_file)
     1889        #fid = NetCDFFile(self.test_MOST_file)
     1890        fid = NetCDFFile(self.test_MOST_file + '_ha.nc')       
    18841891        first_value = fid.variables['HA'][:][0,0,0]
    18851892        fourth_value = fid.variables['HA'][:][0,0,3]
     
    18871894
    18881895        #Call conversion (with zero origin)
    1889         ferret2sww('small', verbose=False,
     1896        #ferret2sww('small', verbose=False,
     1897        #           origin = (56, 0, 0))
     1898        ferret2sww(self.test_MOST_file, verbose=False,
    18901899                   origin = (56, 0, 0))
    1891 
    18921900
    18931901        #Work out the UTM coordinates for first point
     
    18961904
    18971905        #Read output file 'small.sww'
    1898         fid = NetCDFFile('small.sww')
     1906        #fid = NetCDFFile('small.sww')
     1907        fid = NetCDFFile(self.test_MOST_file + '.sww')       
    18991908
    19001909        x = fid.variables['x'][:]
     
    19211930        #Cleanup
    19221931        import os
    1923         os.remove('small.sww')
    1924 
     1932        #os.remove('small.sww')
     1933        os.remove(self.test_MOST_file + '.sww')
    19251934
    19261935
     
    19431952
    19441953        #fid = NetCDFFile('small_ha.nc')
    1945         fid = NetCDFFile(self.test_MOST_file)
     1954        fid = NetCDFFile(self.test_MOST_file + '_ha.nc')
    19461955
    19471956        #Pick a coordinate and a value
     
    19601969
    19611970        #Call conversion (with zero origin)
    1962         ferret2sww('small', verbose=False,
     1971        ferret2sww(self.test_MOST_file, verbose=False,
    19631972                   origin = (56, 0, 0))
    19641973
     
    19681977
    19691978        #Read output file 'small.sww'
    1970         fid = NetCDFFile('small.sww')
     1979        fid = NetCDFFile(self.test_MOST_file + '.sww')
    19711980
    19721981        x = fid.variables['x'][:]
     
    19861995        #Cleanup
    19871996        import os
    1988         os.remove('small.sww')
     1997        os.remove(self.test_MOST_file + '.sww')
    19891998
    19901999
     
    21542163
    21552164        #Call conversion (with nonzero origin)
    2156         ferret2sww('small', verbose=False,
     2165        ferret2sww(self.test_MOST_file, verbose=False,
    21572166                   origin = (56, 100000, 200000))
    21582167
     
    21622171
    21632172        #Read output file 'small.sww'
    2164         fid = NetCDFFile('small.sww', 'r')
    2165 
     2173        #fid = NetCDFFile('small.sww', 'r')
     2174        fid = NetCDFFile(self.test_MOST_file + '.sww')
     2175       
    21662176        x = fid.variables['x'][:]
    21672177        y = fid.variables['y'][:]
     
    21742184
    21752185        #Cleanup
    2176         os.remove('small.sww')
     2186        os.remove(self.test_MOST_file + '.sww')
    21772187
    21782188
Note: See TracChangeset for help on using the changeset viewer.