Changeset 16 in ZHANGProjects


Ignore:
Timestamp:
Aug 25, 2014, 4:10:56 PM (10 years ago)
Author:
zhangj
Message:

Fix the bug to generate random number move_am_() in the Fortran file.

Location:
ICOSIM/CPP/trunk/source
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • ICOSIM/CPP/trunk/source/crystal_dan_FINAL_VERSION_CHvsVR.f

    r15 r16  
    858858     
    859859      kxmcs = DYA*ran_x
    860       kxmcs = DYA*ran_x
     860      kymcs = DYA*ran_y
    861861     
    862862c      kxmcs = DYA*RAN_GAUSS(1.)
    863 c      kymcs = DYA*RAN_GAUSS(1.)
     863c     kymcs = DYA*RAN_GAUSS(1.)
    864864
    865865      XP = xp+kxmcs
    866866      yp = yp+kymcs
    867867     
    868        write(*,*) "ran_x = ", ran_x, "ran_y = ", ran_y
    869 c      write(*,*) "kxmcs = ", kxmcs, "kymcs = ", kymcs
    870 c      write(*,*)'xp + kxmcs:', XP, 'yp + kxmcs', YP
     868c       write(*,*) "ran_x = ", ran_x, "ran_y = ", ran_y
     869      write(*,*) "kxmcs = ", kxmcs, "kymcs = ", kymcs
     870      write(*,*)'xp + kxmcs:', XP, 'yp + kxmcs', YP
    871871     
    872872cc     XP    = XP+DYA*RAN_GAUSS(1.)
    … …  
    938938      endif
    939939
     940      write(*,*) "ichoix = ", ichoix
     941       write(*,*)'xp final:', xp, 'yp final', yp
     942       
    940943!---------- calculate the related kick -----------
    941944

  • ICOSIM/CPP/trunk/source/simcrys.cc

    r15 r16  
    10001000                     
    10011001              crys.Cry_length = length;
     1002             
     1003                  x=part.x;
     1004                  xp=part.xp;
     1005                  z=part.y;
     1006                  zp=part.yp;
     1007                  p=part.PC;
     1008                  s=part.s;
    10021009                   
    10031010               
    … …  
    11781185             
    11791186
     1187                 
     1188                  x=part.x;
     1189                  xp=part.xp;
     1190                  z=part.y;
     1191                  zp=part.yp;
     1192                  p=part.PC;
     1193                   s=part.s;
    11801194             
    11811195                        /*
    … …  
    14351449
    14361450          // =================================================================================FIN DE RETOUR AU COORDONEE DE ICOSIM====================================================================
    1437                 /*
     1451               
    14381452                  part.xp=xp;
    14391453                  part.yp=zp;
    … …  
    14411455                  part.x=x;
    14421456                  part.PC=p;
    1443                 */
     1457               
    14441458          xp1 = part.xp;
    14451459
Note: See TracChangeset for help on using the changeset viewer.