Bug with data type conversion in C++ interface?

[Thomas Orgis <orgis@xxxxxxxxxxxxxxxxxxxx>]

Hi!

I recently started to deal with NetCDF files and first tried to use the C++ 
interface. I reverted to the plain C interface because of what looks like a bug 
to me:

I have a variable of type double and tried to read values into a float array... 
I'll let the code talk (time is of type double and only depends on one 
dimension, namely time):

        NcVar* timo = ncfile->get_var("time");
        cout << "[datagroup] have " << *(timo->edges()) << "times" << endl;
        float* thetimes = new float[*(timo->edges())];
        timo->get(thetimes,*(timo->edges()));
        for(long i = 0; i <= *(timo->edges()); ++i)
        {
                cout << "time " << i << ": " << thetimes[i] << endl;
        }

This results in only zeroes being printed. Without the need of conversion

        double* thetimes = new double[*(timo->edges())];
        timo->get(thetimes,*(timo->edges()));

I see the correct values. It should also be mentioned that the actual values 
fit easily in a float number; in fact they are in the range of 0 to 60 in steps 
of exactly 0.5.

When using the classical nc_get_var_float() I have the correct converted values.

So, I suspect that the automagic conversion via the overloaded methods of the 
NcVar class doesn't really work here. Can one confirm this or is this just an 
oddity of my setup/data/compiler?

-- 
Thomas Orgis, theoretischer Physiker - praktischer Mensch.

Arbeitsgruppe Nichtlineare Dynamik
Institut für Physik der Universität Potsdam