[bd2278d] | 1 | !**************************************************************
|
---|
| 2 | !
|
---|
| 3 | ! This file contains the subroutines: redvar
|
---|
| 4 | !
|
---|
| 5 | ! Copyright 2003-2005 Frank Eisenmenger, U.H.E. Hansmann,
|
---|
| 6 | ! Shura Hayryan, Chin-Ku
|
---|
| 7 | ! Copyright 2007 Frank Eisenmenger, U.H.E. Hansmann,
|
---|
| 8 | ! Jan H. Meinke, Sandipan Mohanty
|
---|
| 9 | !
|
---|
| 10 | ! **************************************************************
|
---|
[e40e335] | 11 |
|
---|
| 12 |
|
---|
| 13 | subroutine redvar
|
---|
| 14 |
|
---|
[bd2278d] | 15 | ! ...................................................................
|
---|
| 16 | !
|
---|
| 17 | ! PURPOSE: Read global parameters for molecules from lines
|
---|
| 18 | !
|
---|
| 19 | ! +--------------------------------------------------+
|
---|
| 20 | ! |@ molecule no. : six floats separated by commas |
|
---|
| 21 | ! +--------------------------------------------------+
|
---|
| 22 | !
|
---|
| 23 | ! NB: 1) if omit field with molecule no. assume: nml=1
|
---|
| 24 | ! 2) last 3 float are angles in deg.
|
---|
| 25 | !
|
---|
| 26 | ! Read and interpret file to SET and FIX internal variables
|
---|
| 27 | ! by commands:
|
---|
| 28 | !
|
---|
| 29 | ! +-----------------------------------------+
|
---|
| 30 | ! | molecule : residue : variable : value |
|
---|
| 31 | ! +-----------------------------------------+
|
---|
| 32 | !
|
---|
| 33 | ! * Lines containing '&' assign FIXED variable(s), they will
|
---|
| 34 | ! not be varied during subsequent minimization etc.
|
---|
| 35 | !
|
---|
| 36 | ! * Empty LINES or lines containing '#' are ignored
|
---|
| 37 | ! * Several commands on same line must be separated by ';'
|
---|
| 38 | ! * Empty COMMANDS, i.e. ' : : ' are ignored
|
---|
| 39 | ! * All spaces are not significant and are therefore ignored
|
---|
| 40 | !
|
---|
| 41 | ! * A command consists of up to 4 (maxfld) fields, separated
|
---|
| 42 | ! by ':'
|
---|
| 43 | !
|
---|
| 44 | ! - last field : value for VARIABLE (REAL)
|
---|
| 45 | ! ! should never be empty
|
---|
| 46 | ! - 1st before last: name(s) (CHAR) or index(ices) of VARIABLE(S)
|
---|
| 47 | ! - 2nd before last: name(s) or index(ices) of RESIDUE(S)
|
---|
| 48 | ! - 3rd before last: name or number(ices) of MOLECULE(S)
|
---|
| 49 | !
|
---|
| 50 | ! * molecules, residues, variables can be identified, either by,
|
---|
| 51 | ! INDICES (zones 'n1-n2' possible) or NAMES
|
---|
| 52 | !
|
---|
| 53 | ! * several identifiers in a field can be separated by ','
|
---|
| 54 | !
|
---|
| 55 | ! * INDICES: for residues - refer to numbering within molecule
|
---|
| 56 | ! : for variables - refer to numbering within residue
|
---|
| 57 | ! * ZONES: '-n2' indicates '1-n2'
|
---|
| 58 | ! 'n1-' indicates 'n1-(all)'
|
---|
| 59 | ! * NAMES or their ends can be indicated by wild-card '*'
|
---|
| 60 | ! are case-sensitive
|
---|
| 61 | !
|
---|
| 62 | ! Example: phi:-65; psi:-45 >set all phi=-65, all psi=-45
|
---|
| 63 | ! om*: 180 & >set all omg, omt ... to 180 & fix them
|
---|
| 64 | ! 5 : x* : -60 >set all xi-angles of residue 5 to 60
|
---|
| 65 | !
|
---|
| 66 | ! CALLS: setvar,extstr,iendst,ibegst,iopfil,iredin,iredrl
|
---|
| 67 | ! ......................................................................
|
---|
[e40e335] | 68 |
|
---|
| 69 |
|
---|
| 70 | include 'INCL.H'
|
---|
| 71 |
|
---|
[cb47b9c] | 72 | ! functions
|
---|
| 73 | integer iopfil, iendst, iredin, iredrl, ibegst
|
---|
[bd2278d] | 74 | ! maxfld: max. # of fields in one command
|
---|
| 75 | ! maxide: max. # of identifiers in a field
|
---|
| 76 | ! maxcmd: max. # of commands to be interpreted
|
---|
| 77 | ! ilrg: a large integer
|
---|
[e40e335] | 78 |
|
---|
[cb47b9c] | 79 | integer maxfld, maxide, maxcmd, ilrg, ifdend, icb, i, ib, ibz
|
---|
| 80 | integer ife, ifb, ide, id, ice, ie, ieh, ieh1, iez, ihz, ifx, ifld
|
---|
| 81 | integer ile, ii, ihy, ilb, iml, it, in ,io, inum, kbz, iv, ity, jb
|
---|
| 82 | integer j, k, l, kez, kk, kv, lez, lbz,ll, ll1, ll2, n, nfi, ncmd
|
---|
| 83 | integer nfld, nide, nml, ntlvr
|
---|
| 84 | double precision vlvrx, rn, vr, val
|
---|
[e40e335] | 85 | parameter (maxfld=4,
|
---|
[bd2278d] | 86 | & maxide=30,
|
---|
| 87 | & maxcmd=5000,
|
---|
| 88 | & ilrg=1000000)
|
---|
[e40e335] | 89 |
|
---|
| 90 | character spcm,spfd,spcc,sphy,cmt,wdc,sfix,blnk, sglp,
|
---|
[bd2278d] | 91 | & line*132,lincmd*132,linfld(maxfld)*132,linide*132,
|
---|
| 92 | & linh*132,strg(6)*17
|
---|
[e40e335] | 93 | dimension ifdend(maxfld),vlvrx(mxvr),rn(6)
|
---|
| 94 | logical fix,did,exa,forml(mxml),forrs(mxrs),forvr(mxvr),
|
---|
[bd2278d] | 95 | & stvr(mxvr)
|
---|
[e40e335] | 96 | data spcm/';'/,spfd/':'/,spcc/','/,sphy/'-'/,cmt/'#'/,wdc/'*'/,
|
---|
[bd2278d] | 97 | & sfix/'&'/,blnk/' '/, sglp/'@'/
|
---|
[e40e335] | 98 |
|
---|
| 99 |
|
---|
[bd2278d] | 100 | ! ___________________________________ Checks
|
---|
[e40e335] | 101 | ntlvr=ivrml1(ntlml)+nvrml(ntlml)-1
|
---|
| 102 | if (ntlvr.eq.0) then
|
---|
[38d77eb] | 103 | write(logString,*)'redvar> No variables defined in molecule(s)'
|
---|
[e40e335] | 104 | return
|
---|
| 105 | endif
|
---|
[bd2278d] | 106 | ! ___________________________________ Initialize
|
---|
[e40e335] | 107 |
|
---|
| 108 | io=iopfil(lunvar,varfil,'old','formatted')
|
---|
| 109 | if (io.eq.0) then
|
---|
[38d77eb] | 110 | write (logString, '(a,/,a,i3,2a)')
|
---|
[bd2278d] | 111 | & ' redvar> ERROR opening file to set variables:',
|
---|
| 112 | & ' LUN=',lunvar,' FILE=',varfil(1:iendst(varfil))
|
---|
[e40e335] | 113 | stop
|
---|
| 114 | elseif (io.eq.-1) then
|
---|
| 115 | return
|
---|
| 116 | endif
|
---|
[bd2278d] | 117 | ! ___________________________________ Initialization
|
---|
[e40e335] | 118 | do i=1,ntlml
|
---|
| 119 | forml(i)=.true.
|
---|
| 120 | do j=irsml1(i),irsml2(i)
|
---|
| 121 | forrs(j)=.true.
|
---|
| 122 | enddo
|
---|
| 123 | enddo
|
---|
| 124 |
|
---|
| 125 | do i=1,ntlvr
|
---|
| 126 | fxvr(i)=.false.
|
---|
| 127 | forvr(i)=.true.
|
---|
| 128 | stvr(i)=.false.
|
---|
| 129 |
|
---|
| 130 | it=ityvr(i) ! var. type
|
---|
| 131 | if (it.eq.3) then ! torsion
|
---|
| 132 | vr=toat(iatvr(i))
|
---|
| 133 | elseif (it.eq.2) then ! b.angle
|
---|
| 134 | vr=baat(iatvr(i))
|
---|
| 135 | elseif (it.eq.1) then ! b.length
|
---|
| 136 | vr=blat(iatvr(i))
|
---|
| 137 | else
|
---|
[38d77eb] | 138 | write(logString,*) 'redvar> unknown variable type: ',it,' !'
|
---|
[e40e335] | 139 | stop
|
---|
| 140 | endif
|
---|
| 141 | vlvrx(i)=vr
|
---|
| 142 | enddo
|
---|
| 143 |
|
---|
| 144 | ncmd=0
|
---|
| 145 |
|
---|
| 146 | 1 read (lunvar,'(a)',end=2) line
|
---|
| 147 | ile=iendst(line)
|
---|
[bd2278d] | 148 | ! _________________________________ ! ignore empty and commentary lines
|
---|
[e40e335] | 149 | if (ile.gt.0.and.index(line(1:ile),cmt).le.0) then
|
---|
| 150 |
|
---|
[bd2278d] | 151 | ! _________________________________________ Global variables
|
---|
[e40e335] | 152 | ilb = index(line(1:ile),sglp)+1
|
---|
| 153 | if (ilb.ge.2) then
|
---|
| 154 |
|
---|
| 155 | if (index(line(ilb:ile),spfd).gt.0) then ! field with mol.#
|
---|
| 156 |
|
---|
| 157 | call extstr(spfd,ilb,ile,line,lincmd,l)
|
---|
| 158 |
|
---|
| 159 | if (iredin(lincmd,nml).le.0.or.
|
---|
[bd2278d] | 160 | & nml.le.0.or.nml.gt.ntlml) then
|
---|
[38d77eb] | 161 | write (logString, *)
|
---|
| 162 | & 'redvar> ','Incorrect molecule number >',
|
---|
[bd2278d] | 163 | & lincmd(1:l),'< Must be in range [1,',
|
---|
| 164 | & ntlml,'] !'
|
---|
[e40e335] | 165 | close(lunvar)
|
---|
| 166 | stop
|
---|
| 167 | endif
|
---|
| 168 |
|
---|
| 169 | else
|
---|
| 170 | nml = 1 ! assume mol. #1
|
---|
| 171 | endif
|
---|
| 172 |
|
---|
| 173 | l=ile-ilb+1
|
---|
| 174 | if (l.le.0) goto 105
|
---|
| 175 | lincmd=blnk
|
---|
| 176 | lincmd(1:l)=line(ilb:ile)
|
---|
| 177 |
|
---|
| 178 | k = 1
|
---|
| 179 | do i = 1,5 ! try to read 5 parameters
|
---|
| 180 | call extstr(spcc,k,l,lincmd,linh,n)
|
---|
| 181 | if (k.gt.l.or.iredrl(linh,rn(i)).le.0) goto 105
|
---|
| 182 | enddo
|
---|
| 183 |
|
---|
| 184 | n=l-k+1 ! try 6th parameter
|
---|
| 185 | if (n.le.0) goto 105
|
---|
| 186 | linh=blnk
|
---|
| 187 | linh(1:n)=lincmd(k:l)
|
---|
| 188 | if (iredrl(linh,rn(6)).le.0) goto 105
|
---|
| 189 |
|
---|
[bd2278d] | 190 | ! ---------------------------------------- check global angles
|
---|
[e40e335] | 191 | if ( abs(rn(4)).gt.(1.8d2+1d-6)
|
---|
[bd2278d] | 192 | & .or. abs(rn(5)).gt.(9d1+1d-6)
|
---|
| 193 | & .or. abs(rn(6)).gt.(1.8d2+1d-6)
|
---|
| 194 | & ) goto 106
|
---|
[e40e335] | 195 |
|
---|
| 196 | do i = 1,3
|
---|
| 197 | gbpr(i,nml) = rn(i)
|
---|
| 198 | enddo
|
---|
| 199 | do i = 4,6
|
---|
| 200 | gbpr(i,nml) = rn(i)*cdr
|
---|
| 201 | enddo
|
---|
| 202 |
|
---|
| 203 | goto 1
|
---|
| 204 |
|
---|
| 205 | endif ! global vars
|
---|
| 206 |
|
---|
| 207 |
|
---|
| 208 | ilb=1
|
---|
| 209 |
|
---|
| 210 | do while (ilb.le.ile) ! ________________________ Commands
|
---|
| 211 | call extstr(spcm,ilb,ile,line,lincmd,ice)
|
---|
| 212 |
|
---|
| 213 | if (ice.gt.0) then ! ignore empty commands
|
---|
| 214 | ncmd=ncmd+1
|
---|
| 215 | if (ncmd.gt.maxcmd) goto 101
|
---|
| 216 |
|
---|
| 217 | ifx=index(lincmd(1:ice),sfix)
|
---|
| 218 | if (ifx.gt.0) then ! check for commands to fix variables
|
---|
| 219 | fix=.true.
|
---|
| 220 | lincmd(ifx:ifx)=blnk
|
---|
| 221 | if (ifx.eq.ice) then
|
---|
| 222 | ice=ice-1
|
---|
| 223 | if (ice.eq.0) then ! fix all
|
---|
| 224 | ice=1
|
---|
| 225 | lincmd(1:1)=wdc
|
---|
| 226 | endif
|
---|
| 227 | endif
|
---|
| 228 | else
|
---|
| 229 | fix=.false.
|
---|
| 230 | endif
|
---|
| 231 |
|
---|
[bd2278d] | 232 | ! _________________________________________ Extract Command Fields
|
---|
[e40e335] | 233 | nfld=0
|
---|
| 234 | icb=1
|
---|
| 235 | do while (icb.le.ice)
|
---|
| 236 | nfld=nfld+1
|
---|
| 237 | if (nfld.gt.maxfld) goto 100
|
---|
| 238 | call extstr(spfd,icb,ice,lincmd,linfld(nfld),ifdend(nfld))
|
---|
| 239 |
|
---|
| 240 | if (ifdend(nfld).le.0) then ! empty field means 'all'
|
---|
| 241 | linfld(nfld)(1:1)=wdc
|
---|
| 242 | ifdend(nfld)=1
|
---|
| 243 | endif
|
---|
| 244 |
|
---|
| 245 | enddo
|
---|
[bd2278d] | 246 | ! _______________________________ Interpret Command Fields (except last)
|
---|
[e40e335] | 247 | do i=1,nfld-1
|
---|
| 248 | ii=i
|
---|
| 249 | ifld=nfld-i
|
---|
| 250 |
|
---|
| 251 | if (ifld.eq.3) then ! Initialize Molecules
|
---|
| 252 | do j=1,ntlml
|
---|
| 253 | forml(j)=.false.
|
---|
| 254 | enddo
|
---|
| 255 | elseif (ifld.eq.2) then ! Initialize Residues
|
---|
| 256 | do j=1,ntlml
|
---|
| 257 | do k=irsml1(j),irsml2(j)
|
---|
| 258 | forrs(k)=.false.
|
---|
| 259 | enddo
|
---|
| 260 | enddo
|
---|
| 261 | elseif (ifld.eq.1) then ! Initialize Variables
|
---|
| 262 | do j=1,ntlvr
|
---|
| 263 | forvr(j)=.false.
|
---|
| 264 | enddo
|
---|
| 265 | endif
|
---|
[bd2278d] | 266 | ! __________________________________ Identifiers in field
|
---|
[e40e335] | 267 | nide=0
|
---|
| 268 | ifb=1
|
---|
| 269 | ife=ifdend(i)
|
---|
| 270 | do while (ifb.le.ife)
|
---|
| 271 | nide=nide+1
|
---|
| 272 | if (nide.gt.maxide) goto 103
|
---|
| 273 | call extstr(spcc,ifb,ife,linfld(ii),linide,ide)
|
---|
| 274 | if (ide.le.0.or.linide(1:1).eq.wdc) then ! ... All
|
---|
| 275 | if (ifld.eq.3) then ! Mol.
|
---|
| 276 | do j=1,ntlml
|
---|
| 277 | forml(j)=.true.
|
---|
| 278 | enddo
|
---|
| 279 | elseif (ifld.eq.2) then ! Res.
|
---|
| 280 | do j=1,ntlml
|
---|
| 281 | if (forml(j)) then
|
---|
| 282 | do k=irsml1(j),irsml2(j)
|
---|
| 283 | forrs(k)=.true.
|
---|
| 284 | enddo
|
---|
| 285 | endif
|
---|
| 286 | enddo
|
---|
| 287 | elseif (ifld.eq.1) then ! Var.
|
---|
| 288 | do j=1,ntlml
|
---|
| 289 | if (forml(j)) then
|
---|
| 290 | do k=irsml1(j),irsml2(j)
|
---|
| 291 | if (forrs(k)) then
|
---|
| 292 | ll=ivrrs1(k)
|
---|
| 293 | do l=ll,ll+nvrrs(k)-1
|
---|
| 294 | forvr(l)=.true.
|
---|
| 295 | enddo
|
---|
| 296 | endif
|
---|
| 297 | enddo
|
---|
| 298 | endif
|
---|
| 299 | enddo
|
---|
| 300 | endif
|
---|
| 301 |
|
---|
| 302 | else ! ...................... Identifier .ne. wdc
|
---|
| 303 |
|
---|
| 304 | ihy=index(linide(1:ide),sphy) ! ? zone of numbers
|
---|
| 305 |
|
---|
| 306 | if (ihy.le.0) then ! _____ No zone
|
---|
| 307 | if (iredin(linide,inum).gt.0) then ! ... number
|
---|
| 308 | if (ifld.eq.3) then ! Mol.
|
---|
| 309 |
|
---|
[bd2278d] | 310 | ! ################### impossible # (inum) of molecule
|
---|
[e40e335] | 311 |
|
---|
| 312 | if (inum.le.0.or.inum.gt.ntlml) then
|
---|
[38d77eb] | 313 | write (logString, *) ' # 1: ',inum
|
---|
[e40e335] | 314 | goto 104
|
---|
| 315 | endif
|
---|
| 316 |
|
---|
| 317 | forml(inum)=.true.
|
---|
| 318 | elseif (ifld.eq.2) then ! Res.
|
---|
| 319 | do j=1,ntlml
|
---|
| 320 | if (forml(j)) then
|
---|
| 321 | nfi=irsml1(j)
|
---|
| 322 | k=inum+nfi-1
|
---|
| 323 |
|
---|
[bd2278d] | 324 | ! ################### impossible # of residue (inum) in molecule
|
---|
[e40e335] | 325 |
|
---|
| 326 | if (k.lt.nfi.or.k.gt.irsml2(j)) then
|
---|
[38d77eb] | 327 | write (logString, *) ' # 2: ',inum
|
---|
[e40e335] | 328 | goto 104
|
---|
| 329 | endif
|
---|
| 330 |
|
---|
| 331 | forrs(k)=.true.
|
---|
| 332 | endif
|
---|
| 333 | enddo
|
---|
| 334 | elseif (ifld.eq.1) then ! Var.
|
---|
| 335 | do j=1,ntlml
|
---|
| 336 | if (forml(j)) then
|
---|
| 337 | do k=irsml1(j),irsml2(j)
|
---|
| 338 | if (forrs(k)) then
|
---|
| 339 | nfi=ivrrs1(k)
|
---|
| 340 | l=inum+nfi-1
|
---|
| 341 |
|
---|
[bd2278d] | 342 | ! ################### impossible # of variable (inum) in residue
|
---|
[e40e335] | 343 |
|
---|
| 344 | if (l.lt.nfi.or.
|
---|
[bd2278d] | 345 | & l.gt.nfi+nvrrs(k)-1) then
|
---|
[38d77eb] | 346 | write (logString, *) ' # 3: ',inum
|
---|
[e40e335] | 347 | goto 104
|
---|
| 348 | endif
|
---|
| 349 |
|
---|
| 350 | forvr(l)=.true.
|
---|
| 351 | endif
|
---|
| 352 | enddo
|
---|
| 353 | endif
|
---|
| 354 | enddo
|
---|
| 355 | endif
|
---|
| 356 |
|
---|
| 357 | else ! ... Name
|
---|
| 358 | if (linide(ide:ide).eq.wdc) then
|
---|
| 359 | id=ide-1
|
---|
| 360 | exa=.false.
|
---|
| 361 | else ! exact match of names
|
---|
| 362 | id=ide
|
---|
| 363 | exa=.true.
|
---|
| 364 | endif
|
---|
| 365 |
|
---|
| 366 | if (ifld.eq.3) then ! Mol.
|
---|
| 367 | do j=1,ntlml
|
---|
| 368 | ib=ibegst(nmml(j))
|
---|
| 369 | if (ib.gt.0) then
|
---|
| 370 | linh=blnk
|
---|
| 371 | ieh=iendst(nmml(j))
|
---|
| 372 | ieh1=ieh-ib+1
|
---|
| 373 | linh(1:ieh1)=nmml(j)(ib:ieh)
|
---|
| 374 | if (((exa.and.ieh1.eq.id).or.
|
---|
[bd2278d] | 375 | & (.not.exa.and.ieh1.ge.id)).and.
|
---|
| 376 | & linh(1:id).eq.linide(1:id))
|
---|
| 377 | & forml(j)=.true.
|
---|
[e40e335] | 378 | endif
|
---|
| 379 | enddo
|
---|
| 380 | elseif (ifld.eq.2) then ! Res.
|
---|
| 381 | do j=1,ntlml
|
---|
| 382 | if (forml(j)) then
|
---|
| 383 | do k=irsml1(j),irsml2(j)
|
---|
| 384 | ib=ibegst(seq(k))
|
---|
| 385 | if (ib.gt.0) then
|
---|
| 386 | linh=blnk
|
---|
| 387 | ieh=iendst(seq(k))
|
---|
| 388 | ieh1=ieh-ib+1
|
---|
| 389 | linh(1:ieh1)=seq(k)(ib:ieh)
|
---|
| 390 | if (((exa.and.ieh1.eq.id).or.
|
---|
[bd2278d] | 391 | & (.not.exa.and.ieh1.ge.id))
|
---|
| 392 | & .and.linh(1:id).eq.linide(1:id))
|
---|
| 393 | & forrs(k)=.true.
|
---|
[e40e335] | 394 | endif
|
---|
| 395 | enddo
|
---|
| 396 | endif
|
---|
| 397 | enddo
|
---|
| 398 | elseif (ifld.eq.1) then ! Var.
|
---|
| 399 | do j=1,ntlml
|
---|
| 400 | if (forml(j)) then
|
---|
| 401 | do k=irsml1(j),irsml2(j)
|
---|
| 402 | if (forrs(k)) then
|
---|
| 403 | ll=ivrrs1(k)
|
---|
| 404 | do l=ll,ll+nvrrs(k)-1
|
---|
| 405 | ib=ibegst(nmvr(l))
|
---|
| 406 | if (ib.gt.0) then
|
---|
| 407 | linh=blnk
|
---|
| 408 | ieh=iendst(nmvr(l))
|
---|
| 409 | ieh1=ieh-ib+1
|
---|
| 410 | linh(1:ieh1)=nmvr(l)(ib:ieh)
|
---|
| 411 | if (((exa.and.ieh1.eq.id)
|
---|
[bd2278d] | 412 | & .or.(.not.exa.and.ieh1.ge.id))
|
---|
| 413 | & .and.linh(1:id).eq.linide(1:id))
|
---|
| 414 | & forvr(l)=.true.
|
---|
[e40e335] | 415 | endif
|
---|
| 416 | enddo
|
---|
| 417 | endif
|
---|
| 418 | enddo
|
---|
| 419 | endif
|
---|
| 420 | enddo
|
---|
| 421 | endif
|
---|
| 422 |
|
---|
| 423 | endif
|
---|
| 424 |
|
---|
| 425 | else ! ___ Zone
|
---|
| 426 |
|
---|
[bd2278d] | 427 | ! ################### impossible zone '-' (without integer)
|
---|
[e40e335] | 428 |
|
---|
| 429 | if (ide.eq.1.and.ihy.eq.ide) then
|
---|
[38d77eb] | 430 | write (logString, *) ' # 4: ',ide
|
---|
[e40e335] | 431 | goto 104
|
---|
| 432 | endif
|
---|
| 433 |
|
---|
| 434 | if (ihy.eq.1) then
|
---|
| 435 | ibz=1
|
---|
| 436 | else
|
---|
| 437 | linh=blnk
|
---|
| 438 | linh=linide(1:ihy-1)
|
---|
| 439 |
|
---|
[bd2278d] | 440 | ! ################### impossible (to read) integer before '-'
|
---|
[e40e335] | 441 |
|
---|
| 442 | if (iredin(linh,ibz).le.0.or.ibz.le.0)
|
---|
[bd2278d] | 443 | & then
|
---|
[38d77eb] | 444 | write (logString, *) ' # 5 '
|
---|
[e40e335] | 445 | goto 104
|
---|
| 446 | endif
|
---|
| 447 |
|
---|
| 448 | endif
|
---|
| 449 | if (ihy.eq.ide) then
|
---|
| 450 | iez=ilrg
|
---|
| 451 | else
|
---|
| 452 | linh=blnk
|
---|
| 453 | linh=linide(ihy+1:ide)
|
---|
| 454 |
|
---|
[bd2278d] | 455 | ! ################### impossible (to read) integer after '-'
|
---|
[e40e335] | 456 |
|
---|
| 457 | if (iredin(linh,iez).le.0.or.iez.le.0.or.
|
---|
[bd2278d] | 458 | & iez.lt.ibz) then
|
---|
[38d77eb] | 459 | write (logString, *) ' # 6 '
|
---|
[e40e335] | 460 | goto 104
|
---|
| 461 | endif
|
---|
| 462 |
|
---|
| 463 | endif
|
---|
| 464 |
|
---|
| 465 | if (ifld.eq.3) then ! Mol.
|
---|
| 466 | if (iez.gt.ntlml) iez=ntlml
|
---|
| 467 | do j=ibz,iez
|
---|
| 468 | forml(j)=.true.
|
---|
| 469 | enddo
|
---|
| 470 | elseif (ifld.eq.2) then ! Res.
|
---|
| 471 | do j=1,ntlml
|
---|
| 472 | if (forml(j)) then
|
---|
| 473 | kbz=irsml1(j)+ibz-1
|
---|
| 474 | kez=irsml1(j)+iez-1
|
---|
| 475 | if (kez.gt.irsml2(j)) then
|
---|
| 476 | kk=irsml2(j)
|
---|
| 477 | else
|
---|
| 478 | kk=kez
|
---|
| 479 | endif
|
---|
| 480 | do k=kbz,kk
|
---|
| 481 | forrs(k)=.true.
|
---|
| 482 | enddo
|
---|
| 483 | endif
|
---|
| 484 | enddo
|
---|
| 485 | elseif (ifld.eq.1) then ! Var.
|
---|
| 486 | do j=1,ntlml
|
---|
| 487 | if (forml(j)) then
|
---|
| 488 | do k=irsml1(j),irsml2(j)
|
---|
| 489 | kv=nvrrs(k)
|
---|
| 490 | if (forrs(k).and.kv.gt.0) then
|
---|
| 491 | ll=ivrrs1(k)
|
---|
| 492 | lbz=ll+ibz-1
|
---|
| 493 | if (iez.gt.kv) then
|
---|
| 494 | lez=ll+kv-1
|
---|
| 495 | else
|
---|
| 496 | lez=ll+iez-1
|
---|
| 497 | endif
|
---|
| 498 | do l=lbz,lez
|
---|
| 499 | forvr(l)=.true.
|
---|
| 500 | enddo
|
---|
| 501 | endif
|
---|
| 502 | enddo
|
---|
| 503 | endif
|
---|
| 504 | enddo
|
---|
| 505 | endif
|
---|
| 506 |
|
---|
| 507 | endif
|
---|
| 508 | endif
|
---|
| 509 |
|
---|
| 510 | enddo ! ... identifiers
|
---|
| 511 | enddo ! ... Fields (excl. value)
|
---|
| 512 |
|
---|
[bd2278d] | 513 | ! _____________________________________________________ Execute Command
|
---|
[e40e335] | 514 |
|
---|
| 515 | if (iredrl(linfld(nfld),val).gt.izero) then ! Read Value
|
---|
| 516 | did=.false.
|
---|
| 517 | do i=1,ntlvr
|
---|
| 518 | if (forvr(i)) then
|
---|
| 519 | did=.true.
|
---|
| 520 | vlvrx(i)=val
|
---|
| 521 |
|
---|
| 522 | fxvr(i)=fix
|
---|
| 523 |
|
---|
| 524 | stvr(i)=.true.
|
---|
| 525 | endif
|
---|
| 526 | enddo
|
---|
[38d77eb] | 527 | if (.not.did) write (logString, '(3a)')
|
---|
[bd2278d] | 528 | & ' redvar> No variables affected by command >',
|
---|
| 529 | & lincmd(1:ice),'<'
|
---|
[e40e335] | 530 | else
|
---|
| 531 |
|
---|
| 532 | ll1=ibegst(linfld(nfld))
|
---|
| 533 | ll2=iendst(linfld(nfld))
|
---|
[38d77eb] | 534 | write (logString, *) 'll1,ll2, linfld(nfld): ',ll1,ll2,
|
---|
[bd2278d] | 535 | & '>',linfld(nfld)(ll1:ll2),'<'
|
---|
[e40e335] | 536 |
|
---|
| 537 | goto 102
|
---|
| 538 | endif
|
---|
| 539 |
|
---|
| 540 | endif
|
---|
| 541 | enddo ! ... Commands at one line
|
---|
| 542 | endif
|
---|
| 543 | goto 1
|
---|
| 544 |
|
---|
| 545 | 2 close(lunvar)
|
---|
[bd2278d] | 546 | ! __________________________ Summary
|
---|
[e40e335] | 547 | iv=0
|
---|
| 548 | do i=1,ntlml
|
---|
| 549 |
|
---|
| 550 | ie=iendst(nmml(i))
|
---|
| 551 |
|
---|
| 552 | do j =1,6
|
---|
| 553 | if (gbpr(j,i).ne.zero) then
|
---|
| 554 |
|
---|
| 555 | do k = 1,3
|
---|
| 556 | write(strg(k),'(f17.6)') gbpr(k,i)
|
---|
| 557 | enddo
|
---|
| 558 | do k = 4,6
|
---|
| 559 | write(strg(k),'(f17.6)') (gbpr(k,i)*crd)
|
---|
| 560 | enddo
|
---|
| 561 |
|
---|
[38d77eb] | 562 | write (logString, '(3a,/,1x,5(a,2x),a)') ' redvar> ',
|
---|
| 563 | & nmml(i)(1:ie),
|
---|
[bd2278d] | 564 | & ' with global parameters:',
|
---|
| 565 | & (strg(k)(ibegst(strg(k)):),k=1,6)
|
---|
[e40e335] | 566 | call setvar(i,vlvrx)
|
---|
| 567 | goto 3
|
---|
| 568 | endif
|
---|
| 569 | enddo
|
---|
| 570 |
|
---|
| 571 | 3 if (nvrml(i).gt.0) then
|
---|
| 572 | iml=i
|
---|
| 573 | did=.false.
|
---|
| 574 | in=0
|
---|
| 575 | jb=irsml1(i)-1
|
---|
| 576 | do j=irsml1(i),irsml2(i)
|
---|
| 577 | kk=ivrrs1(j)
|
---|
| 578 | do k=kk,kk+nvrrs(j)-1
|
---|
| 579 | iv=iv+1
|
---|
| 580 | if (stvr(iv)) then
|
---|
| 581 | did=.true.
|
---|
| 582 | if (fxvr(iv)) then
|
---|
[38d77eb] | 583 | write(logString,'(3a,i4,1x,4a,f10.3,a)') ' redvar> ',
|
---|
[bd2278d] | 584 | & nmml(i)(1:ie),': residue ',j-jb,seq(j),
|
---|
| 585 | & ': ',nmvr(iv),' set ',vlvrx(iv),' Fixed'
|
---|
[e40e335] | 586 | else
|
---|
[38d77eb] | 587 | write (logString, '(3a,i4,1x,4a,f10.3)') ' redvar> ',
|
---|
[bd2278d] | 588 | & nmml(i)(1:ie),': residue ',j-jb,seq(j),
|
---|
| 589 | & ': ',nmvr(iv),' set ',vlvrx(iv)
|
---|
[e40e335] | 590 | endif
|
---|
| 591 | ity=ityvr(iv)
|
---|
| 592 | if (ity.eq.3.or.ity.eq.2)
|
---|
[bd2278d] | 593 | & vlvrx(iv)=vlvrx(iv)*cdr ! angles
|
---|
[e40e335] | 594 |
|
---|
| 595 | else
|
---|
| 596 | in=in+1
|
---|
| 597 | endif
|
---|
| 598 | enddo
|
---|
| 599 | enddo
|
---|
| 600 | if (did) then
|
---|
[38d77eb] | 601 | if (in.gt.0) write (logString, '(3a,i5,a)')
|
---|
[bd2278d] | 602 | & ' redvar> Molecule ',nmml(i)(1:ie),': ',in,
|
---|
| 603 | & ' variable(s) remain unchanged'
|
---|
[e40e335] | 604 | call setvar(iml,vlvrx)
|
---|
| 605 | else
|
---|
[38d77eb] | 606 | write (logString, '(3a)') ' redvar> Molecule ',
|
---|
[bd2278d] | 607 | & nmml(i)(1:ie),': No internal variables changed'
|
---|
[e40e335] | 608 | endif
|
---|
| 609 | endif
|
---|
| 610 | enddo
|
---|
| 611 |
|
---|
| 612 | return
|
---|
[bd2278d] | 613 | ! ____________________________________________________________ Errors
|
---|
[38d77eb] | 614 | 100 write (logString, '(3a)') ' redvar> Cannot interpret command >',
|
---|
[bd2278d] | 615 | & lincmd(1:ice),'<'
|
---|
[e40e335] | 616 | close(lunvar)
|
---|
| 617 | stop
|
---|
[38d77eb] | 618 | 101 write (logString, '(a,i5,a)') ' redvar> Command number ',ncmd,
|
---|
| 619 | & ' reached'
|
---|
[e40e335] | 620 | close(lunvar)
|
---|
| 621 | stop
|
---|
[38d77eb] | 622 | 102 write (logString, '(3a)') ' redvar> Cannot read value from >',
|
---|
[bd2278d] | 623 | & lincmd(1:ice),'<'
|
---|
[e40e335] | 624 | close(lunvar)
|
---|
| 625 | stop
|
---|
[38d77eb] | 626 | 103 write (logString, '(a,i3,3a)') ' redvar> Cannot read >',maxide,
|
---|
[bd2278d] | 627 | & ' identifiers from >',linfld(ii)(1:ife),'<'
|
---|
[e40e335] | 628 | close(lunvar)
|
---|
| 629 | stop
|
---|
[38d77eb] | 630 | 104 write (logString, '(5a)') ' redvar> Error in identifier >',
|
---|
[bd2278d] | 631 | & linide(1:ide),'< of command >',lincmd(1:ice),'<'
|
---|
[e40e335] | 632 | close(lunvar)
|
---|
| 633 | stop
|
---|
[38d77eb] | 634 | 105 write (logString, '(a,/,a,/,2a,/)')
|
---|
| 635 | & ' redvar> line with global paramters:',
|
---|
[bd2278d] | 636 | & line(1:ile),' must contain 6 floating',
|
---|
| 637 | & ' point numbers separated by commas !'
|
---|
[e40e335] | 638 | close(lunvar)
|
---|
| 639 | stop
|
---|
| 640 |
|
---|
[38d77eb] | 641 | 106 write (logString, '(a,/,a,/,2a,/)')
|
---|
| 642 | & ' redvar> line with global paramters:',
|
---|
[bd2278d] | 643 | & line(1:ile),' angles must be inside ',
|
---|
| 644 | &'ranges [-180,180], [-90,90], and [-180,180] Deg., respectively !'
|
---|
[e40e335] | 645 | close(lunvar)
|
---|
| 646 | stop
|
---|
| 647 |
|
---|
| 648 | end
|
---|