20050101
CLS
PI = 3.14159
CR1 = 18
CD1 = 9
PRINT " input 1 2 3 6 OR 9 for standard? or known PARAMETERS "
INPUT " DIAMETER METERS ="; ID
IF ID = 1 THEN GOSUB 910
IF ID = 2 THEN GOSUB 920
IF ID = 3 THEN GOSUB 930
IF ID = 6 THEN GOSUB 960
IF ID = 9 THEN GOSUB 990
IF GA = 60 THEN GOTO 100
LOCATE 2, 1
INPUT " DENSITY Kg/m^3 (3000ish ) = "; IR
LOCATE 2, 1
INPUT " VEOCITY Km/s (20 to 30) = "; IK
LOCATE 2, 1
INPUT " ANGLE DEG (30 to 90) = "; GA
100
PRINT ""
PRINT "for an impactor of "
PRINT "DIAMETER = "; ID; TAB(20); "meters"; TAB(35);
PRINT "DENSITY = "; IR; "Kg/m^3"
PRINT "VELOCITY = "; IK; TAB(20); " Km/s"; TAB(35);
PRINT "entry ANGLE = "; GA; " degrees "
IV = IK * 1000: VI = PI * ID * ID * ID / 6
GF = (SIN(GA / 180 * PI)) ^ .33
MI = IR * VI: KE = .5 * MI * IV * IV
ER = KE * 1E+07
KWHR = KE / 3600000!
KT = KE / 4.2E+12
RH = KE / 1E+17
TP = KE / 1E+23
SF = KE / 1E+29
OB = KE / 1E+35
CD = 2 * CR1 * KT ^ .3 * GF: CA = CD * 1.25
CZ = CD1 * KT ^ .3 * GF: CL = CZ / 1.25
CV = .5 * PI * CD * CD / 4 * CZ: CE = CD * 2.15
PRINT " WITH KE = "; KE; " JOULE"
PRINT ""
PRINT "INPACTOR"; TAB(38); "CREATOR DIAMETER"
PRINT USING " VOLUME ##.# M^3 "; VI; TAB(36);
PRINT USING "ACTUAL ###### M "; CD
PRINT USING " MASS ###.# Kg "; MI; TAB(36);
PRINT USING "APPEARENT ###### M "; CA
PRINT TAB(42); "DEPTH"
PRINT "( 0) ERG ";
PRINT TAB(12); INT(ER); TAB(36);
PRINT USING "ACTUAL ###### M"; CZ
PRINT "( 6) JOULE";
PRINT TAB(12); INT(KE); TAB(36);
PRINT USING "APPEARENT ###### M"; CL
PRINT "(12) KWHR ";
PRINT TAB(12); INT(KWHR)
PRINT "(18) KT ";
PRINT TAB(12); INT(KT)
PRINT "(24) Rh ";
PRINT TAB(12); RH; TAB(28);
PRINT "(24)Cr "; INT(CD / 1000); TAB(46);
IF RH < 250 THEN PRINT "(24)Ri "; INT(KT ^ .3 / 10) ELSE PRINT " "
PRINT "(30) Tp "; INT(TP)
PRINT "(36) Sf "; INT(SF); TAB(42);
PRINT "EJECTRA"
PRINT "(42) Ob "; INT(OB); TAB(36);
PRINT USING "REMOVED ##.# M^3"; CV
PRINT ""
END
910
PRINT "_ 100 MEGATONNE STANDARD 1KM CRATER _"
ID = 100
IK = 25
IR = 2600
GA = 60
RETURN
920
PRINT "Cr_2"
ID = 200
IK = 27
IR = 2700
GA = 60
RETURN
930
PRINT "STANDARD 9 Ri QUAKE?"
ID = 350
IK = 27
IR = 2000
GA = 60
RETURN
960
PRINT "STANDARD 100KM CRATER ~> Copernicus "
ID = 13000
IK = 28
IR = 3300
GA = 60
RETURN
990
PRINT "STANDARD 2BIG 4Moon Craters"
ID = 130000
IK = 28
IR = 3300
GA = 60
RETURN
|