Previous Thread
Next Thread
Print Thread
checking for water transits
#36780 01/22/18 06:54 PM
Joined: Sep 2003
Posts: 8,470
rmv Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 8,470
The attached files are from a workflow I used with CHARMM to count water transits through a bilayer. It assumes that any drift in z has been removed by some alignment procedure that ensures the bilayer is intact and centered at a known position. Briefly,
  • watchk.inp :: identifies waters that visit the bilayer midplane region, creates a log file
  • BldWatEnt.csh :: processes the log file to get a unique list of RESIDs, and creates a stream file of ENTER commands for use with CORREL
  • wtrans.inp :: creates time series of z coord for identified waters using the stream file from BldWatEnt.csh, one file per water in a subdir (env var $ZWAT)
  • wtransits.f90 :: Fortran90 program that analyzes the water z time series produced, and reports on true transits

The idea is to greatly reduce the number of water molecules to consider explicitly by determining which ones actually visit the hydrophobic midplane region in the initial screening (watchk.inp).

The Fortran program expects command line args, and gives a brief help message when none are provided.

(The .txt extensions should be removed before use; make sure the .csh file has execute permission.)



Attached Files
watchk.inp.txt (1.03 KB, 326 downloads)
BldWatEnt.csh.txt (358 Bytes, 301 downloads)
wtrans.inp.txt (663 Bytes, 765 downloads)
wtransits.f90.txt (7.11 KB, 425 downloads)

Rick Venable
computational chemist

Re: checking for water transits
rmv #37182 10/25/18 02:10 AM
Joined: Aug 2009
Posts: 132
L
lqz Offline
Forum Member
Offline
Forum Member
L
Joined: Aug 2009
Posts: 132
The BldWatEnt.csh program could not handle more than 500 water molecules in one run. Is there a way to increase its size limit? Thanks.

Re: checking for water transits
rmv #37185 10/25/18 06:38 PM
Joined: Sep 2003
Posts: 8,470
rmv Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 8,470
There's nothing in the BldWatEnt.csh script that limits the number of waters.


Rick Venable
computational chemist

Re: checking for water transits
rmv #37186 10/25/18 09:20 PM
Joined: Aug 2009
Posts: 132
L
lqz Offline
Forum Member
Offline
Forum Member
L
Joined: Aug 2009
Posts: 132
Hi, Rick, when I run this program, it gave me error message as "Found 1201 unique water RESIDs, foreach: Too many arguments". When I cut the trajectory into two halves, and ran each piece one after another, it worked. That is why I suspected some kind of size limitation involved.

Re: checking for water transits
rmv #37187 10/25/18 09:59 PM
Joined: Sep 2003
Posts: 8,470
rmv Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 8,470
It's possible the 'foreach' command has a limit of 1200 for the number of items inside the parentheses; I hadn't run into that.


Rick Venable
computational chemist


Moderated by  lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.4
(Release build 20200307)
Responsive Width:

PHP: 5.6.33-0+deb8u1 Page Time: 0.014s Queries: 25 (0.007s) Memory: 0.9227 MB (Peak: 1.0163 MB) Data Comp: Off Server Time: 2020-07-12 18:46:58 UTC
Valid HTML 5 and Valid CSS