- z_c=get_rand_reject(d3_l.max_z,ne_max,n_e_loss);
- process_cell(&d3_l,x_c,y_c,z_c,&my_info);
+ // z_c=get_rand_reject(d3_l.max_z,ne_max,n_e_loss);
+ z_c=get_rand(d3_l.max_z);
+ nel_z=(URAND_MAX*n_e_loss[z_c])/ne_max;
+ process_cell(&d3_l,x_c,y_c,z_c,&my_info,nel_z);