implemented saturation value of carbon in amorph cells

[physik/nlsop.git] / nlsop.c

diff --git a/nlsop.c b/nlsop.c
index 0e06c77..07313e4 100644 (file)
--- a/nlsop.c
+++ b/nlsop.c
@@ -80,6 +80,7 @@ int usage(void)
  puts("-P <file> \t specify implantation profile file");
  puts("-N <file> \t specify nuclear energy loss profile file");
  printf("-H <value> \t collisions per ion in simulation window (default %d)\n",CPI);
  puts("-P <file> \t specify implantation profile file");
  puts("-N <file> \t specify nuclear energy loss profile file");
  printf("-H <value> \t collisions per ion in simulation window (default %d)\n",CPI);

+ puts("-m <value> \t specify c->a carbon saturation");

  
  return 1;
 }
  
  return 1;
 }


@@ -150,7 +151,7 @@ int distrib_c(d3_lattice *d3_l,info *my_info,int step,double c_ratio,u32 rj_m,u3
    {
     offset=i+j*d3_l->max_x+k*d3_l->max_x*d3_l->max_y;
     /* case amorph: amorph <- cryst diffusion */
    {
     offset=i+j*d3_l->max_x+k*d3_l->max_x*d3_l->max_y;
     /* case amorph: amorph <- cryst diffusion */

-    if(*(d3_l->status+offset)&AMORPH)
+    if(*(d3_l->status+offset)&AMORPH && *(d3_l->extra+offset)<my_info->c_sat)

     {
      for(c=-1;c<=1;c++)
      {
     {
      for(c=-1;c<=1;c++)
      {


@@ -829,6 +830,7 @@ int main(int argc,char **argv)
  my_info.dr_cc=DR_CC;
  my_info.diff_rate=DIFF_RATE;
  my_info.cpi=CPI;
  my_info.dr_cc=DR_CC;
  my_info.diff_rate=DIFF_RATE;
  my_info.cpi=CPI;

+ my_info.c_sat=C_SAT;

  nowait=0;
  quit=0;
  escape=0;
  nowait=0;
  quit=0;
  escape=0;


@@ -942,6 +944,9 @@ int main(int argc,char **argv)
     case 'H':
      my_info.cpi=atoi(argv[++i]);
      break;
     case 'H':
      my_info.cpi=atoi(argv[++i]);
      break;

+    case 'm':
+     my_info.c_sat=atoi(argv[++i]);
+     break;

     default:
      usage();
      return -1;
     default:
      usage();
      return -1;