 Timestamp:
 Oct 17, 2010 10:32:31 PM (12 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/sandbox/widplt
r953 r1029 563 563 } 564 564 } 565 # assume FWHM add as square roots 565 # Use polynomial of P. Thompson, D.E. Cox & J.B. Hastings, 566 # J. Appl. Cryst.,20,7983, 1987 (GSAS manual 9/26/04 p157) 567 # to compute composite FWHM 566 568 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 }] 568 579 } 569 }580 } 570 581 if $lflag { 571 582 catch { … … 938 949 set tantt [expr tan($tt * $torad ) ] 939 950 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 944 958 proc LFWHM {tt X Y} { 945 959 set pi 3.14159
Note: See TracChangeset
for help on using the changeset viewer.