bandw wrote:
> 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.
Hmm, odd. Anyways, follow my advice: use minimum.reduce() and wrap
results in float() or array() if you really need floats or rank-0 arrays.
--
Robert Kern
[EMAIL PROTECTED]
"In the fields of hell where the grass grows high
Are the graves of dreams allowed to die."
-- Richard Harter
--
http://mail.python.org/mailman/listinfo/python-list
