I have been playing around with SB to try and get used to some of its features, and have been having fun in the process, which I guess is the main objective! I have written some messy code which I am having some problems with. I'll try and explain!
Basically, the code randomly picks a locations for 2 radioactive sources. It then scans from top to bottom of the screen, calculating the activity/dose that would be observed at each point over an area of 400 x 500 pixels, and creates what I am calling a "dose map". An example of this is shown below.
- IMG_DF24B0C7DA1A-1.jpeg (118.26 KiB) Viewed 3489 times
The programme then generates a cursor (rectangle) controlled by 4 movement boxes (left, right, up and down) to allow the user to "walk" over the dose map. The code then calculates the dose at the cursor position and displays the combined dose.
However, there is a problem!
When you move the cursor over the radiation "hot spots", you'd expect to get the highest dose rate. However, this is not the case. You actually get the highest rate at what appears to be twice the x and y position that you'd expect. You can see this by moving the cursor to the right and below a hotspot, assuming its been randomly generated sufficiently close to x=0 and y=0 to allow you to move the cursor far enough away!
I think I have the maths code correct, but wondered if there was a problem with my saving and pasting of the radiation map.
I'm a little confused.
Any thoughts very much appreciated.
All the best,
Paul.
Oh, and to exit the program once its running, tap the far right and bottom of the screen!
graphics
graphics clear 1,1,1
refresh on
set toolbar off
x=0
y=0
source1=rnd(1000000)
source2=rnd(1000000)
' Choose random location for radiation source 1
a1=rnd(320)
b1=rnd(524)
' Choose random location for radiation source 2
a2=rnd(320)
b2=rnd(524)
draw color 0,0,0
loop:
' use Pythag theory to calculate the distance from the x,y position to source 1
distx1=abs(a1-x)
disty1=abs(b1-y)
shinepath1=sqrt((distx1*distx1)+(disty1*disty1))
' use Pythag theory to calculate the distance from the x,y position to source 2
distx2=abs(a2-x)
disty2=abs(b2-y)
shinepath2=sqrt((distx2*distx2)+(disty2*disty2))
' prevents a divide by 0 error
if shinepath1=0 then
shinepath1=1
end if
' prevents a divide by 0 error
if shinepath2=0 then
shinepath2=1
end if
' calculate the dose from source 1 and source 2, based on inverse square law
dose1=(source1*(1/(shinepath1*shinepath1)))
dose2=(source2*(1/(shinepath2*shinepath2)))
' total dose at position x,y is the sum of the dose from each source
dosetotal=dose1+dose2
' Determine what colour to plot the total dose, depending on its magnitude
if dosetotal >0 and dosetotal=<50 then
draw color 0,0,0
endif
if dosetotal >50 and dosetotal=<100 then
draw color 1,0.2,0.2
endif
if dosetotal >100 and dosetotal=<150 then
draw color 1,0.4,0.4
endif
if dosetotal >150 and dosetotal=<200 then
draw color 1,0.6,0.6
endif
if dosetotal>200 and dosetotal=<250 then
draw color 1,0.7,0.7
endif
if dosetotal>250 and dosetotal=<300 then
draw color 1,0.8,0.8
endif
' draw the dose at position x,y
draw pixel x,y
' once 500 lines in the Y axis has been reach move onto the section of code that allows the cursor to be moved around the zone
if y>500 then
goto survey
endif
' increment x until 400 pixels has been reached, then increment Y value
x=x+1
if x>400 then
y=y+1
x=1
end if
' continue scanning the whole area
goto loop
' this section of code allows the user to move a cursor around the area to determine the dose at each point
survey:
' save the "dose map" to an image. This prevents having to rec-calculate the dose map each time the cursor is moved
graphics save 0,0,200,250 to "dosemap.jpg"
' draw 4 buttons to allow the user to move a cursor around the area
fill color 0,0.5,0.5
fill rect 100,250 to 150,300
fill rect 150,330 to 200,380
fill rect 50,330 to 100,380
fill rect 100,410 to 150,460
cursmove:
' draw the dose map to the screen
draw image "dosemap.jpg" at 0,0 scale 1 angle 0
' see which direction button has been pressed by the user
GET TOUCH 0 AS xx,yy
if xx=> 100 and xx=<150 and yy=>250 and yy=<300 then cursy=cursy-1
if xx=>150 and xx=<200 and yy=>330 and yy=<380 then cursx=cursx+1
if xx=>50 and xx=<100 and yy=>330 and yy=<380 then cursx=cursx-1
if xx=>100 and xx=<150 and yy=>410 and yy=<460 then cursy=cursy+1
' Stop the cursor from going off the edge of the dosemap
if cursx>=190 then cursx=190
if cursx<=0 then cursx=0
if cursy>=240 then cursy=240
if cursy<=0 then cursy=0
' exit if screen touch on far right of screen
if xx=>300 and yy=>440 then end
'calculate the dose at the cursor location
adistx1=(abs(a1)-(cursx))
adisty1=(abs(b1)-(cursy))
ashinepath1=sqrt((adistx1*adistx1)+(adisty1*adisty1))
adistx2=(abs(a2)-(cursx))
adisty2=(abs(b2)-(cursy))
ashinepath2=sqrt((adistx2*adistx2)+(adisty2*adisty2))
if ashinepath1=0 then
ashinepath1=1
end if
if ashinepath2=0 then
ashinepath2=1
end if
dosea=(source1*(1/(ashinepath1*ashinepath1)))
doseb=(source2*(1/(ashinepath2*ashinepath2)))
dosetotalatcurs=dosea+doseb
' draw cursor
fill color 0,1,0
fill rect cursx,cursy to cursx+10,cursy+10
refresh on
fill color 1,1,1
fill rect 20,490 to 200,510
draw text dosetotalatcurs at 20,495
goto cursmove
