0 CLS 10 PRINT "SADEYA/INFORTECNICA" PRINT PRINT "Lepant, 422 08025 Barcelona" PRINT "0034 934334243, www.sadeya.org" PRINT PRINT "La valoraci¢n r pida de las consecuencias de" PRINT "un terremoto permite tomar medidas inmediatas" PRINT "que pueden salvar muchas vidas." PRINT "Este programa puede ser £til para anticiparse" PRINT "a la informaci¢n confusa en el inicio." PRINT 1 REM 20 PRINT 30 PRINT "CONSECUENCIAS DE UN TERREMOTO" 40 PRINT 42 INPUT "CAUSAS(6) CONSECUENCIAS(1) RICHTER(2) MULTIPLE(3) DISTANCIA(4) PROBABILIDAD(5)"; UUU CLS 43 IF UUU = 2 THEN 1000 IF UUU = 4 THEN 4000 IF UUU = 5 THEN 10000 IF UUU = 6 THEN 15000 44 PRINT NUMT = 0 TNT = 0 RICHEP = 0 50 IF UUU = 3 THEN INPUT "NUMERO DE SISMOS ="; NUSI IF UUU = 3 THEN INPUT "INTERVALO DE TIEMPO EN UNIDADES T.(A¥OS, O DIAS, O SEG. ETC..)="; unt 55 REM CLS 56 NUMT = NUMT + 1 hipo = 20 57 IF UUU = 3 THEN PRINT "SISMO NUMERO "; NUMT 58 INPUT "INTENSIDAD EN EL EPICENTRO, ESCALA RICHTER ="; RICHEP CLS 60 INPUT "DISTANCIA AL EPICENTRO EN KMS. ="; dist INPUT "PROFUNDIDAD HIPOCENTRO KM. (MEDIO=20)"; hipo dist = dist + hipo 62 IF hipo < .01 THEN hipo = .01 hipe = hipo IF hipe > 30 THEN hipe = 30 63 VALOR = RICHEP * 25 * (hipe / 20) 70 RICHLOC = RICHEP - (dist / VALOR) DIS2 = (5 - RICHEP) * VALOR * -1 80 PRINT "INTENSIDAD ORIGEN Y MEDIA LOCAL (RICHTER) ="; RICHEP; RICHLOC 90 RICH = RICHLOC 100 DURA1 = .25 * 3 ^ (RICHEP - 2.5) DURA = .25 * 3 ^ (RICHLOC - 2.5) IF DURA < 1.0001 THEN DURA = 1 DISG = dist * 90 / 10000 velmg = 4.3 + DISG / 36 TP = DISG / velmg * 60: TPP = TP * 1.25: TS = TP * 1.8: TSS = TP * 2.3 110 PRINT "RETRASO SE¥AL: P,PP,S,SS="; INT(TP); INT(TPP); INT(TS); INT(TSS); "seg." PRINT "DURACION ORIGEN Y LOCAL DENTRO DEL GRADO="; DURA1; DURA; "seg." IF RICH < -5 THEN RICH = -5 NUMER = 1.3 * 10 ^ (8 - RICH) PROB = 1 / NUMER 111 INTEN = 33.3 ^ (RICH - 2) 112 INTENJ = INTEN * 3.57E+09 115 PRINT "ENERGIA LOCAL EN KILOTON (TNT) Y EN JOULS ="; INTEN / 1000; INTENJ 117 REM PRINT "LOCAL: SISMOS/A¥O GLOBAL TIERRA < "; NUMER; " INTERVALO A¥OS > "; PROB IF UUU < 3 THEN 118 TNT = INTEN + TNT TNT2 = TNT / unt RICHLOG = (LOG(TNT) / LOG(10) + LOG(33) / LOG(10) * 2) / (LOG(33) / LOG(10)) RICHLOG2 = (LOG(TNT2) / LOG(10) + LOG(33) / LOG(10) * 2) / (LOG(33) / LOG(10)) IF UUU = 3 THEN PRINT "ENERGIA MEGATONES ACUMULADA ="; TNT / 1000000 IF UUU = 3 THEN PRINT "RICHTER ACUMULADO ="; RICHLOG IF UUU = 3 THEN PRINT "RICHTER MEDIO POR UNIDAD DE TIEMPO="; RICHLOG2 IF NUMT < NUSI THEN 55 RICHEP = RICHLOG: RICH = RICHLOG 118 PRINT IF UUU = 3 THEN INPUT "OBTENER RESULTADOS SOBRE EL TOTAL DE ENERGIA (1)"; NNUO IF UUU = 3 GOTO 1 REM PRINT NUMER = 1.3 * 10 ^ (8 - RICHEP) PROB = 1 / NUMER INTEN = 33.3 ^ (RICHEP - 2) INTENJ = INTEN * 3.57E+09 PRINT "ORIGEN Y CONSECUENCIAS" PRINT "ENERGIA EN ORIGEN EN KILOTON (TNT) Y EN JOULS ="; INTEN / 1000; INTENJ PRINT "SISMOS/A¥O GLOBAL TIERRA < "; NUMER; " INTERVALO A¥OS > "; PROB AREF = DIS2 * DIS2 * 3.1416 IF RICHEP > 5 THEN PRINT "RADIO Y AREA DE AFECTACION DESTRUCTIVA (>5RICHTER) KM.KM2.="; INT(DIS2); INT(AREF) IF RICHEP > 10.5 THEN PRINT "Si epicentro marino, tsunami extremo seguro." IF RICHEP > 10.5 THEN PRINT "Probable causa: impacto infrecuente de gran meteoroide" IF RICHEP > 10.5 THEN PRINT "en cualquier lugar de la Tierra." IF RICHEP > 10.5 THEN PRINT "Actividad volc nica y s¡smica extrema consecuente." IF RICHEP > 10.5 THEN GOTO 120 IF RICHEP > 8.999 THEN PRINT "Si epicentro marino, tsunami muy probable." IF RICHEP > 8.999 THEN PRINT "M£ltiples sismos e hipocentros por colapsos sucesivos de placa" IF RICHEP > 8.999 THEN PRINT "en  reas de subducci¢n. R‚plicas muy intensas y duraderas." IF RICHEP > 8.999 THEN PRINT "Si hay volcanes en el  rea, probable activaci¢n." IF RICHEP > 8.999 THEN GOTO 120 IF RICHEP > 7.999 THEN PRINT "Si epicentro marino, tsunami probable." IF RICHEP > 7.999 THEN PRINT "Origen probable en  rea de subducci¢n." IF RICHEP > 7.999 THEN PRINT "R‚plicas intensas." IF RICHEP > 7.999 THEN GOTO 120 IF RICHEP > 6.999 THEN PRINT "Si epicentro marino, tsunami posible." IF RICHEP > 6.999 THEN PRINT "Origen tect¢nico. R‚plicas muy probables." IF RICHEP > 6.999 THEN GOTO 120 IF RICHEP < 6 THEN 119 IF RICHEP < 7 THEN PRINT "Si epicentro marino, tsunami poco probable." 119 IF RICHEP < 7 THEN PRINT "Origen tect¢nico o volc nico. Terremoto factible fuera" IF RICHEP < 7 THEN PRINT "de  reas de m xima actividad s¡smica." 120 PRINT 121 PRINT "AFECTACION LOCAL" 122 mercalli = RICH * 1.3 125 IF mercalli > 12.999 THEN mercalli = 12 127 IF mercalli < 1 THEN mercalli = 0 mercalli = INT(mercalli) 130 PRINT "ESCALA MERCALLI =" 131 IF mercalli > 11 THEN 1140 IF mercalli > 10 THEN 1130 IF mercalli > 9 THEN 1120 IF mercalli > 8 THEN 1110 IF mercalli > 7 THEN 1100 IF mercalli > 6 THEN 1090 IF mercalli > 5 THEN 1080 IF mercalli > 4 THEN 1070 IF mercalli > 3 THEN 1060 IF mercalli > 2 THEN 1050 IF mercalli > 1 THEN 1040 IF mercalli > .1 THEN 1030 IF mercalli < .1 THEN 1025 138 PRINT 139 INPUT "RESISTENCIA DE MATERIALES Y VICTIMAS (1)"; res IF res > 1 THEN END PRINT 140 PRINT "TIPOS DE CONSTRUCCION MAYORITARIA" 150 PRINT "TIPO A: BAJA CALIDAD: MAMPOSTERIA, BARRO = 0" PRINT "TIPO B: CALIDAD MEDIANA: LADRILLOS, VIGAS = 1" PRINT "TIPO C: ELEVADA RESISTENCIA, HORMIGON, EST.METALICAS = 2" INPUT "INTRODUZCA TIPO"; r PRINT INPUT "POBLACION DEL AREA "; POB IF hipo < 5 THEN hipo = 5 VICTIM = 10 ^ (RICH - (6 + r)) * (20 / hipo) IF VICTIM > 99 THEN VICTIM = 99 IF VICTIM < .01 THEN VICTIM = 0 EDIF = VICTIM * 9 IF EDIF > 99 THEN EDIF = 99 NUMV = VICTIM * POB / 100 PRINT "VICTIMAS SOBRE LA POBLACION ="; "<"; VICTIM; "%, NUM. MAX.="; INT(NUMV) PRINT "EDIFICACIONES DESTRUIDAS ="; (EDIF); "%" INPUT "PARA NUEVO CALCULO = 1 "; NUECA IF NUECA = 1 THEN 0 300 END 1000 INPUT "ENERGIA DEL TERREMOTO EQUIVALENTE EN EXPLOSION TM. DE TNT"; TNT RICHEP = (LOG(TNT) / LOG(10) + LOG(33) / LOG(10) * 2) / (LOG(33) / LOG(10)) PRINT "I. ESCALA RICHTER = "; RICHEP GOTO 60 1020 END 1025 PRINT "Imperceptible: S¢lo se registra en los sism¢grafos." GOTO 138 1030 PRINT "I. Muy d‚bil: Imperceptible para la mayor¡a excepto en condiciones favorables." GOTO 138 1040 PRINT "II. D‚bil: Perceptible s¢lo por algunas personas en reposo," PRINT " aquellas que se encuentran en pisos superiores deedificios." PRINT " Los objetos colgantes suelen oscilar. " GOTO 138 1050 PRINT "III. Leve: Perceptible por algunas personas dentro de los edificios, en pisos altos." PRINT " Muchos no lo reconocen como terremoto." PRINT " Los autom¢viles detenidos se mueven ligeramente." PRINT " Sensaci¢n semejante al paso de un cami¢n peque¤o." GOTO 138 1060 PRINT "IV. Moderado: Perceptible por la mayoria de personas dentro de los edificios y" PRINT " por pocas personas en el exterior durante el d¡a." PRINT " Durante la noche algunas personas pueden despertarse." PRINT " Perturbaci¢n en puertas y ventanas." PRINT " Las paredes suelen hacer ruido." PRINT " Los autom¢viles detenidos se mueven con m s energ¡a." PRINT " Sensaci¢n semejante al paso de un cami¢n grande." GOTO 138 1070 PRINT "V. Poco Fuerte: La mayor¡a de los objetos se caen." GOTO 138 1080 PRINT "VI. Fuerte: Lo perciben todas las personas, muchas asustadas salen" PRINT " al exterior. Paso insostenible. Ventanas, platos y cristaler¡as da¤ados." PRINT " Los objetos se caen de sus lugares, muebles movidos o ca¡dos." PRINT " Revoques da¤ados. Da¤os leves en estructuras." GOTO 138 1090 PRINT "VII. Muy fuerte: Pararse es dif¡cil. Da¤os insignificantes en estructuras" PRINT " de buen dise¤o y construcci¢n. De leves a moderados en estructuras" PRINT " ordinarias bien construidas y considerables en estructuras deficientes." PRINT " Mamposter¡a afectada." PRINT " Perceptible por personas en veh¡culos en movimiento." GOTO 138 1100 PRINT "VIII.Destructivo: Da¤os leves en estructuras t‚cnicas y considerables en" PRINT " estructuras ordinarias bien construidas, algunos colapsos de edificios." PRINT " Da¤o severo en estructuras pobremente construidas." PRINT " Mamposter¡a seriamente da¤ada o destruida." PRINT " Muebles completamente fuera de lugar." GOTO 138 1110 PRINT "IX. Ruinoso: P nico generalizado. Da¤os en estructuras resistentes," PRINT " paredes fuera de plomo. Grandes da¤os en importantes edificios," PRINT " con colapsos parciales. Edificios desplazados fuera de las bases." GOTO 138 1120 PRINT "X. Desastroso: Destruidas algunas estructuras de madera bien fabricadas." PRINT " Caen la mayor¡a de las estructuras de mamposter¡a." PRINT " El marco destruido con sus bases. Rieles doblados." GOTO 138 1130 PRINT "XI. Muy desastroso: Pocas estructuras de mamposter¡a permanecen en pie." PRINT " Puentes destruidos. Rieles curvados en gran medida." GOTO 138 1140 PRINT "XII. Catastr¢fico: Destrucci¢n total con pocos sobrevivientes." PRINT " Los objetos saltan al aire." PRINT " Los niveles y perspectivas quedan distorsionadas." PRINT " Cambios profundos en los relieves." GOTO 138 PRINT "" GOTO 138 4000 PRINT 4005 PRINT "DETERMINACION DE LA DISTANCIA A PARTIR DE LOS RETRASOS ONDAS S-P" 4010 INPUT "ENERGIA LOCAL RICHTER"; eloc PRINT PRINT "TIEMPO S-P OBSERVADO EN SISMOGRAFO SI ORIENTACION OPTIMA, SINO MENOR" INPUT "RETRASO EN SEGUNDO ONDAS S-P"; RET TPM = RET * 1.25 / 60 PRINT "RETRASO SE¥AL EN MINUTOS="; TPM DISG = -4.3 / (1 / 36 - 1 / TPM): disk = DISG * 111.1111 PRINT "DISTANCIA EN GRADOS Y KM.="; DISG; disk ep = (-25 * eloc - SQR(25 * eloc * 25 * eloc + 100 * disk)) / -50 PRINT "RICHTER EPICENTRO MEDIO(20KM.PROF.)="; ep RICHEP = ep dist = disk unt = 1 INPUT "PROSEGUIR CON MAS DATOS (1)"; PTU IF PTU = 1 THEN GOTO 62 END 10000 PRINT PRINT "SADEYA/INFORTECNICA" PRINT "PROBABILIDAD POISSON DE N SUCESOS (TERREMOTOS, ERUPCIONES..." PRINT "A PARTIR DE UN TIEMPO DE RECURRENCIA Y DE UNA MAGNITUD" PRINT INPUT "NUMERO DE SUCESOS ESPERADOS"; N INPUT "PERIODO EN UNIDAD DE TIEMPO (ANS,DIAS,ETC...) DEL ESTUDIO"; T INPUT "EN LA MISMA UNIDAD TIEMPO DE RECURRENCIA DEL SUCESO"; tr INPUT "MAGNITUD ESCALA R=10(RICHTER, VEI...)DEL SUCESO ESPERADO"; MAGP INPUT "MAGNITUD ESCALA R=10 DEL SUCESO OBSERVADO"; MAGR PRINT REM TERMINO TER tr = 10 ^ (MAGP - MAGR) * tr REM ter = (1 / (TR * 10 ^ (magr - MAGP)) * (T)) ter = (T / tr) REM FACTORIAL M = N TOT = M 10030 M = M - 1 TOT = TOT * M IF M > 2 THEN 10030 IF TOT < 2 THEN TOT = 1 REM PRINT TOT AA = 1 / TOT AB = (ter) ^ N AC = 2.71828183# ^ -ter PR = AA * AB * AC REM PR = 1 - PR PRINT "PROB.P-Q DE"; N; "SUCESOS DE MAG."; MAGP; "EN TIEMPO"; T; " ="; PR; 1 - PR REM PR2 = 1 - (1 - (ter)) ^ T PR3 = 1 - 2.71828183# ^ -ter PRINT "PROBABILIDAD PARA UN SUCESO="; PR3 PRINT INPUT "NUEVO CALCULO (1)"; NUC IF NUC = 1 THEN 40 END 2 CLS 15000 PRINT "SADEYA INFORTECNICA" PRINT PRINT "CAUSAS DE UN TERREMOTO" INPUT "METEORITO(1) DESPLOME(2) EXPLOSION-LIBERACION(3) ERUPCION(4)"; TEP PRINT IF TEP = 4 THEN 40000 IF TEP = 2 THEN 30000 IF TEP = 3 THEN 35000 INPUT "VELOCIDAD INFINITA O RELATIVA KM/S"; VEL INPUT "DENSIDAD (AGUA=1): COMETAS DE .3 A 1, ASTEROIDES 2.5 A 3"; DEN INPUT "DIAMETRO EN METROS"; DIA MASA = (4 / 3) * 3.1416 * (DIA / 2) ^ 3 * DEN VOLU = (4 / 3) * 3.1416 * (DIA / 2) ^ 3 VOLUK = VOLU / 1000000000 NUMER = 5000 / VOLUK vef = SQR(VEL * VEL + 11.2 * 11.2) PRINT jul = ((MASA * 1000) / 2) * (vef * 1000 * vef * 1000) PA = .001 * DIA ^ 3: PC = 1 / PA PRINT "METEOROIDE" PRINT "VELOCIDAD DE ENTRADA O IMPACTO KM/S="; vef PRINT "MASA EN TONELADAS="; MASA PRINT "ENERGIA JOULS, KILOT ="; jul; (jul / 850000000000#) / 4.2# TNT = (jul / 850000000000#) / 4.2# 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)) PRINT "I. ESCALA RICHTER = "; RICHEP PRINT GOTO 40 30000 PRINT "DESPLOME TECTONICO" INPUT "LONGITUD, ANCHURA, GROSOR MED. METROS. MATERIAL DESPLOME"; LONGI, ANCH, GROS VOL = LONGI * ANCH * GROS INPUT "DENSIDAD (AGUA=1)"; DEN INPUT "ALTURA DESPLOME METROS"; ESP MASA = VOL * DEN * 1000 TIE = SQR(ESP / (9.8 * .5)) vef = 9.8 * TIE jul = (MASA / 2) * (vef * vef) PRINT "DESPLOME" PRINT "VELOCIDAD DE IMPACTO M/S="; vef PRINT "MASA DESPLOME EN TONELADAS="; MASA PRINT "ENERGIA JOULS, KILOT ="; jul; (jul / 850000000000#) / 4.2# TNT = (jul / 850000000000#) / 4.2# 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)) PRINT "I. ESCALA RICHTER = "; RICHEP PRINT GOTO 40 35000 PRINT "CAUSA INESPECIFICA: EXPLOSION, LIBERACION TECTONICA..." INPUT "EN JOULES(1) EN KTM DE TNT(2)"; TIP IF TIP = 2 THEN 35020 INPUT "ENERGIA EN JOULS"; jul TNT = (jul / 850000000000#) / 4.2# GOTO 35030 35020 INPUT "ENERGIA DE LA EXPLOSION EN KTM. DE TNT"; TNT 35030 TNT = TNT * 1000 RICHEP = (LOG(TNT) / LOG(10) + LOG(33) / LOG(10) * 2) / (LOG(33) / LOG(10)) PRINT "I. ESCALA RICHTER = "; RICHEP PRINT GOTO 40 40000 PRINT "ERUPCION VOLCANICA " INPUT "ENERGIA ERUPCION ESCALA VEI="; escv PRINT ijo0 = 1 / 12 * (1.7 * 1D+18) * 10 ^ (escv - 5) IMEG0 = ijo0 * 239 * (1 / 1D+18) * .5 PRINT "ENERGIA EXPLOSION ERUPTIVA EN JOUL, MEG.="; ijo0; IMEG0 ijo = (1.7 * 1D+18) * 10 ^ (escv - 5) IMEG = ijo * 239 * (1 / 1D+18) * .5 PRINT "ENERGIA GLOBAL LIBERADA JOUL, MEG.="; ijo; IMEG TNT = IMEG * 100000 TNT2 = IMEG0 * 10000 RICHEP = (LOG(TNT) / LOG(10) + LOG(33) / LOG(10) * 2) / (LOG(33) / LOG(10)) richep2 = (LOG(TNT2) / LOG(10) + LOG(33) / LOG(10) * 2) / (LOG(33) / LOG(10)) PRINT "I. ESCALA RICHTER EXPLOSION Y GLOBAL= "; richep2; RICHEP PRINT GOTO 40