db2ed1bb4c9847f60fbc5824cc26cfeb75025fc6
[unres.git] / source / unres / src_MD / src / check_sc_distr.f
1       subroutine check_sc_distr
2       implicit real*8 (a-h,o-z)
3       include 'DIMENSIONS'
4       include 'COMMON.TIME1'
5       include 'COMMON.INTERACT'
6       include 'COMMON.NAMES'
7       include 'COMMON.GEO'
8       include 'COMMON.HEADER'
9       include 'COMMON.CONTROL'
10       logical fail
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
16       print *,'dv=',dv
17       do 10 it=1,1 
18         if (it.eq.10) goto 10 
19         open (20,file=restyp(it)//'_distr.sdc',status='unknown')
20         call gen_side(it,90.0D0*deg2rad,al,om,fail)
21         close (20)
22         goto 10
23         open (20,file=restyp(it)//'_distr1.sdc',status='unknown')
24         do i=0,90
25           do j=0,72
26             prob(j,i)=0.0D0
27           enddo
28         enddo
29         do isample=1,MaxSample
30           call gen_side(it,90.0D0*deg2rad,al,om)
31           indal=rad2deg*al/2
32           indom=(rad2deg*om+180.0D0)/5
33           prob(indom,indal)=prob(indom,indal)+delt
34         enddo
35         do i=45,90
36           do j=0,72 
37             write (20,'(2f10.3,1pd15.5)') 2*i+0.0D0,5*j-180.0D0,
38      &              prob(j,i)/dV
39           enddo
40         enddo
41    10   continue
42       return
43       end