/usr/share/xcrysden/Tcl/forces.tcl is in xcrysden-data 1.5.60-1build3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 | #############################################################################
# Author: #
# ------ #
# Anton Kokalj Email: Tone.Kokalj@ijs.si #
# Department of Physical and Organic Chemistry Phone: x 386 1 477 3523 #
# Jozef Stefan Institute Fax: x 386 1 477 3811 #
# Jamova 39, SI-1000 Ljubljana #
# SLOVENIA #
# #
# Source: $XCRYSDEN_TOPDIR/Tcl/forces.tcl
# ------ #
# Copyright (c) 1996-2003 by Anton Kokalj #
#############################################################################
proc forceVectors can {
global check sInfo
if { ! [info exists sInfo(lforce)] } {
return
}
if { ! $sInfo(lforce) } {
return
}
if { $check(forces) } {
xc_forces $can on
$can render
xcSwapBuffers
} else {
xc_forces $can off
$can render
xcSwapBuffers
}
}
proc forceVectorsSet can {
global check forceVec sInfo
if { ! [info exists sInfo(lforce)] } {
return
}
set c [lindex [split $can .] end]
if { [winfo exists .force$c] } {
return
}
if { ![info exists forceVec(scalefunction)] } {
global mody
set forceVec(scalefunction) linear
set forceVec(threshold) 0.0005
set forceVec(lengthfactor) 200
set forceVec(rod_thickf) [xc_getdefault $mody(D_FORCE_RODTHICKF)]
set forceVec(arr_thickf) [xc_getdefault $mody(D_FORCE_ARRTHICKF)]
set forceVec(arr_lenf) [xc_getdefault $mody(D_FORCE_ARRLENF)]
set forceVec(color) [xc_getdefault $mody(D_FORCE_COLOR)]
puts stderr "forceVec: color == $forceVec(color)"
}
# from here on: widget managing
set t [xcToplevel .force$c "Forces: Settings" "Forces" . -0 0 1]
set f1 [frame $t.f1 -relief ridge -bd 2]
set f2 [frame $t.f2 -relief ridge -bd 2]
set f3 [frame $t.f3]
pack $f1 $f2 $f3 -side top -expand 1 -fill x -padx 3m -pady 3m -ipadx 3m -ipady 3m
pack $f3 -side top -expand 1
#
# FRAME-1: Scale-Function + Threshold + Length-factor
#
set mb [xcMenuButton $f1 \
-labeltext "Scale Function:" \
-labelwidth 15 \
-textvariable forceVec(scalefunction) \
-menu { \
{Linear} {set forceVec(scalefunction) linear} \
{Natural Logarithm} {set forceVec(scalefunction) log} \
{Decadic Logaritm} {set forceVec(scalefunction) log10} \
{Square Root} {set forceVec(scalefunction) sqrt} \
{Cubic Root} {set forceVec(scalefunction) root3} \
{Exponential} {set forceVec(scalefunction) exp} \
{exp(x*x)} {set forceVec(scalefunction) exp2} \
}]
set m1 [message $f1.m1 -justify left -aspect 800 \
-text "Threshold::\ndo not show the forces whose magnitude is smaller then the specified threshold.\n\n\nLength Factor::\nassumed unit for force is Hartree/Angstrom\nforceVector_length == forceVectorSize * LengthFactor"]
pack $mb $m1 -side top -fill x -anchor w -expand 1
FillEntries $f1 {Treshold {Length Factor}} \
{forceVec(threshold) forceVec(lengthfactor)} 13 10 top left
#
# FRAME-2: attributes + Threshold + Length-factor
#
FillEntries $f2 {
"Vectors thickness factor:"
"Thickness factor for arrow-cap:"
"Length factor for arrow-cap:"
} {forceVec(rod_thickf) forceVec(arr_thickf) forceVec(arr_lenf)} 31 10 top left
button $f2.button -text "Set vector's color" -command [list forceVectors_Color $t]
button $f2.reset -text "Reset vector's attributes" -command forceVectors_ResetAttributes
pack $f2.button -side left -padx 5 -pady 5
pack $f2.reset -side right -padx 5 -pady 5
#
# FRAME-3: Close + Update buttons
#
set b1 [button $f3.close -text "Close" -command [list DestroyWid $t]]
set b2 [button $f3.update -text "Update" \
-command [list forceUpdate $can]]
pack $b1 $b2 -side left -padx 3m -pady 3m -expand 1
}
proc forceUpdate can {
global forceVec mody
xc_forces $can scalefunction $forceVec(scalefunction)
xc_forces $can threshold $forceVec(threshold)
xc_forces $can lengthfactor $forceVec(lengthfactor)
xc_newvalue .mesa $mody(R_FORCE_RODTHICKF) $forceVec(rod_thickf)
xc_newvalue .mesa $mody(R_FORCE_ARRTHICKF) $forceVec(arr_thickf)
xc_newvalue .mesa $mody(R_FORCE_ARRLENF) $forceVec(arr_lenf)
eval xc_newvalue .mesa $mody(R_FORCE_COLOR) $forceVec(color)
$can render
xcSwapBuffers
}
proc forceVectors_Color {parent} {
global forceVec
set t [xcToplevel [WidgetName] \
"Set Color of Vectors" "Vector's color" $parent 0 0 1]
set init_color [rgb_f2h $forceVec(color)]
xcModifyColor $t "Set Color of Vectors" $init_color \
groove left left 100 100 70 5 20
set forceVec(colorID) [xcModifyColorGetID]
set ok [DefaultButton [WidgetName $t] -text "OK" \
-command [list forceVectors_ColorOK $t]]
set can [button [WidgetName $t] -text "Cancel" \
-command [list destroy $t]]
pack $ok $can -padx 10 -pady 10 -expand 1
}
proc forceVectors_ColorOK {t} {
global forceVec mody_col mody
set id $forceVec(colorID)
set alpha [lindex $forceVec(color) 3]
set forceVec(color) [list $mody_col($id,red) $mody_col($id,green) $mody_col($id,blue) $alpha]
# update vector's color
eval xc_newvalue .mesa $mody(L_FORCE_COLOR) $forceVec(color)
destroy $t
}
proc forceVectors_ResetAttributes {} {
global forceVec mody
# reset attributes
xc_resetvar .mesa $mody(R_FORCE_RODTHICKF)
xc_resetvar .mesa $mody(R_FORCE_ARRTHICKF)
xc_resetvar .mesa $mody(R_FORCE_ARRLENF)
xc_resetvar .mesa $mody(R_FORCE_COLOR)
# now reload default values to forceVec
set forceVec(rod_thickf) [xc_getdefault $mody(D_FORCE_RODTHICKF)]
set forceVec(arr_thickf) [xc_getdefault $mody(D_FORCE_ARRTHICKF)]
set forceVec(arr_lenf) [xc_getdefault $mody(D_FORCE_ARRLENF)]
set forceVec(color) [xc_getdefault $mody(D_FORCE_COLOR)]
}
|