source/unres/src_MD/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