FONcMap.cc

// FONcMap.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 "config.h"
 
#include <cstring>
 
#include <libdap/Array.h>
 
#include <BESDebug.h>
 
#include "FONcArray.h"
#include "FONcMap.h"
#include "FONcUtils.h"

FONcMap::FONcMap(FONcArray *a, bool ingrid)
        : _arr(a), _ingrid(ingrid), _defined(false), _ref(1)
{
}

FONcMap::~FONcMap()
{
    if (_ingrid) {
        delete _arr;
        _arr = 0;
    }
}

void
FONcMap::decref()
{
    _ref--;
    if (!_ref) delete this;
}

bool FONcMap::compare(libdap::Array *tomap)
{
    bool isequal = true;
 
    libdap::Array *map = _arr->array();
 
    BESDEBUG("fonc", "FONcMap::compare - comparing " << tomap->name()
                                                     << " to " << map->name() << endl);
 
    // compare the name
    if (isequal && tomap->name() != map->name()) {
        isequal = false;
    }
 
    // compare the type
    if (isequal && tomap->var()->type() != map->var()->type()) {
        isequal = false;
    }
 
    // compare the length of the array
    if (isequal && tomap->length() != map->length()) {
        isequal = false;
    }
 
    // compare the number of dimensions
    if (isequal && tomap->dimensions() != map->dimensions()) {
        isequal = false;
    }
 
    // the variable name needs to be the same as the dimension name
    if (isequal && map->dimension_name(map->dim_begin()) != map->name()) {
        isequal = false;
    }
 
    // compare the dimension name
    if (isequal && tomap->dimension_name(tomap->dim_begin()) != map->dimension_name(map->dim_begin())) {
        isequal = false;
    }
 
    // compare the dimension size. Is this the same as the length of the array
    if (isequal && tomap->dimension_size_ll(tomap->dim_begin(), true) != map->dimension_size_ll(map->dim_begin(), true)) {
        isequal = false;
    }
 
    if (isequal) {
        // compare the values of the array
        char *map_buf = map->get_buf();
        char *tomap_buf = tomap->get_buf();
 
        // Only compare the values if both the value vectors are not null.
        // Optimizations in this handler might cause one of these to be
        // null, resulting in a segmenation fault. If one or both are null,
        // assume that by this point in the series of tests, the two maps
        // are equal. jhrg 6/8/21
        if (map_buf && tomap_buf && 0 != memcmp(map_buf, tomap_buf, map->width()))
            isequal = false;
    }
 
    BESDEBUG("fonc",
             "FONcMap::compare - done comparing " << tomap->name() << " to " << map->name() << ": " << isequal << endl);
    return isequal;
}

void
FONcMap::add_grid(const string &name)
{
    _shared_by.push_back(name);
}

void
FONcMap::clear_embedded()
{
    _arr->clear_embedded();
}

void
FONcMap::define(int ncid)
{
    if (!_defined) {
        _arr->define(ncid);
        _defined = true;
    }
}

void
FONcMap::write(int ncid)
{
    _arr->write(ncid);
}

void
FONcMap::dump(ostream &strm) const
{
    strm << BESIndent::LMarg << "FONcMap::dump - ("
         << (void *) this << ")" << endl;
    BESIndent::Indent();
    strm << BESIndent::LMarg << "array:";
    if (_arr) {
        strm << endl;
        BESIndent::Indent();
        _arr->dump(strm);
        BESIndent::UnIndent();
    }
    else {
        strm << " not set" << endl;
    }
    strm << BESIndent::LMarg << "shared by: ";
    vector<string>::const_iterator i = _shared_by.begin();
    vector<string>::const_iterator e = _shared_by.end();
    bool first = true;
    for (; i != e; i++) {
        if (!first) strm << ", ";
        strm << (*i);
        first = false;
    }
    strm << endl;
    BESIndent::UnIndent();
}