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*(1-phi)*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[num-1]) 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 slit-smearing 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 q-range %g to %g is %g ^-3 cm^-1\r\r",qw[0],qw[(num-1)],inv334 Printf "The invariant over the measured q-range %g to %g is %g A^-3 cm^-1\r\r",qw[0],qw[(num-1)],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 power-law fit instead … … 400 400 401 401 nume=numpnts(extr_lqq) 402 Printf "The invariant over the q-range %g to %g is %g ^-3 cm^-1\r\r",extr_lqq[0],extr_lqq[(nume-1)],inv402 Printf "The invariant over the q-range %g to %g is %g A^-3 cm^-1\r\r",extr_lqq[0],extr_lqq[(nume-1)],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 q-range %g to %g is %g ^-3 cm^-1\r\r",extr_hqq[0],extr_hqq[(nume-1)],inv452 Printf "The invariant over the q-range %g to %g is %g A^-3 cm^-1\r\r",extr_hqq[0],extr_hqq[(nume-1)],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.