| Topic: |
Science > Physics |
| User: |
"Douglas Eagleson" |
| Date: |
15 Apr 2007 10:10:12 AM |
| Object: |
A Radiation Shielding Qbasic Calculator |
Here is some source code for gamma ray shielding. It is good for
0.03Mev-3 Mev
It was reviewed for consistency, but not absolute accuracy. So it is
only for evaluation purposes.
The idea was to automate the basic , point, line, disk sources.
Just get the gamma ray energy and yield from a reference and plug in
the inputs!
Microshield the commercial product is the thing to use for real work,
but soon I will compare the two. Mine is free. Mine has 3 geometries,
Microshield has 15. But I once made a set of line sources to simulate
all geometries. And min ehas no fancy graphics, just old time command
line inputs/outputs.
I can email a executable for Windows-xp if anybody is interested. Also
google goofs-up the long lined data sets. So be careful if
reformating for a local qbasic run.
********begin real
code**********************************************************
REM DOSECALC IS A PROGRAM TO ASSIST IN GAMMA RAY SAFETY.
REM IT CALCULATES POINT SOURCE DOSE RATE AND LINE SOURCE DOSE
REM RATE FOR THE SELECTED SHIELDING MATERIALS.
REM VERSION 1.0a
REM BY: DOUGLAS EAGLESON
REM VERSION 1.01a WILL BE COMING OUT SOON.
REM THE ONLY VALIDATION IS CURSORY.
REM SO YOU GET WHAT YOU PAY FOR.
REM I TRIED HARD TO GET ACCURATE TABLES.
REM AND THE BUILDUP COEFFICIENTS TABLES WERE
REM GRAPHICALLY DETERMINED. THE RAD HEALTH
REM HANDBOOK SUPPLIED THE GRAPHS.
REM IT IS HARD TO PROGRAM THIS STUFF
REM ANY CLEARER.
DIM BD(6)
DIM I(6)
DIM T(3, 13)
DIM A(3, 13)
DIM L(3, 13)
DIM F(3, 13)
DIM W(3, 13)
DIM C(3, 13)
DIM D(13)
DIM M(18)
DIM U(1 TO 7, 1 TO 18)
RESTORE DLIST
FOR I = 1 TO 13
READ D(I)
NEXT I
RESTORE MLIST
FOR I = 1 TO 18
READ M(I)
NEXT I
RESTORE BDLIST
FOR I = 1 TO 6
READ BD(I)
NEXT I
RESTORE ILIST
FOR I = 1 TO 6
READ I(I)
NEXT I
RESTORE TLIST
FOR I = 1 TO 13
READ T(1, I)
NEXT I
FOR I = 1 TO 13
READ T(2, I)
NEXT I
FOR I = 1 TO 13
READ T(3, I)
NEXT I
RESTORE ALIST
FOR I = 1 TO 13
READ A(1, I)
NEXT I
FOR I = 1 TO 13
READ A(2, I)
NEXT I
FOR I = 1 TO 13
READ A(3, I)
NEXT I
RESTORE LLIST
FOR I = 1 TO 13
READ L(1, I)
NEXT I
FOR I = 1 TO 13
READ L(2, I)
NEXT I
FOR I = 1 TO 13
READ L(3, I)
NEXT I
RESTORE FLIST
FOR I = 1 TO 13
READ F(1, I)
NEXT I
FOR I = 1 TO 13
READ F(2, I)
NEXT I
FOR I = 1 TO 13
READ F(3, 13)
NEXT I
RESTORE WLIST
FOR I = 1 TO 13
READ W(1, I)
NEXT I
FOR I = 1 TO 13
READ W(2, I)
NEXT I
FOR I = 1 TO 13
READ W(3, I)
NEXT I
RESTORE CLIST
FOR I = 1 TO 13
READ C(1, I)
NEXT I
FOR I = 1 TO 13
READ C(2, I)
NEXT I
FOR I = 1 TO 13
READ C(3, I)
NEXT I
RESTORE ULIST
FOR I = 1 TO 18
READ U(1, I)
NEXT I
FOR I = 1 TO 18
READ U(2, I)
NEXT I
FOR I = 1 TO 18
READ U(3, I)
NEXT I
FOR I = 1 TO 18
READ U(4, I)
NEXT I
FOR I = 1 TO 18
READ U(5, I)
NEXT I
FOR I = 1 TO 18
READ U(6, I)
NEXT I
FOR I = 1 TO 18
READ U(7, I)
NEXT I
P$(1) = "1-LEAD"
P$(2) = "2-CONCRETE"
P$(3) = "3-WATER"
P$(4) = "4-IRON"
P$(5) = "5-TISSUE"
P$(6) = "6-AIR"
P$(7) = "7-NO SHIELD"
GA$(1) = "1= POINT SOURCE"
GA$(2) = "2= LINE SOURCE"
GA$(3) = "3= DISK SOURCE"
GA$(4) = "0= Restart from VERY Beginning"
70 DIM U$
75 DIM H(17, 5)
80 DIM N(19, 7)
85 DIM R$(6)
90 DIM E$(6)
91 DIM E(6)
95 DIM K$(6)
96 DIM K(6)
100 DIM S$(6)
101 DIM S(6)
105 DIM Q$(6)
106 DIM Q(6)
107 DIM GB(3)
108 DIM X$(7)
110 DIM X(7)
115 DIM Y(7)
116 DIM Y$(7)
120 DIM Z(7)
121 DIM Z$(7)
122 DIM O$(6)
123 DIM M$(6)
REM THE FIRST SECTION READS IN THE GAMMA ENERGIES
REM AND GAMMA YIELDS. IT THEN GOES TO THE INTERACTION
REM COEFFICIENT ARRAYS TO LOOKUP THE CORRECT VALUES.
REM THE DATA IS IN DIM 18 ARRAYS FOR THE SELECTED MATERIALS.
124 PRINT "DOSECALC vs1.0"
125 PRINT " "
130 PRINT "ENTER ISOTOPE ID -OR- EXIT"
135 INPUT R$
136 IF R$ = "EXIT" GOTO 2170
137 IF R$ = "exit" GOTO 2170
140 CLS
141 OPEN "C:\DOSECALC\output.txt" FOR APPEND AS #1
145 PRINT "ENTER: GAMMA ENERGY(MEV) INTENSITY(0.0)"
150 PRINT "INPUT FORM- E1 'ENTER' I1 'ENTER'"
155 PRINT "MAXIMUM OF SIX ENERGIES"
156 PRINT "ENTER A ZERO TO END THE INPUT"
157 PRINT " "
160 PRINT "ENERGY(MEV)"
161 PRINT "INTENSITY"
165 FOR I = 1 TO 6
170 INPUT E(I)
171 LET E$(I) = STR$(E(I))
175 IF E(I) = 0 THEN GOTO 225
180 INPUT K(I)
181 LET K$(I) = STR$(K(I))
185 PRINT E$(I); " "; K$(I)
190 NEXT I
195 LET I = 0
REM THE SHIELD MATERIAL IS SELECTED AND ARRAY LOOKUP FOR THE SHIELD
REM IS DONE. TISSUE IS LOOKED UP ALSO.
225 PRINT "SELECT SHIELD MATERIAL, ENTER NUMBER"
230 FOR I = 1 TO 7
235 PRINT " "; P$(I)
240 NEXT I
245 INPUT B
250 IF B = 7 THEN GOTO 295
255 FOR I = 1 TO 6
260 IF E(I) = 0 THEN GOTO 295
265 FOR J = 1 TO 18
270 IF E(I) <= U(7, J) THEN GOTO 280
275 NEXT J
280 IF J = 1 OR J = 19 THEN GOTO 340
REM SHIELD LOOKUP AND INTERPOLATION
285 S(I) = U(B, J - 1) + ((E(I) - U(7, J - 1)) / (U(7, J) - U(7, J -
1))) * (U(B, J) - U(B, J - 1))
290 NEXT I
295 FOR I = 1 TO 6
300 IF E(I) = 0 THEN GOTO 346
305 FOR J = 1 TO 18
310 IF E(I) <= U(7, J) THEN GOTO 320
315 NEXT J
320 IF J = 1 OR J = 18 THEN GOTO 340
REM TISSUE LOOKUP AND INTERPOLATION
325 LET Q(I) = M(J - 1) + ((E(I) - U(7, J - 1)) / (U(7, J) - U(7, J -
1))) * (M(J) - M(J - 1))
330 NEXT I
335 GOTO 346
340 PRINT "ENERGY OUT OF RANGE"
346 PRINT " "
REM THE INITIAL DATA IS PRINTED OUT.
347 PRINT "DOSECALC vs2.10"
348 PRINT " "
350 PRINT "ISOTOPE="; R$
355 PRINT "SHIELD="; P$(B)
360 PRINT "E(MEV) INT U-S U-T"
REM 6092 PRINT " TOTAL="; : PRINT USING " #.#####^^^^"; AL; :
PRINT USING " #.#####^^^^ "; EL
365 FOR I = 1 TO 6
370 IF E(I) = 0 THEN GOTO 385
371 LET E$(I) = STR$(E(I))
372 LET K$(I) = STR$(K(I))
373 LET S$(I) = STR$(S(I))
374 LET Q$(I) = STR$(Q(I))
REM 375 PRINT E$(I); " "; K$(I); " "; S$(I); " "; Q$(I)
375 PRINT USING "#.### "; E(I); : PRINT USING " #.## "; K(I); : PRINT
USING " #.#####^^^^ "; S(I); : PRINT USING " #.#####^^^^ "; Q(I)
380 NEXT I
385 PRINT "U-S=SHIELD LINEAR ATTEN. COEFF. (CM-1)"
390 PRINT "U-T=TISSUE MASS ABSORB. COEFF.(CM2/G)"
391 PRINT " "
410 PRINT "SELECT GEOMETRY- For Selected Shield ID"
REM THE THIRD CHOICE IS TO BE THE DISK SOURCE.
REM I AM STILL TRYING TO GET A PUBLISHED
REM REFERENCE FOR THE EXACT EQUATIONS.
REM MANY REFERENCES ARE UNCOMPREHENSIVE.
415 FOR I = 1 TO 4
420 PRINT GA$(I)
425 NEXT I
430 INPUT C
434 FOR I = 1 TO 6
435 LET X(I) = 0
440 LET Y(I) = 0
445 LET Z(I) = 0
446 NEXT I
450 IF C = 0 THEN GOTO 130
455 IF C = 1 THEN GOTO 475
456 IF C = 2 THEN GOTO 1000
460 IF C = 3 THEN GOTO 2000
470 GOTO 410
475 PRINT "POINT SOURCE"
480 PRINT "ENTER SOURCE ACTIVITY (CURIES)"
485 INPUT SA
505 CLS
510 PRINT "ENTER DISTANCE (CM)"
515 INPUT H
520 CLS
525 LET D = SA * 3.7E+10
530 IF B = 7 THEN GOTO 565
535 FOR I = 1 TO 6
540 IF E(I) = 0 THEN GOTO 585
545 GOSUB 4000
550 LET Z(I) = K(I) * D / (4 * 3.1416 * (H ^ 2)) * ((1 + S(I) * H *
GB(1) + ((S(I) * H) ^ 2) * GB(2) + ((S(I) * H) ^ 3) * GB(3)) * EXP(-
S(I) * H))
555 LET X(I) = K(I) * D / (4 * 3.1416 * (H ^ 2)) * EXP(-S(I) * H)
560 NEXT I
561 GOTO 585
REM THIS SMALL SECTION CALCS. THE NO SHIELD DOSE.
565 FOR I = 1 TO 6
570 IF E(I) = 0 THEN GOTO 585
575 LET X(I) = K(I) * D / (4 * 3.1416 * (H ^ 2))
580 NEXT I
REM 3000 CALCULATES THE DOES VALUES.
REM X(I) IS THE FLUX
REM Y(I) IS THE NO BUILDUP DOSE
REM Z(I) IS THE BUILDUP DOSE
REM 6000 PRINTS THE RESULTS.
585 GOSUB 3000
590 GOSUB 6000
595 GOTO 410
1000 PRINT "LINE SOURCE"
1005 PRINT "ENTER ACTIVTY(CURIES/CM)"
1010 INPUT SA
1015 CLS
1035 PRINT "ENTER DISTANCE (CM)"
1040 INPUT H
1045 CLS
1050 PRINT "ENTER SOURCE DIMENSION (CM)"
1055 INPUT DIMA
1060 LET D = SA * 3.7E+10
1065 LET THE = ATN(DIMA / (2 * H))
1070 IF B = 7 THEN GOTO 1145
1075 FOR I = 1 TO 6
1080 IF E(I) = 0 THEN GOTO 1165
1085 GOSUB 4000
1095 LET N = 1
1100 LET LOW = -THE
1105 LET UP = THE
1110 LET NUM = 1
REM 5000 COMPUTES THE LINE INTEGRAL WITH NO BUILDUP
1115 GOSUB 5000
1120 LET X(I) = (D / (4 * 3.1416 * H)) * PA * K(I)
REM 7000 COMPUTES THE LINE INTEGRAL WITH BUILDUP
1130 GOSUB 7000
1135 LET Z(I) = (D / (4 * 3.1416 * H)) * PA * K(I)
1136 PRINT "PA+"; PA
1140 NEXT I
1142 GOTO 1165
REM NO SHIELD SECTION
1145 FOR I = 1 TO 6
1150 IF E(I) = 0 THEN GOTO 1165
1155 LET X(I) = (D / (4 * 3.1416 * H)) * 2 * THE * K(I)
1160 NEXT I
1165 GOSUB 3000
1170 GOSUB 6000
1175 GOTO 410
REM THIS SECTION IS REM'ED OUT DUE TO A LACK OF
REM COMPREHENSIVE REFERENCE.
2000 PRINT "AREA SOURCE"
2005 PRINT "ENTER ACTIVITY (CURIES/CM SQ.)"
2010 INPUT SA
2030 CLS
2035 PRINT "ENTER DISTANCE (CM)"
2040 INPUT H
2045 CLS
2050 PRINT "ENTER SOURCE DIMENSION-DISK DIAMETER- (CM)"
2055 INPUT DIMA
2060 LET D = SA * 3.7E+10
2065 IF B = 7 GOTO 2140
2075 FOR I = 1 TO 6
2080 IF E(I) = 0 THEN GOTO 2160
2085 GOSUB 4000
2090 LET LOW = H
2095 LET UP = SQR((DIMA / 2) ^ 2 + H ^ 2)
2100 LET NUM = 1
2110 GOSUB 8000
REM PA-RETURNED FROM INTEGRATION
2115 LET X(I) = (K(I) * D) / 2 * PA
PRINT "PANO="; PA
2125 GOSUB 9000
2130 LET Z(I) = (K(I) * D) / 2 * PA
PRINT "PABU="; PA
2135 NEXT I
2136 GOTO 2160
2140 FOR I = 1 TO 6
2145 IF E(I) = 0 THEN GOTO 2160
2146 PRINT "DIMA="; DIMA
2150 LET X(I) = K(I) * (D / (2)) * (LOG(SQR((DIMA / 2) ^ 2 + H ^ 2)) -
LOG(H))
2151 PRINT "X="; X(I)
2155 NEXT I
2160 GOSUB 3000
2165 GOSUB 6000
2166 GOTO 410
2170 STOP
3000 FOR I = 1 TO 6
3005 IF E(I) = 0 THEN GOTO 3035
3010 LET Y(I) = E(I) * X(I) * Q(I) * .0000576
3015 LET Z(I) = Z(I) * E(I) * Q(I) * .0000576
3020 NEXT I
3035 RETURN
REM A TABLE LOOKUP AND INTERPOLOATION SECTION FOR
REM POWER FUNCTION BUILDUP FACTORS.
4000 IF E(I) <= M(I) THEN GOTO 4025
4005 FOR J = 2 TO 13
4010 IF E(I) <= M(J) THEN GOTO 4035
4015 NEXT J
4020 GOTO 4035
4025 LET J = 1
4030 GOTO 4040
4035 LET J = J - 1
4040 FOR A = 1 TO 3
4041 IF B = 1 THEN GOTO 4052
4042 IF B = 2 THEN GOTO 4054
4043 IF B = 3 THEN GOTO 4056
4044 IF B = 4 THEN GOTO 4058
4045 IF B = 5 THEN GOTO 4060
4046 IF B = 6 THEN GOTO 4062
4047 IF B = 7 THEN GOTO 4064
4052 LET GB(A) = L(A, J)
4053 GOTO 4065
4054 LET GB(A) = C(A, J)
4055 GOTO 4065
4056 LET GB(A) = W(A, J)
4057 GOTO 4065
4058 LET GB(A) = F(A, J)
4059 GOTO 4065
4060 LET GB(A) = T(A, J)
4061 GOTO 4065
4062 LET GB(A) = A(A, J)
4063 GOTO 4065
4064 LET GB(A) = A(A, J)
4065 NEXT A
4070 RETURN
REM 5000 AND 7000 IS NUMERIC INTEGRATION FOR THE LINE SOURCE.
REM 5000 IS NO BUILDUP, 7000 IS BUILDUP
5000 LET PA = 0
5005 LET FA = 0
5006 LET OA = 0
REM
5020 LET CA = 0
5025 LET DA = UP
5030 LET FA = 0
5035 FOR K = 1 TO 6
5040 LET XX = (I(K) * (DA - CA) / 2 + ((DA + CA) / 2))
REM LA- INTEGRAL FUNCTION
5045 LET LA = EXP(-S(I) * H / COS(XX))
5050 LET GAA = LA * BD(K)
5055 LET FA = FA + GAA
5060 NEXT K
5065 LET OA = ((DA - CA)) * FA
5070 LET PA = PA + OA
REM
5076 PRINT "S(I)="; S(I)
5077 PRINT "PA="; PA; FA
5080 RETURN
6000 PRINT GA$(C)
6001 K = H
6005 PRINT "ACTIVITY="; SA; "CURIES"
6010 PRINT "SOURCE DIMENSION="; DIMA; "Point:CM-line:CM-diameter"
6012 IF B = 7 THEN K = 0
6015 PRINT "SHIELD THICKNESS="; K; "CM"
6016 PRINT "DISTANCE TO DOSE POINT="; H; "CM"
6017 PRINT " "
6020 PRINT "NB-FLUX DOSE RATE-NB DOSE RATE-B"
6030 LET AL = 0
6035 LET EL = 0
6040 FOR I = 1 TO 6
6060 IF E(I) = 0 THEN GOTO 6085
6065 PRINT USING " #.#####^^^^ "; X(I); Y(I); Z(I)
REM 6065 PRINT X$(I); " "; Y$(I); " "; Z$(I)
6070 LET AL = AL + Y(I)
6075 EL = EL + Z(I)
6080 NEXT I
6085 LET O$ = STR$(AL)
6090 LET M$ = STR$(EL)
6091 PRINT " "
6092 PRINT " TOTAL="; : PRINT USING " #.#####^^^^"; AL; : PRINT
USING " #.#####^^^^ "; EL
6100 PRINT "NB-FLUX (GAMMAS/CM2-SEC)"
6105 PRINT "NB- NO BUILDUP (RAD/HR)"
6110 PRINT "B-WITH BUILDUP (RAD/HR)"
6111 PRINT " "
6115 RETURN
7000 LET PA = 0
7005 LET FA = 0
7010 LET OA = 0
7020 LET CA = LOW
7025 LET DA = UP
7030 LET FA = 0
7035 FOR K = 1 TO 6
7040 LET XX = ((I(K) * (DA - CA) / 2) + ((DA + CA) / 2))
REM LA-INTEGRAL FUNCTION BUILDUP
7045 LET LA = (1 + (S(I) * (H / COS(XX))) * GB(1) + ((S(I) * (H /
COS(XX))) ^ 2) * GB(2) + ((S(I) * (H / COS(XX))) ^ 3) * GB(3)) * EXP((-
S(I) * H) / COS(XX))
7046 PRINT "GB="; GB(1); GB(2); GB(3)
7050 LET GAA = LA * BD(K)
7055 LET FA = FA + GAA
7060 NEXT K
7065 LET OA = ((DA - CA) / 2) * FA
7070 LET PA = PA + OA
7075 PRINT "PA_buildup)="; PA
7080 RETURN
REM 8000 AND 9000 IS NUMERIC INTEGRATION FOR THE DISK SOURCE.
8000 LET PA = 0
8005 LET FA = 0
8010 LET OA = 0
8020 LET CA = LOW
8025 LET DA = UP
8030 LET FA = 0
8035 FOR K = 1 TO 6
8040 LET XX = ((I(K) * (DA - CA) / 2) + ((DA + CA) / 2))
REM NO BUILDUP INTEGRAL FUNCTION-LA
8045 LET LA = (EXP(-XX * S(I))) / XX
8050 LET GAA = LA * BD(K)
8055 LET FA = FA + GAA
8060 NEXT K
8065 LET OA = ((DA - CA) / 2) * FA
8070 LET PA = PA + OA
8080 RETURN
9000 LET PA = 0
9005 LET FA = 0
9010 LET OA = 0
9020 LET CA = LOW
9025 LET DA = UP
9030 LET FA = 0
9035 FOR K = 1 TO 6
9040 LET XX = ((I(K) * (DA - CA) / 2) + (DA + CA) / 2)
REM BUILDUP INTEGRAL ALL WRONG-LA
9045 LET LA = (1 + (S(I) * XX) * GB(1) + ((S(I) * XX) ^ 2) * GB(2) +
((S(I) * XX) ^ 3) * GB(3)) * EXP(-XX * S(I)) / XX
9050 LET GAA = LA * BD(K)
9055 LET FA = FA + GAA
9060 NEXT K
9065 LET OA = ((DA - CA) / 2) * FA
9070 LET PA = PA + OA
9080 RETURN
REM D-BUILDUP ENERGIES MEV
DLIST:
DATA .03,.04,.05,.06,.08,.1,.2,.3,.4,.5,1.0,2.0,3.0
REM M-MASS ENERGY ABSORP. COEFF. TISSUE. CM2/G
MLIST:
DATA .
1540,.0677,.0409,.0312,.0255,.0252,.0276,.0297,.0317,.0325,.0327,.0326,.0318,.0308,.0281,.0257,.0225,.0225
REM GAUSSIAN QUADRATURE COEFFICIENTS B AND I
BDLIST:
DATA .171324492,.360761573,.467913935,.467913935,.360761573,.171324492
ILIST:
DATA -.932469514,-.661209386,-.
238619186,.238619186,.661209386,.932469514
REM TISSUE BUILDUP COEFFICIENTS
TLIST:
DATA
1.277,2.310,2.866,2.906,2.783,2.602,1.595,1.118,1.028,1.001,.955,.810,.810
DATA -.
0634,.1184,.5643,.9967,.1498,.9078,.6541,.6617,.5502,.4488,.1812,.0478,.0478
DATA .0023,-.0005,.0024,.0267,.1288,.1931,.1191,.0413,.0164,.0050,-.
0028,-.0012,-.0012
REM AIR BUILDUP COEFFICIENTS
ALIST:
DATA
1.219,2.251,2.852,2.960,2.719,2.485,1.602,1.117,1.045,1.000,.948,.798,.798
DATA -.
0673,.0905,.5033,.9288,1.1714,.0343,.6458,.6743,.5391,.4492,.1824,.0487,.0487
DATA .0025,-.0002,.0015,.0215,.1095,.1600,.1167,.0366,.0163,.0038,-.
0028,-.0012,-.0012
REM LEAD BUILDUP COEFFICIENTS
LLIST:
DATA .496,.496,.496,.496,.496,.496,.496,.496,.496,.496,.766,.564,.353
DATA -.0522,-.0522,-.0522,-.0522,-.0522,-.0522,-.0522,-.0522,-.0522,-.
0522,-.0522,-.0439,.0080,.0054
DATA .
002,.002,.002,.002,.002,.002,.002,.002,.002,.002,.00115,.0004,.0003
REM IRON BUILDUP COEFFICENTS
FLIST:
DATA
1.592,1.592,1.592,1.592,1.592,1.592,1.592,1.592,1.592,1.592,1.074,.750,.553
DATA .
1629,.1629,.1629,.1629,.1629,.1629,.1629,.1629,.1629,.1629,.1104,.0358,.0269
DATA .008,.008,.008,.008,.008,.008,.008,.008,.008,.008,-.0009,-.0003,-.
0004
REM WATER BUILDUP COEFFICIENTS
WLIST:
DATA
1.320,22.382,2.910,2.877,2.714,2.639,1.532,1.129,1.004,.992,.950,.799,.678
DATA -.
0600,.1512,.6340,1.0688,1.1797,.8770,.6829,.6313,.5595,.4510,.1785,.0481,.0146
DATA .0023,-.0008,.0041,.0354,.1434,.2113,.1168,.0486,.0139,.0038,-.
0026,-.0012,-.0004
REM CONCRETE BUILDUP COEFFICIENTS
CLIST:
DATA 1.11,1.11,1.11,1.11,1.11,1.1,1.11,1.11,1.11,1.11,.95,.74,.61
DATA .15,.15,.15,.15,.15,.15,.15,.15,.15,.15,.09,.04,.02
DATA .
00167,.00167,.00167,.00167,.00167,.00167,.00167,.00167,.00167,.00167,-.
00057,-.00051,-.00016
REM LINEAR ATTEN. COEFF. FOR LEAD,CONCRETE,WATER,IRON,TISSUE,AIR
REM U(7,X) IS THE ENRGY VALUES USED FOR INTERPOLATION-MEV.
ULIST:
DATA
337.1,158.9,88.64,55.27,26.45,61.29,22.36,11.25,4.585,2.622,1.827,1.419,1.005,.8036,.5868,.5164,.4744,.4744
DATA
2.737,1.392,.9016,.6785,.4899,.4117,.3312,.2921,.2484,.2215,.2017,.1863,.1631,.1465,.1194,.1030,.0840,.0840
DATA .
370,.267,.227,.206,.184,.171,.151,.137,.119,.106,.0968,.0896,.0786,.0707,.0575,.0494,.0397,.0397
DATA
63.9,28.45,15.25,9.43,4.677,2.908,1.541,1.148,.8646,.7388,.6602,.6044,.5258,.4708,.3836,.3341,.2845,.2845
DATA .
4737,.3404,.2883,.2604,.2324,.2159,.1892,.1727,.1499,.1334,.1219,.1128,.0989,.0889,.0724,.0621,.0499,.0499
DATA .
000451,.000321,.000269,.000243,.000216,.000199,.000176,.000159,.000138,.000123,.000113,.000104,.000091,.000082,.000067,.000058,.000046,.000046
DATA .03,.04,.05,.06,.08,.10,.15,.2,.3,.4,.5,.6,.8,1.0,1.5,2.0,3.0,3.0
.
|
|
| User: "H. Wabnig .... .-- .- -... -. .. --. @ .- --- -. DOT .- -" |
|
| Title: Re: A Radiation Shielding Qbasic Calculator |
16 Apr 2007 03:20:11 AM |
|
|
On 15 Apr 2007 08:10:12 -0700, "Douglas Eagleson"
<eaglesondouglas@yahoo.com> wrote:
Here is some source code for gamma ray shielding. It is good for
0.03Mev-3 Mev
It was reviewed for consistency, but not absolute accuracy. So it is
only for evaluation purposes.
The idea was to automate the basic , point, line, disk sources.
Just get the gamma ray energy and yield from a reference and plug in
the inputs!
Microshield the commercial product is the thing to use for real work,
but soon I will compare the two. Mine is free. Mine has 3 geometries,
Microshield has 15. But I once made a set of line sources to simulate
all geometries. And min ehas no fancy graphics, just old time command
line inputs/outputs.
I can email a executable for Windows-xp if anybody is interested. Also
google goofs-up the long lined data sets. So be careful if
reformating for a local qbasic run.
********begin real
code**********************************************************
REM DOSECALC IS A PROGRAM TO ASSIST IN GAMMA RAY SAFETY.
REM IT CALCULATES POINT SOURCE DOSE RATE AND LINE SOURCE DOSE
REM RATE FOR THE SELECTED SHIELDING MATERIALS.
REM VERSION 1.0a
REM BY: DOUGLAS EAGLESON
REM VERSION 1.01a WILL BE COMING OUT SOON.
.........
Thanks for your efforts,
but where does Windows XP hide that QBASIC interpreter,
and willl it be available under Windows Vista?
w.
.
|
|
|
| User: "Douglas Eagleson" |
|
| Title: Re: A Radiation Shielding Qbasic Calculator |
16 Apr 2007 08:21:55 AM |
|
|
On Apr 16, 4:20 am, H. Wabnig <.... .-- .- -... -. .. --. @ .-
--- -. DOT .- -> wrote:
On 15 Apr 2007 08:10:12 -0700, "Douglas Eagleson"
<eaglesondoug...@yahoo.com> wrote:
Here is some source code for gamma ray shielding. It is good for
0.03Mev-3 Mev
It was reviewed for consistency, but not absolute accuracy. So it is
only for evaluation purposes.
The idea was to automate the basic , point, line, disk sources.
Just get the gamma ray energy and yield from a reference and plug in
the inputs!
Microshield the commercial product is the thing to use for real work,
but soon I will compare the two. Mine is free. Mine has 3 geometries,
Microshield has 15. But I once made a set of line sources to simulate
all geometries. And min ehas no fancy graphics, just old time command
line inputs/outputs.
I can email a executable for Windows-xp if anybody is interested. Also
google goofs-up the long lined data sets. So be careful if
reformating for a local qbasic run.
********begin real
code**********************************************************
REM DOSECALC IS A PROGRAM TO ASSIST IN GAMMA RAY SAFETY.
REM IT CALCULATES POINT SOURCE DOSE RATE AND LINE SOURCE DOSE
REM RATE FOR THE SELECTED SHIELDING MATERIALS.
REM VERSION 1.0a
REM BY: DOUGLAS EAGLESON
REM VERSION 1.01a WILL BE COMING OUT SOON.
.........
Thanks for your efforts,
but where does Windows XP hide that QBASIC interpreter,
and willl it be available under Windows Vista?
w.- Hide quoted text -
- Show quoted text -
Qbasic 4.5 makes executables that run under windows xp. A small
command window encapsulates the basic.
I will try the code on a windows vista machine over at COMP USA.
.
|
|
|
| User: "Douglas Eagleson" |
|
| Title: Re: A Radiation Shielding Qbasic Calculator |
16 Apr 2007 03:56:58 PM |
|
|
On Apr 16, 9:21 am, "Douglas Eagleson" <eaglesondoug...@yahoo.com>
wrote:
On Apr 16, 4:20 am, H. Wabnig <.... .-- .- -... -. .. --. @ .-
--- -. DOT .- -> wrote:
On 15 Apr 2007 08:10:12 -0700, "Douglas Eagleson"
<eaglesondoug...@yahoo.com> wrote:
Here is some source code for gamma ray shielding. It is good for
0.03Mev-3 Mev
It was reviewed for consistency, but not absolute accuracy. So it is
only for evaluation purposes.
The idea was to automate the basic , point, line, disk sources.
Just get the gamma ray energy and yield from a reference and plug in
the inputs!
Microshield the commercial product is the thing to use for real work,
but soon I will compare the two. Mine is free. Mine has 3 geometries,
Microshield has 15. But I once made a set of line sources to simulate
all geometries. And min ehas no fancy graphics, just old time command
line inputs/outputs.
I can email a executable for Windows-xp if anybody is interested. Also
google goofs-up the long lined data sets. So be careful if
reformating for a local qbasic run.
********begin real
code**********************************************************
REM DOSECALC IS A PROGRAM TO ASSIST IN GAMMA RAY SAFETY.
REM IT CALCULATES POINT SOURCE DOSE RATE AND LINE SOURCE DOSE
REM RATE FOR THE SELECTED SHIELDING MATERIALS.
REM VERSION 1.0a
REM BY: DOUGLAS EAGLESON
REM VERSION 1.01a WILL BE COMING OUT SOON.
.........
Thanks for your efforts,
but where does Windows XP hide that QBASIC interpreter,
and willl it be available under Windows Vista?
w.- Hide quoted text -
- Show quoted text -
Qbasic 4.5 makes executables that run under windows xp. A small
command window encapsulates the basic.
I will try the code on a windows vista machine over at COMP USA.- Hide quoted text -
- Show quoted text -
I tried it on VISTA. Qbasic 4.5 works there and makes a VISTA
runnable executable application file.
The only problem was a file reference in the posted source file.. It
needs removal so the bad path does not crash the application.
.
|
|
|
|
|
|
|
Related Articles |
|
|
