oups, c conc will never increase. fixed ! :)

[physik/nlsop.git] / nlsop_client.c

diff --git a/nlsop_client.c b/nlsop_client.c
index 770ecb9..beeda9f 100644 (file)
--- a/nlsop_client.c
+++ b/nlsop_client.c
@@ -205,7 +205,7 @@ int distrib_c(d3_lattice *d3_l,info *my_info,int step,u32 rj_m,u32 *rj_g)
  int carry;
 
  /* put one c ion somewhere in the lattice */
- if(c_ratio*step<my_info->cc) {
+ if(c_ratio*step<=my_info->cc) {
    x=get_rand(d3_l->max_x);
    y=get_rand(d3_l->max_y);
    z=get_rand_reject(d3_l->max_z,rj_m,rj_g);
@@ -623,6 +623,7 @@ int get_data_and_calc(t_event *event,void *allineed) {
   /* maybe we need a ratio! */
   if(!do_sputter) {
     get_c_ratio(&c_ratio,p_file,&my_info,&d3_l);
+    printf("calculated c ratio: %f\n",c_ratio);
   }
 
   /*