source: XIOS/trunk/src/io/onetcdf4_impl.hpp @ 1050

Last change on this file since 1050 was 1050, checked in by ymipsl, 7 years ago
  • Add label attribute for axis
  • if label attribute is present, then only only string label for each level will be output in netcdf file, ie value and bounds will not be output.

YM

  • Property copyright set to
    Software name : XIOS (Xml I/O Server)
    http://forge.ipsl.jussieu.fr/ioserver
    Creation date : January 2009
    Licence : CeCCIL version2
    see license file in root directory : Licence_CeCILL_V2-en.txt
    or http://www.cecill.info/licences/Licence_CeCILL_V2-en.html
    Holder : CEA/LSCE (Laboratoire des Sciences du CLimat et de l'Environnement)
    CNRS/IPSL (Institut Pierre Simon Laplace)
    Project Manager : Yann Meurdesoif
    yann.meurdesoif@cea.fr
File size: 3.5 KB
Line 
1#ifndef __ONETCDF4_IMPL_HPP__
2#define __ONETCDF4_IMPL_HPP__
3
4#include "onetcdf4.hpp"
5#include "netCdfInterface.hpp"
6
7namespace xios
8{
9  template <class T, int ndim>
10  void CONetCDF4::writeData(const CArray<T, ndim>& data, const StdString & name,
11                            bool collective, StdSize record,
12                            const std::vector<StdSize> * start,
13                            const std::vector<StdSize> * count)
14  {
15    int grpid = this->getCurrentGroup();
16    int varid = this->getVariable(name);
17    StdSize array_size = 1;
18    std::vector<StdSize> sstart, scount;
19
20    if (this->wmpi && collective)
21    CNetCdfInterface::varParAccess(grpid, varid, NC_COLLECTIVE);
22    if (this->wmpi && !collective)
23    CNetCdfInterface::varParAccess(grpid, varid, NC_INDEPENDENT);
24
25    this->getWriteDataInfos
26    (name, record, array_size,  sstart, scount, start, count);
27    if (data.numElements() != array_size)
28    {
29      ERROR("CONetCDF4::writeData(...)",
30      << "[ input array size = "  << data.numElements()
31      << ", intern array size = " << array_size
32      << " ] Invalid input data !" );
33    }
34
35    this->writeData_(grpid, varid, sstart, scount, data.dataFirst());
36  }
37
38  template <>
39  void CONetCDF4::writeData(const CArray<StdString, 1>& data, const StdString & name,
40                            bool collective, StdSize record,
41                            const std::vector<StdSize> * start,
42                            const std::vector<StdSize> * count)
43  {
44    int grpid = this->getCurrentGroup();
45    int varid = this->getVariable(name);
46    StdSize array_size = 1;
47    std::vector<StdSize> sstart, scount;
48
49    if (this->wmpi && collective)
50    CNetCdfInterface::varParAccess(grpid, varid, NC_COLLECTIVE);
51    if (this->wmpi && !collective)
52    CNetCdfInterface::varParAccess(grpid, varid, NC_INDEPENDENT);
53
54    this->getWriteDataInfos
55    (name, record, array_size,  sstart, scount, start, count);
56    if (data.numElements()*stringArrayLen != array_size)
57    {
58      ERROR("CONetCDF4::writeData(...)",
59      << "[ input array size = "  << data.numElements()
60      << ", intern array size = " << array_size
61      << " ] Invalid input data !" );
62    }
63    char* ArrayStr ;
64    char *PtrArrayStr ;
65    PtrArrayStr=ArrayStr=new char[data.numElements()*stringArrayLen] ;
66    Array<StdString,1>::const_iterator it, itb=data.begin(), ite=data.end() ;
67    for(it=itb;it!=ite;++it,PtrArrayStr+=stringArrayLen)
68    {
69      it->copy(PtrArrayStr,it->size()) ;
70      PtrArrayStr[it->size()]='\0' ;
71    }
72    this->writeData_(grpid, varid, sstart, scount, ArrayStr);
73    delete [] ArrayStr ;
74  }
75
76//----------------------------------------------------------------
77
78  template <class T>
79  void CONetCDF4::setDefaultValue(const StdString & varname, const T * value)
80  {
81    int grpid = this->getCurrentGroup();
82    int varid = this->getVariable(varname);
83
84    if (value != NULL)
85    {
86      // nc_def_var_fill will automatically set the _FillValue attribute when
87      // using the NetCDF 4 format but we need to do it manually otherwise
88      if (useClassicFormat)
89        this->addAttribute(StdString("_FillValue"), *value, &varname);
90      else
91        CNetCdfInterface::defVarFill(grpid, varid, 0, (void*)value);
92      this->addAttribute(StdString("missing_value"), *value, &varname);
93    }
94    else if (!useClassicFormat)
95    {
96      // The "no-fill mode" is set globally for the classic NetCDF format
97      CNetCdfInterface::defVarFill(grpid, varid, 1, NULL);
98    }
99  }
100
101  ///---------------------------------------------------------------
102
103}
104
105
106
107#endif
Note: See TracBrowser for help on using the repository browser.