source: trunk/expgui @ 969

Last change on this file since 969 was 953, checked in by toby, 12 years ago

remove NIST refs; review space in dir warn; use PGPLOT_DIR on win-NT, etc.; doc updates; tkcon problems on windows, add open console; move update menu from help

  • Property svn:keywords set to Author Date Revision Id
File size: 158.1 KB
RevLine 
[670]1#!/bin/sh
2# the next line restarts this script using wish found in the path\
3exec wish "$0" "$@"
4# If this does not work, change the #!/usr/bin/wish line below
5# to reflect the actual wish location and delete all preceeding lines
6#
7# (delete here and above)
8#!/usr/bin/wish
[82]9# $Id: expgui 953 2010-03-11 20:00:14Z toby $
[10]10set expgui(Revision) {$Revision: 953 $ $Date: 2010-03-11 20:00:14 +0000 (Thu, 11 Mar 2010) $}
11
[801]12package require Tk
13
[10]14# to do:
15#
[30]16# need to change heading and button label depending on where getExpFileName
[139]17# is called from?
[28]18#
[10]19# global background editing & profile work differently: should both
20# start out blank with a "load from option"?
21#
22# idea:
23# a scroll list for all histogram refinement flags ; click on takes you to the
24# appropriate menu.
25#
[93]26# idea:
27#   change cell parameters to labels and have a edit cell button
28#   that enforces metric symmetry
29#
[20]30# to allow "global" access on phase page
[10]31#   change buttons from radio to multiple
32#   -- or display all 9 cell flag/damps and all atoms
33#   make editMultipleRecords work with multiple phases, also cell flag/damp
34#   blank cell entries
35#   add phase to atom number in listing
[539]36#   DisplayAllAtoms needs to loop over phases
[139]37#
38# idea: load more than one bank from a multi-bank .RAW file
39#
[16]40if {$tcl_version < 8.0} {
41    tk_dialog .expFileErrorMsg "Version Error" \
[327]42            "EXPGUI requires Tcl/Tk version 8.0 or higher" error 0 "Exit"
[16]43    exit
44}
45
[306]46# the shell variable tells expgui to act as a shell for GSAS
47# as well as edit the EXP file. In no-shell mode, it can
48# only be used to edit the .EXP file as a callable program
49set expgui(shell) 1
50catch {if $env(EXPGUI_NOSHELL) {set expgui(shell) 0}}
51
[10]52if {$argv != ""} {
[306]53    if {[string match *noshell* [string tolower $argv]]} {
[539]54        #puts noshell
[306]55        set expgui(shell) 0
56        set expgui(expfile) [lindex $argv 1]
57    } else {
58        set expgui(expfile) [lindex $argv 0]
59    }
[809]60    if {[string match -nocase {.o[0-9a-f][0-9a-f]} \
61             [file extension $expgui(expfile)]]} {
62        # this is an archived file -- need to handle this later
63        tk_dialog .expFileErrorMsg "No archived files yet" \
64            "At present you must use open (expnam) to open archived files" warning 0 "Continue"
65        set expgui(expfile) {}
66    } elseif {[string toupper [file extension $expgui(expfile)]] != ".EXP"} {
[20]67        append expgui(expfile) ".EXP"
[10]68    }
69} else {
[30]70    set expgui(expfile) {}
[10]71}
72
[199]73set expgui(curhist) {}
74set expmap(powderlist) {}
[306]75set expgui(bkgcolor1) #fdf
[199]76
[10]77set expgui(debug) 0
78catch {if $env(DEBUG) {set expgui(debug) 1}}
79#set expgui(debug) 1
80
[395]81# location for web pages, if not found locally
[953]82set expgui(website) 11bm.xor.aps.anl.gov/expguidoc/
[93]83# default for archive mode = on
[10]84set expgui(archive) 1
[93]85# default for autoexec load = off
86set expgui(autoexpload) 0
[908]87# default for execprompt = on
88set expgui(execprompt) 1
[670]89# default for autostart GRWND = off
90set expgui(autoGRWND) 0
[359]91# by default expgui is iconified while GENLES, etc runs
92set expgui(autoiconify) 1
[323]93# default for show EXPTOOL output = off
94set expgui(showexptool) 0
[20]95# save the name of the wish executable
96set wishshell [info nameofexecutable] 
97# misc constants
[181]98set txtvw(font) "Courier"
[419]99set expgui(font) 14
[20]100set liveplot(hst) 1
101set liveplot(legend) 1
[30]102set expgui(filesort) 1
[118]103set expgui(initstring) {}
[139]104# use a separate window for DISAGL (default)
105set expgui(disaglSeparateBox) 1
[234]106set expgui(DefaultPeakType) 0
[267]107# default: keep current atoms when replacing a phase
108set expgui(DeleteAllAtoms) 0
[670]109# flags for running POWPREF
110set expgui(needpowpref) 0
111set expgui(needpowpref_why) ""
[774]112# on Mac associate a app with .EXP file (on by default)
113set expgui(MacAssignApp) 1
[10]114#=============================================================================
[20]115#----------------------------------------------------------------
116# where are we?
117set expgui(script) [info script]
118# translate links -- go six levels deep
119foreach i {1 2 3 4 5 6} {
120    if {[file type $expgui(script)] == "link"} {
121        set link [file readlink $expgui(script)]
122        if { [file  pathtype  $link] == "absolute" } {
[48]123            set expgui(script) $link
[20]124        } {
125            set expgui(script) [file dirname $expgui(script)]/$link
126        }
127    } else {
128        break
129    }
130}
[21]131# fixup relative paths
132if {[file pathtype $expgui(script)] == "relative"} {
133    set expgui(script) [file join [pwd] $expgui(script)]
134}
[20]135set expgui(scriptdir) [file dirname $expgui(script) ]
136set expgui(gsasdir) [file dirname $expgui(scriptdir)]
[21]137set expgui(gsasexe) [file join $expgui(gsasdir) exe]
[323]138set expgui(docdir) [file join $expgui(scriptdir) doc]
[20]139#----------------------------------------------------------------
[823]140# use EXPGUI directory for packages
[133]141lappend auto_path $expgui(scriptdir)
142#----------------------------------------------------------------
[419]143source [file join $expgui(scriptdir) opts.tcl]
[20]144# fetch EXP file processing routines
145source [file join $expgui(scriptdir) readexp.tcl]
146# commands for running GSAS programs
147source [file join $expgui(scriptdir) gsascmds.tcl]
148# contents of GSAS menus
149source [file join $expgui(scriptdir) gsasmenu.tcl]
[93]150# commands for adding phases, histograms & atoms
151source [file join $expgui(scriptdir) addcmds.tcl]
[199]152# commands for preferred orientation
153source [file join $expgui(scriptdir) orient.tcl]
[478]154# setting data range/excluded regions
155source [file join $expgui(scriptdir) exclinit.tcl]
[20]156#---------------------------------------------------------------------------
157# override options with locally defined values
[696]158set filelist [file join $expgui(scriptdir) localconfig]
159if {$tcl_platform(platform) == "windows"} {
160    lappend filelist "c:/gsas.config"
161} else {
162    lappend filelist [file join ~ .gsas_config]
163}
[645]164if {[catch {
[670]165    foreach file $filelist {
[660]166        if [file exists $file] {source $file}
[645]167    }
168} errmsg]} {
[660]169    set msg "Error reading file $file (aka [file nativename $file]): $errmsg"
[645]170    MyMessageBox -parent . -title "Customize warning" \
[660]171        -message $msg -icon warning -type Ignore -default ignore \
172        -helplink "expguierr.html Customizewarning"
[20]173}
[823]174if [catch {package require BWidget}] {
175    set msg "Error loading the BWidget package: This should not happen if you use the version of Tcl/Tk (tcl84+.exe) distributed with EXPGUI"
176    MyMessageBox -parent . -title "BWidget load error" \
177        -message $msg -icon error -type Error -default error \
178        -helplink "debug_tcltk.html"
179    destroy .
180}
[419]181SetTkDefaultOptions $expgui(font)
[327]182#---------------------------------------------------------------------------
[838]183set expgui(resize) 0
[327]184# platform-specific code
[359]185if {$tcl_platform(platform) == "windows" \
186        && $tcl_platform(os) == "Windows 95"} {
[809]187    if {[catch {package require winexec}] && \
188            [catch {package require winutils}]} {
[327]189        MyMessageBox -parent . -title "WINEXEC Error" \
[809]190            -message "Error -- Unable to load the WINEXEC or the WINUTILS package. This is needed in Win-95/-98/-me machines" \
191            -icon error -type Quit -default quit \
192            -helplink "expgui_Win_readme.html Winexec"
[327]193        destroy .
194    }
[359]195}
196if {$tcl_platform(platform) == "windows"} {
197    # check the path -- can DOS use it?
198    if {[string first {\\} $expgui(script) ] != -1} {
[605]199        MyMessageBox -parent . -title "Networked Path" \
200                -message "Note -- You may have problems running EXPGUI/GSAS from a network drive. If you have errors, map $expgui(gsasdir) to a \"drive\" (like F:). (Use \"Map network drive\" to do this.)" \
201                -icon error -type {"Be brave"} -default "be brave" \
[359]202                -helplink "expgui_Win_readme.html NetPath"
203    }
[884]204    # test for error in COMSPEC environment variable setting
205    if {[catch {set env(COMSPEC)}]} { 
206        set env(COMSPEC) {c:\windows\system32\cmd.exe}
207    }
208    if {[auto_execok $env(COMSPEC)] == ""} {
209        set env(COMSPEC) {c:\windows\system32\cmd.exe}
210    }
211    if {[auto_execok $env(COMSPEC)] == ""} {
212        MyMessageBox -parent . -title "cmd.exe not found " \
213            -message "The Windows COMSPEC environment variable is not set correctly. It should point to a copy of cmd.exe (or command.com on -95/-98/-ME). EXPGUI also did not find the file in the most common location (c:\\windows\\system32\\cmd.exe). Please fix: EXPGUI cannot run GSAS programs without this." \
214            -icon "error" -type Sorry -default sorry
215#           -helplink "expguierr.html Customizewarning"
216        return
217    }
[359]218    set expgui(exptool) [file join $expgui(gsasexe) exptool.exe]
[93]219} else {
220    set expgui(exptool) [file join $expgui(gsasexe) exptool]
[746]221    if {$tcl_platform(os) != "Darwin"} {
222        if [catch {set env(GSASBACKSPACE)}] {set env(GSASBACKSPACE) 1}
223    }
[93]224}
[935]225# find the pgplot directory
226if {[file exists [file join $expgui(gsasdir) pgl]]} {
227    set expgui(pgplotdir) [file join $expgui(gsasdir) pgl]
228} elseif {[file exists [file join $expgui(gsasdir) pgplot]]} {
229    set expgui(pgplotdir) [file join $expgui(gsasdir) pgplot]
230} else {
231    MyMessageBox -parent . -title "PGPLOT Error" \
232            -message "Warning -- Unable to find PGPLOT directory. GSAS graphics will not work without this. Is GSAS correctly installed?" \
233            -icon warning -type {"Limp Ahead"} -default "Limp Ahead" \
234            -helplink "expguierr.html NoPGPLOT"
235    set expgui(resize) 1
236    set expgui(pgplotdir) ""
[395]237}
[935]238
239# old code used to fix up bad distributions
240# # do we have a PGPLOT fonts file?
241# # if it is in the "wrong" place/name -- make it "right"
242# if {![file exists [file join $expgui(pgplotdir) grfont.dat]] && \
243#       [file exists [file join $expgui(gsasdir) fonts grfont.dat]]} {
244#     catch {file mkdir $expgui(pgplotdir)}
245#     file copy [file join $expgui(gsasdir) fonts grfont.dat] \
246#           [file join $expgui(pgplotdir)  grfont.dat]
247# }
248# # do we have a PGPLOT fonts file?
249# if {![file exists [file join $expgui(pgplotdir)  grfont.dat]] && \
250#       [file exists [file join $expgui(gsasdir) fonts pgfont.dat]]} {
251#     catch {file mkdir $expgui(pgplotdir)}
252#     file copy [file join $expgui(gsasdir) fonts pgfont.dat] \
253#           [file join $expgui(pgplotdir)  grfont.dat]
254# }
255# # do we have a PGPLOT fonts file?
256# if {![file exists [file join $expgui(pgplotdir)  grfont.dat]] && \
257#       [file exists [file join $expgui(pgplotdir)  pgfont.dat]]} {
258#     file copy [file join $expgui(pgplotdir)  pgfont.dat] \
259#           [file join $expgui(pgplotdir)  grfont.dat]
260# }
261
262if { $expgui(pgplotdir) != "" && \
263         ![file exists [file join $expgui(pgplotdir)  grfont.dat]]} {
[390]264    MyMessageBox -parent . -title "PGPLOT Error" \
[935]265            -message "Warning -- Unable to find file GRFONT.DAT in $expgui(pgplotdir). GSAS graphics will not work. Is GSAS correctly installed?" \
[390]266            -icon warning -type {"Limp Ahead"} -default "Limp Ahead" \
267            -helplink "expguierr.html NoPGPLOT"
[838]268    set expgui(resize) 1
[390]269}
[20]270#---------------------------------------------------------------------------
[30]271if {$expgui(expfile) != ""} {
[670]272   # is there a space in the EXP name?
273    if {[string first " " [file tail $expgui(expfile)]] != -1} {
[30]274        update
[605]275        MyMessageBox -parent . -title "File Name Error" \
[670]276                -message "File name \"$expgui(expfile)\" is invalid -- EXPGUI cannot process experiment files with spaces in the name" \
[605]277                -icon warning -type Continue -default continue
278#               -helplink "expguierr.html OpenErr"
279        set expgui(expfile) {}
[838]280        set expgui(resize) 1
[670]281    } elseif {[string first " " [file dirname $expgui(expfile)]] != -1} {
282        update
283        MyMessageBox -parent . -title "Good luck..." \
[908]284            -message "You are using a directory with a space in the name ([file dirname $expgui(expfile)]) -- You may encounter bugs in EXPGUI. Please e-mail them to Brian.Toby@ANL.gov so they can be fixed." \
[670]285                -icon warning -type Continue -default continue
286#               -helplink "expguierr.html OpenErr"
[838]287        set expgui(resize) 1
[605]288    } elseif ![file exists $expgui(expfile)] {
289        update
[327]290        set ans [
291        MyMessageBox -parent . -title "File Open Error" \
292                -message "File [file tail $expgui(expfile)] does not exist in [file dirname $expgui(expfile)]. OK to create?" \
293                -icon question -type {"Select other" "Create"} -default "select other" \
294                -helplink "expguierr.html OpenErr"
295        ]
296        if {[string tolower $ans] != "create"} {set expgui(expfile) {}}
[838]297        set expgui(resize) 1
[20]298    }
299}
[30]300if {$expgui(expfile) == ""} {
[935]301    # place the parent window because the getExpFileName window will be centered above it.
302    set xpos 50
303    catch {set xpos $expgui(xpos)}
304    set ypos 50
305    catch {set ypos $expgui(ypos)}
306    set poscntr "nsew"
[944]307    catch {set poscntr $expgui(poscenter)}
[935]308    LocateWindow "." $xpos $ypos $poscntr
[74]309    # windows needed this update before when using tk_getOpenFile.
310    # I am not sure it is still needed.
[30]311    update
[74]312    #
[746]313    set expgui(expfile) [getExpFileName ""]
314    set expgui(resize) 1
[30]315}
316if {$expgui(expfile) == ""} exit
317# you've been warned this .EXP does not exist!
318if ![file exists $expgui(expfile)] {
319    # create an "empty" exp file
320    createexp $expgui(expfile) \
321            [getstring "title for experiment $expgui(expfile)" 60 0]
322}
323catch {cd [string trim [file dirname $expgui(expfile)]]}
[20]324
[10]325#
326# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
327# <<<<<<<<<<    BEGINNING OF MAIN: GLOBAL AREA FOR DATA EXTRACTION >>>>>>>>>>>
328# >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
329# load exp file and set up dialogs
330proc loadexp {expfile} {
[133]331    global expgui expmap entryvar entrycmd tcl_platform
[390]332    # is this a compressed archive file?
[466]333    if {[string match {*.O[0-9A-F][0-9A-F]} $expfile]} {
334        set expgui(expfile) {}
335        set expnam [file rootname $expfile]
[390]336        set ans [MyMessageBox -parent . -title "Load Archived File" \
[466]337                -message "Loading archived version of $expnam. Do you want to continue using the same experiment name or work with the archived version under a new name?" \
338                -icon question -type "{Use New Name} {Continue with current}" \
[390]339                -default {Use New Name} \
340                -helplink "expguierr.html LoadArchived"
341        ]
[466]342        # archive the current .EXP file
343        if {$ans != "use new name" && [file exists $expfile]} {
344            # get the last archived version
[478]345            set lastf [lindex [lsort [glob -nocomplain $expnam.{O\[0-9A-F\]\[0-9A-F\]}]] end]
[466]346            if {$lastf == ""} {
347                set num 01
348            } else {
349                regexp {.*\.O([0-9A-F][0-9A-F])$} $lastf a num
350                scan $num %x num
351                if {$num >= 255} {
352                    set num FF
353                } else {
354                    set num [string toupper [format %.2x [incr num]]]
355                }
356            }
357            catch {
358                set newfile $expnam.O$num
359                file rename -force $expnam.EXP $newfile
[484]360                set fp [open $expnam.LST a+]
361                puts $fp "\n----------------------------------------------"
362                puts $fp "     Regressing to archive file [file tail $expfile]"
363                puts $fp "     but first archiving [file tail $expnam.EXP] as [file tail $newfile]"
364                puts $fp "----------------------------------------------\n"
365                close $fp
[466]366            }
[746]367            file copy -force $expfile $expnam.EXP
368            set expfile $expnam.EXP
[390]369        }
370        if {$ans == "use new name"} {
[746]371            set newexpfile [getExpFileName new]
372            if {$newexpfile == ""} return 
373            file copy -force $expfile $newexpfile
374            catch {cd [string trim [file dirname $expgui(expfile)]]}
375            set expfile [file tail $newexpfile]
[750]376            set expgui(needpowpref) 2
377            set expgui(needpowpref_why) "\tA new .EXP file was created\n" 
[390]378        }
[746]379    }
380    set expgui(expfile) $expfile
[774]381
382    # change the icon and assign an app to this .EXP file
383    global tcl_platform
384    if {$tcl_platform(os) == "Darwin" && $expgui(MacAssignApp)} {
385        MacSetResourceFork $expfile
386    }
387
388    # read in the .EXP file
[746]389    set fmt [expload $expfile]
390    # if the file was not in the correct format, force a rewrite before use
391    if {$fmt < 0} {
392        # read error
393        return
394    } elseif {$fmt == 1} {
395        set expgui(changed) 0
[385]396    } else {
[746]397        set expgui(changed) 1
[133]398    }
399    # force exp files to be upper case, force save if name changes
400    set filetail [file tail $expfile]
401    set filetailcaps [string toupper $filetail]
402    if {$tcl_platform(platform) == "unix" && $filetail != $filetailcaps} {
403        set expgui(changed) 1
404    }
405    if {[file dirname $expfile] == "."} {
406        set expgui(expfile) $filetailcaps
407    } else {
408        set expgui(expfile) [file join \
409                [file dirname $expfile] $filetailcaps]
410    }
411
412    mapexp
[908]413    if {$expgui(MacroRunning)} {
414        set expgui(MacroStatus) "?"
415        catch {
416            set expgui(MacroStatus) \
417                "[expinfo cyclesrun] cycles; Chi2=$expgui(LastChi2); $expgui(LastnVars) vars"
418        }
419    }
420
[133]421    set expgui(expModifiedLast) 0
[10]422    catch {
[133]423        set expgui(expModifiedLast) [file mtime $expgui(expfile)]
[10]424    }
[20]425    set expgui(last_History) [string range [string trim [lindex [exphistory last] 1]] 0 50 ]
[10]426    # set the window/icon title
[866]427    wm title . "EXPGUI interface to GSAS: $expfile"
[16]428    set expgui(titleunchanged) 1
[10]429    wm iconname . [file tail $expfile]
430
[908]431    # load current information into EXPGUI pages
432    PaintEXPGUIpages
433
434    # start checking for external changes
435    afterawhile
436}
437
438# [re]load all screens with current state of EXPGUI file
439proc PaintEXPGUIpages {} {
440    global expgui expmap entryvar entrycmd
[539]441    # reset the phase buttons
442    set expgui(curPhase) ""
[10]443    # set the number of phases on the phase page
444    setphases
445
[234]446    # disable the "global options" that don't make sense based on
447    # which histograms present
[10]448    foreach num {1 2 3 4 5} {
449        set flag($num) 0
450    }
[234]451    # save a list of the allowed modes, too
452    set expgui(AllowedHistSelectModes) {0 6}
[10]453    foreach h $expmap(powderlist) {
454        if {[string range $expmap(htype_$h) 2 2] == "T"} {set flag(1) 1}
455        if {[string range $expmap(htype_$h) 1 2] == "NC"} {set flag(2) 1}
456        if {[string range $expmap(htype_$h) 1 2] == "XC" && \
457                [histinfo $h lam2] != 0.0} {set flag(3) 1}
458        if {[string range $expmap(htype_$h) 1 2] == "XC" && \
459                [histinfo $h lam2] == 0.0} {set flag(4) 1}
460        if {[string range $expmap(htype_$h) 1 2] == "XE"} {set flag(5) 1}
461    }
462    foreach num {1 2 3 4 5} \
463            lbl {TOF "CW Neutron" "Alpha12 Xray" "Monochromatic Xray" \
464            "Energy Disp Xray"} {
465        if $flag($num) {
[16]466            $expgui(fm).option.menu.editmode entryconfigure $lbl -state normal
[234]467            lappend expgui(AllowedHistSelectModes) $num
[10]468        } else {
[16]469            $expgui(fm).option.menu.editmode entryconfigure $lbl -state disabled
[10]470        }
471    }
[12]472    # disable traces on entryvar until we are ready
473    set entrycmd(trace) 0
[117]474    trace vdelete entryvar w entvartrace
[12]475
[28]476    # propogate changes on the least squares page
[10]477    set entryvar(cycles) [expinfo cycles]
478    set entrycmd(cycles) "expinfo cycles"
[838]479    set entryvar(mbw) [expinfo mbw]
480    set entrycmd(mbw) "expinfo mbw"
[112]481    # set expgui(globalmode) 0
[10]482    set expgui(printopt) "Print Options ([expinfo print])"
[20]483    set entryvar(title) [expinfo title]
[10]484    global printopts
485    foreach num [array names printopts] {
486        set entrycmd(printopt$num) "printsetting $num"
487        set entryvar(printopt$num) [printsetting $num]
488    }
[12]489    # enable traces on entryvar
490    set entrycmd(trace) 1
[117]491    trace variable entryvar w entvartrace
492
[133]493    # set fo extraction on LS page
[20]494    SetupExtractHist
[399]495    # set convergence criterion
496    InitLSvars
[28]497
[238]498    # update the histogram list & update the page
[28]499    sethistlist
[10]500}
501
[93]502# called to reread the .EXP file
[10]503proc rereadexp {expfile} {
504    global expgui
505    if $expgui(changed) {
506        set decision [tk_dialog .instrSaveData {Save .EXP changes} \
507                {You have made changes to the Experiment. Rereading will cause the changes to be lost. Select an option:} \
508                {} 0 "Save and reread" "Reread without Save" "Cancel reread command"]
509        switch $decision {
[20]510            0 { savearchiveexp }
511            1 { }
512            2 { return }
[10]513        }
514    }
[20]515    loadexp $expgui(expfile)
[10]516}
517
518proc SaveAsFile {} {
[20]519    global expgui
[30]520    set newexpfile [getExpFileName new]
[10]521    if {$newexpfile == ""} return 
[71]522    expwrite $newexpfile
[20]523    set expgui(expfile) $newexpfile
[21]524    catch {cd [string trim [file dirname $expgui(expfile)]]}
[774]525    # change the icon and assign an app to this .EXP file
526    global tcl_platform
527    if {$tcl_platform(os) == "Darwin" && $expgui(MacAssignApp)} {
528        MacSetResourceFork $expgui(expfile)
529    }
[10]530    set expgui(changed) 0
[20]531    set expgui(expModifiedLast) [file mtime $expgui(expfile)]
532    set expgui(last_History) [string range [string trim [lindex [exphistory last] 1]] 0 50 ]
[10]533    # set the window/icon title
[20]534    wm title . $expgui(expfile)
[16]535    set expgui(titleunchanged) 1
[20]536    wm iconname . [file tail $expgui(expfile)]
[399]537    # set convergence criterion
538    InitLSvars
[750]539    set expgui(needpowpref) 2
540    set expgui(needpowpref_why) "\tA new .EXP file was created\n" 
[10]541}
542
543# called to read a different .EXP file
544proc readnewexp {} {
[30]545    global expgui expmap
[10]546    if $expgui(changed) {
547        set decision [tk_dialog .instrSaveData {Save .EXP changes} \
[466]548                {You have made changes to the Experiment. Reading a different file without first saving will cause the changes to be lost. Select an option:} \
[10]549                {} 0 "Save and read" "Read without Save" "Cancel read command"]
550        switch $decision {
[20]551            0 { savearchiveexp }
552            1 {                }
553            2 { return }
[10]554        }
555    }
[30]556    set newexpfile [getExpFileName old]
[10]557    if {$newexpfile == ""} return 
[71]558
[152]559    # switch to the 1st page
560    RaisePage lsFrame
561
[71]562    if ![file exists $newexpfile] {
563        # you've been warned this .EXP does not exist!
564        # create an "empty" exp file
565        createexp $newexpfile \
566                [getstring "title for experiment $newexpfile" 60 0]
567    }
[20]568    set expgui(expfile) $newexpfile
[21]569    catch {cd [string trim [file dirname $expgui(expfile)]]}
[112]570    set expgui(globalmode) 0
[20]571    loadexp $expgui(expfile)
[30]572
[539]573    # reset the phase selection
574    set expgui(curPhase) {}
575
[30]576    # select the first histogram in the list by default (if there are any)
[71]577    if {[llength $expmap(histlistboxcontents)] > 0} {
578        set expgui(curhist) 0
579    } else {
580        set expgui(curhist) {}
581    }
[117]582    if {[CountHistory] > 100} {
583        DeleteHistoryRecords "This .EXP file has [CountHistory] history records\nErasing most will speed EXPGUI"
584    }
[10]585}
586
587#------------- set up data read/write layer ----------------------
588# trace routine on entryvar
589proc entvartrace {array elem action} {
[303]590    global expgui entrycmd entryvar entrybox
[10]591    if !$entrycmd(trace) return
592   
593    catch {
[20]594        if {$entrycmd($elem) == ""} return
[10]595        incr expgui(changed)
596        if $expgui(debug) {puts "$entrycmd($elem)  set $entryvar($elem) "}
597        if {$entrycmd($elem) == ""} return
598        if [catch {
[303]599            set result [eval $entrycmd($elem) set [list $entryvar($elem)]]
600            if {!$result} {
601                if $expgui(debug) {puts "error with $entrycmd($elem)"}
602                catch {$entrybox($elem) config -fg red}
603            } else {
604                catch {$entrybox($elem) config -fg black}
[908]605                RecordMacroEntry "incr expgui(changed); $entrycmd($elem) set [list $entryvar($elem)]" 1
[303]606            }
[539]607            if {[string match "*atominfo" [lindex $entrycmd($elem) 0]]} {
608                after idle "UpdateAtomLine \
609                        [list [lindex $entrycmd($elem) 2]] \
610                        [lindex $entrycmd($elem) 1]"
[10]611            }
[267]612        } errmsg] {error $errmsg}       
[10]613    }
614}
615
616# disable traces on entryvar until we are ready
617set entrycmd(trace) 0
618trace variable entryvar w entvartrace
619
620#
621#
622#
623##############################################################################
624#####                    #####################################################
625##### PROCEDURES SECTION #####################################################
626#####                    #####################################################
627##############################################################################
628
[670]629# save some of the global options in ~/.gsas_config or ~/gsas.config in Windows
[28]630proc SaveOptions {} {
[478]631    global expgui env tcl_platform graph peakinfo
[670]632    if {$tcl_platform(platform) == "windows"} {
[696]633        set fp [open c:/gsas.config a]
[670]634    } else {
635        set fp [open [file join ~ .gsas_config] a]
636    }
637
[776]638    puts $fp "# EXPGUI saved options from [clock format [clock seconds]]"
[774]639    set itemlist {archive asorttype hsorttype filesort disaglSeparateBox \
[908]640        font autoexpload autoiconify autotick execprompt ShowGENLES}
[774]641    if {$tcl_platform(os) == "Darwin"} {
642        lappend itemlist MacAssignApp
643    }
644    if {$tcl_platform(platform) == "windows" && \
645            $tcl_platform(os) == "Windows 95"} {
646        lappend itemlist autoGRWND
647    }
648    foreach item $itemlist {
[660]649        puts $fp "set expgui($item) [list $expgui($item)]"
[139]650    }
[28]651    if {$tcl_platform(platform) != "windows"} {
[660]652        puts $fp "set env(GSASBACKSPACE) [list $env(GSASBACKSPACE)]"
[28]653    }
[478]654    foreach v {printout legend outname outcmd autoraise color_excl \
655            color_obs color_calc} {
[660]656        puts $fp "set graph($v) [list $graph($v)]"
[478]657    }
658    foreach v {obssym obssize exclsym exclsize} {
[660]659        puts $fp "set peakinfo($v) [list $peakinfo($v)]"
[478]660    }
[28]661    close $fp
662}
663
[944]664# save some of the global options in ~/.gsas_config or ~/gsas.config in Windows
665proc SavePosition {} {
666    global expgui env tcl_platform graph peakinfo
667    if {$tcl_platform(platform) == "windows"} {
668        set fp [open c:/gsas.config a]
669    } else {
670        set fp [open [file join ~ .gsas_config] a]
671    }
672
673    puts $fp "# EXPGUI saved screen position on [clock format [clock seconds]]"
674    puts $fp "set expgui(xpos) [expr {[winfo x .]*100./[winfo screenwidth .]}]"
675    puts $fp "set expgui(ypos) [expr {[winfo y .]*100./[winfo screenheight .]}]"
676    puts $fp "set expgui(poscenter) NW"
677    close $fp
678}
679
[10]680proc About { } {
681    global expgui expmap
[952]682    set version "Revision [lindex $expgui(Revision) 1] (readexp.tcl [lindex $expmap(Revision) 1])"
683    catch {
684        set version "[GetSVNVersion $::expgui(gsasdir)]\ntrack: $::expgui(SVNversion)"
685    }
[10]686    tk_dialog .about {About...} \
687"EXPGUI\n\
[838]688Brian Toby\n\
[908]689APS, Argonne National Laboratory\n\n\
[838]690Not subject to copyright\n\n\
[952]691$version\n\n\
[449]692Cite: B. H. Toby, EXPGUI, a graphical\n\
693user interface for GSAS, J. Appl. Cryst.\n\
69434, 210-21 (2001).
695\n\n\
696Generalized Structure Analysis System\n(GSAS)\n\
[953]697R. B. Von Dreele,\n
698APS, Argonne National Laboratory\n
699and A. C. Larson, Los Alamos (retired)\n\n\
[20]700" \
[10]701        info 0 OK
702}
[449]703proc Cite { } {
704    global expgui expmap
705    tk_dialog .about {Citations...} \
706"If you use EXPGUI, please cite\n\n\
707B.H. Toby, EXPGUI, a graphical\n\
708user interface for GSAS, J. Appl. Cryst.\n\
70934, 210-21 (2001).\n\n\
710as well as\n\n\
711A.C. Larson and R.B. Von Dreele,\n\
712\"General Structure Analysis System (GSAS)\",\n\
713Los Alamos National Laboratory Report\n\
[838]714LAUR 86-748 (2004)." \
[449]715        info 0 OK
716}
[10]717
[774]718# this proc is no longer called, but I am leaving it here as it may
719# be of use in the future
[760]720proc MakeAppleScript {} {
721    global wishshell expgui
722    # create a local script directory, if it does not exist
723    if {![file exists ~/Library/Scripts]} {file mkdir  ~/Library/Scripts}
[761]724    set tmpfile [file nativename ~/tmpscriptfile]
725    set startdir [tk_chooseDirectory -initialdir [pwd] -mustexist 1 \
726                      -title "Choose GSAS starting directory"]
727    if {$startdir == ""} {set startdir "~"}
[760]728    set dir [file nativename  ~/Library/Scripts]
[761]729    if {[set xterm [auto_execok xterm]] == ""} {
730        MyMessageBox -parent . -title "xterm not found " \
731            -message "The AppleScript could not be created because the X11 xterm application was not found. Please correct your path and try again." \
732            -icon "error" -type Sorry -default sorry
733#           -helplink "expguierr.html Customizewarning"
734        return
735    }
736    set file [tk_getSaveFile -initialdir $dir -initialfile EXPGUI.scpt \
737                      -title "Choose location to save script"]
738    set path {$PATH:}
739    append path [file dirname $xterm]
[760]740    set fp [open $tmpfile w]
741    # the applescript starts here
742    puts $fp {on run}
743    puts $fp {  tell application "Finder"}
744    puts $fp {     launch application "X11"}
745    puts $fp {   end tell}
[761]746    puts $fp "  set results to do shell script \"cd $startdir; DISPLAY=:0.0 PATH=$path $wishshell $expgui(script)  > /dev/null 2>&1 &\""
[760]747    puts $fp {end run}
[763]748    # drag & drop
749    puts $fp {on open these_files}
750    puts $fp {  tell application "Finder"}
751    puts $fp {     launch application "X11"}
752    puts $fp {   end tell}
753    puts $fp {  repeat with this_file in these_files}
754    puts $fp "  set results to do shell script \"cd $startdir; DISPLAY=:0.0 PATH=$path $wishshell $expgui(script) \" & the quoted form of the POSIX path of this_file & \" > /dev/null 2>&1 &\""
755    puts $fp {  end repeat}
756    puts $fp {end open}
[760]757    close $fp
758    if {[catch {
759        exec osacompile -l AppleScript -o $file  $tmpfile
760        file delete -force $tmpfile
761        MyMessageBox -parent . -title "AppleScript created" \
762            -message "Script $file created & compiled. You may wish to use the Script Editor to save this as an application" \
[763]763            -icon "info" -type OK -default ok \
764            -helplink "osx.html CompileAppleScript"
[760]765    } errmsg]} {
766        MyMessageBox -parent . -title "AppleScript warning" \
[763]767        -message "An error occurred while attempting to create the script. Please report this bug, including these details:\n$errmsg"\
[760]768            -icon warning -type Ignore -default ignore
769#       -helplink "expguierr.html Customizewarning"
770    }
771}
[774]772
[539]773# this proc gets called when the export coordinate button is pressed
774# it loads export
775proc BuildCoordExpMenu {menu} {
776    global expgui_cmdlist expgui
777    # do this only once
778    $menu config -postcommand {}
779    $menu delete 1 end
780    # add the cif export routine
781    set cmd gsas2cif
782    set action {}
783    catch {set action [lindex $expgui_cmdlist($cmd) 0]}
784    if {$action != "" && $action != "-"} {
785        $menu add command -label $cmd -command [subst $action]
786    }
787    # get a list of files to read
788    set filelist [glob -nocomplain [file join $expgui(scriptdir) export_*.tcl]]
789    foreach file $filelist {
790        source $file
791        $menu add command -label $label -command $action
792    }
793}
794
795# utility export routines for the export_*.tcl files:
796# make a box for export
797proc MakeExportBox {win title webref} {
798    global expmap expgui
799    catch {destroy $win}
800    toplevel $win
801    wm title $win "Export coordinates"
802    if {$webref != ""} {
803        bind $win <Key-F1> $webref
804    }
805    pack [label $win.lbl -text $title] -side top -anchor center
806    pack [frame $win.ps] -side top -anchor w
807    pack [label $win.ps.lbl -text "Select phase: "] -side left
808    foreach num $expmap(phaselist) {
809        pack [button $win.ps.$num -text $num \
810                    -command "SetExportPhase $num $win"] -side left
811    }
812    # leave a place for format-specific items
813    pack [frame $win.special] -side top
814    pack [frame $win.but] -side top -fill x -expand yes
815    pack [button $win.but.1 -text Write -command "destroy $win"] -side left
816    SetExportPhase [lindex $expmap(phaselist) 0] $win
[928]817    pack [button $win.but.2 -text Cancel \
[539]818            -command "set expgui(export_phase) 0;destroy $win"] -side left
819    pack [button $win.but.help -text Help -bg yellow \
820            -command "MakeWWWHelp expgui.html ExportMSI"] \
821            -side right
822}
823
824# set the phase in response to the button
825proc SetExportPhase {num win} {
826    global expmap expgui
827    foreach n $expmap(phaselist) { 
828        if {$n == $num} { 
829            $win.ps.$n config -relief sunken; 
830            set expgui(export_phase) $num 
831        } else { 
832            $win.ps.$n config -relief raised
833        }
834    }
835}
836
[16]837# wait until idle
838proc afterawhile {} {
839    # cancel any other instances of this loop
840    after cancel afterawhile
841    after cancel whenidle
842    after cancel whenidle
843    after idle whenidle
844}
845
[133]846# This is called every 2 seconds to check for changes to the .EXP file
[16]847proc whenidle {} {
[359]848    global expgui tcl_platform
[16]849    if $expgui(titleunchanged) {
850        if {$expgui(changed) != 0} {
[866]851            wm title . "EXPGUI interface to GSAS: $expgui(expfile) (modified)"
[16]852            set expgui(titleunchanged) 0
853        }
854    }
[133]855    if {$expgui(expModifiedLast) == 0} {
856        after 2000 afterawhile
857        return
858    }
[466]859    if {![file exists $expgui(expfile)]} {
860        after 2000 afterawhile
861        return
862    }
[133]863    if {[file mtime $expgui(expfile)] != $expgui(expModifiedLast)} {
[359]864        # we are "locked". Note that whenidle loop will be restarted later
865        if {$tcl_platform(platform) == "windows" && [file exists expgui.lck]} {
866            return
867        }
[133]868        set ans [ReloadExpMsg [file tail $expgui(expfile)] $expgui(changed)]
[359]869
[16]870        if {$ans == 0} {
[20]871            loadexp $expgui(expfile)
[16]872        } elseif {$ans == 1} {
873            # reset the time to the next version
[20]874            set expgui(expModifiedLast) [file mtime $expgui(expfile)]
[16]875        } elseif {$ans == 2} {
[323]876            SaveAsFile
[16]877        }
878    }
879    after 2000 afterawhile
880}
[133]881
882# place a message about changes over the main window
883proc ReloadExpMsg {file changes} {
884    global expgui tcl_platform
885    set msg "File $file has been modified by another program"
886    if {$changes == 1} {
887        append msg " and you have made a change to this version.\n"
888    } elseif {$changes > 0} {
889        append msg " and you have made $changes changes to this version.\n"
890    } else {
891        append msg ".\n"
892    }
[199]893    append msg "Do you want to use the newer (modified) version or continue with the older (previous) version of the file?"
894
[133]895    set w .ask
896    catch {destroy $w}
897    toplevel $w -class Dialog
898    wm title $w "Reload?"
899    wm iconname $w "Reload?"
900    wm protocol $w WM_DELETE_WINDOW { }
901    wm transient $w .
[323]902    bind $w <Key-F1> "MakeWWWHelp expguierr.html Overwrite"
903    pack [button $w.help -text Help -bg yellow \
904            -command "MakeWWWHelp expguierr.html Overwrite"] \
905            -side top -anchor e
[133]906    frame $w.bot
907    pack $w.bot -side bottom
[419]908    frame $w.top -class FixedFont
[133]909    pack $w.top -side top -fill both -expand 1
910    label $w.top.msg -justify left \
911            -wraplength 5i -font {Times 18} \
912            -text $msg 
913    if {$tcl_platform(platform) == "windows"} {
914        $w.top.msg config -font {Times 14}
915    }
916    pack $w.top.msg  -side right -expand 1 -fill both -padx 3m -pady 3m
917    pack [button $w.bot.1 -text "Load new" \
918            -default active -command "set expgui(dialogbutton) 0" \
919            ] -side left -expand 1 -padx 3m -pady 2m
[199]920    pack [button $w.bot.2 -text "Continue with old" \
[133]921            -command "set expgui(dialogbutton) 1"] \
922            -side left -expand 1 -padx 3m -pady 2m
923    if {$changes > 0} {
924        pack [button $w.bot.3 -text "Save edited version" \
925            -command "set expgui(dialogbutton) 2"] \
926            -side left -expand 1 -padx 3m -pady 2m
927    }
928    # Create a binding for <Return> on the dialog
[746]929    bind $w <Return> "$w.bot.1 invoke"
[133]930    wm withdraw $w
931    update idletasks
932
[809]933    # compute widths with a catch, since this has been showing an error
934    # at least for one Mac user
935    set askwid 40
936    catch {set askwid [winfo reqwidth $w]}
937    set askhgt 80
938    catch {set askhgt [winfo reqheight $w]}
[908]939    catch {
[133]940        # for now, always center the message over the main window
941        # center the new window in the middle of the parent
942        set x [expr [winfo x .] + [winfo width .]/2 - \
[809]943                $askwid/2 - [winfo vrootx .]]
[133]944        set y [expr [winfo y .] + [winfo height .]/2 - \
[809]945                $askhgt/2 - [winfo vrooty .]]
[133]946        wm geom $w +$x+$y
[908]947    }
[133]948    wm deiconify $w
949
950    # Grab the focus
951    set oldFocus [focus]
952    set oldGrab [grab current $w]
953    if {[string compare $oldGrab ""]} {
954        set grabStatus [grab status $oldGrab]
955    }
[605]956    catch {grab $w}
[133]957    focus $w.bot.1
[359]958    # for windows rearrange window stacking
959    # -- Removed since this will normally happen after the GSAS program
960    # has finished
961    #if {$tcl_platform(platform) == "windows"} {
962        #lower .
963        #raise $w .
964    #}
[133]965    update idletasks
966
967    tkwait variable expgui(dialogbutton)
968    catch {focus $oldFocus}
969    destroy $w
970    if {[string compare $oldGrab ""]} {
[289]971        if {![string compare $grabStatus "global"]} {
[605]972            catch {grab -global $oldGrab}
[133]973        } else {
[605]974            catch {grab $oldGrab}
[133]975        }
976    }
[289]977    # for windows rearrange window stacking
[359]978    #if {$tcl_platform(platform) == "windows"} {
979        #raise .
980    #}
[133]981    return $expgui(dialogbutton)
982}
983
[10]984# --------  called to confirm before exiting
985proc catchQuit {} {
986    if {[confirmBeforeSave] == "Continue"} {
987        destroy .
988    }
989}
990# save the .EXP file before exiting?
991proc confirmBeforeSave {} {
[20]992    global expgui
[10]993    if !$expgui(changed) {
994        return "Continue"
995    }
996    set decision [tk_dialog .instrSaveData {Save .EXP changes} \
997            {You have made changes to the Experiment, but the changes are not saved. Select an option:} \
998            {} 0 "Save and Exit" "Exit without Save" "Cancel exit command"]
999    switch $decision {
[20]1000        0 { savearchiveexp;  return "Continue" }
1001        1 {                  return "Continue" }
1002        2 {                  return "Cancel"   }
[10]1003    }
1004}
1005
[133]1006# setup buttons for each phase on the phase page
[10]1007proc setphases {} {
1008    global expgui expmap
[152]1009    eval destroy [winfo children $expgui(phaseFrame).top.ps]
[255]1010    pack [label $expgui(phaseFrame).top.ps.0 -text Phase:] -side left
[539]1011    foreach num $expmap(phaselist) {
[10]1012        pack [button $expgui(phaseFrame).top.ps.$num -text $num \
[255]1013                -command "SelectOnePhase $num" -padx 1.5m] -side left
[10]1014    }
[255]1015    if {[file executable $expgui(exptool)] && \
1016            [llength $expmap(phaselist)]} {
[323]1017        pack [button $expgui(phaseFrame).top.ps.10 \
1018                -text "Replace" -command MakeReplacePhaseBox \
1019                ] -side left
[255]1020    }
[10]1021}
1022
1023# Procedure to respond to changes the phase.
1024#  This loads the "phases" widgets with data corresponding to the selected phase.
1025proc SelectOnePhase {num} {
[303]1026    global entryvar entrycmd entrybox expmap expgui
[539]1027    # if no phase has been selected, select the first one
1028    if {$num == ""} {set num [lindex $expmap(phaselist) 0]}
1029
[133]1030    set crsPhase {}
[255]1031    $expgui(atomxform) config -text "Xform Atoms" -state disabled
[133]1032    foreach n $expmap(phaselist) type $expmap(phasetype) {
[539]1033        if {$n == $num} {
[238]1034            catch {$expgui(phaseFrame).top.ps.$num config -relief sunken}
[133]1035            set crsPhase $num
1036            if {$type == 3} {
1037                set expgui(phasetype) "Magnetic\nOnly"
1038            } elseif {$type == 2} {
1039                set expgui(phasetype) "Magnetic\n& Nuclear"
1040            } elseif {$type == 4} {
1041                set expgui(phasetype) "Macromolecular"
[838]1042            } elseif {$type == 10} {
1043                set expgui(phasetype) "Pawley"
[133]1044            } else {
1045                set expgui(phasetype) ""
1046            }
[10]1047        } else { 
[238]1048            catch {$expgui(phaseFrame).top.ps.$n config -relief raised}
[10]1049        }
[20]1050    }
1051    # no phase is selected
[71]1052    if {$crsPhase == "" || [llength $expmap(phaselist)] == 0} {
[20]1053        # disable traces on entryvar
1054        set entrycmd(trace) 0
1055        set entrycmd(phasename) ""
1056        set entryvar(phasename) ""
[670]1057        foreach ent {a b c alpha beta gamma} {
1058            set entryvar($ent) ""
1059        }
1060        foreach ent {cellref celldamp} {
[20]1061            set entrycmd($ent) ""
1062            set entryvar($ent) ""
1063        }
[71]1064        set expgui(curPhase) {}
[20]1065        # enable traces on entryvar
1066        set entrycmd(trace) 1
[97]1067        $expgui(EditingAtoms) config -text ""
1068        DisplayAtom 0 0
1069        DisplayU 0 0
1070        DisplayRefFlags 0 0
1071        $expgui(atomlistbox) delete 0 end
[20]1072        return
1073    }
1074
[539]1075    # don't reload the last displayed phase
1076    if {$expgui(curPhase) == $crsPhase} return
[10]1077
1078    ##########################################################
[539]1079    # load and display a phase
[10]1080    ##########################################################
[539]1081    # disable traces on entryvar while loading
1082    set entrycmd(trace) 0
[10]1083    # phase title
[20]1084    set entrycmd(phasename) "phaseinfo $crsPhase name"
[10]1085    set entryvar(phasename) [phaseinfo $crsPhase name]
1086    # cell parameters & flags
[670]1087    foreach ent {a b c alpha beta gamma} {
1088        set entryvar($ent) [phaseinfo $crsPhase $ent]
1089    }
1090    foreach ent {cellref celldamp} {
[10]1091        set entrycmd($ent) "phaseinfo $crsPhase $ent"
1092        set entryvar($ent) [phaseinfo $crsPhase $ent]
1093    }
1094
[539]1095    # initialize atoms display & disable
[10]1096    DisplayAtom 0 0
1097    DisplayU 0 0
1098    DisplayRefFlags 0 0
1099    $expgui(EditingAtoms) config -text ""
1100
[539]1101    DisplayAllAtoms $crsPhase
[10]1102
1103    # enable traces on entryvar now
1104    set entrycmd(trace) 1
1105}
1106
1107set expgui(noreenterDisplayAllAtoms) 0
[539]1108# Populate expgui(atomlistbox) (a ScrolledListBox) with atoms
1109# from the selected phase.
1110proc DisplayAllAtoms {curPhase "mode reset"} {
[10]1111    global entryvar entrycmd expmap expgui
[838]1112    # make sure that atomlistboxcontents element exists
1113    if {[catch {set expmap(atomlistboxcontents)}]} {
1114        set expmap(atomlistboxcontents) {}
1115    }
[539]1116    # if it does not show, we don't have a phase or we are already displaying
1117    # don't bother
[10]1118    if {$expgui(pagenow) != "phaseFrame"} return
[539]1119    if {$curPhase == ""} return
[10]1120    if $expgui(noreenterDisplayAllAtoms) return
[539]1121    # prevent reentry
[10]1122    set expgui(noreenterDisplayAllAtoms) 1
[539]1123    # set the current phase
1124    set expgui(curPhase) $curPhase
[10]1125    if {$mode != "reset"} {
[63]1126        # save the scrolled position
[10]1127        set pos [lindex [$expgui(atomlistbox) yview] 0]
[63]1128    } else {
[539]1129        # for reset, do not keep the previously selected atoms
[63]1130        set expgui(selectedatomlist) {}
[10]1131    }
1132    $expgui(atomlistbox) delete 0 end
[539]1133    # displaying a macromolecular phase?
1134    if {[lindex $expmap(phasetype) [expr {$expgui(curPhase) - 1}]] == 4} {
1135        set mm 1
1136        $expgui(phaseFrame).top.ps.10 config -state disabled
1137        $expgui(AddAtomBut)  config -state disabled
1138        pleasewait "loading atoms..."
1139    } else {
1140        set mm 0
1141        if {[file executable $expgui(exptool)]} {
1142            $expgui(phaseFrame).top.ps.10 config -state normal
1143            $expgui(AddAtomBut) config -state normal
1144        }
1145    }
1146
1147    # prepare header info
[10]1148    set maxline I
1149    set phase $expgui(curPhase)
1150    set atomlist {}
[234]1151    set typehead "type  "
1152    set namehead "  name  "
1153    set multhead "Mult"
1154    set coordhead "   "
[539]1155    if {$mm} {
1156        set cmd mmatominfo
1157        set frachead "Occ."
1158    } else {
1159        set cmd atominfo
1160        set frachead "Occupancy"
1161    }
1162    set reshead "res/grp/#"
1163    # sort the atoms, as requested
[10]1164    if  {$expgui(asorttype) == "type"} {
1165        # sort on atom type
[234]1166        set typehead "type* "
[10]1167        foreach atom $expmap(atomlist_$phase) {
[539]1168            lappend atomlist "$atom [$cmd $phase $atom type] $phase"
[10]1169        }
1170        set expmap(atomlistboxcontents) [lsort -ascii -index 1 $atomlist]
1171    } elseif {$expgui(asorttype) == "number"} {
1172        # sort on atom number
[234]1173        set namehead "* name  "
[10]1174        foreach atom $expmap(atomlist_$phase) {
1175            lappend atomlist "$atom $atom $phase"
1176        }
1177        set expmap(atomlistboxcontents) [lsort -integer -index 1 $atomlist]
[234]1178    } elseif {$expgui(asorttype) == "mult"} {
[539]1179        if {$mm} {
1180            set reshead "res*/grp/#"
1181            foreach atom $expmap(atomlist_$phase) {
1182                lappend atomlist "$atom [mmatominfo $phase $atom residue] $phase"
1183            }
1184            set expmap(atomlistboxcontents) [lsort -ascii -index 1 $atomlist]
1185        } else {
1186            # sort on atom number
1187            set multhead "Mlt*"
1188            foreach atom $expmap(atomlist_$phase) {
1189                lappend atomlist "$atom [atominfo $phase $atom mult] $phase"
1190            }
1191            set expmap(atomlistboxcontents) [lsort -integer -decreasing -index 1 $atomlist]
[234]1192        }
[255]1193    } elseif {$expgui(asorttype) == "occupancy"} {
1194        # sort on atom number
[539]1195        if {$mm} {
1196            set frachead "  Occ* "
1197        } else {
1198            set frachead "  Occup* "
1199        }
[255]1200        foreach atom $expmap(atomlist_$phase) {
[539]1201            lappend atomlist "$atom [$cmd $phase $atom frac] $phase"
[255]1202        }
1203        set expmap(atomlistboxcontents) [lsort -real -decreasing -index 1 $atomlist]
[10]1204    } elseif {$expgui(asorttype) == "x"} {
1205        # sort on x
[234]1206        set coordhead "(x*)"
[10]1207        foreach atom $expmap(atomlist_$phase) {
[539]1208            lappend atomlist "$atom [$cmd $phase $atom x] $phase"
[10]1209        }
1210        set expmap(atomlistboxcontents) [lsort -real -index 1 $atomlist]
1211    } elseif {$expgui(asorttype) == "y"} {
1212        # sort on y
[234]1213        set coordhead "(y*)"
[10]1214        foreach atom $expmap(atomlist_$phase) {
[539]1215            lappend atomlist "$atom [$cmd $phase $atom y] $phase"
[10]1216        }
1217        set expmap(atomlistboxcontents) [lsort -real -index 1 $atomlist]
1218    } elseif {$expgui(asorttype) == "z"} {
1219        # sort on z
[234]1220        set coordhead "(z*)"
[10]1221        foreach atom $expmap(atomlist_$phase) {
[539]1222            lappend atomlist "$atom [$cmd $phase $atom z] $phase"
[10]1223        }
1224        set expmap(atomlistboxcontents) [lsort -real -index 1 $atomlist]
1225    } else {
[234]1226        error "Bad expgui(asorttype) = $expgui(asorttype)"
[10]1227    }
1228
[539]1229    set expgui(atomlistboxline) {}
1230    # loop over atoms
[10]1231    foreach tuple $expmap(atomlistboxcontents) {
1232        set atom [lindex $tuple 0]
1233        set phase [lindex $tuple 2]
[539]1234        lappend expgui(atomlistboxline) $atom
1235        $expgui(atomlistbox) insert end \
1236                [FormatAtomLine $atom $phase maxline]
[10]1237    }
1238    $expgui(atomtitle) delete 0 end
[539]1239
1240    # create the header
1241    if {$mm} {
1242        $expgui(atomtitle) insert end [format "%12s %9s  %6s %8s%29s %4s  %s" \
1243                $namehead $reshead $typehead "ref/damp  " \
1244                "fractional coordinates$coordhead" \
1245                "$frachead" \
1246                " Uiso"]
1247        donewait
1248    } elseif {$maxline == "A"} {
[234]1249        $expgui(atomtitle) insert end [format "%10s %6s %8s%29s %9s  %s" \
1250                $namehead $typehead "ref/damp  " \
1251                "fractional coordinates$coordhead" \
[255]1252                "$multhead $frachead" \
1253                " Uiso/Uij                                            "]
[10]1254    } else {
[234]1255        $expgui(atomtitle) insert end [format "%10s %6s %8s%29s %9s  %s" \
1256                $namehead $typehead "ref/damp  " \
1257                "fractional coordinates$coordhead" \
[255]1258                "$multhead $frachead" \
1259                " Uiso"]
[10]1260    }
1261    if {$mode != "reset"} {
[63]1262        # restore the selected items
1263        foreach i $expgui(selectedatomlist) {
[10]1264            $expgui(atomlistbox) selection set $i
1265        }
[63]1266        # restore the last scrolled position
1267        $expgui(atomlistbox) yview moveto $pos
[10]1268    }
[539]1269    # clear the reentry flag
[10]1270    set expgui(noreenterDisplayAllAtoms) 0
1271}
1272
[539]1273proc FormatAtomLine {atom phase maxline_var} {
1274    global expmap
1275    if {[lindex $expmap(phasetype) [expr {$phase - 1}]] == 4} {
1276        foreach type {x u f} {
1277            if {[mmatominfo $phase $atom ${type}ref]} {
1278                append refflag "[string toupper $type][mmatominfo $phase $atom ${type}damp] "
1279            } else {
1280                append refflag " [mmatominfo $phase $atom ${type}damp] "
1281            }   
1282        }
1283        set line [format \
1284                "%5d %-6s %-3s%-2s%4d  %-6s %8s %9.5f%9.5f%9.5f%8.4f %7.4f" \
1285                $atom \
1286                [mmatominfo $phase $atom label] \
1287                [mmatominfo $phase $atom residue] \
1288                [mmatominfo $phase $atom group] \
1289                [mmatominfo $phase $atom resnum] \
1290                [mmatominfo $phase $atom type] \
1291                $refflag \
1292                [mmatominfo $phase $atom x] \
1293                [mmatominfo $phase $atom y] \
1294                [mmatominfo $phase $atom z] \
1295                [mmatominfo $phase $atom frac] \
1296                [mmatominfo $phase $atom Uiso] 
1297        ]
1298    } elseif {[atominfo $phase $atom temptype] == "A"} {
1299        foreach type {x u f} {
1300            if {[atominfo $phase $atom ${type}ref]} {
1301                append refflag "[string toupper $type][atominfo $phase $atom ${type}damp] "
1302            } else {
1303                append refflag " [atominfo $phase $atom ${type}damp] "
1304            }
1305        }
1306        # want to set maxline in parent
1307        upvar $maxline_var maxline
1308        set maxline A
1309        # aniso
1310        set line [format "%3d %-6s %-6s %8s %10.6f%10.6f%10.6f%4d%9.4f" \
1311                $atom \
1312                [atominfo $phase $atom label] \
1313                [atominfo $phase $atom type] \
1314                $refflag \
1315                [atominfo $phase $atom x] \
1316                [atominfo $phase $atom y] \
1317                [atominfo $phase $atom z] \
1318                [atominfo $phase $atom mult] \
1319                [atominfo $phase $atom frac] 
1320        ]
1321        append line [format "  %9.5f%9.5f%9.5f%9.5f%9.5f%9.5f" \
1322                [atominfo $phase $atom U11] \
1323                [atominfo $phase $atom U22] \
1324                [atominfo $phase $atom U33] \
1325                [atominfo $phase $atom U12] \
1326                [atominfo $phase $atom U23] \
1327                [atominfo $phase $atom U13] 
1328        ]
1329    } else {
1330        foreach type {x u f} {
1331            if {[atominfo $phase $atom ${type}ref]} {
1332                append refflag "[string toupper $type][atominfo $phase $atom ${type}damp] "
1333            } else {
1334                append refflag " [atominfo $phase $atom ${type}damp] "
1335            }   
1336        }
1337        set line [format \
1338                "%3d %-6s %-6s %8s %10.6f%10.6f%10.6f%4d%9.4f  %9.5f" \
1339                $atom \
1340                [atominfo $phase $atom label] \
1341                [atominfo $phase $atom type] \
1342                $refflag \
1343                [atominfo $phase $atom x] \
1344                [atominfo $phase $atom y] \
1345                [atominfo $phase $atom z] \
1346                [atominfo $phase $atom mult] \
1347                [atominfo $phase $atom frac] \
1348                [atominfo $phase $atom Uiso] 
1349        ]
1350    }
1351    return $line
1352}
1353
1354# update the display of atom as they are changed
1355proc UpdateAtomLine {atomlist phase} {
1356    global expgui
1357    # for lots of atoms, it is faster to repaint the listbox
1358    if {[llength $atomlist] > 25} {
1359        DisplayAllAtoms $expgui(curPhase) noreset
1360        return
1361    }
1362    foreach atom $atomlist {
1363        set linenum [lsearch -exact $expgui(atomlistboxline) $atom]
1364        $expgui(atomlistbox) delete $linenum
1365        $expgui(atomlistbox) insert $linenum \
1366                [FormatAtomLine $atom $phase maxline]
1367    }
1368    # restore the selected items
1369    foreach i $expgui(selectedatomlist) {
1370        $expgui(atomlistbox) selection set $i
1371    }
1372}
1373
[10]1374# Procedure to select all atoms in response to a right-click
1375proc SelectAllAtoms {} {
1376    global expgui
1377    $expgui(atomlistbox) selection set 0 end
[63]1378    # call editRecord in case trace was called before the selection was made
[10]1379    editRecord
1380}
1381
1382# Procedure to respond to left mouse release in the atoms Pane
1383proc editRecord { args } {
1384    global entrycmd expgui
[63]1385    set expgui(selectedatomlist) [$expgui(atomlistbox) curselection]
[10]1386    # disable traces on entryvar for right now
1387    set entrycmd(trace) 0
1388
[63]1389    if {[llength $expgui(selectedatomlist)] == 0} {
[255]1390        if $expgui(debug) {error "Attempt display non-existent atoms"}
[63]1391    } elseif {[llength $expgui(selectedatomlist)] == 1} {
1392        editOneRecord $expgui(selectedatomlist)
[10]1393    } else {
[63]1394        editMultipleRecords $expgui(selectedatomlist)
[10]1395    }
1396    # reenable traces on entryvar
1397    set entrycmd(trace) 1
1398    # repaint the atoms box in case anything was changed
[63]1399    #    DisplayAllAtoms noreset
[10]1400}
1401
1402proc editOneRecord { AtomIndex } {
1403    global expmap expgui
[838]1404    # make sure that atomlistboxcontents element exists
1405    if {[catch {set expmap(atomlistboxcontents)}]} return
1406
[10]1407    # get atom number & phase
1408    set tuple [lindex $expmap(atomlistboxcontents) $AtomIndex]
1409    set atomnum [lindex $tuple 0]
1410    set p [lindex $tuple 2]
1411    DisplayU $atomnum $p
1412    DisplayAtom $atomnum $p
1413    DisplayRefFlags $atomnum $p
1414    $expgui(EditingAtoms) config -text "Editing atom #$atomnum -- [atominfo $p $atomnum label]"
[255]1415    $expgui(atomxform) config -text "Xform Atom" -state normal
[10]1416}
1417
1418# this will not work for a multi-phase list of atoms (yet)
1419proc editMultipleRecords { AtomIndexList } {
1420    global expmap expgui
[838]1421    # make sure that atomlistboxcontents element exists
1422    if {[catch {set expmap(atomlistboxcontents)}]} return
1423
[10]1424    set numberList {}
1425    # current phase
1426    set p $expgui(curPhase)
1427    foreach AtomIndex $AtomIndexList {
1428        # get atom number & phase
1429        set tuple [lindex $expmap(atomlistboxcontents) $AtomIndex]
1430        lappend numberList [lindex $tuple 0]
1431#       set p [lindex $tuple 2]
1432    }
1433    # this needs to track by phase
1434    $expgui(EditingAtoms) config -text \
1435            "Set refinement options: atoms [CompressList $numberList]"
1436    DisplayU 0 0
1437    DisplayAtom 0 0
1438    # this needs to track by phase
1439    DisplayRefFlags $numberList $p
[255]1440    $expgui(atomxform) config -text "Xform Atoms" -state normal
[10]1441}
1442
1443# format a string of numbers to save space, e.g. "1 2 3 4 6 7 19 13 14 15"
1444# becomes "1-4,6,7,13-15,19"
[539]1445proc CompressList {numberList "max 9999"} {
[10]1446    # format the number list to save space
1447    set lastnum -99
1448    set flist {}
1449    set count 0
[539]1450    set length 0
[199]1451    if [catch {set sortlist [lsort -integer $numberList]}] {return $numberList}
1452    foreach num $sortlist {
[539]1453        set next [expr $lastnum+1]
1454        if {$num != $next} {
1455            if {$count == 0 && $flist != ""} {
1456                if {[string length $flist] - $length > $max} {
1457                    set length [string length $flist]
1458                    append flist ",\n$num"
1459                } else {
[199]1460                    append flist ",$num"
[539]1461                }
1462            } elseif {$count == 1 && $flist != ""} {
1463                if {[string length $flist] - $length > $max} {
1464                    set length [string length $flist]
1465                    append flist ",$lastnum,\n$num"
1466                } else {
[199]1467                    append flist ",$lastnum,$num"
[539]1468                }
1469            } elseif {$flist != ""} {
1470                if {[string length $flist] - $length > $max} {
1471                    set length [string length $flist]
1472                    append flist "-$lastnum,\n$num"
1473                } else {
[199]1474                    append flist "-$lastnum,$num"
[539]1475                }
[10]1476            } else {
1477                append flist "$num"
1478            }
1479            set lastnum $num
1480            set count 0
1481        } else {
1482            incr count
1483            incr lastnum
1484        }
1485    }
1486    if {$count == 1 && $flist != ""} {
1487        append flist ",$lastnum"
1488    } elseif {$flist != "" && $count > 1} {
1489        append flist "-$lastnum"
1490    }
1491    return $flist
1492}
1493
1494# Procedure to display Isotropic or Anisotropic temperature factors
1495#  Changes the display to one entry widget for Isotropic motion OR
1496#   6 entry widgets for Anisotropic motion in Frame3.
1497#   or disables the widet entirly if atom = 0
1498proc DisplayU { atomnum p} {
[539]1499    global expgui entryvar entrycmd expmap
1500    set mm 0
[10]1501    if {$atomnum == 0} {
1502        set iOrA disable
[539]1503    } elseif {[lindex $expmap(phasetype) 0] == 4} {
1504        set mm 1
1505        set iOrA I
[10]1506    } else {
1507        set iOrA [atominfo $p $atomnum temptype]
1508    }
1509
1510    set firstbox [lindex $expgui(anisolabels) 0]
1511    if { $iOrA == "A" } {
1512        $firstbox config -text "U11 "
1513        foreach item $expgui(anisolabels) {
1514            $item config -fg black
1515        }
1516        foreach item $expgui(anisoentry) var {U11 U22 U33 U12 U13 U23} {
1517            set entrycmd($var) "atominfo $p $atomnum $var"
1518            set entryvar($var) [eval $entrycmd($var)]
1519            $item config -fg black -state normal  -bg white
1520        }
1521    } elseif { $iOrA == "I" || $iOrA == "disable"} {
1522        foreach item $expgui(anisolabels) {
[306]1523            $item config -fg $expgui(bkgcolor1)
[10]1524        }
1525        foreach item [lrange $expgui(anisoentry) 1 end] \
1526                var {U22 U33 U12 U13 U23} {
1527            set entrycmd($var) ""
1528            set entryvar($var) ""
[306]1529            $item config -fg $expgui(bkgcolor1) -bg $expgui(bkgcolor1) \
1530                    -state disabled
[10]1531        }
1532        if { $iOrA == "disable"} {
1533            set entrycmd($var) ""
1534            set entryvar($var) ""
[306]1535            [lindex $expgui(anisoentry) 0] config \
1536                    -fg $expgui(bkgcolor1) -bg $expgui(bkgcolor1) \
1537                    -state disabled
[539]1538        } elseif {$mm} {
1539            set entrycmd(U11) "mmatominfo $p $atomnum Uiso"
1540            set entryvar(U11) [eval $entrycmd(U11)]
1541            $firstbox config -text Uiso -fg black
1542            [lindex $expgui(anisoentry) 0] config -fg black -bg white -state normal
[10]1543        } else {
1544            set entrycmd(U11) "atominfo $p $atomnum Uiso"
1545            set entryvar(U11) [eval $entrycmd(U11)]
1546            $firstbox config -text Uiso -fg black
1547            [lindex $expgui(anisoentry) 0] config -fg black -bg white -state normal
1548        }
1549    }
1550}
1551
1552# need to think about multiple phases
1553
1554# Procedure to display refinement flags
1555proc DisplayRefFlags { atomnum p} {
[539]1556    global expgui entryvar entrycmd expmap
[10]1557    if {$atomnum == 0} {
1558        foreach label $expgui(atomreflbl) {
[306]1559            $label config -fg $expgui(bkgcolor1)
[10]1560        }
1561        foreach entry $expgui(atomref) {
[306]1562            $entry config -state disabled \
1563                    -fg $expgui(bkgcolor1) -bg $expgui(bkgcolor1)
1564            # turn off checkbuttons
1565            catch {$entry deselect}
1566
[10]1567        }
1568        return
1569    }
1570    foreach label $expgui(atomreflbl) {
1571        $label config -fg black
1572    }
1573    foreach entry $expgui(atomref) {
[306]1574        $entry config -state normal -fg black -bg $expgui(bkgcolor1)
[10]1575    }
[539]1576    if {[lindex $expmap(phasetype) 0] == 4} {
1577        foreach var {xref uref fref xdamp udamp fdamp}  {
1578            set entrycmd($var) "mmatominfo $p [list $atomnum] $var"
1579            set entryvar($var) [eval $entrycmd($var)]
1580        }
1581    } else {
1582        foreach var {xref uref fref xdamp udamp fdamp}  {
1583            set entrycmd($var) "atominfo $p [list $atomnum] $var"
1584            set entryvar($var) [eval $entrycmd($var)]
1585        }
[10]1586    }
1587}
1588
1589# Procedure to display an atom in the atom edit boxes
1590proc DisplayAtom { atomnum p} {
[539]1591    global expgui entryvar entrycmd expmap
[10]1592    if {$atomnum == 0} {
1593        foreach label $expgui(atomlabels) {
[306]1594            $label config -fg $expgui(bkgcolor1)
[10]1595        }
1596        foreach entry $expgui(atomentry) {
[306]1597            $entry config -state disabled \
1598                    -fg $expgui(bkgcolor1) -bg $expgui(bkgcolor1)
[10]1599        }
1600        return
1601    }
1602    foreach label $expgui(atomlabels) {
1603        $label config -fg black
1604    }
1605    foreach entry $expgui(atomentry) {
1606        $entry config -state normal -fg black -bg white
1607    }
[539]1608    if {[lindex $expmap(phasetype) 0] == 4} {
1609        foreach var {x y z label frac } {
1610            set entrycmd($var) "mmatominfo $p $atomnum $var"
1611            set entryvar($var) [eval $entrycmd($var)]
1612        }
1613    } else {
1614        foreach var {x y z label frac } {
1615            set entrycmd($var) "atominfo $p $atomnum $var"
1616            set entryvar($var) [eval $entrycmd($var)]
1617        }
[10]1618    }
1619}
1620
[234]1621# make a histogram box; used in MakeHistPane,
1622proc MakeHistBox {frm} {
1623    global expgui
1624    grid [label $frm.mode -text "Select a Histogram" \
[866]1625            -bg beige -anchor center -bd 2 -relief raised] \
[234]1626            -row 0 -column 0 -columnspan 2 -sticky ew
1627    bind $frm.mode <Button-1> {
1628        set i [lsearch $expgui(AllowedHistSelectModes) $expgui(globalmode)]
1629        set expgui(globalmode) [lindex \
1630                "$expgui(AllowedHistSelectModes) \
1631                $expgui(AllowedHistSelectModes)" [incr i]]
1632        sethistlist
1633    }
1634    bind $frm.mode <Button-3> {set expgui(globalmode) 0; sethistlist}
1635    grid [listbox $frm.title -height 1 -relief flat \
[419]1636            -exportselection 0 ] -row 1 -column 0 -sticky ew
[234]1637    grid [listbox $frm.lbox -height 10 -width 25 \
1638            -exportselection 0 \
1639            -xscrollcommand "$frm.x set" \
1640            -yscrollcommand "$frm.y set" \
1641            ] -row 2 -column 0 -sticky news
1642    lappend expgui(HistSelectList) $frm
1643    grid [scrollbar $frm.x -orient horizontal \
[449]1644            -command "move2boxesX \" $frm.title $frm.lbox \" " 
[234]1645    ] -row 3 -column 0 -sticky ew
1646    grid [scrollbar $frm.y \
1647            -command "$frm.lbox yview"] \
1648            -row 2 -column 1 -sticky ns
1649    grid columnconfigure $frm 0 -weight 1
1650    grid rowconfigure $frm 2 -weight 1
1651}
1652
[10]1653# update the histogram list
1654# to do: show histogram ref flags?
1655proc sethistlist {} {
1656    global expgui expmap
[234]1657    array set lbl {
1658        1 "Select 1 or more\nTOF Histograms"
1659        2 "Select 1 or more\nCW Neutron Histograms"
1660        3 "Select 1 or more\nAlpha 1,2 X-ray Histograms"
1661        4 "Select 1 or more\nmonochromatic X-ray Histograms"
[478]1662        5 "Select 1 or more Energy\nDispersive X-ray Histograms"
[234]1663        6 "Select 1 or more of\n any type Histograms"
1664    }
[10]1665    foreach lbox $expgui(HistSelectList) {
1666        $lbox.title delete 0 end
1667        $lbox.lbox delete 0 end
1668        if {$expgui(globalmode) != 0} {
1669            $lbox.lbox config -selectmode extended
[234]1670            $lbox.mode config -text $lbl($expgui(globalmode)) -bg yellow
[10]1671        } else {
1672            $lbox.lbox config -selectmode browse
[234]1673            $lbox.mode config -text "Select a histogram" -bg beige
[10]1674        }
1675    }
1676    # disable the unallowed pages in all mode
1677    if {$expgui(globalmode) == 6} {
1678        foreach pair $expgui(GlobalModeAllDisable) {
1679            if {$expgui(pagenow) == [lindex $pair 0]} {
[118]1680                RaisePage lsFrame
[10]1681            }
1682            eval [lindex $pair 1] -state disabled
1683        }
1684    } else {
1685        foreach pair $expgui(GlobalModeAllDisable) {
1686            eval [lindex $pair 1] -state normal
1687        }
1688    }
[20]1689    set histlist {}
[10]1690    if  {$expgui(hsorttype) == "type"} {
1691        # sort on histogram type
1692        foreach h [lsort -integer -increasing $expmap(powderlist)] {
1693            lappend histlist "$h [string range $expmap(htype_$h) 1 2]"
1694        }
1695        set expmap(histlistboxcontents) [lsort -ascii -index 1 $histlist]
1696    } elseif {$expgui(hsorttype) == "number"} {
1697        # sort on histogram number
1698        foreach h [lsort -integer -increasing $expmap(powderlist)] {
1699            lappend histlist "$h $h"
1700        }
1701        set expmap(histlistboxcontents) [lsort -integer -index 1 $histlist]
1702    } elseif {$expgui(hsorttype) == "bank"} {
1703        # sort on original bank number
1704        foreach h [lsort -integer -increasing $expmap(powderlist)] {
1705            lappend histlist "$h [histinfo $h bank]"
1706        }
1707        set expmap(histlistboxcontents) [lsort -integer -index 1 $histlist]
1708    } elseif {$expgui(hsorttype) == "angle"} {
1709        # sort on wavelength (CW) or angle (E disp.)
1710        foreach h [lsort -integer -increasing $expmap(powderlist)] {
1711            if {[string range $expmap(htype_$h) 2 2] == "T"} {
1712                set det [format %8.2f [histinfo $h tofangle]]
1713            } elseif {[string range $expmap(htype_$h) 2 2] == "C"} {
1714                set det [format %8.5f [histinfo $h lam1]]
1715            } elseif {[string range $expmap(htype_$h) 2 2] == "E"} {
1716                set det [format %8.2f [histinfo $h lam1]]
1717            } else {
1718                set det {}
1719            }
1720            lappend histlist "$h $det"
1721        }
1722        set expmap(histlistboxcontents) [lsort -real -index 1 $histlist]
1723    }
[20]1724
[10]1725    # title field needs to match longest title
1726    foreach lbox $expgui(HistSelectList) {
1727        $lbox.title insert end [format "%2s %s %4s %8s  %-67s" \
1728                "h#" \
1729                type \
1730                bank \
1731                "ang/wave" \
1732                "    title" \
1733                ]
1734    }
1735    foreach tuple $expmap(histlistboxcontents) {
1736        set h [lindex $tuple 0]
1737
1738        if {$expgui(globalmode) == 1} {
1739            if {[string range $expmap(htype_$h) 2 2] != "T"} continue
1740        } elseif {$expgui(globalmode) == 2} {
1741            if {[string range $expmap(htype_$h) 1 2] != "NC"} continue
1742        } elseif {$expgui(globalmode) == 3} {
1743            if {[string range $expmap(htype_$h) 1 2] != "XC" || \
1744                    [histinfo $h lam2] == 0.0} continue
1745        } elseif {$expgui(globalmode) == 4} {
1746            if {[string range $expmap(htype_$h) 1 2] != "XC" || \
1747                    [histinfo $h lam2] != 0.0} continue
1748        } elseif {$expgui(globalmode) == 5} {
1749            if {[string range $expmap(htype_$h) 1 2] != "XE"} continue
1750        }
1751
1752        if {[string range $expmap(htype_$h) 2 2] == "T"} {
1753            set det [format %8.2f [histinfo $h tofangle]]
1754        } elseif {[string range $expmap(htype_$h) 2 2] == "C"} {
1755            set det [format %8.5f [histinfo $h lam1]]
1756        } elseif {[string range $expmap(htype_$h) 2 2] == "E"} {
1757            set det [format %8.2f [histinfo $h lam1]]
1758        } else {
1759            set det {}
1760        }
1761        foreach lbox $expgui(HistSelectList) {
[449]1762            $lbox.lbox insert end [format "%2d  %s %4d %8s  %-67s" \
[10]1763                    $h \
[449]1764                    [string range $expmap(htype_$h) 1 3] \
[10]1765                    [histinfo $h bank] \
1766                    $det \
1767                    [string range [histinfo $h title] 0 66] \
1768                    ]
1769        }
1770    }
[238]1771    UpdateCurrentPage
1772}
1773
1774proc UpdateCurrentPage {} {
1775    global expgui
[10]1776    foreach set $expgui(frameactionlist) {
[238]1777        if {$expgui(pagenow) == [lindex $set 0]} {catch [lindex $set 1]}
[10]1778    }
1779}
1780
1781#-----------------------------------------------------------------------
1782# ----------- draw Histogram page
1783#-----------------------------------------------------------------------
1784proc DisplayHistogram {} {
[303]1785    global expgui entrycmd entryvar entrybox expmap
[10]1786
[112]1787    # trap if more than one histogram is selected unless global mode
1788    if {$expgui(globalmode) == 0 && [llength $expgui(curhist)] > 1} {
1789        set expgui(curhist) [lindex $expgui(curhist) 0] 
1790    }
1791
[327]1792    # disable the add histogram button if no phases are present
[390]1793    catch {
[457]1794        foreach c [winfo children $expgui(histFrame).bb] {
1795            if {[llength $expmap(phaselist)] == 0} {
1796                $c configure -state disabled
1797            } else {
1798                $c configure -state normal
1799            } 
[390]1800        }
[327]1801    }
1802
[10]1803    # display the selected histograms
1804    $expgui(histFrame).hs.lbox selection clear 0 end
1805    foreach h $expgui(curhist) {
1806        $expgui(histFrame).hs.lbox selection set $h
1807    }
1808
1809    # disable traces on entryvar for right now
1810    set entrycmd(trace) 0
1811
1812    # get histogram list
1813    set histlist {}
1814    foreach item $expgui(curhist) {
1815        lappend histlist [lindex $expmap(powderlist) $item]
1816    }
[20]1817    # must have at least one histogram selected here
1818    if {[llength $histlist] == 0} {
1819        set expgui(backtermlbl) ""
1820        set expgui(backtypelbl) ""
[557]1821        foreach var {bref bdamp absref absdamp} {
[20]1822            set entrycmd($var) ""
1823            set entryvar($var) ""
1824        }
1825        $expgui(histFrame).top.txt config -text "No Selected Histograms"
1826        grid $expgui(histFrame).top -column 1 -row 0 -sticky nsew       
1827        set expgui(bkglbl) ""
[557]1828        set expgui(abslbl) ""
[152]1829        eval destroy [winfo children $expgui(diffBox)]
[20]1830        set entrycmd(trace) 1
1831        return
1832    }
1833
[10]1834    if {$expgui(globalmode) != 0} {
1835        set expgui(backtermlbl) ""
1836        set expgui(backtypelbl) ""
[557]1837        foreach var {bref bdamp absref absdamp} {
[10]1838            set entrycmd($var) "histinfo [list $histlist] $var"
1839            set entryvar($var) [histinfo [lindex $histlist 0] $var]
1840        }
1841    } else {
1842        set hist $histlist
1843        set terms [histinfo $hist backterms]
1844        set expgui(backtermlbl) "($terms terms)"
[720]1845        # background type 3 & 9 have gone away
1846        if {[histinfo $hist backtype] == 3 || [histinfo $hist backtype] == 9} {
1847            MyMessageBox -parent . -title "Background Change" \
1848                -type ok -default ok \
1849                -icon warning \
1850                -message "Background function #[histinfo $hist backtype] is no longer supported -- the function will now be changed to type #1 & the values reset"
1851            histinfo $histlist backtype set 1
[908]1852            RecordMacroEntry "histinfo [list $histlist] backtype set 1" 0
[720]1853            incr expgui(changed)
1854            for {set num 1 } { $num <= $terms } { incr num } {
1855                set var "bterm$num"
1856                histinfo $histlist $var set 0
[908]1857                RecordMacroEntry "histinfo [list $histlist] $var set 0" 0
[720]1858                incr expgui(changed)
1859            }
[908]1860            RecordMacroEntry "incr expgui(changed)" 0
[720]1861        }
[10]1862        set expgui(backtypelbl) "Function type [histinfo $hist backtype]"
[557]1863        foreach var {bref bdamp absref absdamp} {
[10]1864            set entrycmd($var) "histinfo $hist $var"
1865            set entryvar($var) [eval $entrycmd($var)]
1866        }
1867    }
[14]1868    # Top box
[823]1869    catch {destroy $expgui(histFrame).pflag}
[10]1870    if {$expgui(globalmode) != 0} {
1871        $expgui(histFrame).top.txt config \
1872                -text "Selected Histograms: [CompressList $histlist]"
1873        grid $expgui(histFrame).top -column 1 -row 0 -sticky nsew       
1874        set expgui(bkglbl) "Globally Edit Background"
[557]1875        set expgui(abslbl) "Globally Edit Absorption"
[10]1876    } else {
1877        grid forget $expgui(histFrame).top
1878        set expgui(bkglbl) "Edit Background"
[557]1879        set expgui(abslbl) "Edit Abs./Refl."
[823]1880        if {[llength $expmap(phaselist)] > 1} {
[250]1881            TitleFrame $expgui(histFrame).pflag  \
1882                    -borderwidth 4 -side left -relief groove \
1883                    -text "Phase Flags"
1884            set expgui(pflag) [$expgui(histFrame).pflag getframe]
1885            grid $expgui(histFrame).pflag -column 1 -row 1 -sticky nsew
1886            grid rowconfigure $expgui(histFrame) 2 -minsize 35
1887            foreach p $expmap(phaselist) {
1888                pack [checkbutton $expgui(pflag).$p \
1889                        -command "GetPhaseFlags $hist" \
1890                        -variable expgui(pflag$p) -text $p] -side left
1891                if {[lsearch $expmap(phaselist_$hist) $p] == -1} {
1892                    set expgui(pflag$p) 0
1893                } else {
1894                    set expgui(pflag$p) 1
1895                }
1896            }
1897        }
[10]1898    }
1899
1900    # diffractometer constants
1901    foreach var {lam1 lam2 kratio pola ipola ddamp zero \
1902            wref pref dcref daref ratref ttref zref } {
1903        set entrycmd($var) "histinfo [list $histlist] $var"
1904        set entryvar($var) [histinfo [lindex $histlist 0] $var]
1905    }
1906
[152]1907    eval destroy [winfo children $expgui(diffBox)]
[10]1908    if {$expgui(globalmode) == 0} {
1909        if {[string range $expmap(htype_$hist) 2 2] == "T"} {
1910        #------
1911        # TOF |
1912        #------
1913            grid [ label $expgui(diffBox).lDCrc -text "Refine DIFC" ] \
1914                    -column 1 -row 1
1915            grid [ checkbutton $expgui(diffBox).rfDCrc -variable entryvar(dcref) ] \
1916                    -column 2 -row 1
1917            grid [ label $expgui(diffBox).lDCdifc -text DIFC ] \
1918                    -column 3 -row 1 -sticky w
1919            grid [ entry $expgui(diffBox).eDCdifc -textvariable entryvar(lam1) \
1920                    -width 15 ] -column 4 -row 1
[303]1921            set entrybox(lam1) $expgui(diffBox).eDCdifc
[10]1922            #
1923            grid [ label $expgui(diffBox).lDCra -text "Refine DIFA" ] \
1924                    -column 1 -row 2
1925            grid [ checkbutton $expgui(diffBox).rfDCra -variable entryvar(daref) ] \
1926                    -column 2 -row 2
1927            grid [ label $expgui(diffBox).lDCdifa -text DIFA ] \
1928                    -column 3 -row 2
1929            grid [ entry $expgui(diffBox).eDCdifa -textvariable entryvar(lam2) \
1930                    -width 15 ] -column 4 -row 2
[303]1931            set entrybox(lam2) $expgui(diffBox).eDCdifa
[10]1932            #
1933            grid [ label $expgui(diffBox).lDCzero -text "Zero"] \
1934                    -column 3 -row 3
1935            grid [ entry $expgui(diffBox).eDCzero -textvariable entryvar(zero) \
1936                    -width 15 ] -column 4 -row 3
[303]1937            set entrybox(zero) $expgui(diffBox).eDCzero
[10]1938            grid [ label $expgui(diffBox).lDCzref -text "Refine zero" ] \
1939                    -column 1 -row 3 -sticky w
1940            grid [ checkbutton $expgui(diffBox).rfDCzref \
1941                    -variable entryvar(zref) ] -column 2 -row 3
1942        } elseif {[string range $expmap(htype_$hist) 1 2] == "NC"} { 
1943        #---------------
1944        # CW - neutron |
1945        #---------------
1946            grid [ label $expgui(diffBox).lDC1 -text "Refine wave" ] \
1947                    -column 1 -row 1
1948            grid [ checkbutton $expgui(diffBox).rfDC1 -variable entryvar(wref) ] \
1949                    -column 2 -row 1
1950            grid [ label $expgui(diffBox).lDCdifc -text wave ] \
1951                    -column 3 -row 1 -sticky w
1952            grid [ entry $expgui(diffBox).eDCdifc -textvariable entryvar(lam1) \
1953                    -width 15 ] -column 4 -row 1
[303]1954            set entrybox(lam1) $expgui(diffBox).eDCdifc
[10]1955            #
1956            grid [ label $expgui(diffBox).lDCzref -text "Refine zero" ] \
1957                    -column 1 -row 3 -sticky w
1958            grid [ checkbutton $expgui(diffBox).rfDCzref \
1959                    -variable entryvar(zref) ] -column 2 -row 3
1960            grid [ label $expgui(diffBox).lDCzero -text "Zero"] \
1961                    -column 3 -row 3
1962            grid [ entry $expgui(diffBox).eDCzero -textvariable entryvar(zero) \
1963                    -width 15 ] -column 4 -row 3
[303]1964            set entrybox(zero) $expgui(diffBox).eDCzero
[10]1965        } elseif {[string range $expmap(htype_$hist) 1 2] == "XC" && \
1966                [histinfo $hist lam2] == 0.0} {
1967        #--------------------------
1968        # CW - x-ray 1 wavelength |
1969        #--------------------------
1970            grid [ label $expgui(diffBox).lDC1 -text "Refine wave" ] \
1971                    -column 1 -row 1
1972            grid [ checkbutton $expgui(diffBox).rfDC1 -variable entryvar(wref) ] \
1973                    -column 2 -row 1
1974            grid [ label $expgui(diffBox).lDCdifc -text wave ] \
1975                    -column 3 -row 1 -sticky w
1976            grid [ entry $expgui(diffBox).eDCdifc -textvariable entryvar(lam1) \
1977                    -width 15 ] -column 4 -row 1
[303]1978            set entrybox(lam1) $expgui(diffBox).eDCdifc
[10]1979            #
1980            grid [ label $expgui(diffBox).lDCzref -text "Refine zero" ] \
1981                    -column 1 -row 3 -sticky w
1982            grid [ checkbutton $expgui(diffBox).rfDCzref \
1983                    -variable entryvar(zref) ] -column 2 -row 3
1984            grid [ label $expgui(diffBox).lDCzero -text "Zero"] \
1985                    -column 3 -row 3
1986            grid [ entry $expgui(diffBox).eDCzero -textvariable entryvar(zero) \
1987                    -width 15 ] -column 4 -row 3
[303]1988            set entrybox(zero) $expgui(diffBox).eDCzero
[10]1989            #
1990            grid [ label $expgui(diffBox).lDCpref -text "Refine POLA" ] \
1991                    -column 1 -row 4 -sticky w
1992            grid [ checkbutton $expgui(diffBox).rfDCpref \
1993                    -variable entryvar(pref) ] -column 2 -row 4
1994            grid [ label $expgui(diffBox).lDCpola -text POLA ] \
1995                    -column 3 -row 4
1996            grid [ entry $expgui(diffBox).eDCpola \
1997                    -textvariable entryvar(pola) -width 15 ] -column 4 -row 4
[303]1998            set entrybox(pola) $expgui(diffBox).eDCpola
[10]1999            grid [ label $expgui(diffBox).lDCipola -text "IPOLA" ] \
2000                    -column 5 -row 4
2001            grid [ entry $expgui(diffBox).eDCipola -width 2 \
2002                    -textvariable entryvar(ipola)] -column 6 -row 4
[303]2003            set entrybox(ipola) $expgui(diffBox).eDCipola
[10]2004        } elseif {[string range $expmap(htype_$hist) 1 2] == "XC"} {
2005        #---------------------------
2006        # CW - x-ray 2 wavelengths |
2007        #---------------------------
2008            grid [ label $expgui(diffBox).lDCdifc -text wavelengths ] \
2009                    -column 3 -row 1 -sticky w
2010            grid [ entry $expgui(diffBox).eDCdifc -textvariable entryvar(lam1) \
2011                    -width 15 ] -column 4 -row 1
[303]2012            set entrybox(lam1) $expgui(diffBox).eDCdifc
[10]2013            grid [ entry $expgui(diffBox).eDCdifa -textvariable entryvar(lam2) \
2014                    -width 15 ] -column 5 -row 1
[303]2015            set entrybox(lam2) $expgui(diffBox).eDCdifa
[10]2016            #
2017            grid [ label $expgui(diffBox).lDCrref -text "Refine ratio" ] \
2018                    -column 1 -row 2 -sticky w
2019            grid [ checkbutton $expgui(diffBox).rfDCrref \
2020                    -variable entryvar(ratref) ] -column 2 -row 2
2021            grid [ label $expgui(diffBox).lDCratio -text Ratio ] \
2022                    -column 3 -row 2
2023            grid [ entry $expgui(diffBox).eDCkratio \
2024                    -textvariable entryvar(kratio) \
2025                    -width 15 ] -column 4 -row 2
[303]2026            set entrybox(kratio) $expgui(diffBox).eDCkratio
[10]2027            #
2028            grid [ label $expgui(diffBox).lDCzero -text "Zero"] \
2029                    -column 3 -row 3
2030            grid [ entry $expgui(diffBox).eDCzero -textvariable entryvar(zero) \
2031                    -width 15 ] -column 4 -row 3
2032            grid [ label $expgui(diffBox).lDCzref -text "Refine zero" ] \
2033                    -column 1 -row 3 -sticky w
[303]2034            set entrybox(zero) $expgui(diffBox).eDCzero
[10]2035            grid [ checkbutton $expgui(diffBox).rfDCzref \
2036                    -variable entryvar(zref) ] -column 2 -row 3
2037            grid [ label $expgui(diffBox).lDCpref -text "Refine POLA" ] \
2038                    -column 1 -row 4 -sticky w
2039            grid [ checkbutton $expgui(diffBox).rfDCpref \
2040                    -variable entryvar(pref) ] -column 2 -row 4
2041            grid [ label $expgui(diffBox).lDCpola -text POLA ] \
2042                    -column 3 -row 4
2043            grid [ entry $expgui(diffBox).eDCpola \
2044                    -textvariable entryvar(pola) -width 15 ] -column 4 -row 4
[303]2045            set entrybox(pola) $expgui(diffBox).eDCpola
[10]2046            grid [ label $expgui(diffBox).lDCipola -text "IPOLA" ] \
2047                    -column 5 -row 4
2048            grid [ entry $expgui(diffBox).eDCipola -width 2 \
2049                    -textvariable entryvar(ipola)] -column 6 -row 4
[303]2050            set entrybox(ipola) $expgui(diffBox).eDCipola
[10]2051        } elseif {[string range $expmap(htype_$hist) 1 2] == "XE"} {
2052        #-------------
2053        # ED - x-ray |
2054        #-------------
2055            grid [ label $expgui(diffBox).lDC1 -text "Refine 2theta" ] \
2056                    -column 1 -row 1
2057            grid [ checkbutton $expgui(diffBox).rfDC1 -variable entryvar(ttref) ] \
2058                    -column 2 -row 1
2059            grid [ label $expgui(diffBox).lDCdifc -text 2Theta ] \
2060                    -column 3 -row 1 -sticky w
2061            grid [ entry $expgui(diffBox).eDCdifc -textvariable entryvar(lam1) \
2062                    -width 15 ] -column 4 -row 1
[303]2063            set entrybox(lam1) $expgui(diffBox).eDCdifc
[10]2064            #
2065            grid [ label $expgui(diffBox).lDCpref -text "Refine POLA" ] \
2066                    -column 1 -row 4 -sticky w
2067            grid [ checkbutton $expgui(diffBox).rfDCpref \
2068                    -variable entryvar(pref) ] -column 2 -row 4
2069            grid [ label $expgui(diffBox).lDCpola -text POLA ] \
2070                    -column 3 -row 4
2071            grid [ entry $expgui(diffBox).eDCpola \
2072                    -textvariable entryvar(pola) -width 15 ] -column 4 -row 4
[303]2073            set entrybox(pola) $expgui(diffBox).eDCpola
[10]2074            grid [ label $expgui(diffBox).lDCipola -text "IPOLA" ] \
2075                    -column 5 -row 4
2076            grid [ entry $expgui(diffBox).eDCipola -width 2 \
2077                    -textvariable entryvar(ipola)] -column 6 -row 4
[303]2078            set entrybox(ipola) $expgui(diffBox).eDCipola
[10]2079        }
2080    } elseif {$expgui(globalmode) == 1} {
2081        #-------------
2082        # Global TOF |
2083        #-------------
2084        grid [ label $expgui(diffBox).lDCrc -text "Refine DIFC" ] \
2085                -column 1 -row 1
2086        grid [ checkbutton $expgui(diffBox).rfDCrc -variable entryvar(dcref) ] \
2087                -column 2 -row 1
2088        grid [button $expgui(diffBox).bDCdifc -text "Set DIFC Globally" \
2089                -command "editglobalparm histinfo difc {DIFC}"] -column 3 -row 1
2090        #
2091        grid [ label $expgui(diffBox).lDCra -text "Refine DIFA" ] \
2092                -column 1 -row 2
2093        grid [ checkbutton $expgui(diffBox).rfDCra -variable entryvar(daref) ] \
2094                -column 2 -row 2
2095        grid [ button $expgui(diffBox).bDCdifa -text "Set DIFA Globally" \
2096                -command "editglobalparm histinfo difa {DIFA}"] -column 3 -row 2
2097        #
2098        grid [ label $expgui(diffBox).lDCzref -text "Refine zero" ] \
2099                -column 1 -row 3 -sticky w
2100        grid [ checkbutton $expgui(diffBox).rfDCzref \
2101                -variable entryvar(zref) ] -column 2 -row 3
2102        grid [ button $expgui(diffBox).bDCzero -text "Set ZERO Globally" \
2103                -command "editglobalparm histinfo zero {Zero}"] -column 3 -row 3
2104    } elseif {$expgui(globalmode) == 2} {
2105        #--------------------
2106        # Global CW neutron |
2107        #--------------------
2108        grid [ label $expgui(diffBox).lDC1 -text "Refine wave" ] \
2109                -column 1 -row 1
2110        grid [ checkbutton $expgui(diffBox).rfDC1 -variable entryvar(wref) ] \
2111                -column 2 -row 1
2112        grid [button $expgui(diffBox).bDCdifc -text "Set Wave Globally" \
2113                -command "editglobalparm histinfo lam1 Wavelength"] \
2114                -column 3 -row 1
2115        #
2116        grid [ label $expgui(diffBox).lDCzref -text "Refine zero" ] \
2117                -column 1 -row 3 -sticky w
2118        grid [ checkbutton $expgui(diffBox).rfDCzref \
2119                -variable entryvar(zref) ] -column 2 -row 3
2120        grid [button $expgui(diffBox).bDCzero -text "Set Zero Globally" \
2121                -command "editglobalparm histinfo zero Zero"] -column 3 -row 3
2122    } elseif {$expgui(globalmode) == 4} {
2123        #----------------------
2124        # Global CW mono xray |
2125        #----------------------
2126        grid [ label $expgui(diffBox).lDC1 -text "Refine wave" ] \
2127                -column 1 -row 1
2128        grid [ checkbutton $expgui(diffBox).rfDC1 -variable entryvar(wref) ] \
2129                -column 2 -row 1
2130        grid [button $expgui(diffBox).bDCdifc -text "Set Wave Globally" \
2131                -command "editglobalparm histinfo lam1 Wavelength"] \
2132                -column 3 -row 1
2133        #
2134        grid [ label $expgui(diffBox).lDCzref -text "Refine zero" ] \
2135                -column 1 -row 3 -sticky w
2136        grid [ checkbutton $expgui(diffBox).rfDCzref \
2137                -variable entryvar(zref) ] -column 2 -row 3
2138        grid [button $expgui(diffBox).bDCzero -text "Set Zero Globally" \
2139                -command "editglobalparm histinfo zero Zero"] -column 3 -row 3
2140        #
2141        grid [ label $expgui(diffBox).lDCpref -text "Refine POLA" ] \
2142                -column 1 -row 4 -sticky w
2143        grid [ checkbutton $expgui(diffBox).rfDCpref \
2144                -variable entryvar(pref) ] -column 2 -row 4
2145        grid [button $expgui(diffBox).bDCpola -text "Set POLA Globally" \
2146                -command "editglobalparm histinfo pola POLA"] -column 3 -row 4
2147        grid [button $expgui(diffBox).bDCipola -text "Set IPOLA Globally" \
2148                -command "editglobalparm histinfo ipola IPOLA"] -column 4 -row 4
2149    } elseif {$expgui(globalmode) == 3} {
2150        #------------------------
2151        # Global alpha 1,2 xray |
2152        #------------------------
2153        grid [button $expgui(diffBox).bDCl1 -text "Set Wave1 Globally" \
2154                -command "editglobalparm histinfo lam1 {Wavelength 1}"] \
2155                -column 3 -row 1
2156        grid [button $expgui(diffBox).bDCl2 -text "Set Wave2 Globally" \
2157                -command "editglobalparm histinfo lam2 {Wavelength 2}"] \
2158                -column 4 -row 1
2159        #
2160        grid [ label $expgui(diffBox).lDCratref -text "Refine Ratio" ] \
[557]2161                -column 1 -row 2 -sticky w
[10]2162        grid [ checkbutton $expgui(diffBox).rfDCratref \
[557]2163                -variable entryvar(ratref) ] -column 2 -row 2
[10]2164        grid [button $expgui(diffBox).bDCrrat -text "Set Ratio Globally" \
2165                -command "editglobalparm histinfo ratio {Wavelength Ratio}"] \
[557]2166                -column 3 -row 2
[10]2167        #
2168        grid [ label $expgui(diffBox).lDCzref -text "Refine zero" ] \
2169                -column 1 -row 3 -sticky w
2170        grid [ checkbutton $expgui(diffBox).rfDCzref \
2171                -variable entryvar(zref) ] -column 2 -row 3
2172        grid [button $expgui(diffBox).bDCzero -text "Set Zero Globally" \
2173                -command "editglobalparm histinfo zero Zero"] -column 3 -row 3
2174        #
2175        grid [ label $expgui(diffBox).lDCpref -text "Refine POLA" ] \
2176                -column 1 -row 4 -sticky w
2177        grid [ checkbutton $expgui(diffBox).rfDCpref \
2178                -variable entryvar(pref) ] -column 2 -row 4
2179        grid [button $expgui(diffBox).bDCpola -text "Set POLA Globally" \
2180                -command "editglobalparm histinfo pola POLA"] -column 3 -row 4
2181        grid [button $expgui(diffBox).bDCipola -text "Set IPOLA Globally" \
2182                -command "editglobalparm histinfo ipola IPOLA"] -column 4 -row 4
2183    } elseif {$expgui(globalmode) == 5} {
2184        #-----------------
2185        # Global ED xray |
2186        #-----------------
2187        grid [ label $expgui(diffBox).lDC1 -text "Refine 2theta" ] \
2188                -column 1 -row 1
2189        grid [ checkbutton $expgui(diffBox).rfDC1 -variable entryvar(ttref) ] \
2190                -column 2 -row 1
2191        grid [button $expgui(diffBox).bDCdifc -text "Set 2Theta Globally" \
2192                -command "editglobalparm histinfo ratio {Fixed 2Theta}"] \
2193                -column 3 -row 1
2194        #
2195        grid [ label $expgui(diffBox).lDCpref -text "Refine POLA" ] \
2196                -column 1 -row 4 -sticky w
2197        grid [ checkbutton $expgui(diffBox).rfDCpref \
2198                -variable entryvar(pref) ] -column 2 -row 4
2199        grid [button $expgui(diffBox).bDCpola -text "Set POLA Globally" \
2200                -command "editglobalparm histinfo pola POLA"] -column 3 -row 4
2201        grid [button $expgui(diffBox).bDCipola -text "Set IPOLA Globally" \
2202                -command "editglobalparm histinfo ipola IPOLA"] -column 4 -row 4
2203    }
2204    if {$expgui(globalmode) == 0} {
[557]2205        grid [frame $expgui(diffBox).d] -column 5 -row 1 -rowspan 3 \
[10]2206                -columnspan 2 -sticky e
2207    } else {
[557]2208        grid [frame $expgui(diffBox).d] -column 4 -row 2 -rowspan 2 \
[10]2209                -columnspan 2 -sticky e
2210    }
2211    grid [label $expgui(diffBox).d.lDamp -text "Damping  "] \
2212            -column 1 -row 1 
2213    tk_optionMenu $expgui(diffBox).d.om entryvar(ddamp) 0 1 2 3 4 5 6 7 8 9
2214    grid $expgui(diffBox).d.om -column 2 -row 1 
2215    grid columnconfigure $expgui(diffBox) 9  -weight 1
2216    grid columnconfigure $expgui(diffBox) 0  -weight 1
2217    update idletasks
2218    # enable traces on entryvar now
2219    set entrycmd(trace) 1
2220}
2221
[255]2222# this gets the phase flags as set in the expgui(pflag*) elements
2223# (linked to phase flag checkbuttons) and the sets the "HST xx NPHAS" flags
2224# accordingly using SetPhaseFlag
[250]2225proc GetPhaseFlags {hist} {
2226    global expmap expgui
2227    set plist {}
2228    foreach p $expmap(phaselist) {
2229        if {$expgui(pflag$p)} {lappend plist $p}
2230    }
2231    SetPhaseFlag $hist $plist
[908]2232    RecordMacroEntry "SetPhaseFlag [list $hist] [list $plist]" 0
[250]2233    incr expgui(changed)
[908]2234    RecordMacroEntry "incr expgui(changed)" 0
[670]2235    # set the powpref warning (1 = suggested)
2236    set expgui(needpowpref) 2
2237    set msg "Phase flags" 
2238    if {[string first $msg $expgui(needpowpref_why)] == -1} {
2239        append expgui(needpowpref_why) "\t$msg were changed\n"
2240    }
[250]2241    mapexp
[838]2242    # reset the phase selection
2243    set expgui(curPhase) {}
[250]2244}
2245
[10]2246#-----------------------------------------------------------------------
[14]2247# populate the Scaling page
[10]2248#-----------------------------------------------------------------------
2249proc DisplayFrac {} {
[303]2250    global expgui entrycmd entryvar entrybox expmap
[10]2251
[112]2252    # trap if more than one histogram is selected unless global mode
2253    if {$expgui(globalmode) == 0 && [llength $expgui(curhist)] > 1} {
2254        set expgui(curhist) [lindex $expgui(curhist) 0] 
2255    }
2256
[10]2257    # display the selected histograms
2258    $expgui(fracFrame).hs.lbox selection clear 0 end
2259    foreach h $expgui(curhist) {
2260        $expgui(fracFrame).hs.lbox selection set $h
2261    }
2262
2263    # disable traces on entryvar
2264    set entrycmd(trace) 0
2265
2266    # get histogram list
2267    set histlist {}
2268    foreach item $expgui(curhist) {
2269        lappend histlist [lindex $expmap(powderlist) $item]
2270    }
2271
[20]2272    # must have at least one histogram selected here
2273    if {[llength $histlist] == 0} {
2274        foreach var {scale sref sdamp} {
2275            set entrycmd($var) ""
2276            set entryvar($var) ""
2277        }
2278        set parm [grid info $expgui(scaleBox).but1]
2279        if {$parm != ""} {
2280            grid forget  $expgui(scaleBox).but1
2281            eval grid $expgui(scaleBox).ent1 $parm
2282        }
2283        # destroy the contents of the frame
2284        set phaseFractf1 $expgui(FracBox).f
[152]2285        eval destroy [winfo children $phaseFractf1]
[20]2286        # reenable traces on entryvar
2287        set entrycmd(trace) 1
2288        return
2289    }
2290
[14]2291    #--------------
2292    # Scale factor
2293    #--------------
2294    if {$expgui(globalmode) != 0} {
2295        foreach var {scale sref sdamp} {
2296            set entrycmd($var) "histinfo [list $histlist] $var"
2297            set entryvar($var) [histinfo [lindex $histlist 0] $var]
2298        }
[303]2299        # reset scale to black
2300        catch {$entrybox(scale) config -fg black}
[14]2301        set parm [grid info $expgui(scaleBox).ent1]
2302        if {$parm != ""} {
2303            grid forget  $expgui(scaleBox).ent1
2304            eval grid $expgui(scaleBox).but1 $parm
2305        }
2306    } else {
2307        set hist $histlist
2308        foreach var {scale sref sdamp} {
2309            set entrycmd($var) "histinfo $hist $var"
2310            set entryvar($var) [eval $entrycmd($var)]
2311        }
[303]2312        # reset scale to black
2313        catch {$entrybox(scale) config -fg black}
[14]2314        set parm [grid info $expgui(scaleBox).but1]
2315        if {$parm != ""} {
2316            grid forget  $expgui(scaleBox).but1
2317            eval grid $expgui(scaleBox).ent1 $parm
2318        }
2319    }
2320
2321    #----------------
2322    # Phase Fractions
2323    #----------------
2324    set phaseFractf1 $expgui(FracBox).f
[10]2325    # destroy the contents of the frame
[152]2326    eval destroy [winfo children $phaseFractf1]
[10]2327    if {$expgui(globalmode) != 0} {
[133]2328        set txt "Phase Fractions for Histograms: [CompressList $histlist]"
[14]2329    } else {
[133]2330        set txt "Phase Fractions"
[10]2331    }
[823]2332    $expgui(fracFrame).f1.phaseFrac configure -text $txt
[10]2333    # Create the frame inside the canvas, One frame for each Phase.
2334    foreach i {1 2 3 4 5 6 7 8 9} {set phasehistlist($i) ""}
2335    foreach hist $histlist {
2336        foreach i $expmap(phaselist_$hist) {
2337            lappend phasehistlist($i) $hist
2338        }
2339    }
2340    foreach i {1 2 3 4 5 6 7 8 9} {
2341        if {[llength $phasehistlist($i)] == 0} continue
2342        set framePF [frame $phaseFractf1.pF$i -relief groove  -bd 4]
[14]2343        grid $framePF -column 0 -row $i -sticky ew
[10]2344        # Label Heading for each phase.
2345        if {$expgui(globalmode) != 0} {
2346            grid [label $framePF.l1 \
[133]2347                    -text "Phase $i Hist: [CompressList $phasehistlist($i)]"] \
[10]2348                    -column 0 -row 0 -sticky nws
2349            grid [button $framePF.but1 -text "Set Globally" \
2350                    -command "editglobalparm hapinfo frac \"Phase $i Fraction\" \
2351                    [list $phasehistlist($i)] $i" \
2352                    ] -column 1 -row 0
[866]2353            set entrycmd(frref$i) "hapinfo [list $histlist] $i frref"
2354            set entryvar(frref$i) [hapinfo $hist $i frref]
2355            set entrycmd(frdamp$i) "hapinfo [list $histlist] $i frdamp"
[10]2356        } else {
[133]2357            grid [label $framePF.l1  -text "Phase $i"] \
[10]2358                    -column 0 -row 0 -sticky nws
2359            grid [entry $framePF.ent -textvariable entryvar(frac$i) -width 15]\
2360                    -column 1 -row 0
[303]2361            set entrybox(frac$i) $framePF.ent
[866]2362            set entrycmd(frref$i) "hapinfo $hist $i frref"
2363            set entryvar(frref$i) [hapinfo $hist $i frref]
2364            set entrycmd(frdamp$i) "hapinfo $hist $i frdamp"
[10]2365        }
2366        set entrycmd(frac$i) "hapinfo $hist $i frac"
2367        set entryvar(frac$i) [hapinfo $hist $i frac]
2368        grid [label $framePF.l2  -text "  Refine"] \
2369                -column 2 -row 0 -sticky nws
2370        grid [checkbutton $framePF.cb -variable entryvar(frref$i)] \
2371                -column 3 -row 0 -sticky nws
2372        grid [label $framePF.l3  -text "  Damping"] \
2373                -column 4 -row 0 -sticky nws
2374        tk_optionMenu $framePF.tkOptDamp entryvar(frdamp$i) \
2375                0 1 2 3 4 5 6 7 8 9     
2376        set entryvar(frdamp$i) [hapinfo $hist $i frdamp]
2377        grid $framePF.tkOptDamp -row 0 -sticky nsw -column 5
2378    }
[14]2379    # resize the scroll window to match the actual
[10]2380    update idletasks
[14]2381    $expgui(FracBox) config -scrollregion [grid bbox $expgui(FracBox).f]
2382    $expgui(FracBox) config -width [lindex [grid bbox $expgui(FracBox).f] 2]
2383    update idletasks
[10]2384    # enable traces on entryvar now
2385    set entrycmd(trace) 1
2386}
2387
2388#-----------------------------------------------------------------------
2389# display the profile page
2390#-----------------------------------------------------------------------
2391proc DisplayProfile {} {
[303]2392    global expgui entrycmd entryvar entrybox expmap
[112]2393
2394    # trap if more than one histogram is selected unless global mode
2395    if {$expgui(globalmode) == 0 && [llength $expgui(curhist)] > 1} {
2396        set expgui(curhist) [lindex $expgui(curhist) 0] 
2397    }
2398    # display the selected histograms
[10]2399    $expgui(profFrame).hs.lbox selection clear 0 end
2400    foreach h $expgui(curhist) {
2401        $expgui(profFrame).hs.lbox selection set $h
2402    }
2403
2404    # destroy the contents of the frame
[152]2405    eval destroy [winfo children $expgui(ProfileBox).f]
[327]2406    # since the next steps can take a while, do a screen update
2407    update idletasks
[10]2408
2409    if {$expgui(globalmode) == 0} {
[20]2410        # must have at least one histogram selected here
2411        if {[llength $expgui(curhist)] == 0} return
[21]2412        # disable traces on entryvar for right now
2413        set entrycmd(trace) 0
[10]2414        set hist [lindex $expmap(powderlist) $expgui(curhist)]
[97]2415        # no defined histograms?
2416        if {$hist == ""} return
[10]2417        # Create one frame for each Phase.
2418        set ind -1
2419        set htype [string range $expmap(htype_$hist) 2 2]
[838]2420        set zflag 0
2421        if {$htype == "C"} {
2422            set zflag [histinfo $hist zref]
2423        }
[10]2424        foreach i $expmap(phaselist_$hist) {
2425            incr ind
2426            # Label Heading for each phase.
2427            set ptype [string trim [hapinfo $hist $i proftype]]
[823]2428            grid [TitleFrame $expgui(ProfileBox).f.$i \
2429                      -text "Hist $hist -- Phase $i (type $ptype)" \
2430                      -relief groove -bd 2] \
2431                -column 0 -row $ind -sticky ew
2432            set ProfileFrame [$expgui(ProfileBox).f.$i getframe]
2433            grid [frame $ProfileFrame.1] \
2434                -column 0 -row 0 -columnspan 10
2435            pack [label $ProfileFrame.1.l  \
2436                      -text Damping]\
2437                -side left
[133]2438            tk_optionMenu $ProfileFrame.1.tkOptDamp entryvar(pdamp_$i) \
[10]2439                    0 1 2 3 4 5 6 7 8 9
2440            set entrycmd(pdamp_$i) "hapinfo $hist $i pdamp"
2441            set entryvar(pdamp_$i) [hapinfo $hist $i pdamp]
[133]2442            pack $ProfileFrame.1.tkOptDamp -side left
2443            pack [label $ProfileFrame.1.l1 \
[10]2444                    -text "  Peak cutoff"]\
2445                    -side left
[133]2446            pack [entry $ProfileFrame.1.e1  \
[10]2447                    -width 10 -textvariable entryvar(pcut_$i)]\
2448                    -side left
[303]2449            set entrybox(pcut_$i) $ProfileFrame.1.e1
[10]2450            set entrycmd(pcut_$i) "hapinfo $hist $i pcut"
2451            set entryvar(pcut_$i) [hapinfo $hist $i pcut]
[234]2452
2453            pack [button $ProfileFrame.1.b1  \
2454                    -text "Change Type" \
2455                    -command "ChangeProfileType $hist $i"]\
2456                    -side left
[10]2457           
2458            set col -1
2459            set row 1
2460            set nterms [hapinfo $hist $i profterms]
[234]2461            set lbls "dummy [GetProfileTerms $i $hist $ptype]"
[10]2462            for { set num 1 } { $num <= $nterms } { incr num } {
2463                set term {}
2464                catch {set term [lindex $lbls $num]}
2465                if {$term == ""} {set term $num}
2466                incr col
[133]2467                grid [label $ProfileFrame.l${num}_${i} -text "$term"] \
[10]2468                        -row $row -column $col
2469                incr col
[133]2470                grid [checkbutton $ProfileFrame.ref${num}_${i} \
[10]2471                        -variable entryvar(pref${num}_$i)] -row $row -column $col
2472                set entrycmd(pref${num}_$i) "hapinfo $hist $i pref$num"
2473                set entryvar(pref${num}_$i) [hapinfo $hist $i pref$num]
2474                incr col
[133]2475                grid [entry $ProfileFrame.ent${num}_${i} \
[10]2476                        -textvariable entryvar(pterm${num}_$i)\
[28]2477                        -width 12] -row $row -column $col
[303]2478                set entrybox(pterm${num}_$i) $ProfileFrame.ent${num}_${i}
[10]2479                set entrycmd(pterm${num}_$i) "hapinfo $hist $i pterm$num"
2480                set entryvar(pterm${num}_$i) [hapinfo $hist $i pterm$num]
[838]2481                # disable trns & shft when zero is refined
2482                if {$zflag && ($term == "trns" || $term == "shft")} {
2483                    if {$entryvar(pref${num}_$i)} {
2484                        incr expgui(changed)
2485                        set entryvar(pref${num}_$i) 0
2486                    }
2487                    $ProfileFrame.l${num}_${i} config -fg gray
2488                    $ProfileFrame.ref${num}_${i} config -state disabled
2489                    $ProfileFrame.ent${num}_${i} config -fg gray
2490                }
[10]2491                if {$col > 6} {set col -1; incr row}
2492            }
2493        }
2494        grid columnconfigure $expgui(ProfileBox).f 0 -weight 1
2495    } else {
2496        # get histogram list
2497        set histlist {}
2498        foreach item $expgui(curhist) {
2499            lappend histlist [lindex $expmap(powderlist) $item]
2500        }
[20]2501        # must have at least one histogram selected here
2502        if {[llength $histlist] == 0} return
[21]2503        # disable traces on entryvar for right now
2504        set entrycmd(trace) 0
[112]2505        # loop through histograms & phases, set up an array by phase & profile type
2506        catch {unset prtyparray histarray phasearray}
[10]2507        foreach hist $histlist {
[112]2508            foreach phase $expmap(phaselist_$hist) {
2509                set prtyp [string trim [hapinfo $hist $phase proftype]]
2510                set key ${prtyp}_$phase
2511                lappend prtyparray($key) $hist
2512                lappend histarray($key) $hist
2513                lappend phasearray($key) $phase
[10]2514            }
2515        }
2516       
2517        set ptype ""
[112]2518        set i -1
2519        # loop over all combined phases and profile types, sorted 1st by profile number
2520        foreach key [lsort [array names prtyparray]] {
[10]2521            # split key
[112]2522            scan $key %d_%d prftyp p
2523
2524            if {$ptype != $prftyp || !$expgui(globalphasemode)} {
2525                set ptype $prftyp
2526                set curhistlist $histarray($key)
2527                set curphaslist $phasearray($key)
2528               
2529                set hist1 [lindex $curhistlist 0]
2530                set phase1 [lindex $curphaslist 0]
2531                set nterms [hapinfo $hist1 $phase1 profterms]
2532                set htype [string range $expmap(htype_$hist1) 2 2]
[234]2533                set lbls "dummy [GetProfileTerms $phase1 $hist1 $ptype]"
[112]2534                # Create a frame for this type
2535                incr i
2536                set boxtitle "Phase $p, hist [CompressList $histarray($key)]"
[823]2537                grid [TitleFrame $expgui(ProfileBox).f.$i \
2538                          -text "(type $ptype)" \
2539                          -relief groove -bd 2] \
2540                    -column 0 -row $i -sticky ew
2541                set ProfileFrame [$expgui(ProfileBox).f.$i getframe]
2542                grid [frame $ProfileFrame.0] \
2543                    -column 0 -row 0 -columnspan 20 -sticky ew
[133]2544                grid [label $ProfileFrame.0.1  \
[112]2545                        -anchor w] -row 0 -column 1
[133]2546                grid [frame $ProfileFrame.1] \
[10]2547                        -column 0 -row 1 -columnspan 20 -sticky ew
[133]2548                grid [label $ProfileFrame.1.2  \
[10]2549                        -text "Damping"] -row 0 -column 2
[133]2550                tk_optionMenu $ProfileFrame.1.tkOptDamp \
[112]2551                        entryvar(pdamp_$i) 0 1 2 3 4 5 6 7 8 9
[133]2552                grid $ProfileFrame.1.tkOptDamp -row 0 -column 3
2553                grid [button $ProfileFrame.1.edit \
[112]2554                        -text "Global Edit"] -row 0 -column 4 -sticky w
2555                set entryvar(pdamp_$i) [hapinfo $hist $phase pdamp]
[234]2556                grid [button $ProfileFrame.1.b1 -text "Change Type"] \
2557                        -row 0 -column 5 -sticky w
[10]2558                set col -1
2559                set row 2
2560                for { set num 1 } { $num <= $nterms } { incr num } {
2561                    set term {}
2562                    catch {set term [lindex $lbls $num]}
2563                    if {$term == ""} {set term $num}
2564                    incr col
[133]2565                    grid [label $ProfileFrame.l${num}_${i} \
[112]2566                            -text "$term"] -row $row -column $col
[10]2567                    incr col
[133]2568                    grid [checkbutton $ProfileFrame.ref${num}_${i} \
[10]2569                            -variable entryvar(pref${num}_$i)] \
2570                            -row $row -column $col
2571                    set entryvar(pref${num}_$i) [hapinfo $hist $phase pref$num]
2572                    if {$col > 10} {set col -1; incr row}
2573                }
2574                grid columnconfigure $expgui(ProfileBox).f 0 -weight 1
[112]2575            } else {
2576                # add to the current entry
2577                eval lappend curhistlist $histarray($key)
2578                eval lappend curphaslist $phasearray($key)
2579                append boxtitle "\nPhase $p, hist [CompressList $histarray($key)]"
[10]2580            }
[133]2581            $ProfileFrame.0.1 config -text $boxtitle
2582            $ProfileFrame.1.edit config -command "\
[112]2583                    EditProfile \"\n$boxtitle\" \
2584                    [list $curhistlist] \
2585                    [list $curphaslist]"
[234]2586            $ProfileFrame.1.b1 config -command "ChangeProfileType \
2587                    [list $curhistlist] [list $curphaslist]" 
[112]2588            set entrycmd(pdamp_$i) "hapinfo \
2589                    [list $curhistlist] \
2590                    [list $curphaslist] pdamp"
2591            for { set num 1 } { $num <= $nterms } { incr num } {
2592                set entrycmd(pref${num}_$i) "hapinfo \
2593                        [list $curhistlist] \
2594                        [list $curphaslist] pref$num"
2595            }
[10]2596        }
2597    }
2598   
2599    # resize the scroll window to match the actual
2600    update idletasks
2601    $expgui(ProfileBox) config -scrollregion [grid bbox $expgui(ProfileBox).f]
2602    $expgui(ProfileBox) config -width [lindex [grid bbox $expgui(ProfileBox).f] 2]
2603    update idletasks
[133]2604    ResizeNotebook
[10]2605    # enable traces on entryvar now
2606    set entrycmd(trace) 1
2607}
2608
2609# process the bit settings in the print options
2610#   bitnum -- the number of the bit to be tested/set starting at 0 for the LSBit
2611proc printsetting {bitnum "action get" "value {}"} {
2612    global entryvar expgui
2613    if {$action == "get"} {
2614        return [expr ([expinfo print] & int(pow(2,$bitnum))) != 0]
2615    } elseif $value {
2616        set newval [expr ([expinfo print] | int(pow(2,$bitnum)))]
2617    } else {
2618        set newval [expr ([expinfo print] & ~int(pow(2,$bitnum)))]
2619    }
2620    expinfo print set $newval
[908]2621    RecordMacroEntry "incr expgui(changed); expinfo print set $newval" 1
[10]2622    set expgui(printopt) "Print Options ([expinfo print])"
2623}
[20]2624
2625# need to respond to mouse presses -- control variable associated with extract Fobs
2626# and set the LeBail extraction flags
2627proc SetupExtractHist {} {
2628    global expgui entrycmd entryvar expmap
2629
2630    # display the selected histograms
2631    $expgui(lsFrame).hs.lbox selection clear 0 end
2632    foreach h $expgui(curhist) {
2633        $expgui(lsFrame).hs.lbox selection set $h
2634    }
2635
2636    # get histogram list
2637    set histlist {}
2638    foreach item $expgui(curhist) {
[30]2639        set hist [lindex $expmap(powderlist) $item]
2640        if {$hist != ""} {lappend histlist $hist}
[20]2641    }
2642    set entrycmd(fobsextract) "histinfo [list $histlist] foextract"
[30]2643    if {[llength $histlist] == 0 || [string trim $histlist] == ""} {
[399]2644        set entrycmd(LBdamp) ""
[20]2645        foreach phase {1 2 3 4 5 6 7 8 9} {
[133]2646            $expgui(FobsExtractFrame).l$phase config -fg grey
[20]2647            set expgui(Fextract$phase) {}
[28]2648            foreach item $expgui(ExtractSettingsRadiobuttons) {
2649                ${item}$phase config -state disabled -bd 1
[20]2650            }
2651        }
2652    } elseif {[llength $histlist] == 1} {
[21]2653        # disable traces on entryvar
2654        set entrycmd(trace) 0
[20]2655        set entryvar(fobsextract) [histinfo $histlist foextract]
[399]2656        set entrycmd(LBdamp) "histinfo $histlist LBdamp"
2657        set entryvar(LBdamp) [histinfo $histlist LBdamp]
[20]2658        foreach phase {1 2 3 4 5 6 7 8 9} {
2659            # is the phase present?
2660            if {[lsearch -exact $expmap(phaselist_$histlist) $phase] == -1} {
[133]2661                $expgui(FobsExtractFrame).l$phase config -fg grey
[20]2662                set expgui(Fextract$phase) {}
[28]2663                foreach item $expgui(ExtractSettingsRadiobuttons) {
2664                    ${item}$phase config -state disabled -bd 1
[20]2665                }
2666            } else {
[133]2667                $expgui(FobsExtractFrame).l$phase config -fg black
[28]2668                foreach item $expgui(ExtractSettingsRadiobuttons) {
2669                    ${item}$phase config -state normal -bd 2
[20]2670                }
2671                set expgui(Fextract$phase) [hapinfo $histlist $phase extmeth]
2672            }
2673        }
2674    } elseif {[llength $histlist] > 1} {
[21]2675        # disable traces on entryvar
[399]2676        set entrycmd(LBdamp) "histinfo [list $histlist] LBdamp"
2677        set entryvar(LBdamp) [histinfo [lindex $histlist 0] LBdamp]
[21]2678        set entrycmd(trace) 0
[20]2679        # multiple histograms need phases in any histogram
2680        foreach phase {1 2 3 4 5 6 7 8 9} {
2681            set gotphase($phase) 0
2682        }           
2683        foreach hist $histlist {
2684            foreach phase $expmap(phaselist_$hist) {
2685                set gotphase($phase) 1
2686            }
2687        }
2688        foreach phase {1 2 3 4 5 6 7 8 9} {
2689            set expgui(Fextract$phase) {}
2690            if $gotphase($phase) {
[133]2691                $expgui(FobsExtractFrame).l$phase config -fg black
[28]2692                foreach item $expgui(ExtractSettingsRadiobuttons) {
2693                    ${item}$phase config -state normal -bd 2
[20]2694                }
2695            } else {
[133]2696                $expgui(FobsExtractFrame).l$phase config -fg grey
[28]2697                foreach item $expgui(ExtractSettingsRadiobuttons) {
2698                    ${item}$phase config -state disabled -bd 1
[20]2699                }
2700            }
2701        }
2702    }
2703    # reenable traces
2704    set entrycmd(trace) 1
2705}
2706# respond to a change in the fobs extraction method for a phase
2707# force the main extraction flag on, if fobs extraction is selected for any phase
2708proc HistExtractSet {phase} {
2709    global expgui entryvar expmap
2710    foreach item $expgui(curhist) {
2711        lappend histlist [lindex $expmap(powderlist) $item]
2712    }
2713    hapinfo $histlist $phase extmeth set $expgui(Fextract$phase)
[908]2714    RecordMacroEntry "incr expgui(changed); hapinfo [list $histlist] $phase extmeth set $expgui(Fextract$phase)" 1
[48]2715    incr expgui(changed)
[20]2716    if {$expgui(Fextract$phase) != 0} {set entryvar(fobsextract) 1}
2717}
[10]2718#---------------------------- Global Edit Functions ------------------------
2719proc editbackground {} {
2720    global expgui expmap entrycmd
2721    set histlist {}
2722    foreach n $expgui(curhist) {
2723        lappend histlist [lindex $expmap(powderlist) $n]
2724    }
[20]2725    if {[llength $histlist] == 0} return
2726
2727    set w .back
2728    catch {destroy $w}
2729    toplevel $w -bg beige
[10]2730    if {$expgui(globalmode) != 0} {
[250]2731        wm title $w "Global Edit Background" 
2732    } else {
[10]2733        wm title $w "Edit Background" 
2734    }
2735   
2736    pack [frame $w.0 -bd 6 -relief groove  -bg beige \
2737            ] -side top -expand yes -fill both
2738    if {[llength $histlist] > 1} {
2739        grid [label $w.0.a \
2740            -text "Setting background terms for histograms [CompressList $histlist]" \
2741            -bg beige] -row 0 -column 0 -columnspan 10
2742    } else {
2743        grid [label $w.0.a \
2744            -text "Setting background terms for histogram $histlist" \
[399]2745            -bg beige] -row 0 -column 0 -columnspan 4
2746        grid [button $w.0.bkg -text "Fit Background\nGraphically" \
2747                -command "QuitEditBackground $w; bkgedit $histlist"] \
2748                -row 0 -column 4 -rowspan 2
2749        grid columnconfig $w.0 0 -weight 1
2750        grid columnconfig $w.0 4 -weight 1
[10]2751    }
2752    set hist [lindex $histlist 0]
[399]2753    grid [label $w.0.b -text "Function type" -bg beige]  -row 1 -column 0 -sticky e
[10]2754
2755    # disable traces on  expgui(backtype) & expgui(backterms) now
2756    set entrycmd(trace) 0
2757
2758    # number of terms
2759    set expgui(backtype) [histinfo $hist backtype]
2760    set expgui(orig_backtype) $expgui(backtype)
2761    set expgui(prev_backtype) $expgui(backtype)
[250]2762    set typemenu [tk_optionMenu $w.0.type expgui(backtype) null]
2763    $typemenu delete 0 end
2764    foreach item {
[399]2765        "1 - Shifted Chebyschev"
[250]2766        "2 - Cosine Fourier series"
2767        "4 - Power series in Q**2n/n!"
2768        "5 - Power series in n!/Q**2n"
2769        "6 - Power series in Q**2n/n! and n!/Q**2n"
2770        "7 - Linear interpolation function"
2771        "8 - Reciprocal interpolation function"
2772    } {
2773        set val [lindex $item 0]
2774        $typemenu insert end radiobutton -variable expgui(backtype) \
2775                -label $item -value $val
2776    }
[720]2777# removed
2778#       "3 - Radial distribution peaks"
[250]2779
[10]2780    grid $w.0.type   -row 1 -column 1
2781    grid [label $w.0.c -text "  Number of terms"  -bg beige] -row 1 -column 2
2782
2783    # function type
2784    set expgui(backterms) [histinfo $hist backterms]
2785    set expgui(orig_backterms) $expgui(backterms) 
2786    set list {}; for {set i 1} {$i <= 36} {incr i} {lappend list $i}
2787    eval tk_optionMenu $w.0.terms expgui(backterms) $list
2788    grid $w.0.terms   -row 1 -column 3 
2789    # enable traces on  expgui(backtype) & expgui(backterms) now
2790    set entrycmd(trace) 1
2791
2792    #set background terms
2793    for {set num 1 } { $num <= 36 } { incr num } {
2794        set var "bterm$num"
2795        set expgui($var) {}
2796        set expgui(orig_$var) {}
2797    }
2798    if {[llength $histlist] == 1} {
2799        for {set num 1 } { $num <= $expgui(backterms) } { incr num } {
2800            set var "bterm$num"
2801            set expgui($var) [histinfo $histlist $var]
2802            set expgui(orig_$var) $expgui($var)
2803        }
2804    }
2805
2806    pack [frame $w.1 -bd 6 -relief groove  -bg beige] -side top \
2807            -expand yes -fill both
2808    ShowBackTerms $w.1
2809
2810    set expgui(temp) {}
[323]2811    pack [frame $w.b -bg beige] -fill x -expand yes -side top
[928]2812    grid [button $w.b.2 -text Continue -command "destroy $w"] -row 0 -column 1
2813    grid [button $w.b.3 -text Cancel \
[323]2814            -command "QuitEditBackground $w"] -row 0 -column 2
2815    grid [button $w.b.help -text Help -bg yellow \
2816            -command "MakeWWWHelp expgui3.html EditBackground"] \
2817            -row 0 -column 4
2818    grid columnconfig $w.b 0 -weight 1
2819    grid columnconfig $w.b 3 -weight 1
2820    bind $w <Key-F1> "MakeWWWHelp expgui3.html EditBackground"
[93]2821    bind $w <Return> "destroy $w"
2822
[10]2823    # force the window to stay on top
[100]2824    putontop $w
[10]2825
2826    focus $w.b.2
2827    tkwait window $w
[93]2828    afterputontop
[10]2829
2830    if {$expgui(temp) != ""} return
2831
2832    if {$expgui(orig_backtype) != $expgui(backtype)} {
2833        histinfo $histlist backtype set $expgui(backtype)
[908]2834        RecordMacroEntry "histinfo [list $histlist] backtype set $expgui(backtype)" 0
[48]2835        incr expgui(changed)
[10]2836    }
2837    if {$expgui(orig_backterms) != $expgui(backterms)} {
2838        histinfo $histlist backterms set $expgui(backterms)
[908]2839        RecordMacroEntry "histinfo [list $histlist] backterms set $expgui(backterms)" 0
[48]2840        incr expgui(changed)
[10]2841    }
2842    for {set num 1 } { $num <= $expgui(backterms) } { incr num } {
2843        set var "bterm$num"
2844        if {$expgui(orig_$var) != $expgui($var)} {
2845            histinfo $histlist $var set $expgui($var)
[908]2846            RecordMacroEntry "histinfo [list $histlist] $var set $expgui($var)" 0
[48]2847            incr expgui(changed)
[10]2848        }
2849    }
[908]2850    RecordMacroEntry "incr expgui(changed)" 0
[10]2851
2852    if {$expgui(globalmode) == 0} {
2853        set expgui(backtypelbl) "Function type [histinfo $hist backtype]"
2854        set expgui(backtermlbl) "([histinfo $hist backterms] terms)"
2855    }
2856}
2857
2858trace variable expgui(backterms) w ChangeBackTerms
2859proc ChangeBackTerms {a b c} {
2860    global entrycmd expgui
2861    if !$entrycmd(trace) return
2862    ShowBackTerms .back.1
2863}
2864
2865trace variable expgui(backtype) w ChangeBackType
2866# reset the terms to 1, 0, 0... when the number of terms increase
2867proc ChangeBackType {a b c} {
2868    global entrycmd expgui
2869    if !$entrycmd(trace) return
2870    if {$expgui(prev_backtype) == $expgui(backtype)} return
2871    set expgui(prev_backtype) $expgui(backtype)
2872    for {set num 1 } { $num <= $expgui(backterms) } { incr num } {
2873        set var "bterm$num"
2874        if {$num == 1} {
2875            set expgui($var) 1.0
2876        } else {
2877            set expgui($var) 0.0
2878        }
2879    }
2880}
2881
2882proc ShowBackTerms {w } {
2883    global expgui expmap
2884    # destroy the contents of the frame
[152]2885    eval destroy [winfo children $w]
[10]2886    set histlist {}
2887    foreach n $expgui(curhist) {
2888        lappend histlist [lindex $expmap(powderlist) $n]
2889    }
2890    set widgetsPerRow 4
2891    for {set rows 2; set num 1 } { $num <= $expgui(backterms) } { incr rows } {
2892        for {set cols 0} { (2*$widgetsPerRow > $cols) && ($num <= $expgui(backterms)) }  { incr num }  {
2893            set var "bterm$num"
2894            grid [label $w.l$num -text $num -bg beige]  \
2895                    -row $rows -column $cols -sticky nes
2896            incr cols
2897            grid [entry $w.e$num -width 15 -textvariable expgui($var) \
2898                    ] -row $rows  -column $cols  -sticky news
2899            incr cols
2900        }
2901    }
2902}
2903
2904proc QuitEditBackground {w} {
2905    global expgui
2906    # lets find out if anything changed
2907    set changed 0
2908    if {$expgui(orig_backtype) != $expgui(backtype)} {
2909        set changed 1
2910    }
2911    if {$expgui(orig_backterms) != $expgui(backterms)} {
2912        set changed 1
2913    }
2914    for {set num 1 } { $num <= $expgui(backterms) } { incr num } {
2915        set var "bterm$num"
2916        if {$expgui(orig_$var) != $expgui($var)} {
2917            set changed 1
2918            break
2919        }
2920    }
2921    if $changed {
2922        set decision [tk_dialog .changes "Abandon Changes" \
2923                "You have made changes to the background. Ok to abandon changes?" \
2924                warning 0 Abandon Keep]
2925        if !$decision {
2926            set expgui(temp) "Quit"
2927            destroy $w
2928        }
2929    } else {
2930        set expgui(temp) "Quit"
2931        destroy $w
2932    }
2933}
2934
2935# this probably needs work
2936proc editglobalparm {cmd variable title "histlist {}" "phase {}"} {
2937    global expgui expmap
2938    set w .global
2939    catch {destroy $w}
2940    toplevel $w -bg beige
2941    wm title $w "Edit Global Parameter"
2942    set expgui(temp) {}
2943    if {[llength $histlist] == 0} {
2944        set hist {}
2945        foreach n $expgui(curhist) {
2946            lappend hist [lindex $expmap(powderlist) $n]
2947        }
2948    } else {
2949        set hist $histlist
2950    }
2951    pack [frame $w.0 -bd 6 -relief groove -bg beige] \
2952            -side top -expand yes -fill both
2953    grid [label $w.0.a -text "Setting $title for histograms [CompressList $hist]"\
2954            -bg beige] \
2955            -row 0 -column 0 -columnspan 10
2956    grid [entry $w.0.b -textvariable expgui(temp)] \
2957            -row 1 -column 0 
2958
2959
[323]2960    pack [frame $w.b -bg beige] -fill x -expand yes -side top
[928]2961    pack [button $w.b.2 -text Continue -command "destroy $w"] -side left
2962    pack [button $w.b.3 -text Cancel -command "set expgui(temp) {}; destroy $w"] -side left
[323]2963    pack [button $w.b.help -text Help -bg yellow \
2964            -command "MakeWWWHelp expgui3.html EditParm"] -side right
2965    bind $w <Key-F1> "MakeWWWHelp expgui3.html EditParm"
[10]2966    bind $w <Return> "destroy $w"
2967
[93]2968    # force the window to stay on top
2969    putontop $w
[10]2970    focus $w.b.2
2971    tkwait window $w
[93]2972    afterputontop
2973
[10]2974    if {$expgui(temp) != ""} {
2975        foreach h $hist {
2976            if {$cmd == "histinfo"} {
2977                histinfo $h $variable set $expgui(temp)
[908]2978                RecordMacroEntry "histinfo $h $variable set $expgui(temp)" 1
[48]2979                incr expgui(changed)
[10]2980                if $expgui(debug) {
2981                    puts "histinfo $h $variable set $expgui(temp)"
2982                }
2983            } elseif {$cmd == "hapinfo"} {
2984                hapinfo $h $phase $variable set $expgui(temp)
[908]2985                RecordMacroEntry "hapinfo $h [list $phase] $variable set $expgui(temp)" 1
[48]2986                incr expgui(changed)
[10]2987                if $expgui(debug) {
[908]2988                    puts "hapinfo $h $phase $variable set $expgui(temp)"
[10]2989                }
2990            } else {
2991                error "$cmd unimplemented"
2992            }
2993        }
[908]2994        RecordMacroEntry "incr expgui(changed)" 0
[10]2995    }
2996}
2997
2998proc EditProfile {title histlist phaselist} {
2999    global expgui expmap entrycmd
3000    set w .back
3001    catch {destroy $w}
3002    toplevel $w -bg beige
3003    wm title $w "Global Edit Profile"
3004    set hist [lindex $histlist 0]
3005    set phase [lindex $phaselist 0]
3006    set ptype [string trim [hapinfo $hist $phase proftype]]
3007    set htype [string range $expmap(htype_$hist) 2 2]
3008    set nterms [hapinfo $hist $phase profterms]
3009   
3010    pack [frame $w.0 -bd 6 -relief groove  -bg beige \
3011            ] -side top -expand yes -fill both
3012    grid [label $w.0.a \
3013            -text "Setting profile terms: $title" \
3014            -bg beige] -row 0 -column 0 -columnspan 10
3015    grid [label $w.0.b -text "Function type $ptype"  -bg beige]  -row 1 -column 0
3016    grid [label $w.0.c -text "  Peak cutoff" -bg beige] -row 1 -column 3 
3017    grid [entry $w.0.d -width 10 ]  -row 1 -column 4
3018    set entrylist {}
3019    lappend entrylist "pcut $w.0.d"
3020
3021    set col -1
3022    set row 1
[234]3023    set lbls "dummy [GetProfileTerms $phase $hist $ptype]"
[10]3024    pack [frame $w.1 -bd 6 -relief groove  -bg beige \
3025            ] -side top -expand yes -fill both
3026    for { set num 1 } { $num <= $nterms } { incr num } {
3027        set term {}
3028        catch {set term [lindex $lbls $num]}
3029        if {$term == ""} {set term $num}
3030        incr col
3031        grid [label $w.1.l${num} -text "$term" -bg beige] \
3032                -row $row -column $col
3033        incr col
3034        grid [entry $w.1.ent${num} \
3035                -width 14] -row $row -column $col
3036        lappend entrylist "pterm$num $w.1.ent${num}"   
3037        if {$col > 6} {set col -1; incr row}
3038    }
[323]3039    pack [frame $w.b -bg beige] -fill x -expand yes -side top
[928]3040    grid [button $w.b.2 -text Continue \
[10]3041            -command "SetEditProfile [list $entrylist] [list $phaselist] \
[323]3042            [list $histlist] $w"] -row 0 -column 1
[928]3043    grid [button $w.b.3 -text Cancel \
[323]3044            -command "QuitEditProfile $w [list $entrylist]"] -row 0 -column 2
3045    grid [button $w.b.help -text Help -bg yellow \
3046            -command "MakeWWWHelp expgui5.html GlobalEdit"] \
3047            -row 0 -column 4
3048    grid columnconfig $w.b 0 -weight 1
3049    grid columnconfig $w.b 3 -weight 1
3050    bind $w <Key-F1> "MakeWWWHelp expgui5.html GlobalEdit"
3051    bind $w <Return> "QuitEditProfile $w [list $entrylist]"
[10]3052
[93]3053    # force the window to stay on top
3054    putontop $w
[10]3055    focus $w.b.2
3056    tkwait window $w
[93]3057    afterputontop
[10]3058}
3059
3060proc SetEditProfile {entrylist phaselist histlist w} {
3061    global expgui
3062    foreach item $entrylist {
3063        set value [ [lindex $item 1] get ]
3064        if {$value != ""} {
3065            hapinfo $histlist $phaselist [lindex $item 0] set $value
[908]3066            RecordMacroEntry "hapinfo [list $histlist] [list $phaselist] [lindex $item 0] set $value" 1
[48]3067            incr expgui(changed)
[908]3068            RecordMacroEntry "incr expgui(changed)" 0
[10]3069            if $expgui(debug) {
3070                puts "hapinfo [list $phaselist] [list $histlist] [lindex $item 0] set $value"
3071            }
3072        }
3073    }
3074    destroy $w
3075}
3076
3077proc QuitEditProfile {w entrylist} {
3078    global expgui
3079    # lets find out if anything changed
3080    set changed 0
3081    foreach item $entrylist {
3082        if {[ [lindex $item 1] get ] != ""} {set changed 1; break}
3083    }
3084    if $changed {
3085        set decision [tk_dialog .changes "Abandon Changes" \
3086                "You have made changes to the Profile. Ok to abandon changes?" \
3087                warning 0 Abandon Keep]
3088        if !$decision {destroy $w}
3089    } else {
3090        destroy $w
3091    }
3092}
3093
[557]3094# this is called to change the absorption correction mode and to
3095# change the absorption correction model.
3096proc editabsorption {} {
3097    global expgui expmap
3098    set histlist {}
3099    foreach n $expgui(curhist) {
3100        lappend histlist [lindex $expmap(powderlist) $n]
3101    }
3102    if {[llength $histlist] == 0} return
3103
3104    set w .abs
3105    catch {destroy $w}
3106    toplevel $w -bg beige
3107    if {$expgui(globalmode) != 0} {
3108        wm title $w "Global Edit Absorption/Reflectivity" 
3109    } else {
3110        wm title $w "Edit Absorption/Reflectivity"
3111    }
3112   
3113    pack [frame $w.0 -bd 6 -relief groove  -bg beige \
3114            ] -side top -expand yes -fill both
3115    if {[llength $histlist] > 1} {
3116        grid [label $w.0.a \
3117            -text "Changing settings for histograms [CompressList $histlist]" \
3118            -bg beige] -row 0 -column 0 -columnspan 10
3119    } else {
3120        grid [label $w.0.a \
3121            -text "Changing settings for histogram $histlist" \
3122            -bg beige] -row 0 -column 0 -columnspan 4
3123        #grid columnconfig $w.0 4 -weight 1
3124    }
3125    grid rowconfig $w.0 1 -min 10
3126    set hist [lindex $histlist 0]
3127
3128    grid [label $w.0.lb1 -text "Absorption Coefficient(s)" -bg beige] \
3129            -row 2 -column 1  -columnspan 2
3130    grid [label $w.0.lb1a -text "1" -bg beige] -row 3 -column 1
3131    set expgui(abs2box1) $w.0.lb2a
3132    grid [label $w.0.lb2a -text "2" -bg beige] -row 3 -column 2
3133    grid [label $w.0.lb3 -text Absorption\nFunction -bg beige] \
3134            -row 2 -column 6 -rowspan 2 -columnspan 2
3135    grid [entry $w.0.ent1 -textvariable expgui(abscor1) -width 15] \
3136            -row 4 -column 1
3137    set expgui(abs2box2) $w.0.ent2
3138    grid [entry $w.0.ent2 -textvariable expgui(abscor2) -width 15] \
3139            -row 4 -column 2 
3140    trace vdelete expgui(abstype) w AbsSetoptmsg
3141    eval tk_optionMenu $w.0.m1 expgui(abstype) 0 1 2 3 4
3142    trace variable expgui(abstype) w AbsSetoptmsg
3143    grid $w.0.m1 -row 4 -column 6 -columnspan 2
3144    grid [label $w.0.lb8 -textvariable expgui(opttxt) -bg beige \
3145          -wrap 300 -justify left] -row 5 -column 1  -sticky ne -columnspan 7
3146    grid rowconfig $w.0 5 -min 100
3147    # set the values, note the trace on abstype
3148    foreach var {abscor1 abscor2 abstype} {
3149        set expgui($var) [histinfo $hist $var]
3150    }
3151
3152    pack [frame $w.b -bg beige] -fill x -expand yes -side top
[928]3153    grid [button $w.b.2 -text Continue -command "AbsSaveEdit $w [list $histlist]"] \
[557]3154            -row 0 -column 1
[928]3155    grid [button $w.b.3 -text Cancel \
[557]3156            -command "destroy $w"] -row 0 -column 2
3157    grid [button $w.b.help -text Help -bg yellow \
3158            -command "MakeWWWHelp expgui3.html EditAbsorption"] \
3159            -row 0 -column 4
3160    grid columnconfig $w.b 0 -weight 1
3161    grid columnconfig $w.b 3 -weight 1
3162    bind $w <Key-F1> "MakeWWWHelp expgui3.html EditAbsorption"
3163    bind $w <Return> "destroy $w"
3164
3165    # force the window to stay on top
3166    putontop $w
3167
3168    focus $w.b.2
3169    tkwait window $w
3170    afterputontop
3171}
3172
3173proc AbsSetoptmsg {args} {
3174    global expgui
3175    array set opttxt {
3176        0 "Correction for cylindrical samples [Lobanov & Alte da Veiga]. OK for all data types, but not for Bragg-Brentano flat-plate geometry. Set term 1 to mu*R/lambda (TOF: mu*R for lambda=1). For CW x-ray/neutron, do not refine!"
3177        1 "Wavelength-dependent correction for container penetration. Use with TOF & Energy Disp x-ray only."
3178        2 "Surface roughness correction [Pitschke, Hermann & Muttern]. Use with flat-plate reflection geometry (usually Bragg-Brentano) only."
3179        3 "Surface roughness correction, [Suortti]. Use with flat-plate reflection geometry (usually Bragg-Brentano) only."
3180        4 "Flat plate samples in transmission mode. OK for all data types, but not Bragg-Brentano geometry. Term 2 is angle w/r to beam (usually 0). For CW, do not refine."
3181    }
3182    set expgui(opttxt) ""
3183    catch {set expgui(opttxt) [set opttxt($expgui(abstype))]}
3184    switch $expgui(abstype) {
3185        0 -
3186        1 {
3187            $expgui(abs2box1) config -fg gray
3188            $expgui(abs2box2) config -state disabled -fg gray
3189        } 
3190        2 -
3191        3 -
3192        4 {
3193            $expgui(abs2box1) config -fg black
3194            $expgui(abs2box2) config -state normal -fg black
3195        }
3196        default {
3197            set expgui(opttxt) "Please select an absorption function"
3198        }
3199    }
3200}
3201proc AbsSaveEdit {top histlist} {
3202    global expgui expmap
3203    # sanity check: look at the histogram type
3204    set h [lindex $histlist 0]
3205    if {[string range $expmap(htype_$h) 2 2] == "T"} {set flag 1}
3206    if {[string range $expmap(htype_$h) 1 2] == "NC"} {set flag 2}
3207    if {[string range $expmap(htype_$h) 1 2] == "XC" && \
3208            [histinfo $h lam2] != 0.0} {set flag 3}
3209    if {[string range $expmap(htype_$h) 1 2] == "XC" && \
3210            [histinfo $h lam2] == 0.0} {set flag 4}
3211    if {[string range $expmap(htype_$h) 1 2] == "XE"} {set flag 5}
3212
3213    set msg {}
3214    if {$expgui(abstype) == 0 && ($flag == 3 || $flag == 4)} {
3215        set msg "Mode 0 is appropriate for cylindrical (Debye-Scherrer) geometry only"
3216    } elseif {$expgui(abstype) == 1 && ($flag != 1 && $flag != 5)} {
3217        set msg "Mode 1 is appropriate for wavelength-dispersive (TOF/E.D. X-ray) data only"
3218    } elseif {($expgui(abstype) == 2 || $expgui(abstype) == 3) \
3219            && $flag != 3 && $flag != 4} {
3220        set msg "Mode 1 is appropriate for reflection geometry flat-plate (typically Bragg-Brentano) data only"
3221    } elseif {$expgui(abstype) == 4 && $flag <= 3} {
3222        set msg "Mode 4 is appropriate for flat-plate samples in transmission"
3223    }
3224    if {$msg != ""} {
3225        set result [\
3226                MyMessageBox -parent $top -title "Sanity check" \
3227                -type okcancel -default cancel \
3228                -icon warning -helplink "expgui3.html AbsorptionSanity" \
3229                -message "$msg  -- are you sure you want to do this?"]
3230        if {$result == "cancel"} return
3231    }
3232
3233    # validate abscor1 & abscor2 (if needed)
3234    set msg {}
3235    if {![validreal expgui(abscor1) 15 8]} {
3236        set msg "Term 1 is invalid"
3237    }
3238    if {$expgui(abstype) > 1} {
3239        if {![validreal expgui(abscor2) 15 8]} {
3240            if {$msg != ""} {append msg "\n"}
3241            append msg "Term 2 is invalid"
3242        }
3243    }
3244    if {$msg != ""} {
3245        MyMessageBox -parent $top -title "Entry error" \
3246                -type ok -default ok \
3247                -icon warning -helplink "" \
3248                -message "Invalid data entered. Please correct.\n$msg"
3249        return
3250    }
3251   
3252    histinfo $histlist abstype set $expgui(abstype)
[908]3253    RecordMacroEntry "histinfo [list $histlist] abstype set $expgui(abstype)" 0
[557]3254    histinfo $histlist abscor1 set $expgui(abscor1)