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