This file is indexed.

/usr/share/xcrysden/Tcl/difIsoSurf.tcl is in xcrysden-data 1.5.53-1.

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
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
#############################################################################
# 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/difIsoSurf.tcl                                       
# ------                                                                    #
# Copyright (c) 1996-2003 by Anton Kokalj                                   #
#############################################################################

proc DiffIsoSurf_Widget {{dim 3D}} {
    global prop dif_isosurf
    
    # if window already exist destroy it and create new one
    catch {destroy .isodif}
    set t [xcToplevel .isodif "Difference Maps" "Diff. Maps" . 100 50]

    set top   [frame $t.top       -relief raised -bd 2]
    frame $t.mid
    set left  [frame $t.mid.left  -relief raised -bd 2]
    set right [frame $t.mid.right -relief raised -bd 2]
    set bot   [frame $t.bot       -relief raised -bd 2]

    pack $top $t.mid $bot -side top -expand 1 -fill both
    pack $left $right -side left -expand 1 -fill both

    set l  [label $top.l -text "Difference Map == Map A  -  Map B"]
    set ll [label $left.ll  -text "Map A:" -relief ridge -bd 2]
    set lr [label $right.lr -text "Map B:" -relief ridge -bd 2]
    pack $l -padx 10 -pady 10 -expand 1
    pack $ll $lr -side top -padx 10 -pady 5 -ipadx 20 -ipady 5

    ########################################
    # make everything that goes in LEFT FRAME
    frame $left.f1 -relief groove -bd 2
    set cb [checkbutton $left.f1.cb -text "" -anchor w]
    $cb invoke
    $cb config -state disabled

    set l [label $left.f1.l -textvariable dif_isosurf(file_textvar_A) \
	    -relief flat -anchor w]
    pack $cb $l -side top -fill x -padx 5 -pady 5
    frame $left.f2 -relief groove -bd 2

    xcMenuEntry $left.f2 "Specify Property:" 30  \
	    dif_isosurf(prop_A) $prop(dif_prop3D_list) \
	    -labelwidth 16
    pack $left.f1 $left.f2 -padx 2 -pady 10 -ipady 3 -fill x

    RadioButVarCmd $left "Density matrix to use:" \
	    dif_isosurf(denmat_A) DiffIsoSurf_LeftCmd \
	    top top 1 1 \
	    "SCF density matrix" \
	    "Density matrix as superposition of atomic densities"

    checkbutton $left.f.cb \
	    -text "Give new basis set/electron configuration" \
	    -variable prop(PATO_newbasis_A) \
	    -onvalue 1 \
	    -offvalue 0 \
	    -anchor e
    $left.f.cb config -state disabled

    pack $left.f.cb
    DiffIsoSurf_LeftCmd $dif_isosurf(denmat_A)

    ########################################
    # make everything that goes in RIGHT FRAME
    frame $right.f1 -relief groove -bd 2

    set cb [checkbutton $right.f1.cb \
	    -text "Load Crystal95's unit 9 from another file" \
	    -variable dif_isosurf(load_another_unit9) \
	    -onvalue 1 -offvalue 0 \
	    -command LoadAnotherUnit9 \
	    -anchor w]

    set l [label $right.f1.l -textvariable dif_isosurf(file_textvar_B) \
	    -relief flat -anchor w]
    pack $cb $l -side top -fill x -padx 5 -pady 5
    frame $right.f2 -relief groove -bd 2

    xcMenuEntry $right.f2 "Specify Property:" 30  \
	    dif_isosurf(prop_B) $prop(dif_prop3D_list) \
	    -labelwidth 16
    pack $right.f1 $right.f2 -side top -padx 2 -pady 10 -ipady 3 -fill x

    RadioButVarCmd $right "Density matrix to use:" \
	    dif_isosurf(denmat_B) DiffIsoSurf_RightCmd \
	    top top 1 1 \
	    "SCF density matrix" \
	    "Density matrix as superposition of atomic densities"

    checkbutton $right.f.cb \
	    -text "Give new basis set/electron configuration" \
	    -variable prop(PATO_newbasis_B) \
	    -onvalue 1 \
	    -offvalue 0 \
	    -anchor e
    $right.f.cb config -state disabled
    pack $right.f.cb
    DiffIsoSurf_RightCmd $dif_isosurf(denmat_B)

    ########################################
    # bottom frame
    set can [button $bot.can -text "Cancel" -command [list CancelProc $t]]
    set ok  [button $bot.ok  -text "OK" \
	    -command [list DiffIsoSurf_WidgetOK $t $dim]]
    pack $can $ok -side left -expand 1 -padx 10 -pady 10
}



proc DiffIsoSurf_LeftCmd item {
    global prop

    if { $item == "SCF density matrix" } {
	#.isodif.mid.left.f.cb configure -state disabled
	set prop(PATO_newbasis_A) 0
    } else {
	#.isodif.mid.left.f.cb configure -state normal
    }
}


proc DiffIsoSurf_RightCmd item {
    global prop 

    if { $item == "SCF density matrix" } {
	#.isodif.mid.right.f.cb configure -state disabled
	set prop(PATO_newbasis_B) 0
    } else {
	#.isodif.mid.right.f.cb configure -state normal
    }
}


proc LoadAnotherUnit9 {} {
    global dif_isosurf fileselect

    # if we turn off the "load another Crystal95's unit 9"
    if { $dif_isosurf(load_another_unit9) == 0 } { 
	set dif_isosurf(file_textvar) "loaded UNIT 9:\nsame as in Map A"
	return 
    } 

    fileselect "Load another Crystal95's unit 9"

    # maybe CANCEL button was pressed
    if { $fileselect(path) == "" } {
	set dif_isosurf(file_textvar) "loaded UNIT 9:\nsame as in Map A"
	set dif_isosurf(load_another_unit9) 0
	return
    }
    
    set dif_isosurf(file_textvar_B) "loaded Unit 9:\n$fileselect(path)"
    set dif_isosurf(unit9_B) $fileselect(path)
}


proc DiffIsoSurf_WidgetOK {t {dim 3D}} {
    global prop dif_isosurf system
    
    set match_A 0
    set match_B 0
    foreach item $prop(dif_prop3D_list) {
	if { $item == [string toupper $dif_isosurf(prop_A)] } {
	    set match_A 1
	}
	if { $item == [string toupper $dif_isosurf(prop_B)] } {
	    set match_B 1
	}
    }

    # replace {,} with " in prop(dif_prop3D_lis) and write to $p_list
    regsub -all -- \{ $prop(dif_prop3D_list) \" p_list
    regsub -all -- \} $p_list \" p_list
    
    if { $match_A == 0 } {
	tk_dialog [WidgetName] ERROR "ERROR: unknown property_A \"$dif_isosurf(prop_A)\", should be of $p_list" error 0 OK
	focus .isodif.mid.left.f2.e
	return
    }
    if { $match_B == 0 } {
	tk_dialog [WidgetName] ERROR "ERROR: unknown property_B \"$dif_isosurf(prop_A)\", should be of $p_list" error 0 OK
	focus .isodif.mid.right.f2.e
	return
    }
	
    # we come so far, so everything is OK
    set dif_isosurf(prop_A) [string toupper $dif_isosurf(prop_A)]
    set dif_isosurf(prop_B) [string toupper $dif_isosurf(prop_B)]

    if [winfo exist $t] {destroy $t}

    # now map from $dif_isosurf(prop_X) to {ECHG ECHD POTM}
    set propA [C95Name2Com $dif_isosurf(prop_A)]
    set propB [C95Name2Com $dif_isosurf(prop_B)]
    
    # now space-selection must be specefied 
    SetIsoGridSpace "Difference 3D Map - Grid Specifications" "Grid" \
	    [list $propA $propB] $dim
}



proc C95Name2Com name {
    # thiw proc map from property-name to C95's propertiy-command

    if { $name == "CHARGE DENSITY" } {
	return ECHG\n0
    } elseif { $name == "CHARGE DENSITY GRADIENT" } {
	return ECHG\n1
    } elseif { $name == "ELECTROSTATIC POTENTIAL" } {
	return POTM
    }
    return 0
}