Changeset 1028 for trunk


Ignore:
Timestamp:
Oct 17, 2010 10:32:05 PM (10 years ago)
Author:
toby
Message:

fix FWHM calc

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/widplt

    r953 r1028  
    563563                            }
    564564                        }
    565                         # assume FWHM add as square roots
     565                        # Use polynomial of P. Thompson, D.E. Cox & J.B. Hastings, 
     566                        # J. Appl. Cryst.,20,79-83, 1987 (GSAS manual 9/26/04 p157)
     567                        # to compute composite FWHM
    566568                        lappend tfwhmlist \
    567                                 [expr sqrt($gfwhm*$gfwhm + $lfwhm*$lfwhm)]
     569                            [expr {
     570                                   pow(
     571                                       pow($gfwhm,5) +
     572                                       2.69269*pow($gfwhm,4)*$lfwhm +
     573                                       2.42843*pow($gfwhm,3)*pow($lfwhm,2) +
     574                                       4.47163*pow($gfwhm,2)*pow($lfwhm,3) +
     575                                       0.07842*$gfwhm*pow($lfwhm,4) +
     576                                       pow($lfwhm,5)
     577                                       , 0.2)
     578                               }]
    568579                    }
    569                 }
     580                        }
    570581                if $lflag {
    571582                    catch {
     
    938949    set tantt [expr tan($tt * $torad ) ]
    939950    set costt [expr cos($tt * $torad ) ]
    940     return [expr sqrt \
    941             (8.* log(2) * ($U * $tantt * $tantt + $V * $tantt + $W \
    942             + $P / ($costt * $costt))) / 100.]
    943 }
     951    return [expr {
     952                  sqrt (8.* log(2) *
     953                        ($U * $tantt * $tantt + $V * $tantt + $W +
     954                         $P / ($costt * $costt))
     955                        ) / 100.}]
     956}
     957
    944958proc LFWHM {tt X Y} {
    945959    set pi 3.14159
Note: See TracChangeset for help on using the changeset viewer.