/** @file cdr_advect_diffu_wrap.c

*   @brief Wrapper function for the cdr_advect_diffu function. 

*

*   This function makes it possible to call the Fortran routine cdr_advect_vec

*   a c function. The c routine is named arcos_cdr_advect_vec.

*

*   @author Margreet Nool - 2009

*/

#include "stdlib.h"

#include "stdio.h"

#include "math.h"

int number_cdr=0;

/* cdr_advect: vectorized fortran90 routine; see f90/src/cdr_advect_vec.f90 */

void cdr_advect_vec_(double *mass,double *charge,double *dens_array, 

                 double *d_dens_array,double *er_array,double *ez_array,

		 double *diffusion_coeff,double *dr,double *dz,

		 int *sprite_module,int *r0,int *r1,int *z0,int *z1);

/* Calls hstcrt and checks for errors. If an error occurs, prints the

 corresponding message and exists. */

void 

arcos_cdr_advect_vec(double mass,double charge,double *dens_array, 

                 double *d_dens_array,double *er_array,int er_len,

		 double *ez_array,int ez_len,double diffusion_coeff,

                 double dr,double dz,int sprite_module,

                 int r0,int r1,int z0,int z1)

{

  if ( number_cdr % 1000 == 0 ) {

     fprintf(stdout, "wrap:# calls cdr_advec_vec       =%d\n",number_cdr);

  }

  cdr_advect_vec_(&mass,&charge,dens_array,d_dens_array,er_array,ez_array,

              &diffusion_coeff,&dr,&dz,&sprite_module,&r0,&r1,&z0,&z1);

  number_cdr += 1;

}