moving opening a netcdf file from reading_netcdf to the main function - it caused segfault on Sylwester's computer (the problem is already solved at github.com/scipy/scipy)