2 CLS 5 PRINT "SADEYA INFORTECNICA" 6 PRINT 10 PRINT "ENERGIA DE UN METEOROIDE Y CONSECUENCIA DE LA CAIDA" 15 PRINT "A PARTIR DEL RADIANTE" 20 PRINT INPUT "CONSECUENCIAS(1) AREA(2)"; TEPI 23 IF TEPI = 2 THEN 5000 30 REM INPUT "DECLINACION DEL RADIANTE GRADOS"; DEC 31 INPUT "VELOCIDAD INFINITA O RELATIVA KM/S"; VEL 32 INPUT "DENSIDAD (AGUA=1): COMETAS DE .3 A 1, ASTEROIDES 2.5 A 3"; DEN 33 INPUT "DIAMETRO EN METROS"; DIA INPUT "DISTANCIA DEL LUGAR AL EPIPUNTO KM."; EPP 36 MASA = (4 / 3) * 3.1416 * (DIA / 2) ^ 3 * DEN VOLU = (4 / 3) * 3.1416 * (DIA / 2) ^ 3 VOLUK = VOLU / 1000000000 NUMER = 5000 / VOLUK 37 vef = SQR(VEL * VEL + 11.2 * 11.2) 55 PRINT 77 JUL = ((MASA * 1000) / 2) * (vef * 1000 * vef * 1000) 78 REM LAT = ((DEC / ABS(DEC)) * -90) + DEC 79 PA = .001 * DIA ^ 3: PC = 1 / PA PRINT "METEOROIDE" 80 REM PRINT "LATITUDES LIMITES DE IMPACTO="; 90 * (DEC / ABS(DEC)); LAT 85 PRINT "VELOCIDAD DE ENTRADA O IMPACTO KM/S="; vef 86 PRINT "MASA EN TONELADAS="; MASA 87 PRINT "ENERGIA JOULS, KILOT ="; JUL; (JUL / 850000000000#) / 4.2# TNT = (JUL / 850000000000#) / 4.2# RAYO = JUL / 1000000000 PRINT "ENER. (RAYO=1)="; RAYO altdt = 25.5 - 15 * LOG(DIA) / LOG(10) altd = 42.5 - 25 * LOG(DIA) / LOG(10) IF altd < 0 THEN altd = altd / 10 * DIA / 500 IF altd > 0 THEN PRINT "ALTURA: TRANSF., DESINTEG. KM. ="; altd; altdt IF altd < 0 THEN PRINT "IMPACTO SUP. DISIPACION ENRG. PROF. KM. ="; altd LUMZ = JUL / ((altd * altd) ^ 1) MAGZ = LUMZ / (LOG(2.512) / LOG(10)) hip = SQR((altd * altd) + (EPP * EPP)) LUMZ2 = JUL / (hip * hip) ^ 1 MAGZ2 = LUMZ2 / (LOG(2.512) / LOG(10)) CON = .000000000001# lumzz = LUMZ * CON lumzz2 = LUMZ2 * CON mag = 5 - 2.66 * LOG(MAGZ) / LOG(10) mag2 = 5 - 2.66 * LOG(MAGZ2) / LOG(10) PRINT "MAG.MINIMA VERT. Y DEL LUGAR="; mag; mag2 PRINT "INT RADIACION W/M2s VERT. Y DEL LUGAR ="; LUMZ / 255400000; LUMZ2 / 255400000 PRINT "J. VERT. Y DEL LUGAR ="; LUMZ / 1000000; LUMZ2 / 1000000 IF DIA < 500 THEN PRINT "DISTA.EQUIV.RAYO EN METROS: VERT. Y LUGAR ="; 1 / SQR(LUMZ / 1000000000000000#); 1 / SQR(LUMZ2 / 1000000000000000#) SOND = altd / .34 / 60 SOND2 = (SQR((altd * altd) ^ 1 + (EPP * EPP) ^ 1) / .34 / 60) db1 = ((mag * -1) - 18) * 15 db2 = ((mag2 * -1) - 18) * 15 IF db1 < 0 THEN db1 = 0 IF db2 < 0 THEN db2 = 0 REM PRINT "ONDA P.INT.SONIDO VERT, Y DEL LUGAR W/M2,="; lumzz; lumzz2 PRINT "ONDA P.INT.SONIDO VERT, Y DEL LUGAR dB.="; db1; db2 PRINT "RETRASO SONIDO VERT, Y DEL LUGAR MIN.="; SOND; SOND2 REM INPUT "ENERGIA DEL TERREMOTO EQUIVALENTE EN EXPLOSION TM. DE TNT"; TNT TNT = TNT * 1000 richep = (LOG(TNT) / LOG(10) + LOG(33) / LOG(10) * 2) / (LOG(33) / LOG(10)) IF DIA > 20 THEN PRINT "I. ESCALA RICHTER = "; richep 90 PRINT "PROBABILIDAD EN UN A¥O ="; PC 92 PRINT "PERIODO MEDIO CAIDA A¥OS: TIERRA, CONTIN. ="; 1 / PC; (1 / PC) * 1 / .3 IF DIA < 500 THEN PRINT "PERIODO A¥OS EN AREA: POBLADA, URBANA ="; (1 / PC) * 1 / .05; (1 / PC) * 1 / .005 areat = 6400 ^ 2 * 3.1416 * 4 pk2 = (PC) / areat IF DIA < 1 THEN PRINT "CAIDOS EN LA TIERRA EN UN A¥O POR: KM2, M2="; pk2; pk2 / 1000000 PRINT "NUMERO DE NEAS PROBABLE ="; NUMER MEG = (JUL / 850000000000000#) / 4.2# RAD = SQR(MEG * 200) AREA = RAD ^ 2 * 3.1416 PRINT INPUT PP 93 IF DIA < 20 THEN 98 PRINT "CONSECUENCIAS GEOFISICAS" 94 PRINT "ENERGIA TRANSFORMADA MEGATONES= "; MEG: IF RAD > 20000 THEN RAD = 20000 PRINT "RADIO AFECTACION TERMICA KM.="; RAD: IF RAD > 3000 THEN PRINT "AFECTACION CONTINENTAL " IF RAD < 500 THEN PRINT "RADIO AFECTACION ONDA P.KM.0"; RAD * 10 IF AREA > 100000000 THEN AREA = 100000000 PRINT "AREA AFECTACION TERMICA KM2.="; AREA: IF AREA > 70000000 THEN PRINT "AFECTACION GLOBAL" 95 IF DIA < 62 THEN 98 96 REM IF DIA > 10000 THEN 98 DIAM = 18 * DIA + (DIA / 10) ^ 2 VOLUTS = 100 * (DIA ^ 3) IF DIA > 1000 THEN VOLUTS = 100 * (1000 ^ 3) REM ALTS=SQR(VOLUTS/(6.28*DIMP)) ALTS1 = SQR(VOLUTS / (6.28 * 1000)) ALTS2 = SQR(VOLUTS / (6.28 * 100000)) ALTS3 = SQR(VOLUTS / (6.28 * 10000000)) IF DIA > 50 THEN PRINT "SI IMPACTO MARINO, TSUNAMI ALT/m.:1KM,100KM,10000KM = "; INT(ALTS1); INT(ALTS2); INT(ALTS3) PRINT "SI IMPACTO TERRESTRE, CRATER METROS DIAMETRO= "; DIAM IF DIA > 1000 THEN PRINT "SI IMPACTO EN FONDO MARINO CRATER METROS="; DIAM IF DIAM > 250000 THEN PRINT "PERO DESAPARECE POR FUSION Y EXPANSION DE PLACA" IF DIAM > 300000 THEN PRINT "CON PLEGAMIENTO CORTICAL Y DESPLAZAMIENTOS CONTINENTALES" IF DIAM > 300000 THEN PRINT "POR FORMACION DE EXODORSAL Y TUNEL FLUIDO CON EL NUCLEO" CALT = (JUL / (6 * 10 ^ 27)) MASAFK = JUL / 1000000 MASAFT = MASAFK / 1000 VOLFM = MASAFT / 3 RADFM = VOLFM ^ (1 / 3) PRINT "RADIO FUNDIDO METROS="; RADFM PRINT "INCREMENTO TEMPERATURA CONJUNTO MASA TIERRA GRD.="; CALT Prof = 3 * DIA IF Prof > 1000 THEN Prof = 1000 IF DIAM > 250000 THEN 97 PRINT " PROFUNDIDAD METROS="; Prof IF Prof = 1000 THEN PRINT " FUSION DE FONDOS, CRATER NO SIMPLE" IF DIAM > 400000 THEN PRINT "CON REFRIGERACION POR EL MANTO" 97 KM = 1000 * (DIA / 1000) ^ 3 vol = 1000 * (DIA ^ 3): MASS = 2 * vol * 1000: sedk = MASS / (4 * 3.142 * 6500 ^ 2) PRINT " VOLUMEN MATERIALES DISPERSADOS KM3="; KM REM CORRECION DIFER METEORITO ERUPCION KM = 4 * KM PRINT " DEPOSITO SEDIMENTOS KG/Km2, KG/m2 ="; sedk; sedk / 1000000 DMAT = SQR(sedk / 1000000 * 6) ALTSM = (sedk / 1000000) / 2000: ALTSMM = ALTSM * 1000 PRINT " ALTURA SEDIMENTOS MEDIOS Y DISTANTES mm="; ALTSMM; ALTSMM / 10 PRINT INPUT PP IF DIA > 801 THEN PRINT "INFLUENCIA CLIMATICA" IF DIA < 801 THEN PRINT "INFLUENCIA CLIMATICA SI IMPACTO TERRESTRE" IF DIAM > 2000 THEN 200 98 PRAD = ALTSMM * 20 / 4: PRED1 = SQR(PRAD * 100): IF PRED1 > 99 THEN PRED1 = 99 PRED2 = SQR(PRAD * 100) * .3: IF PRED2 > 99 THEN PRED2 = 99 PRED3 = SQR(PRAD * 100) * .03: IF PRED3 > 99 THEN PRED3 = 99 PRINT "RED RADIANCIA:INMEDIATA,+1 A¥O,+10A¥S %= "; PRED1; PRED2; PRED3 IF PRED1 > 1 THEN PRINT "CON CAMBIO CLIMATICO, ENFRIAMIENTO POR AUMENTO DE ALBEDO" IF PRED2 > 2 THEN PRINT "Y GLACIACION INTENSA" IF DIA > 2000 THEN PRINT "PERO ANTECEDIDO DE CALENTAMIENTO INMEDIATO POR CALCINACION" CIC = KM * 2.35 DES = SQR(PRED1) * 3 ENE = KM * 15 * 7.2 MIN = KM * .38 IF DIA < 100 THEN 99 IF DIA < 100 THEN 99 PRINT "DESCENSO INMEDIATO GRADOS ATMOSFERA="; DES PRINT "DESCENSO PROLONGADO GRADOS OCEANOS="; DES / 5 PRINT "DURACION DEL CICLO HASTA LA RECUPERACION TEMP. ATMOSFERICA A¥OS = "; MIN PRINT "DURACION DEL CICLO HASTA LA RECUPERACION TEMP. OCEANICA A¥OS = "; CIC PRINT PRINT "CONSECUENCIAS BIOLOGICAS" IF RAD > 200 THEN PRINT "AFECTACION Y EXTINCION LOCAL" IF RAD > 3000 THEN PRINT "EXTINCION CONTINENTAL EN LA ZONA DE IMPACTO " IF RAD > 7000 THEN PRINT "EXTINCION GLOBAL CON PERDIDA DE BIOMASA" IF PRED2 > 10 THEN PRINT "EN UN PORCENTAJE DE BIOMASA DEL "; INT(PRED2); " POR CIENTO" IF PRED3 > 2 THEN PRINT "Y UN PORCENTAJE DEL "; INT(PRED3); " POR CIENTO DE ESPECIES." IF DIA > 3000 THEN PRINT "SI OXIGENO, CALCINACION GLOBAL CONTINENTES" IF DIA > 4900 THEN PRINT "EXTINCION CONTINENTAL ESPECIES DE MAS DE 20 KILOS" IF DIA > 6900 THEN PRINT "EXTINCION EXTREMA INCLUSO MARINA" IF DIA > 9990 THEN PRINT "EXTINCION GLOBAL ORGANISMOS PLURICELULARES" IF DIA > 25000 THEN PRINT "SUPERVIVENCIA SOLO MICROORGANISMOS EXTREMOFILOS" 99 PRINT INPUT "PARA NUEVO CALCULO 1"; NC 100 IF NC = 1 THEN 2 102 END 200 PRINT PRINT "FUENTE DE CALOR INTENSA IRRADIADA DE LA ZONA DE IMPACTO" IF DIA > 400 THEN PRINT "CON GRAN ACTIVIDAD VOLCANICA, ADEMAS DE UNA EXTREMADA" 210 PRINT " EMANACION DE CO2, GASES SULFUROSOS Y OTROS" 220 PRINT " DA¥OS GENERALIZADOS DE ESPECIES TERRESTRES" PRINT 230 GOTO 98 5000 PRINT 50010 PRINT "CENTRO APROXIMADO DE AREA DE IMPACTO DE UN METEORITO" PRINT INPUT "T.U. ENTRADA ATMOSFERA H.FRAC="; tui INPUT "A.R. RADIANTE H.FRAC="; ARR INPUT "DEC. RADIANTE GRAD.="; der INPUT "MES.FRACCION EFEMERIDE ="; MES ARS = (MES - 3.7) * 2 dar = ARR - ARS car = 12 + tui - dar car = (car / 24) + 4 car = (car - INT(car)) * 24 car = INT(car * 10) / 10 car2 = 12 - tui - dar LMAX = car + 6 IF LMAX > 24 THEN LMAX = LMAX - 24 LMIN = car - 6 IF LMIN < 0 THEN LMIN = LMIN + 24 LAMAX = der + 90 IF LAMAX > 90 THEN LAMAX = 90 LAMIN = der - 90 IF LAMIN < -90 THEN LAMIN = -90 PRINT PRINT "LONGITUD CENTRO AREA HH="; car PRINT "LONGITUD MAXIMA Y MINIMA HH="; LMAX; LMIN PRINT PRINT "LATITUD CENTRO AREA GRD="; der PRINT "LATITUD MAXIMA Y MINIMA GRD="; LAMAX; LAMIN PRINT INPUT "NUEVO CALCULO =1"; NC IF NC = 1 THEN 5000 END