Changeset 379 for sans/Dev/trunk/NCNR_User_Procedures/SANS/Analysis/Models/Packages/Invariant/Invariant_v40.ipf
 Timestamp:
 Jun 10, 2008 4:00:21 PM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

sans/Dev/trunk/NCNR_User_Procedures/SANS/Analysis/Models/Packages/Invariant/Invariant_v40.ipf
r318 r379 49 49 inv = 2*pi*pi*delta_rho*delta_rho*phi*(1phi)*1e8 50 50 51 Printf "The model invariant is %g ^3 cm^1\r\r",inv51 Printf "The model invariant is %g A^3 cm^1\r\r",inv 52 52 End 53 53 … … 65 65 invar = areaXY(qw,integrand,qw[0],qw[num1]) 66 66 67 return(invar) //units of ^3 cm^167 return(invar) //units of A^3 cm^1 68 68 End 69 69 … … 85 85 invar *= dQv //correct for the effects of slitsmearing 86 86 87 return(invar) //units of ^3 cm^187 return(invar) //units of A^3 cm^1 88 88 End 89 89 … … 140 140 141 141 Printf "I(q=0) = %g (1/cm)\r",G_coef[0] 142 Printf "Rg = %g ( )\r",sqrt(3*G_coef[1])142 Printf "Rg = %g (A)\r",sqrt(3*G_coef[1]) 143 143 144 144 Make/O/D P_coef={0,1,4} //input … … 191 191 ModifyGraph standoff=0 192 192 Label left "Intensity (1/cm)" 193 Label bottom "q (1/ )"193 Label bottom "q (1/A)" 194 194 ModifyGraph rgb($NameofWave(iw))=(0,0,0) 195 195 ModifyGraph opaque($NameofWave(iw))=1 … … 251 251 DrawText 158,211,"Extrapolation" 252 252 SetDrawEnv fstyle= 1 253 DrawText 58,349,"Units are [ ^3 cm^1]"253 DrawText 58,349,"Units are [A^3 cm^1]" 254 254 SetDrawEnv fstyle= 1 255 255 DrawText 61,16,"Calculate the Invariant" … … 332 332 333 333 num=numpnts(qw) 334 Printf "The invariant over the measured qrange %g to %g is %g ^3 cm^1\r\r",qw[0],qw[(num1)],inv334 Printf "The invariant over the measured qrange %g to %g is %g A^3 cm^1\r\r",qw[0],qw[(num1)],inv 335 335 336 336 // update the global display on the panel (there is a dependency for the total) … … 378 378 379 379 Printf "I(q=0) = %g (1/cm)\r",G_coef[0] 380 Printf "Rg = %g ( )\r",sqrt(3*G_coef[1])380 Printf "Rg = %g (A)\r",sqrt(3*G_coef[1]) 381 381 else 382 382 //do a powerlaw fit instead … … 400 400 401 401 nume=numpnts(extr_lqq) 402 Printf "The invariant over the qrange %g to %g is %g ^3 cm^1\r\r",extr_lqq[0],extr_lqq[(nume1)],inv402 Printf "The invariant over the qrange %g to %g is %g A^3 cm^1\r\r",extr_lqq[0],extr_lqq[(nume1)],inv 403 403 404 404 // update the global display on the panel (there is a dependency for the total) … … 450 450 451 451 nume=numpnts(extr_hqq) 452 Printf "The invariant over the qrange %g to %g is %g ^3 cm^1\r\r",extr_hqq[0],extr_hqq[(nume1)],inv452 Printf "The invariant over the qrange %g to %g is %g A^3 cm^1\r\r",extr_hqq[0],extr_hqq[(nume1)],inv 453 453 454 454 // update the global display on the panel (there is a dependency for the total)
Note: See TracChangeset
for help on using the changeset viewer.