Ignore:
Timestamp:
Feb 29, 2012 11:36:57 AM (11 years ago)
Author:
srkline
Message:

fixing the shape filling functions to explicitly take fill as a parameter, and make the default fill value = 10 == 10 e-7 A-2

More additions to the RealSpeca? help file

File:
1 edited

Legend:

Unmodified
Added
Removed
  • sans/Dev/trunk/NCNR_User_Procedures/Analysis/Alpha/Tinker/FFT_Cubes_Includes.ipf

    r838 r840  
    9393End 
    9494 
    95 Proc PutXAxisCoreShellCyl_HexGrid(w,rad,len,sep) 
     95Proc PutXAxisCoreShellCyl_HexGrid(w,rad,len,sep,fill) 
    9696        String w="mat" 
    97         Variable rad=30,len=300,sep=100 
     97        Variable rad=30,len=300,sep=100,fill=10 
    9898        Prompt w,"matrix" 
    9999        prompt rad,"radius of cylinders" 
    100100        prompt len,"length of cylinders" 
    101101        prompt sep,"center-to-center separation of cylinders" 
     102        Prompt fill,"fill SLD value" 
    102103         
    103         $w=0 
    104         X_CoreShellCylinderHexGrid($w,rad,len,sep) 
     104//      $w=0 
     105        X_CoreShellCylinderHexGrid($w,rad,len,sep,fill) 
    105106         
    106107        NumberOfPoints() 
    107108end 
    108109 
    109 Function X_CoreShellCylinderHexGrid(mat,rad,len,sep) 
     110Function X_CoreShellCylinderHexGrid(mat,rad,len,sep,fill) 
    110111        Wave mat 
    111112        variable rad,len                //length of cylinders 
    112113        Variable sep                    // EDGE separation, in same units as cylinder 
     114        Variable fill 
    113115         
     116        NVAR    solventSLD = root:FFT_SolventSLD 
    114117         
    115118        NVAR grid=root:FFT_T 
     
    119122        // fill a 2D plane with points 
    120123        Make/O/B/N=(np,np) plane 
    121         plane = 0 
     124        plane = solventSLD 
    122125 
    123126        spacing = round(sep/grid)               // so it's an integer 
    124         FillPlaneHexagonal(plane,spacing) 
     127        FillPlaneHexagonal(plane,spacing,fill) 
    125128         
    126129        // put it in the proper plane of the matrix 
     
    136139         
    137140        for(ii=0;ii<num;ii+=1) 
    138                 FillXCylinder(mat,grid,rad,x3d[ii],y3d[ii],z3d[ii],len,1)               //cylinder 1 
     141                FillXCylinder(mat,grid,rad,x3d[ii],y3d[ii],z3d[ii],len,fill)            //cylinder 1 
    139142        endfor 
    140143 
    141144// makes a crude core-shell cylinder     
    142145        for(ii=0;ii<num;ii+=1) 
    143                 FillXCylinder(mat,grid,rad-20,x3d[ii],y3d[ii],z3d[ii],len,3)            //cylinder 1 
     146                FillXCylinder(mat,grid,rad-20,x3d[ii],y3d[ii],z3d[ii],len,3*fill)               //cylinder 1 
    144147        endfor 
    145148         
Note: See TracChangeset for help on using the changeset viewer.