Changeset 17 in ZHANGProjects for ICOSIM/CPP/trunk/source/crystal_dan_FINAL_VERSION_CHvsVR.f

Timestamp:

Sep 30, 2014, 11:05:33 AM (10 years ago)

Author:

zhangj

Message:

Added comments...

File:

• ICOSIM/CPP/trunk/source/crystal_dan_FINAL_VERSION_CHvsVR.f (modified) (10 diffs)

ICOSIM/CPP/trunk/source/crystal_dan_FINAL_VERSION_CHvsVR.f

@@ -453 +453 @@
           
 c          write(*,*) "y = ", y
-           write(*,*) "before enter VR: x = ", x, "  xp = ",xp
-           write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
+c           write(*,*) "before enter VR: x = ", x, "  xp = ",xp
+c           write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
             
           xp=xp+0.45*(xp/xpcrit+1)*Ang_avr
@@ -465 +465 @@
      +      xp ,yp,PC)
      
-          write(*,*) "after enter VR: x = ", x, "  xp = ",xp
-          write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
+c          write(*,*) "after enter VR: x = ", x, "  xp = ",xp
+c          write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
           
           x = x + 0.5*s_length * xp
@@ -501 +501 @@
 !     +          DLRi(IS),xp ,yp,PC)
 
-            write(*,*) "before enter VR: x = ", x, "  xp = ",xp
-            write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
+c            write(*,*) "before enter VR: x = ", x, "  xp = ",xp
+c            write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
            
             CALL CALC_ION_LOSS(IS,PC,s_length-Srefl,DESt)
             
-            write(*,*) "before enter VR: x = ", x, "  xp = ",xp
-            write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
-            write(*,*) "IS = ", IS, " NAM = ", NAM
-            write(*,*)  "s_length-Srefl = ",s_length-Srefl
-            write(*,*) "DESt = ", DESt, " DLYi(IS) = ", DLYi(IS)
-            write(*,*) "DLRi(IS) = ", DLRi(IS)
+c            write(*,*) "before enter VR: x = ", x, "  xp = ",xp
+c            write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
+c            write(*,*) "IS = ", IS, " NAM = ", NAM
+c            write(*,*)  "s_length-Srefl = ",s_length-Srefl
+c            write(*,*) "DESt = ", DESt, " DLYi(IS) = ", DLYi(IS)
+c            write(*,*) "DLRi(IS) = ", DLRi(IS)
             
             CALL MOVE_AM_(IS,NAM,s_length-Srefl,DESt,DLYi(IS),
      +          DLRi(IS),xp ,yp,PC)
      
-            write(*,*) "after enter VR: x = ", x, "  xp = ",xp
-            write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
+c            write(*,*) "after enter VR: x = ", x, "  xp = ",xp
+c            write(*,*) "y = ",y, "yp = ", yp, "PC = ", PC
            
             x = x + 0.5 * xp * (s_length - Srefl)
@@ -618 +618 @@
 c      endif
 c      
-      WRITE(*,*)'Crystal process: ',PROC
+c      WRITE(*,*)'Crystal process: ',PROC
 c      write(*,*)'xp_final :', xp
 c      WRITE(*,*)'Crystal process: ',PROC,'Chann Angle',Ch_angle/1000,   +
@@ -838 +838 @@
 
      
-      write(*,*) "In Move_AM(): "
+c      write(*,*) "In Move_AM(): "
 
 
       xp=xp*1000
       yp=yp*1000
-      write(*,*)'xp initial:', xp, 'yp initial', yp
+c      write(*,*)'xp initial:', xp, 'yp initial', yp
 C. DEI - dE/dx stoping energy
       PC    = PC - DEI*DZ    ! energy lost because of ionization process[GeV]
@@ -849 +849 @@
 C. DYA - rms of coloumb scattering
       DYA = (13.6/PC)*SQRT(DZ/DLr)             !MCS (mrad)
-      write(*,*)'dya= ',dya
+c      write(*,*)'dya= ',dya
       
       ran_x = RAN_GAUSS(1.)
-      write(*,*) "ran_x = ", ran_x
+c      write(*,*) "ran_x = ", ran_x
       
       ran_y = RAN_GAUSS(1.)
-      write(*,*) "ran_y = ", ran_y
+c      write(*,*) "ran_y = ", ran_y
       
       kxmcs = DYA*ran_x
@@ -867 +867 @@
       
 c       write(*,*) "ran_x = ", ran_x, "ran_y = ", ran_y
-      write(*,*) "kxmcs = ", kxmcs, "kymcs = ", kymcs
-      write(*,*)'xp + kxmcs:', XP, 'yp + kxmcs', YP
+c      write(*,*) "kxmcs = ", kxmcs, "kymcs = ", kymcs
+c      write(*,*)'xp + kxmcs:', XP, 'yp + kxmcs', YP
       
 cc     XP    = XP+DYA*RAN_GAUSS(1.)
@@ -938 +938 @@
       endif
 
-      write(*,*) "ichoix = ", ichoix
-       write(*,*)'xp final:', xp, 'yp final', yp
+c      write(*,*) "ichoix = ", ichoix
+c       write(*,*)'xp final:', xp, 'yp final', yp
        
 !---------- calculate the related kick -----------
@@ -980 +980 @@
 
 
-      write(*,*) "tx = ", tx, " tz = ", tz
+c      write(*,*) "tx = ", tx, " tz = ", tz
       
       XP = XP + tx
@@ -990 +990 @@
       yp=yp/1000
 
-      write(*,*)'xp final:', xp, 'yp final', yp
+c      write(*,*)'xp final:', xp, 'yp final', yp
 !-----------------------------------------------
 ! print out the energy loss and process experienced