FONcDouble.cc

// FONcDouble.cc
 
// This file is part of BES Netcdf File Out Module
 
// Copyright (c) 2004,2005 University Corporation for Atmospheric Research
// Author: Patrick West <pwest@ucar.edu> and Jose Garcia <jgarcia@ucar.edu>
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
// License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
//
// This library is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
// Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
//
// You can contact University Corporation for Atmospheric Research at
// 3080 Center Green Drive, Boulder, CO 80301
 
// (c) COPYRIGHT University Corporation for Atmospheric Research 2004-2005
// Please read the full copyright statement in the file COPYRIGHT_UCAR.
//
// Authors:
//      pwest       Patrick West <pwest@ucar.edu>
//      jgarcia     Jose Garcia <jgarcia@ucar.edu>
 
#include <BESInternalError.h>
#include <BESDebug.h>
#include <libdap/Float64.h>
 
#include "FONcDouble.h"
#include "FONcUtils.h"
#include "FONcAttributes.h"

FONcDouble::FONcDouble(BaseType *b)
        : FONcBaseType(), _f(0) {
    _f = dynamic_cast<Float64 *>(b);
    if (!_f) {
        string s = (string) "File out netcdf, FONcDouble was passed a "
                   + "variable that is not a DAP Float64";
        throw BESInternalError(s, __FILE__, __LINE__);
    }
}

void
FONcDouble::define(int ncid) {
    FONcBaseType::define(ncid);
 
    if (!d_defined) {
 
        if (d_is_dap4) {
            D4Attributes *d4_attrs = _f->attributes();
            updateD4AttrType(d4_attrs, NC_DOUBLE);
        }
        else {
            AttrTable &attrs = _f->get_attr_table();
            updateAttrType(attrs, NC_DOUBLE);
        }
 
        FONcAttributes::add_variable_attributes(ncid, d_varid, _f, isNetCDF4_ENHANCED(), d_is_dap4);
        FONcAttributes::add_original_name(ncid, d_varid, d_varname, d_orig_varname);
 
        d_defined = true;
    }
}

void
FONcDouble::write(int ncid) {
    BESDEBUG("fonc", "FONcDouble::write for var " << d_varname << endl);
 
    if (d_is_dap4)
        _f->intern_data();
    else
        _f->intern_data(*get_eval(), *get_dds());
 
    double data = _f->value();
    size_t var_index[] = {0};
    int stax = nc_put_var1_double(ncid, d_varid, var_index, &data);
    if (stax != NC_NOERR) {
        string err = (string) "fileout.netcdf - "
                     + "Failed to write double data for "
                     + d_varname;
        FONcUtils::handle_error(stax, err, __FILE__, __LINE__);
    }
 
    BESDEBUG("fonc", "FONcDouble::done write for var " << d_varname << endl);
}

void
FONcDouble::dump(ostream &strm) const {
    strm << BESIndent::LMarg << "FONcDouble::dump - ("
         << (void *) this << ")" << endl;
    BESIndent::Indent();
    strm << BESIndent::LMarg << "name = " << _f->name() << endl;
    BESIndent::UnIndent();
}