/[PAMELA software]/gpamela/gpspe/gpudiffusion.F
ViewVC logotype

Diff of /gpamela/gpspe/gpudiffusion.F

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3.1 by bottai, Tue May 2 11:56:23 2006 UTC revision 3.2 by bottai, Tue May 8 15:02:39 2007 UTC
# Line 26  Line 26 
26  #include "gpstripspe.inc"  #include "gpstripspe.inc"
27  #include "gpgeo.inc"  #include "gpgeo.inc"
28  #include "gpgene.inc"  #include "gpgene.inc"
29    #include "gpkey.inc"
30    
31          INTEGER IACT,NUMVOL(20)          INTEGER IACT,NUMVOL(20)
32          REAL DELOSS, TRAPAR(7),xyzspa(3),VPOS(3),xyzspac(3)          REAL DELOSS, TRAPAR(7),xyzspa(3),VPOS(3),xyzspac(3)
# Line 53  Line 54 
54          zup=TSPA(3)-xyzspa(3)          zup=TSPA(3)-xyzspa(3)
55          zdown=TSPA(3)+xyzspa(3)          zdown=TSPA(3)+xyzspa(3)
56    
57    
58          nearstripx=nearstx(xyzspa(1),xyzspa(2))          nearstripx=nearstx(xyzspa(1),xyzspa(2))
59          if(nearstripx.ne.0) then          if(nearstripx.ne.0) then
60    
# Line 69  Line 71 
71  *        v is along -Z; B is along -Y --> shift is along -X  *        v is along -Z; B is along -Y --> shift is along -X
72  *  *
73  *****************************************************************************  *****************************************************************************
74             IF(FFIELD.NE.0) THEN
75            CALL GUFLD(VPOS,BMAGNET)            CALL GUFLD(VPOS,BMAGNET)
76    
77  c  c
78  c to be checked  c to be checked
79  c  c
80            xshift=xyzspa(1)-zdown*hallmob*1.e-4*BMAGNET(2)/10.            xshift=xyzspa(1)+zdown*hallmob*1.e-4*BMAGNET(2)/10.
81            IF(NSPEPLANE.EQ.6) xshift=xyzspa(1)+            IF(NSPEPLANE.EQ.6) xshift=xyzspa(1)-
82       +     zdown*hallmob*1.e-4*BMAGNET(2)/10.       +     zdown*hallmob*1.e-4*BMAGNET(2)/10.
83    
84                      else
85               xshift=xyzspa(1)
86              endif
87                                    
88  *  *
89  *        Now widths of Gaussian functions can be calculated by means of    *        Now widths of Gaussian functions can be calculated by means of  
Legend:
Removed from v.3.1  
changed lines
  Added in v.3.2
  ViewVC Help
Powered by ViewVC 1.1.23