Changeset 227


Ignore:
Timestamp:
Dec 4, 2009 5:02:31 PM (11 years ago)
Author:
toby
Message:

# on 2000/07/03 21:50:35, toby did:
fix aniso atoms in exp2xtl

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/gsascmds.tcl

    • Property rcs:date changed from 2000/06/09 03:52:16 to 2000/07/03 21:50:35
    • Property rcs:lines changed from +7 -1 to +16 -39
    • Property rcs:rev changed from 1.19 to 1.20
    r201 r227  
    402402        }
    403403       
    404         # are there anisotropic atoms?
    405         set aniso 0
     404        puts $fp "ATOMS"
     405        puts $fp "NAME       X          Y          Z    UISO      OCCUP"
    406406        foreach atom $expmap(atomlist_$phase) {
    407             if {[atominfo $phase $atom temptype] == "A"} {set aniso 1}
    408         }
    409         puts $fp "ATOMS"
    410         if $aniso {
    411             puts $fp "NAME       X          Y          Z    OCCUP U11 U22 U33 U12 U13 U23"
    412             foreach atom $expmap(atomlist_$phase) {
    413                 set label [atominfo $phase $atom label]
    414                 # remove () characters
    415                 if {[atominfo $phase $atom temptype] == "A"} {
    416                     puts $fp "$label [atominfo $phase $atom x] \
    417                             [atominfo $phase $atom y] [atominfo $phase $atom z] \
    418                             [atominfo $phase $atom frac] \
    419                             [atominfo $phase $atom U11] \
    420                             [atominfo $phase $atom U22] \
    421                             [atominfo $phase $atom U33] \
    422                             [atominfo $phase $atom U12] \
    423                             [atominfo $phase $atom U13] \
    424                             [atominfo $phase $atom U23]"
    425                 } else {
    426                     puts $fp "$label [atominfo $phase $atom x] \
    427                             [atominfo $phase $atom y] [atominfo $phase $atom z] \
    428                             [atominfo $phase $atom frac] \
    429                             [atominfo $phase $atom Uiso] \
    430                             [atominfo $phase $atom Uiso] \
    431                             [atominfo $phase $atom Uiso] \
    432                             0 0 0 "
    433                 }
    434             }
    435         } else {
    436             puts $fp "NAME       X          Y          Z    UISO      OCCUP"
    437             foreach atom $expmap(atomlist_$phase) {
    438                 set label [atominfo $phase $atom label]
    439                 # remove () characters
    440                 regsub -all "\[()\]" $label "" label
    441                 puts $fp "$label [atominfo $phase $atom x] \
     407            set label [atominfo $phase $atom label]
     408            # remove () characters
     409            regsub -all "\[()\]" $label "" label
     410            # are there anisotropic atoms?
     411            if {[atominfo $phase $atom temptype] == "A"} {
     412                set uiso [expr \
     413                        ([atominfo $phase $atom U11] + \
     414                        [atominfo $phase $atom U22] + \
     415                        [atominfo $phase $atom U33]) / 3.]
     416            } else {
     417                set uiso [atominfo $phase $atom Uiso]
     418            }
     419            puts $fp "$label [atominfo $phase $atom x] \
    442420                        [atominfo $phase $atom y] [atominfo $phase $atom z] \
    443                         [atominfo $phase $atom Uiso]  [atominfo $phase $atom frac]"
    444             }
     421                        $uiso  [atominfo $phase $atom frac]"
    445422        }
    446423    } errmsg] {
Note: See TracChangeset for help on using the changeset viewer.