projects
/
physik
/
nlsop.git
/ commitdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
| commitdiff |
tree
raw
|
patch
|
inline
| side by side (from parent 1:
1bdb2f4
)
implemented saturation value of carbon in amorph cells
author
hackbard
<hackbard>
Wed, 21 Jul 2004 12:06:25 +0000
(12:06 +0000)
committer
hackbard
<hackbard>
Wed, 21 Jul 2004 12:06:25 +0000
(12:06 +0000)
nlsop.c
patch
|
blob
|
history
nlsop.h
patch
|
blob
|
history
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;
diff --git
a/nlsop.h
b/nlsop.h
index
8b4b63a
..
80ce82b
100644
(file)
--- a/
nlsop.h
+++ b/
nlsop.h
@@
-12,7
+12,7
@@
typedef unsigned int u32;
typedef struct __info
{
typedef struct __info
{
- int cc,steps,range,diff_rate,cpi;
+ int cc,steps,range,diff_rate,cpi
,c_sat
;
double a_el,b_el,a_cd,b_cd,a_ap,b_ap,a_cp,dr_ac,dr_cc;
char z_diff,c_diff;
} info;
double a_el,b_el,a_cd,b_cd,a_ap,b_ap,a_cp,dr_ac,dr_cc;
char z_diff,c_diff;
} info;
@@
-56,6
+56,8
@@
typedef struct __info
#define CPI 20
#define CPI 20
+#define C_SAT (6*6*6)
+
#define FFT_HEIGHT 64
#endif /* NLSOP_H */
#define FFT_HEIGHT 64
#endif /* NLSOP_H */