+++ /dev/null
- subroutine check_sc_distr
- implicit real*8 (a-h,o-z)
- include 'DIMENSIONS'
- include 'COMMON.TIME1'
- include 'COMMON.INTERACT'
- include 'COMMON.NAMES'
- include 'COMMON.GEO'
- include 'COMMON.HEADER'
- include 'COMMON.CONTROL'
- logical fail
- double precision varia(maxvar)
- double precision hrtime,mintime,sectime
- parameter (MaxSample=10000000,delt=1.0D0/MaxSample)
- dimension prob(0:72,0:90)
- dV=2.0D0*5.0D0*deg2rad*deg2rad
- print *,'dv=',dv
- do 10 it=1,1
- if (it.eq.10) goto 10
- open (20,file=restyp(it)//'_distr.sdc',status='unknown')
- call gen_side(it,90.0D0*deg2rad,al,om,fail)
- close (20)
- goto 10
- open (20,file=restyp(it)//'_distr1.sdc',status='unknown')
- do i=0,90
- do j=0,72
- prob(j,i)=0.0D0
- enddo
- enddo
- do isample=1,MaxSample
- call gen_side(it,90.0D0*deg2rad,al,om)
- indal=rad2deg*al/2
- indom=(rad2deg*om+180.0D0)/5
- prob(indom,indal)=prob(indom,indal)+delt
- enddo
- do i=45,90
- do j=0,72
- write (20,'(2f10.3,1pd15.5)') 2*i+0.0D0,5*j-180.0D0,
- & prob(j,i)/dV
- enddo
- enddo
- 10 continue
- return
- end