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
