Line | |
---|
1 | ! Subroutine initlund: Initializes data structures used frequently
|
---|
2 | ! in connection with Biased Gaussian Steps and the energy functions
|
---|
3 | ! from Anders Irback's protein folding model. Calls: none.
|
---|
4 | !
|
---|
5 | subroutine init_lund
|
---|
6 | include 'INCL.H'
|
---|
7 | include 'incl_lund.h'
|
---|
8 | integer i, npprs, j, k
|
---|
9 |
|
---|
10 | logical bgsposs
|
---|
11 | do i=1,mxrs
|
---|
12 | iN(i)=-1
|
---|
13 | iCa(i)=-1
|
---|
14 | iC(i)=-1
|
---|
15 | iphi(i)=-34
|
---|
16 | ipsi(i)=-35
|
---|
17 | enddo
|
---|
18 | ! print *,'total number of variables = ',nvr
|
---|
19 |
|
---|
20 | do i=1,ntlml
|
---|
21 | npprs=1
|
---|
22 | do j=ivrml1(i),ivrml1(i)+nvrml(i)-1
|
---|
23 | mlvr(j)=i
|
---|
24 | if (nmvr(j).eq.'phi') then
|
---|
25 | iphi(npprs)=j
|
---|
26 | ! Now if the residue is a proline, there is no phi angle in the variable
|
---|
27 | ! list in SMMP, and so iphi(j) will remain at the initial value.
|
---|
28 | ! So, make sure you never use iphi(i) for proline.
|
---|
29 | endif
|
---|
30 | if (nmvr(j).eq.'psi'.or.nmvr(j).eq.'pst') then
|
---|
31 | ipsi(npprs)=j
|
---|
32 | npprs=npprs+1
|
---|
33 | endif
|
---|
34 | enddo
|
---|
35 | do j=irsml1(i),irsml2(i)
|
---|
36 | iN(j)=iatrs1(j)
|
---|
37 | do k=iatrs1(j),iatrs2(j)
|
---|
38 | if (nmat(k)(1:2).eq.'ca') then
|
---|
39 | iCa(j)=k
|
---|
40 | endif
|
---|
41 | if (nmat(k)(1:1).eq.'c') then
|
---|
42 | iC(j)=k
|
---|
43 | endif
|
---|
44 | enddo
|
---|
45 | ! print *,'determined phi,psi serial for residue',j,' as '
|
---|
46 | ! # ,iphi(j),ipsi(j)
|
---|
47 | enddo
|
---|
48 | enddo
|
---|
49 | abgs=300.0
|
---|
50 | bbgs=10.0
|
---|
51 | bgsnvar=0
|
---|
52 | ! JHM: Took the following lines out to avoid another dependency.
|
---|
53 | ! do i=1,nvr
|
---|
54 | ! if (bgsposs(i)) then
|
---|
55 | ! bgsnvar=bgsnvar+1
|
---|
56 | ! bgsvar(bgsnvar)=i
|
---|
57 | ! endif
|
---|
58 | ! enddo
|
---|
59 | end subroutine init_lund
|
---|
60 |
|
---|
Note:
See
TracBrowser
for help on using the repository browser.