source: XIOS/dev/dev_trunk_omp/src/io/onetcdf4_impl.hpp @ 1666

Last change on this file since 1666 was 1666, checked in by yushan, 2 years ago

MARK: branch merged with trunk @1663. boost.tar.gz updated to 1652 for visualization.

  • 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: 4.3 KB
Line 
1#ifndef __ONETCDF4_IMPL_HPP__
2#define __ONETCDF4_IMPL_HPP__
3
4#include "onetcdf4.hpp"
5#include "netCdfInterface.hpp"
6#include "timer.hpp"
7
8namespace xios
9{
10  template <class T, int ndim>
11  void CONetCDF4::writeData(const CArray<T, ndim>& data, const StdString & name,
12                            bool collective, StdSize record,
13                            const std::vector<StdSize> * start,
14                            const std::vector<StdSize> * count)
15  {
16    int grpid = this->getCurrentGroup();
17    int varid = this->getVariable(name);
18    StdSize array_size = 1;
19    std::vector<StdSize> sstart, scount;
20
21    if (this->wmpi && collective)
22    CNetCdfInterface::varParAccess(grpid, varid, NC_COLLECTIVE);
23    if (this->wmpi && !collective)
24    CNetCdfInterface::varParAccess(grpid, varid, NC_INDEPENDENT);
25
26    CTimer::get("Files : get data infos").resume();
27    this->getWriteDataInfos
28    (name, record, array_size,  sstart, scount, start, count);
29    CTimer::get("Files : get data infos").suspend();
30
31    if (data.numElements() != array_size)
32    {
33      ERROR("CONetCDF4::writeData(...)",
34      << "[ input array size = "  << data.numElements()
35      << ", intern array size = " << array_size
36      << " ] Invalid input data !" );
37    }
38
39    this->writeData_(grpid, varid, sstart, scount, data.dataFirst());
40  }
41
42  template <>
43  void CONetCDF4::writeData(const CArray<StdString, 1>& data, const StdString & name,
44                            bool collective, StdSize record,
45                            const std::vector<StdSize> * start,
46                            const std::vector<StdSize> * count)
47  {
48    int grpid = this->getCurrentGroup();
49    int varid = this->getVariable(name);
50    StdSize array_size = 1;
51    std::vector<StdSize> sstart, scount;
52
53    if (this->wmpi && collective)
54    CNetCdfInterface::varParAccess(grpid, varid, NC_COLLECTIVE);
55    if (this->wmpi && !collective)
56    CNetCdfInterface::varParAccess(grpid, varid, NC_INDEPENDENT);
57
58    CTimer::get("CONetCDF4::writeData getWriteDataInfos").resume();
59    this->getWriteDataInfos(name, record, array_size,  sstart, scount, start, count);
60    CTimer::get("CONetCDF4::writeData getWriteDataInfos").suspend();
61 
62    int dimArrayLen;
63    for (int i=0; i<scount.size(); i++)
64      if (scount[i] == stringArrayLen)
65        dimArrayLen = i;
66
67    if (data.numElements()*stringArrayLen != array_size)
68    {
69      ERROR("CONetCDF4::writeData(...)",
70      << "[ input array size = "  << data.numElements()*stringArrayLen
71      << ", intern array size = " << array_size
72      << " ] Invalid input data !" );
73    }
74    char *PtrArrayStr ;
75    PtrArrayStr=new char[stringArrayLen*data.numElements()] ;
76    memset (PtrArrayStr,' ',stringArrayLen*data.numElements());
77    size_t offset=0 ;
78/*    Array<StdString,1>::const_iterator it, itb=data.begin(), ite=data.end() ;
79    for(it=itb;it!=ite;++it, offset+=stringArrayLen)
80    {
81      it->copy(PtrArrayStr+offset,it->size()) ;
82      PtrArrayStr[offset+it->size()]='\0' ;
83    }
84*/
85    for(int i=0;i<data.numElements();i++,offset+=stringArrayLen)
86    {
87      data(i).copy(PtrArrayStr+offset,data(i).size()) ;
88      PtrArrayStr[offset+data(i).size()]='\0' ;
89    } 
90
91    CTimer::get("CONetCDF4::writeData writeData_").resume();
92    this->writeData_(grpid, varid, sstart, scount, PtrArrayStr);
93    CTimer::get("CONetCDF4::writeData writeData_").suspend();
94
95    delete []  PtrArrayStr;
96  }
97
98//----------------------------------------------------------------
99
100  template <class T>
101  void CONetCDF4::setDefaultValue(const StdString & varname, const T * value)
102  {
103    int grpid = this->getCurrentGroup();
104    int varid = this->getVariable(varname);
105
106    if (value != NULL)
107    {
108      // nc_def_var_fill will automatically set the _FillValue attribute when
109      // using the NetCDF 4 format but we need to do it manually otherwise
110      if (useClassicFormat)
111        this->addAttribute(StdString("_FillValue"), *value, &varname);
112      else
113        CNetCdfInterface::defVarFill(grpid, varid, 0, (void*)value);
114      this->addAttribute(StdString("missing_value"), *value, &varname);
115    }
116    else if (!useClassicFormat)
117    {
118      // The "no-fill mode" is set globally for the classic NetCDF format
119      CNetCdfInterface::defVarFill(grpid, varid, 1, NULL);
120    }
121  }
122
123  ///---------------------------------------------------------------
124
125}
126
127
128
129#endif
Note: See TracBrowser for help on using the repository browser.