1 subroutine check_sc_distr
2 implicit real*8 (a-h,o-z)
5 include 'COMMON.INTERACT'
8 include 'COMMON.HEADER'
9 include 'COMMON.CONTROL'
11 double precision varia(maxvar)
12 double precision hrtime,mintime,sectime
13 parameter (MaxSample=10000000,delt=1.0D0/MaxSample)
14 dimension prob(0:72,0:90)
15 dV=2.0D0*5.0D0*deg2rad*deg2rad
19 open (20,file=restyp(it)//'_distr.sdc',status='unknown')
20 call gen_side(it,90.0D0*deg2rad,al,om,fail)
23 open (20,file=restyp(it)//'_distr1.sdc',status='unknown')
29 do isample=1,MaxSample
30 call gen_side(it,90.0D0*deg2rad,al,om)
32 indom=(rad2deg*om+180.0D0)/5
33 prob(indom,indal)=prob(indom,indal)+delt
37 write (20,'(2f10.3,1pd15.5)') 2*i+0.0D0,5*j-180.0D0,