#define NX 312
#define NY 312
#define NZ 90
#define BC_WIDTH 2
#define I1 BC_WIDTH
#define I2 I1+NX-1
#define J1 BC_WIDTH
#define J2 J1+NY-1
#define K1 BC_WIDTH
#define K2 K1+NZ-1
#define NXDIM NX+2*BC_WIDTH
#define NYDIM NY+2*BC_WIDTH
#define NZDIM NZ+2*BC_WIDTH
#define MAXSTEP 1000
#define WORKSIZE (2*NX*MAXSTEP+NX+MAXSTEP)
#define PI 4.0*atan(1.0)
#define nstep 8000
#define rx1 3500.0
#define ry1 999999.0
#define rz1 3500.0
#define rx2 1250.0
#define ry2 1250.0
#define rz2 3500.0