#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)

#define INDICES1 ij,l,iq,iedge,edge,ivertex,vertex,ij_left,ij_right
#define INDICES2 ij_up,ij_down,itrisk,edge_trisk,kup,kdown
#define EDGES edge1,edge2,edge3,edge4,edge5,edge6
#define VERTICES vertex1,vertex2,vertex3,vertex4,vertex5,vertex6
#define SIGNS sign1,sign2,sign3,sign4,sign5,sign6
#define EDGE_ENDS ij_up1,ij_up2,ij_up3,ij_up4,ij_up5,ij_up6,ij_down1,ij_down2,ij_down3,ij_down4,ij_down5,ij_down6
#define LENGTHS le_de1,le_de2,le_de3,le_de4,le_de5,le_de6
#define DECLARE_INDICES INTEGER INDICES1,INDICES2
#define DECLARE_EDGES NUM SIGNS,LENGTHS ; INTEGER EDGES, EDGE_ENDS
#define DECLARE_VERTICES INTEGER VERTICES

#define START_TRACE(id,nprimal,ndual,nedge) CALL enter_trace(id, 8*llm*((nprimal)*primal_num+(ndual)*dual_num+(nedge)*edge_num) )
#define STOP_TRACE CALL exit_trace()