for(j=0;j<my_info.cpi;j++) {
x_c=get_rand(d3_l.max_x);
y_c=get_rand(d3_l.max_y);
- // z_c=get_rand_reject(d3_l.max_z,ne_max,n_e_loss);
+#ifdef NEL_Z
+ z_c=get_rand_reject(d3_l.max_z,ne_max,n_e_loss);
+#else
z_c=get_rand(d3_l.max_z);
+#endif
process_cell(&d3_l,x_c,y_c,z_c,&my_info,nel_z[z_c]);
}
distrib_c(&d3_l,&my_info,i,ip_max,c_profile);