cdr_advect_diffu_wrap.c

Go to the documentation of this file.

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

}