10 CLEAR
20 CLS
30 PRINT TAB(21); "************ EFEMERIDES ***********"
40 PRINT
50 PRINT "                    ___________________________________"
60 PRINT
70 PRINT "                           2007 SADEYA/INFORTECNICA      "
80 PRINT "                            LEPANT 422    "
90 PRINT "                            BARCELONA  08025    "
95  PRINT "                           Tel.93 433 42 43"
100 PRINT "                    www.sadeya.org www.infortecnica.com"
110 PRINT "                    ___________________________________  "
115 PRINT
116 PRINT "                       SOFTWARE SERIE ECOLOGIA    "
117 PRINT "   USO PUBLICO PARA UN MEJOR CONOCIMIENTO DEL CAMBIO CLIMATICO "
120 PRINT
130 PRINT "ESTE PROGRAMA MODELIZA UNA ERUPCION VOLCANICA Y DETERMINA "
132 PRINT "LAS CONSECUENCIAS SOBRE EL CLIMA"
140 PRINT
150 PRINT
142 INPUT "ERUPCIONES SUPERFICIALES (1) ABISALES FISULARES (2)"; TIPU
    CLS
    IF TIPU = 2 THEN 10000
160 INPUT "ESCALA VEI (1), VOLUMEN ERUPCION (2)"; ENT
    IF ENT = 2 THEN 164
    INPUT "ESCALA VEI (0-8 O, SUPERIOR)"; escv
    KM = 10 ^ (escv - 5)
PRINT "VOLUMEN KM3. MATERIALES EYECTADOS="; KM

    GOTO 170

164 INPUT "VOLUMEN DE MATERIALES EYECTADOS EN M3 (1) KM3 (2)"; TKM
    INPUT "ENTRAR VOLUMEN"; KM
    IF TKM = 1 THEN KM = KM / 1000000000

170 PRINT
    escv = LOG(KM) / LOG(10) + 5
    IF escv < .000001 THEN escv = 0

    KJ = KM
    KL = (KJ / 500) ^ 2
    IF KJ > 500 THEN KJ = KJ * KL
    PERID = KJ * 6
    PLUMAS = escv * escv
    PLUMA = PLUMAS
    IF PLUMA > 25 THEN PLUMA = 25
    SUPER = 6300 ^ 2 * 4 * 3.1416
    TMAT = KM * 2000000000
    TVAP = TMAT * .1
    TCO2 = TVAP * .1
    TSOLF = TCO2 * .3
    GSOLF = TSOLF * 1000 / SUPER
    OGA = TSOLF * .2
    SULF = KM * 2000000
175 TEM = (KM ^ (1 / 3)) / 3
177 radia = 20 * TEM * .6
178 ALBED = radia / 1360 * 100
    PRINT "ESCALA VEI="; escv
    ijo = (1.7 * 1D+18) * 10 ^ (escv - 5)

    imeg = ijo * 239 * (1 / 1D+18) * .5
    PRINT "ENERGIA GLOBAL LIBERADA JOUL, MEG.="; ijo; imeg
    PRINT "ENERGIA MAX.EXPLOSION JOUL, MEG.="; ijo / 25; imeg / 25
    PRINT "PERIODICIDAD A„OS = "; PERID
    PRINT "PLUMA EFECTIVA KM. ALTURA, TEORICA="; PLUMA; PLUMAS
    PRINT "MASA GLOBAL DE MATERIALES EXPULSADA TM="; TMAT
    PRINT "MASA DE VAPOR DE AGUA Y CO2 TM ="; TVAP; TCO2
    PRINT "MASA DE SO2 Y OTROS GASES TM="; TSOLF; OGA
    PRINT "KM.SO2/KM2 LOCAL, POLOS="; GSOLF; GSOLF / 100
    PRINT
    PRINT "TIPO MEDIO ERUPCION"
   
    IF escv > 8 THEN 9000
    IF escv > 7 THEN 8000
    IF escv > 6 THEN 7000
    IF escv > 5 THEN 6000
    IF escv > 4 THEN 5000
    IF escv > 3 THEN 4000
    IF escv > 1.7 THEN 3000
    IF escv > .5 THEN 2000
    IF escv > 0 THEN 1000
179 REM
    IF KM < 1.8 THEN 235

    REM PRINT
    PRINT "ERUPCION INTENSA CON CONSECUENCIAS GLOBALES"
180 PRINT "DESCENSO TEMPERATURA MEDIA DE LA TIERRA GRD. DURANTE PRIMER A„O="; TEM
    IF TEM > 3 THEN PRINT "CIRCUNSTANCIAS DE GLACIACION"
181 IF radia > 165 THEN radia = 165
    PRINT "DESCENSO RADIANCIA MEDIA DE LA TIERRA EN W/m2 PRIMER A„O="; radia
182 PRINT "INCREMENTO ALBEDO MEDIO DE LA TIERRA DURANTE EL PRIMER A„O %="; ALBED
190 CIC = KM * 1.35
200 DES = KM / 5
210 ENE = KM * 15
205 MIN = KM * .48
220 PRINT "PERIODO DESCENSO TEMPERATURA EN SUPERFICIE CONTINENTAL EN A„OS="; DES
225 REM PRINT "PERIODO EN A„OS PARA ALCANZAR LA MINIMA TEMPERATURA="; MIN
230 PRINT "DURACION DEL CICLO MAR-ATMOSFERA HASTA RECUPERACION EN A„OS="; CIC
    PRINT "DESC.TEMP.MAR CICLO (CUARTO 1§,2§,3§,4§)="; TEM / 2; TEM / 4; TEM / 6; TEM / 10

235 PRINT
240 INPUT "NUEVO CALCULO SI = 1"; NC
250 IF NC = 1 THEN 140
260 END
1000 PRINT "Erupciones no explosivas, volcanes en actividad continua"
     PRINT "Tipo hawuaiano, Kilauea."
     GOTO 179
2000 PRINT "Erupciones poco explosivas, ligeras de lava fluida"
     PRINT "Tipo principal stromboliano."
     GOTO 179
3000 PRINT "Erupciones explosivas, volcanes activos"
     PRINT "Tipo stromboliano-vulcaniano."
     GOTO 179
4000 PRINT "Erupciones violentas, volcanes menos activos, pirocl sticos"
     PRINT "Tipo vulcaniano, sub-pliniana."
     GOTO 179
5000 PRINT "Erupciones cataclismicas plinianas, volcanes poco activos,"
     PRINT "De lava densa, en general con invasi¢n de agua en la c mara."
     GOTO 179
6000 PRINT "Erupciones paroxisticas, superplinianas."
     PRINT "Grandes c maras magm ticas invadidas por el agua."
     GOTO 179
7000 PRINT "Erupciones colosales, ultra-plinianas."
     PRINT "Gigantescas c maras magm ticas invadidas por el agua."
     GOTO 179
8000 PRINT "Supervolcanes, extensas  reas de la corteza afectadas por."
     PRINT "grandes camaras magm ticos subterr neas."
     GOTO 179
9000 PRINT "Erupciones mega-colosales, en general producidas por choques"
     PRINT "de grandes asteroides y fusi¢n de la corteza."
     GOTO 179
     GOTO 235
10000 PRINT "ERUPCIONES FISULARES ABISALES"
PRINT "Erupciones poco perceptibles pero muy intensas y continuas."
PRINT "Se manifiestan por s£bitas elevaciones de la temperatura oce nica"
PRINT "principalmente en el Pac”fico. Pueden alcanzar por su intensidad"
PRINT "de energ”a valores equivalentes y superiores a VEI=8."
PRINT
10012 INPUT "AREA DEL OCEANO ESTUDIADA COMPARADA CON LA SUPERFICIE DE LA TIERRA"; AREA
SUPTC = 4 * 3.1416 * 637000000 ^ 2
SUPES = AREA * SUPTC
VOLC = SUPES * 370000
INPUT "INCREMENTO TEMPERATURA EN GRADOS C"; INCT
CALOR = INCT * VOLC
jouls = CALOR / 4.1668
MEGAT = 1
MEG = jouls / (4.18 * 1000000000000000#)
vei5 = (LOG(jouls / (1.7 * 10 ^ 18))) / LOG(10)
vei = vei5 + 5
PRINT
PRINT "ENERGIA LIBERADA CALORIAS="; CALOR
PRINT "ENERGIA LIBERADA JOULS="; jouls
PRINT "ENERGIA LIBERADA MEG="; MEG
PRINT "GRADO VEI="; vei
PRINT
INPUT "NUEVO CALCULO (1), INICIO(2)"; NUCA
IF NUCA = 1 THEN 10012
IF NUCA = 2 THEN 10

END