Here you go!
function loadCosmo(fname) {
print "reading " + fname
Var fdata = load(fname).split(" DATA")
Var atoms = fdata[2].lines[3][-2]
Var segments = fdata[3].lines[3][0]
Var n = atoms.length + segments.length
fdata = "" + n + "\ndata are from " + fname + "\n"
print "nAtoms = " + atoms.length + " nSegments = " + segments.length
for (Var i = 1; i <=atoms.length; i++) {
Var radius = (atoms[i])[52][63]
fdata += (atoms[i])[6][51] + " - " + radius + "\n"
}
for (Var i = 1; i <=segments.length; i++) {
Var sigma = (segments[i])[77][88]
fdata += "Xx " + (segments[i])[17][52] + sigma + "\n"
}
load inline @fdata
isosurface select {not _XX} ignore {_XX} solvent 1.2 map select {_XX}
ignore {not _XX} property partialcharge colorscheme "rwb"
delete _Xx
}
loadCosmo "xxx.cos"
------------------------------------------------------------------------------
The ultimate all-in-one performance toolkit: Intel(R) Parallel Studio XE:
Pinpoint memory and threading errors before they happen.
Find and fix more than 250 security defects in the development cycle.
Locate bottlenecks in serial and parallel code that limit performance.
http://p.sf.net/sfu/intel-dev2devfeb
_______________________________________________
Jmol-users mailing list
Jmol-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/jmol-users
