1 CLS
PRINT "SADEYA/INFORTECNICA"

PRINT "CONSTANTES ESTELARES Y RADIACION RECIBIDA POR UN PLANETA"
PRINT

INPUT "EXPONENTE MASA^EXP=LUM SER.PRINC=3.5 GIG=4-5..."; EPT
2 INPUT "Distancia estrella-planeta en Km (1) en u.a. (2) a.luz(3) R.Sol(4)"; td
IF td = 1 THEN 8
IF td = 3 THEN 2000
IF td = 4 THEN 3000

INPUT "Distancia en u.a."; di
dist = di * 149600000
GOTO 9
8  INPUT "Distancia Km"; dist
di = dist / 149600000

9 INPUT "Estrella: entrar masa (1) o luminosidad (2)"; entr
IF entr = 1 THEN 10
INPUT "Luminosidad estrella (Sol=1)"; lum
masa = lum ^ (1 / EPT)
GOTO 11
10 INPUT "Masa de ls estrella (Sol=1)"; masa
lum = masa ^ EPT

11 PRINT "Masa y luminosidad estrella (Sol=1)="; masa; lum
temp = (LOG(lum) / LOG(10) + 6) * 1000 - 220
IF temp < 0 THEN temp = 1
mav = (LOG(lum) / LOG(2.512) - 4.8) * -1
dp = 206264.806248# / (dist / 149600000)
map = mav - 5 - 5 * LOG(dp) / LOG(10)

PRINT "Magnitud absoluta estrella, aparente="; mav; map


PRINT "Temperatura media estrella serie principal k="; temp
INPUT "(Si temp.med.=0). Temperatura superficie k (Sol=5780)"; tems
IF tems = 0 THEN tems = temp
IF tems < 0 THEN tems = 0
lum = masa ^ EPT
cl = (tems ^ 4) / (5780 ^ 4)
rae = SQR(lum / cl) * 348000 * 2
atnn = rae / dist
rap = ATN(atnn) * 180 / 3.1416 * 60 * 2
rasp = rap / 15.99373 / 2
PRINT "Radio estelar te¢rico y aparente ="; rae; "Km."; rasp; "Sol=1"; rap; "d.min."
INPUT "( No modificar radio=0). Radio estrella en Km. (Sol=696000)"; rade
IF rade = 0 THEN rade = rae

clp = (temsp ^ 4) / (286 ^ 4)
REM temp = temsp * SQR((radep * 2) / (distp * 2))
REM PRINT "Temperatura planeta al n£cleo, clp="; temp; clp
INPUT "Incremento radiaci¢n (entre -5 y +5 w/m2) en min/max solares (0 sin)"; incr
tem = tems * SQR((rade * 2) / (dist * 2))
cel = tem - 273
20 PRINT "Temperatura estelar en la zona del planeta: K, celsius ="; tem; cel
REM wat = tem * 4.85
lum = masa ^ 3
watl = 1366 * lum
kw = (watl + incr) / ((dist / 149600000) ^ 2)

PRINT "K: Energ¡a estelar incidente en el planeta watios/m2 ="; kw
kt = kw / tem
REM PRINT "K/K0, K/tem="; kw / 1366; kt
PRINT
INPUT "Calor interno del planeta por teor¡a fisiocinetica Si=1 No=2"; teor
IF teor = 1 THEN 27
GOTO 29


27 INPUT "Masa del planeta (Tierra=1)"; masap
INPUT "Exponente fision*cinet (est ndar=2 Venus=2.5)"; expo

cinet = (masap * masa) / (di ^ 3)
cinet2 = (masap * masa) / (di ^ expo)

INPUT "Radio del planeta en Km. (Tierra=6370)"; radep
rap = radep / 6370

tempp = cinet2 * 278 * 1.13
ct1 = tempp - 273
tems = tempp / (rap ^ 2)
ct2 = tems - 273
tempp2 = cinet2 * tem * 1.13
ct12 = tempp2 - 273
tems2 = tempp2 / (rap ^ 2)
ct22 = tems2 - 273
PRINT "Energ¡a cin‚tica del planeta (Tierra=1)"; cinet2
REM print "Tem.planet(K,C)="; tempp; ct1; "Superf/gas="; tems; ct2
IF cinet2 < 1 THEN 29
REM print "Astro exot‚rmico";" Tem.planet(K,C)="; tempp2; ct12; "Superf/gas="; tems2; ct22 

PRINT "Teor¡a fisiocin‚tica: Astro exot‚rmico."; "t.superf-ambiente= "; tems2 - tem
PRINT "Tem.planet(K,C) RadioTierra="; tempp2; ct12; "Superf alb0, gas="; tems2; ct22

29 PRINT
INPUT "Si seleccion de albedo=2. Albedo estandar del planeta="; albe
IF albe = 2 THEN 50
30 INPUT "Si selecci¢n (=200). Reducci¢n radiancia (x/100) por erupci¢n o impacto="; redu
IF redu = 200 THEN 100
32 albe2 = albe + (redu / 100)
REM INPUT "Radio del planeta Km"; rad
enerintercep = kw * 3.1416 * rad ^ 2
enerabsor = kw * 3.1416 * rad ^ 2 * (1 - albe)
radref = kw * albe
radref2 = kw * albe2
radabso = kw - radref

radabso2 = kw - radref2
ro = 5.67 * 10 ^ -8 / (kw / 1366)
te = (1366 * (1 - albe2)) / (4 * ro)
te = SQR(te)
tet = SQR(te)
rads = radabso * .25
rads2 = radabso2 * .25
INPUT "Calor superficial interno o retenido planeta C. (Tierra=35)"; cop
REM tef = rads * cop * 1.258
PRINT "Albedo circunstancial ="; albe2
PRINT "Radiaci¢n reflejada limbo iluminado w/m2="; radref2
PRINT "Radiaci¢n absorbida (retenida) limbo iluminado w/m2="; radabso2
it = radabso / tem
PRINT "Radiaci¢n retenida/unidad de superficie planeta w/m2="; rads2
REM PRINT "radabs/tem="; it
kw2 = kw * (1 - albe2)
tem2 = kw2 / kt
tef = tet + cop
PRINT "Temperatura albedo: k,cent.="; tet; tet - 273
REM PRINT "Temperatura - albedo: K,cent.="; tem2; tem2 - 273
PRINT "Temperatura efectiva: K, cent. ="; tef; tef - 273
PRINT
40 GOTO 2
50 PRINT "Nubes visuales o de albedo: En la Tierra 20% para 15§C. de"
PRINT "temperatura media. Pulsar 200 para otra temperatura"
INPUT "Porcentaje de nubes de albedo, limbo iluminado"; nubes
IF nubes = 200 THEN 1000

51 INPUT "Si hubiera, porcentaje superficie de nubes volc nicas"; nube2
nubes = nubes + nube2 * .66
INPUT "Porcentaje superficie de hielo y nieve"; nieve
INPUT "Porcentaje resto de continentes despejados"; tierra
mar = 100 - nieve - nubes - tierra
PRINT "Porcentaje de superficie de oceano despejado="; mar
albe = nubes / 100 * .75 + nieve / 100 * .65 + tierra / 100 * .23 + mar / 100 * .15
mar2 = (mar / (tierra + mar)) * nubes
REM PRINT "mar2"; mar2
radam = (mar * .85) / (100 * (1 - albe))
radam2 = (mar2 * .25) / (100 * (1 - albe))
PRINT "Radiaci¢n solar directa y a trav‚s de nubes retenida por mares %="; (radam + radam2) * 100
PRINT "Albedo medio ="; albe
GOTO 30

100 INPUT "Erupci¢n 1, Impacto meteorito 2."; tip
IF tip = 1 THEN 140
IF tip = 2 THEN 200
140 INPUT "Escala VEI (1-8) erupci¢n "; vei
 wats = 1.85 * 2.16 ^ (vei - 4)
 redu = wats / kw * 100
150 PRINT "Reducci¢n radiancia wats/m2="; wats; " % ="; redu
 GOTO 32
200 INPUT "Diametro del asteroide en metros"; diame
    redu = diame / 100 * 1.2
    wats = redu / 100 * kw
    GOTO 150

1000 INPUT "Temperatura media asumida "; temp
nubes = .026666666# * temp ^ 2 + .93333333# * temp
IF nubes > 98 THEN nubes = 98
PRINT "Porcentaje de nubes de albedo="; nubes

GOTO 51

2000 INPUT "Distancia en A.l."; dl
dist = dl * 9460730472581#
di = dist / 149600000
GOTO 9
3000 INPUT "Distancia en R.Sol"; ds
dist = ds * 696000#
di = dist / 149600000
GOTO 9