Cnetcdf.h

Go to the documentation of this file.

/**********************************************************************
 * Copyright (C) 2012, The YaDICs Project Developers. 
 * See the COPYRIGHT file at the top-level directory of this distribution ./COPYRIGHT.
 * See ./COPYING file for copying and redistribution conditions.
 * 
 * This file is part of YaDICs.
 * YaDICs is free software: you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation, either version 3 of the License, or
 * (at your option) any later version.
 * 
 * YaDICs 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 General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with YaDICs.  If not, see <http://www.gnu.org/licenses/>.
 * 
 * Information about how to use the software are provided at http://yadic.univ-lille1.fr/
 **********************************************************************/


#ifndef CNETCDF
#define CNETCDF

//--------------------------------------------------------------------
#include "netcdf.h"
#include <mpi.h>
#include <assert.h>
#include "hdf5.h"
#include <string.h>
#include <stdlib.h>

#include <iostream>
#include <cmath>
#include <vector>
#include <fstream>
#include <ctime>
#include <cstdlib>
#include <stdio.h>
#include <sstream>
#include <unistd.h>
    #define GetCurrentDir getcwd

// #include "../CImg/CImg.h"
// using namespace cimg_library;
//--------------------------------------------------------------------


template<typename Timg> 
class Cnetcdf
{
  
  public:
    
  Cimage()
  {// constructor
  
    class_name = "Netcdf_loader";
    
  }  


  bool read_var()
  {
    
    bool error;

    nc_open_par(const char *path, NC_NOCLOBBER || NC_NETCDF4, MPI_Comm comm, MPI_Info info,int *ncidp);
    
    return error;
    
  }
  
  subroutine read_var3d(file_name,var_name,var,mpid,inter)
! -----------------------------------------------------------------------
! io : read 3d variable in a netcdf file
! -----------------------------------------------------------------------
! Matthieu Marquillie
! 04/2013
!
     use mpi
     use class_md
     implicit none
     character(len=*),intent(in) :: file_name,var_name
     real(rk) :: var(:,:,:)
     integer(ik) :: varid(1),i
     integer(ik) :: ncid

     integer(ik),parameter :: ndim=3
     type(mpi_data),optional :: mpid
     integer(ik) :: dim_len(ndim),dimid(ndim),dim_len_check
     integer(ik) :: dimt(3),coord(3,2)
     integer(ik) :: startv(3),countv(3)
     character(*),optional :: inter

     !-> open file
     if (present(mpid)) then
        call io_check(nf90_open(path=file_name,&
!               mode=IOR(NF90_WRITE,NF90_MPIPOSIX),ncid=ncid,&
             mode=IOR(NF90_WRITE,NF90_MPIIO),ncid=ncid,&
             comm=mpid%comm,info=MPI_INFO_NULL))
     else
        call io_check(nf90_open(path=file_name,mode=nf90_write,ncid=ncid))
     endif

     !-> get variable id
     call io_check(nf90_inq_varid(ncid,var_name,varid(1)))

     !-> recompute dimensions, start and count if mpi
     if (present(mpid)) then
        if (present(inter)) then
           call md_mpi_global_coord(mpid,dimt,coord,inter=inter)
        else
           call md_mpi_global_coord(mpid,dimt,coord)
        endif
        startv=(/1,1,1/)
        countv=(/1,1,1/)
        do i=1,3
           !if (dim_len(i)>1) then
              dim_len(i)=dimt(i)
              startv(i)=coord(i,1)
              countv(i)=coord(i,2)
           !endif
        enddo
     else
        startv=(/1,1,1/)
        countv=get_dim_size(var)
     endif

     !-> read field variable
!    call io_check(nf90_get_var(ncid,varid(1),var))
     call
io_check(nf90_get_var(ncid,varid(1),var,start=startv,count=countv))

     !-> close file
     call io_check(nf90_close(ncid))

   end subroutine read_var3d