Robert, Thanks for your reply. However, I am still having problems. Sometimes I get a scalar return and sometimes I get an array. For example, using the netCDF file:
netcdf simple { dimensions: num = 3 ; variables: float temp0(num) ; int temp1(num) ; data: temp0 = 1., 2., 3. ; temp1 = 1, 2, 3 ; } and running: # import Numeric print Numeric.__version__ from Scientific.IO.NetCDF import NetCDFFile cdf_file1 = NetCDFFile("simple.nc","r") var1 = cdf_file1.variables["temp0"][:] var2 = cdf_file1.variables["temp1"][:] min1 = reduce(Numeric.minimum,var1) min2 = reduce(Numeric.minimum,var2) print "Types of var1, min(var1), min1:",type(var1), type(min(var1)), type(min1) print "Types of var2, min(var2), min2:",type(var2), type(min(var2)), type(min2) I get: 24.0b2 Types of var1, min(var1), min1: <type 'array'> <type 'array'> <type 'array'> Types of var2, min(var2), min2: <type 'array'> <type 'int'> <type 'int'> Even something like: >>> import Numeric >>> a = Numeric.array([1.,2.]) >>> print type(a),type(min(a)) <type 'array'> <type 'float'> does not produce an array. Any comments woud be appreciated. Fred Clare -- http://mail.python.org/mailman/listinfo/python-list