from astropy import units as u atec = exocounts.InstClass() atec.lamb = 1.7 * u.micron #micron R = 100000.0 atec.dlam = atec.lamb / R #micron #atec.dtel = 9.1 #telescope diameter m atec.dtel = 2.0 * u.m #telescope diameter m atec.dstel = 0.00 * u.m #secondary telescope diameter m atec.throughput = 0.05 atec.ndark = 0.17 / u.s #dark current atec.nread = 14.0 #nr atec.fullwell = 80000. target = exocounts.TargetClass() target.name = "HR8799bcde" target.teff = 1000.0 * u.K #K target.rstar = 0.1 * const.R_sun #Rsolar target.d = 39 * u.pc #pc obs = exocounts.ObsClass(atec, target) obs.texposure = 1.0 * 24.0 * 1.0 * u.h #1 d obs.tframe = 7.1 * u.s #time for one frame [sec] obs.napix = 15 # number of the pixels in aperture obs.mu = 1 S = 1.8 * 1.8 * np.pi #core size obs.effnpix = S / 3.0 #3 is an approx. increment factor of PSF obs.mu = 1
import exocounts
import convmag
import nstar
import numpy as np
from astropy import constants as const
from astropy import units as u
magdict = convmag.get_magdict()
alphacenA = exocounts.TargetClass()
alphacenA.name = "alpha cen A"
alphacenA.teff = 5790.0 * u.K #K ## optimistic (no cloud)
alphacenA.rstar = 1.2234 * const.R_sun
alphacenA.d = 1.34 * u.pc #pc
fluxA = nstar.getflux(alphacenA, 10.47 * u.micron)
print(alphacenA.name)
print(fluxA)
print("N mag=", convmag.get_mag("N", fluxA, magdict))
alphacenB = exocounts.TargetClass()
alphacenB.name = "alpha cen B"
alphacenB.teff = 5260.0 * u.K #K ## optimistic (no cloud)
alphacenB.rstar = 0.8632 * const.R_sun
alphacenB.d = 1.34 * u.pc #pc
fluxB = nstar.getflux(alphacenB, 10.47 * u.micron)
print(alphacenB.name)
print(fluxB)
import exocounts
import convmag
import nstar
import numpy as np
import planet
from astropy import constants as const
from astropy import units as u
lamMIR = 10.47 * u.micron
lamNIR = 1.5 * u.micron
magdict = convmag.get_magdict()
print("=========================")
tauceti = exocounts.TargetClass()
tauceti.name = "tau ceti"
print(tauceti.name)
print("-Emission")
tauceti.teff = 255.0 * u.K #K ## optimistic (no cloud)
tauceti.rstar = 1.0 * const.R_earth
tauceti.d = 3.65 * u.pc #pc
fluxB = nstar.getflux(tauceti, lamMIR)
fluxBj = nstar.getfluxJy(tauceti, lamMIR)
fluxBp = nstar.getfluxph(tauceti, lamMIR)
print(fluxB)
print(fluxBj)
print(fluxBp)
print("cts / J", fluxBp.value / fluxB.value)
