fixed nlsop for packaging ..

[physik/nlsop.git] / nlsop_client.c

diff --git a/nlsop_client.c b/nlsop_client.c
index ad1c84f..beeda9f 100644 (file)
--- a/nlsop_client.c
+++ b/nlsop_client.c
@@ -67,6 +67,7 @@ info *gmy_info;
 int *gi;
 unsigned char dc;
 unsigned char shut_down;
+double c_ratio;
 
 int get_data_and_calc(t_event *event,void *allineed);
 int nop(t_event *event,void *allineed);
@@ -104,6 +105,47 @@ int sputter(d3_lattice *d3_l)
  return 1;
 }
 
+int get_c_ratio(double *c_ratio,char *pfile,info *my_info,d3_lattice *d3_l)
+{
+ double all,a,b,d;
+ int i,k;
+ int p_fd;
+ unsigned char buf[32];
+ char *p;
+ unsigned char c;
+
+ if((p_fd=open(pfile,O_RDONLY))<0)
+ {
+  puts("cannot open profile file");
+  return -1;
+ }
+ k=1;
+ d=0;
+ all=0;
+ while(k)
+ {
+  for(i=0;i<32;i++)
+  {
+   k=read(p_fd,&c,1);
+   buf[i]=c;
+   if(c=='\n') break;
+  }
+  if(k)
+  {
+   p=strtok(buf," ");
+   a=atof(p)/10; /* nm */
+   p=strtok(NULL," ");
+   b=atof(p);
+   if(a<=d3_l->max_z*CELL_LENGTH) d+=b;
+   all+=b;
+  }
+ }
+ *c_ratio=d/all;
+ close(p_fd);
+
+ return 1;
+}
+
 int process_cell(d3_lattice *d3_l,u32 x,u32 y,u32 z,info *my_info,u32 nel_z)
 {
  unsigned char *thiz;
@@ -163,11 +205,13 @@ 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 */
- 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);
- *(d3_l->extra+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)+=1;
- (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);
+   *(d3_l->extra+x+y*d3_l->max_x+z*d3_l->max_x*d3_l->max_y)+=1;
+   (my_info->cc)++;
+ }
 
  if(step%my_info->diff_rate==0)
  {
@@ -452,6 +496,7 @@ int get_data_and_calc(t_event *event,void *allineed) {
   ne_max=0;
   ip_max=0;
   do_sputter=1;
+  c_ratio=0;
 
   net=gnet;
   gd3_l=&d3_l;
@@ -575,6 +620,12 @@ int get_data_and_calc(t_event *event,void *allineed) {
     printf("warning: max_z (%d) too small, there may be amorphous volumes\n",
            d3_l.max_z);
 
+  /* 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);
+  }
+
   /*
    * start of simulation
    */