creating vector using netcdf into array
Question: I'm fairly new to python and have found stack overflow one of the best
resources out there, now I'm hoping someone can help me with what I believe is
a fairly basic question.

I'm looking to create a land mask from a list of lats and lons and rainfall
data extracted from a netCDF file. I need to get the data from the netcdf file
to line up so I can remove rows which have a rainfall value of '-9999.'
(indicating no data because its over the ocean). I can access the file, I can
create a mesh grid, but when it comes to inserting the rainfall data for the
final check I'm getting odd shapes and no luck with the logical test. Can
someone have a look at this code and let me know what you think?

    
    
    from netCDF4 import Dataset
    import numpy as np
    
    f=Dataset('/Testing/Ensemble_grid/1970_2012_eMAST_ANUClimate_mon_evap_v1m0_197001.nc')
    
    lat = f.variables['latitude'][:]
    lon = f.variables['longitude'][:]
    rainfall = np.array(f.variables['lwe_thickness_of_precipitation_amount'])
    lons, lats = np.meshgrid(lon,lat)
    full_ary = np.array((lats,lons))
    full_lats_lons = np.swapaxes(full_ary,0,2)
    rain_data = np.squeeze(rainfall,axis=(0,))
    grid = np.array((full_lats_lons,rain_data))
    full_grid = np.expand_dims(grid,axis=1)
    full_grid_col = np.swapaxes(full_grid,0,1)
    land_grid = np.logical_not(full_grid_col[:,1]==-9999.)
    


Answer: Ian, you need to put a repeatable example up here...

I suspect what you need is something like this;

    
    
    x = np.array([[1, 2, 3], [4, 5, 6]], np.int32)
    x.flat
    

