Previous Thread
Next Thread
Print Thread
#36780 01/22/18 06:54 PM
Joined: Sep 2003
Posts: 8,606
Likes: 24
rmv Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 8,606
Likes: 24
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 Images
watchk.inp.txt (1.03 KB, 562 downloads)
BldWatEnt.csh.txt (358 Bytes, 549 downloads)
wtrans.inp.txt (663 Bytes, 1485 downloads)
wtransits.f90.txt (7.11 KB, 737 downloads)

Rick Venable
computational chemist

Joined: Aug 2009
Posts: 134
L
lqz Offline
Forum Member
Offline
Forum Member
L
Joined: Aug 2009
Posts: 134
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.

Joined: Sep 2003
Posts: 8,606
Likes: 24
rmv Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 8,606
Likes: 24
There's nothing in the BldWatEnt.csh script that limits the number of waters.


Rick Venable
computational chemist

Joined: Aug 2009
Posts: 134
L
lqz Offline
Forum Member
Offline
Forum Member
L
Joined: Aug 2009
Posts: 134
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.

Joined: Sep 2003
Posts: 8,606
Likes: 24
rmv Online Content OP
Forum Member
OP Online Content
Forum Member
Joined: Sep 2003
Posts: 8,606
Likes: 24
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.5
(Release build 20201027)
Responsive Width:

PHP: 7.3.31-1~deb10u1 Page Time: 0.005s Queries: 25 (0.002s) Memory: 0.7465 MB (Peak: 0.8178 MB) Data Comp: Off Server Time: 2022-10-06 07:57:13 UTC
Valid HTML 5 and Valid CSS