radnce.f

Go to the documentation of this file.

      subroutine radnce(pi,conv,bmu,the,gio,ii,jpart,jphi)              
c                                                                       
c***********************************************************************
c.....declaration statements............................................
      implicit real*8 (a-h,o-z)
      real*8 pi,conv,bmu(51)                                            
      real*8 gio(4,50,46),the(50)                                        
      integer*4 jphi(91)                                                
c***********************************************************************
      jpartm=jpart-1
c     jpartm=1
      do 97 j=1,jpart,jpartm
      qp=0.0                                                            
      plr=0.0                                                           
      bett=0.0                                                          
      ft=gio(1,ii,j)+gio(2,ii,j)                                      
      ftd=abs(gio(1,ii,j)-gio(2,ii,j) )                               
      if(ft .le.1.0e-9) go to 107                                      
      qp=(gio(1,ii,j)-gio(2,ii,j))**2+gio(3,ii,j)**2+gio(4,ii,j)**2 
      sqp=sqrt(qp)                                                      
      plr=sqp/ft                                                        
      bett=.5*dasin(gio(4,ii,j)/sqp)                                    
107   continue                                                          
      tb=dsin(bett)/dcos(bett)                                            
       ftpi=ft*pi                                                       
      write(6,364)the(ii),bmu(ii),jphi(j),(gio(n,ii,j),n=1,3),ftpi,
     1  plr*100.0
97    continue                                                          
      write(6,46)                                                       
c***********************************************************************
46     format(1x,'    '  )                                              
364    format(1x,f5.1,1x,f7.4,1x,i3,1x,1p4e11.3,0pf7.2)                        
c***********************************************************************
      return                                                            
      end                                                               
c***********************************************************************