#define DBL REAL(C_DOUBLE)
#define FLT REAL(C_FLOAT)
#define TIME DBL

#ifdef CPP_MIXED_PREC
#define NUM FLT
#define MPI_NUM MPI_REAL4
#else
#define NUM DBL
#define MPI_NUM MPI_REAL8
#endif

! NUM and MPI_NUM are used for physical quantities (scalar or array)
! DBL and MPI_REAL8 are used for variables that always require double precision

#define NUM1(m) NUM, DIMENSION(m)
#define NUM2(m,n) NUM, DIMENSION(m,n)
#define NUM3(m,n,p) NUM, DIMENSION(m,n,p)
#define NUM4(m,n,p,q) NUM, DIMENSION(m,n,p,q)
#define INDEX INTEGER(C_INT)

#define FIELD_PS      NUM1(primal_num)
#define FIELD_MASS    NUM2(llm, primal_num)
#define FIELD_Z       NUM2(llm, dual_num)
#define FIELD_U       NUM2(llm, edge_num)
#define FIELD_UL      NUM2(llm+1, edge_num)
#define FIELD_THETA   NUM3(llm, primal_num, nqdyn)
#define FIELD_GEOPOT  NUM2(llm+1, primal_num)

#define BINDC_(thename) BIND(C, name=#thename)
#define BINDC(thename) BINDC_(dynamico_ ## thename)