[½»Á÷] ¡¾×ªÌû¡¿¶¯Á¦Ñ§ÃÉ¿¨·½·¨£¨KMC£©Ä£ÄⱡĤÉú³¤µÄ»ù±¾ÖªÊ¶ ÒÑÓÐ11È˲ÎÓë

Originally posted by zyj8119 at 2010-09-11 11:17:05:
Ò»¡¢±¡Ä¤Éú³¤Ä£ÄâÄ£ÐÍ
Ò»ÀàÊÇÀ©É¢¿ØÖƵijɺËÉú³¤Ä£ÐÍ£¬ÕâÀàÄ£ÐͼÙÉ裬À©É¢µÄÔ­×ÓÒ»µ©µ½´ïÍÅ´Ø»ò¾Û¼¯ÌåµÄ±ßÔµÔò²¢Èë¸ÃÍÅ´Ø¡£
1¡¢À©É¢ÏÞÖƾۼ¯(DLA)Ä£ÐÍ£¬ÕâÊÇ×î³õµÄÄ£ÐÍÆäÖ÷Ҫ˼ÏëΪ£ºÔÚ¾§¸ñÖÐÐÄ´¦·ÅÒ»¸öÖÖÁ£×Ó ...

subroutine  diffusion_conditions ()  
  use variables   
  implicit none

  call Potential_energy_equation() 1  
  if (direction(1)==1) then
  E_xr=v(1)-partical_energy
    if (E_xr<=0) then   
    call Periodic_boundary_conditions ()   
    x=x+1   
    else
    call random_number(num_3)  
      if(num_3       call Periodic_boundary_conditions ()  
      x=x+1   
      else   
      a(x,y)=1
      end if  
    end if
  else if (direction(1)==2) then 
  E_xl=v(2)-partical_energy  
    if (E_xr<=0) then   
    call Periodic_boundary_conditions ()
    x=x-1  
    else   
    call random_number(num_4)   
      if(num_4       call Periodic_boundary_conditions ()  
      x=x-1  
      else
      a(x,y)=1   
      end if   
    end if
  else if (direction(1)==3) then   
  E_yr=v(3)-partical_energy   
    if (E_yr<=0) then  
    call Periodic_boundary_conditions ()  
    y=y+1
    else
    call random_number(num_5)  
      if(num_5       call Periodic_boundary_conditions ()
      y=y+1  
      else   
      a(x,y)=1
      end if h  
    end if
  else if (direction(1)==4) then
  E_yl=v(4)-partical_energy
    if (E_yl<=0) then     
    call Periodic_boundary_conditions ()  
    y=y-1
    else   
    call random_number(num_6)
      if(num_6       call Periodic_boundary_conditions ()  
      y=y-1   
      else
      a(x,y)=1  
      end if
    end if   
  end if
end