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 | logical bgsposs
|
---|
9 | do i=1,mxrs
|
---|
10 | iN(i)=-1
|
---|
11 | iCa(i)=-1
|
---|
12 | iC(i)=-1
|
---|
13 | iphi(i)=-34
|
---|
14 | ipsi(i)=-35
|
---|
15 | enddo
|
---|
16 | ! print *,'total number of variables = ',nvr
|
---|
17 |
|
---|
18 | do i=1,ntlml
|
---|
19 | npprs=1
|
---|
20 | do j=ivrml1(i),ivrml1(i)+nvrml(i)-1
|
---|
21 | mlvr(j)=i
|
---|
22 | if (nmvr(j).eq.'phi') then
|
---|
23 | iphi(npprs)=j
|
---|
24 | ! Now if the residue is a proline, there is no phi angle in the variable
|
---|
25 | ! list in SMMP, and so iphi(j) will remain at the initial value.
|
---|
26 | ! So, make sure you never use iphi(i) for proline.
|
---|
27 | endif
|
---|
28 | if (nmvr(j).eq.'psi'.or.nmvr(j).eq.'pst') then
|
---|
29 | ipsi(npprs)=j
|
---|
30 | npprs=npprs+1
|
---|
31 | endif
|
---|
32 | enddo
|
---|
33 | do j=irsml1(i),irsml2(i)
|
---|
34 | iN(j)=iatrs1(j)
|
---|
35 | do k=iatrs1(j),iatrs2(j)
|
---|
36 | if (nmat(k)(1:2).eq.'ca') then
|
---|
37 | iCa(j)=k
|
---|
38 | endif
|
---|
39 | if (nmat(k)(1:1).eq.'c') then
|
---|
40 | iC(j)=k
|
---|
41 | endif
|
---|
42 | enddo
|
---|
43 | ! print *,'determined phi,psi serial for residue',j,' as '
|
---|
44 | ! # ,iphi(j),ipsi(j)
|
---|
45 | enddo
|
---|
46 | enddo
|
---|
47 | abgs=300.0
|
---|
48 | bbgs=10.0
|
---|
49 | bgsnvar=0
|
---|
50 | ! JHM: Took the following lines out to avoid another dependency.
|
---|
51 | ! do i=1,nvr
|
---|
52 | ! if (bgsposs(i)) then
|
---|
53 | ! bgsnvar=bgsnvar+1
|
---|
54 | ! bgsvar(bgsnvar)=i
|
---|
55 | ! endif
|
---|
56 | ! enddo
|
---|
57 | end subroutine init_lund
|
---|
58 |
|
---|
Note:
See
TracBrowser
for help on using the repository browser.