wham and cluster_wham Adam's new constr_dist multichain
[unres.git] / source / unres / src_MD-M / checkvar.f
1       logical function check_var(var,info)
2       implicit real*8 (a-h,o-z)
3       include 'DIMENSIONS'
4       include 'COMMON.VAR'
5       include 'COMMON.IOUNITS'
6       include 'COMMON.GEO'
7       include 'COMMON.SETUP'
8       dimension var(maxvar)
9       dimension info(3)
10 C AL -------
11        check_var=.false.
12        do i=nphi+ntheta+1,nphi+ntheta+nside
13 ! Check the side chain "valence" angles alpha
14          if (var(i).lt.1.0d-7) then 
15            write (iout,*) 'CHUJ NASTAPIL ABSOLUTNY!!!!!!!!!!!!'
16            write (iout,*) 'Processor',me,'received bad variables!!!!'
17            write (iout,*) 'Variables'
18            write (iout,'(8f10.4)') (rad2deg*var(j),j=1,nvar)
19            write (iout,*) 'Continuing calculations at this point',
20      & ' could destroy the results obtained so far... ABORTING!!!!!!'
21            write (iout,'(a19,i5,f10.4,a4,2i4,a3,i3)') 
22      &     'valence angle alpha',i-nphi-ntheta,var(i),
23      &     'n it',info(1),info(2),'mv ',info(3)
24            write (*,*) 'CHUJ NASTAPIL ABSOLUTNY!!!!!!!!!!!!'
25            write (*,*) 'Processor',me,'received bad variables!!!!'
26            write (*,*) 'Variables'
27            write (*,'(8f10.4)') (rad2deg*var(j),j=1,nvar)
28            write (*,*) 'Continuing calculations at this point',
29      & ' could destroy the results obtained so far... ABORTING!!!!!!'
30            write (*,'(a19,i5,f10.4,a4,2i4,a3,i3)') 
31      &     'valence angle alpha',i-nphi-ntheta,var(i),
32      &     'n it',info(1),info(2),'mv ',info(3)
33            check_var=.true.
34            return
35          endif
36        enddo
37 ! Check the backbone "valence" angles theta
38        do i=nphi+1,nphi+ntheta
39          if (var(i).lt.1.0d-7) then
40            write (iout,*) 'CHUJ NASTAPIL ABSOLUTNY!!!!!!!!!!!!'
41            write (iout,*) 'Processor',me,'received bad variables!!!!'
42            write (iout,*) 'Variables'
43            write (iout,'(8f10.4)') (rad2deg*var(j),j=1,nvar)
44            write (iout,*) 'Continuing calculations at this point',
45      & ' could destroy the results obtained so far... ABORTING!!!!!!'
46            write (iout,'(a19,i5,f10.4,a4,2i4,a3,i3)') 
47      &     'valence angle theta',i-nphi,var(i),
48      &     'n it',info(1),info(2),'mv ',info(3)
49            write (*,*) 'CHUJ NASTAPIL ABSOLUTNY!!!!!!!!!!!!'
50            write (*,*) 'Processor',me,'received bad variables!!!!'
51            write (*,*) 'Variables'
52            write (*,'(8f10.4)') (rad2deg*var(j),j=1,nvar)
53            write (*,*) 'Continuing calculations at this point',
54      & ' could destroy the results obtained so far... ABORTING!!!!!!'
55            write (*,'(a19,i5,f10.4,a4,2i4,a3,i3)') 
56      &     'valence angle theta',i-nphi,var(i),
57      &     'n it',info(1),info(2),'mv ',info(3)
58            check_var=.true.
59            return
60          endif
61        enddo
62        return
63        end