From: Cosseddu, Salvatore (S.M.Cosseddu_at_warwick.ac.uk)
Date: Wed Mar 12 2014 - 07:05:43 CDT
Dear Norman and Thomas,
It works also within tclforces. You must be ready to a slight loss of performance for the additional I/O work. Anyway the following procedure has been tested and have worked fine for me since version 2.8. Hope it will help
------------------------------------
proc readcell {xscfile} {
if {[catch {open "${xscfile}" r} xsc]} {
puts stderr "ERROR opening ${xscfile}
CHECK!!!"
exit 1
}
set in [read $xsc]
set bcell [list [lindex $in 28] [lindex $in 32] [lindex $in 36]]
print "readcell: read cell, at step [lindex $in 27] was $bcell"
close $xsc
return $bcell
}
-------------------------------------
simply usage:
set boxcell [readcell <xsc filename>]
Best regards
Salvatore M Cosseddu
________________________________
From: owner-namd-l_at_ks.uiuc.edu <owner-namd-l_at_ks.uiuc.edu> on behalf of Norman Geist <norman.geist_at_uni-greifswald.de>
Sent: 12 March 2014 07:29
To: 'Thomas Evangelidis'
Cc: Namd Mailing List
Subject: AW: namd-l: how to get cell vectors in tcl forces script
Hi,
I was also searching for such a command, but the only way you are currently able to do it, seems to be using the "output" command to write a set of restarts and read the box info out of the xsc file. No idea if this is also possible within TclForces, but its working from jobscript context. Also, just to mention it, the only way currently to "set" the box during simulation is by loading dcd frames with the "coorfile" command.
Norman Geist.
Von: owner-namd-l_at_ks.uiuc.edu [mailto:owner-namd-l_at_ks.uiuc.edu] Im Auftrag von Thomas Evangelidis
Gesendet: Dienstag, 11. M?rz 2014 21:14
An: namd-l
Betreff: namd-l: how to get cell vectors in tcl forces script
Dear NAMD community,
Is it possible to get the current unit cell vectors on the fly from within a tcl forces script?
-- ====================================================================== Thomas Evangelidis PhD student University of Athens Faculty of Pharmacy Department of Pharmaceutical Chemistry Panepistimioupoli-Zografou 157 71 Athens GREECE email: tevang_at_pharm.uoa.gr<mailto:tevang_at_pharm.uoa.gr> tevang3_at_gmail.com<mailto:tevang3_at_gmail.com> website: https://sites.google.com/site/thomasevangelidishomepage/ ________________________________ [http://static.avast.com/emails/avast-mail-stamp.png] <http://www.avast.com/> Diese E-Mail ist frei von Viren und Malware, denn der avast! Antivirus<http://www.avast.com/> Schutz ist aktiv.
This archive was generated by hypermail 2.1.6 : Thu Dec 31 2015 - 23:20:36 CST