+ this_cell=cell_p+x+y*x_max+z*x_max*y_max;
+ pressure=p0*URAND_2BYTE_MAX;
+
+ for(i=-range;i<=range;i++)
+ {
+ for(j=-range;j<=range;j++)
+ {
+ // z relaxation, no pressure in z direction
+ // for(k=-range;k<=range;k++)
+ // {
+ if(!(i==0 && j==0)) // && k==0))
+ {
+ // if((cell_p+((x+x_max+i)%x_max)+((y+j+y_max)%y_max)*x_max+((z+k+z_max)%z_max)*x_max*y_max)->status&AMORPH) pressure+=faktor*URAND_2BYTE_MAX/(i*i+j*j+k*k);
+ if((cell_p+((x+x_max+i)%x_max)+((y+j+y_max)%y_max)*x_max+z*x_max*y_max)->status&AMORPH) pressure+=faktor*URAND_2BYTE_MAX/(i*i+j*j);
+ }
+ // }
+ }
+ }
+ pressure*=this_cell->conc;
+ if(this_cell->status&AMORPH)
+ {
+ /* wrong probability! just test by now ...*/
+ if(rand_get(URAND_2BYTE_MAX)>pressure)
+ {
+ make_cryst(this_cell);
+ *c_conc=*c_conc+1+this_cell->conc;
+ } else *c_conc+=1;
+ } else
+ {
+ if(rand_get(URAND_2BYTE_MAX)<=pressure)
+ {
+ make_amorph(this_cell);
+ *c_conc=*c_conc+1-this_cell->conc;
+ } else *c_conc+=1;
+ }