AstroToolkit 1.0.0__py3-none-any.whl

AstroToolkit/Figures/HRD.py

@@ -0,0 +1,82 @@
+from astropy.table import Table
+import matplotlib.pyplot as plt
+import numpy as np
+from scipy import stats
+import math
+import cmasher as cmr
+from bokeh.plotting import figure,show
+import os
+from importlib_resources import files
+
+from ..Surveys.Gaia import GaiaQueryDesignation
+
+fits_file = files('WDPlanetsToolkit.Figures').joinpath('backdrop_hrd_allmags.fits')
+
+def get_plot(source=None,sources=None):
+	background=Table.read(fits_file).to_pandas()
+	background_100pc=background.loc[background['type']==b'100pc']
+	background_wd=background.loc[background['type']==b'WD']
+	background_cv=background.loc[background['type']==b'CV']
+	background_wd_dm=background.loc[background['type']==b'WD+dM']
+	background_sd=background.loc[background['type']==b'SD']
+	
+	background_split=[background_100pc,background_wd,background_wd_dm,background_cv,background_sd]
+	
+	x_arr=[]
+	y_arr=[]
+	colour_arr=[]
+	
+	for sample in background_split:
+		x=sample['phot_bp_mean_mag'].values-sample['phot_rp_mean_mag'].values
+		y=sample['phot_g_mean_mag'].values+5*np.log10(sample['parallax'].values/1000)+5
+		
+		x_arr.append(x)
+		y_arr.append(y)
+	
+	colour_maps=['Greens','Purples','Reds','Blues','YlOrBr']
+	
+	for i in range(0,len(colour_maps)):
+		colour_maps[i]=cmr.get_sub_cmap(colour_maps[i],0.4,1.0)
+
+	for i in range(0,len(x_arr)):
+		values=np.vstack([x_arr[i],y_arr[i]])
+		kernel=stats.gaussian_kde(values)
+		weights=kernel(values)
+		weights=weights/weights.max()
+		
+		colour=["#%02x%02x%02x" % (int(r), int(g), int(b)) for r, g, b, _ in 255*colour_maps[i](weights)]
+		colour_arr.append(colour)
+
+	plot=figure(width=400,height=400,title='Gaia HRD',x_axis_label=r'\[\text{bp-rp}\]',y_axis_label=r'\[\text{abs g}\]')
+	plot.y_range.flipped = True
+
+	labels=['100pc','WDs','CVs','WD+dM','SDs']
+	for i in range(0,len(x_arr)):
+		plot.scatter(x_arr[i],y_arr[i],size=3,color=colour_arr[i],legend_label=labels[i])
+
+	if source!=None:
+		data=GaiaQueryDesignation(Designation=source)
+		gmag,bpmag,rpmag,parallax=data['phot_g_mean_mag'].values,data['phot_bp_mean_mag'].values,data['phot_rp_mean_mag'].values,data['parallax'].values
+		x=bpmag-rpmag
+		y=gmag+5*np.log10(parallax/1000)+5
+		
+		if not math.isnan(gmag) and not math.isnan(bpmag) and not math.isnan(rpmag) and not math.isnan(parallax):
+			plot.scatter(x=x,y=y,marker='square_dot',line_color='black',fill_color=None,size=20,line_width=2,legend_label='Source')
+		else:
+			return None
+
+	elif sources!=None:
+		for i in range(0,len(sources)):
+			data=GaiaQueryDesignation(Designation=sources[i])
+			gmag,bpmag,rpmag,parallax=data['phot_g_mean_mag'].values,data['phot_bp_mean_mag'].values,data['phot_rp_mean_mag'].values,data['parallax'].values
+			x=bpmag-rpmag
+			y=gmag+5*np.log10(parallax/1000)+5
+			
+			if not math.isnan(gmag) and not math.isnan(bpmag) and not math.isnan(rpmag) and not math.isnan(parallax):
+				plot.scatter(x=x,y=y,marker='square_dot',line_color='black',fill_color=None,size=20,line_width=2,legend_label='Source')
+			else:
+				pass
+
+	plot.legend.click_policy="hide"
+
+	return plot

AstroToolkit/Figures/SED.py

@@ -0,0 +1,364 @@
+import astropy.units as u
+import math
+import pandas as pd
+import matplotlib
+import numpy as np
+from bokeh.plotting import figure,show
+from bokeh.models import Whisker, ColumnDataSource
+
+from ..Surveys.Gaia import GaiaQueryDesignation
+from ..Miscellaneous.ProperMotionCorrection import PMCorrection as CorrectPM
+from ..Surveys.SDSS import get_photometry as get_phot_sdss
+from ..Surveys.WISE import get_photometry as get_phot_wise
+from ..Surveys.TWOMASS import get_photometry as get_phot_twomass
+from ..Surveys.Gaia import GaiaGetPhotometryCoords
+from ..Surveys.GALEX import GALEXGetPhotometryCoords
+from ..Surveys.PanSTARRS import PanSTARRSGetPhotometryCoords
+from ..Surveys.SkyMapper import SkyMapperGetPhotometryCoords
+
+def zeropoint(mag,zp,wl):
+	flux=zp*10**(-0.4*mag)
+	c=2.988*10**18
+	fnl=1*10**(-23)
+	flux=flux/((fnl*c)/wl**2)*1000
+	return flux
+
+def get_plot(source=None,pos=None,radius=3):
+	if source!=None and pos!=None:
+		raise Exception('simulatenous source and pos input detected')
+	
+	if source!=None:
+		gaia_data=GaiaQueryDesignation(Designation=source)
+		ra,dec,pmra,pmdec=gaia_data['ra'].values[0],gaia_data['dec'].values[0],gaia_data['pmra'].values[0],gaia_data['pmdec'].values[0]
+		
+		pos_corrected=CorrectPM([2016,0],[2012,0],ra,dec,pmra,pmdec)
+		ra_panstarrs,dec_panstarrs=pos_corrected[0],pos_corrected[1]
+		
+		pos_corrected=CorrectPM([2016,0],[2016,0],ra,dec,pmra,pmdec)
+		ra_skymapper,dec_skymapper=pos_corrected[0],pos_corrected[1]
+		
+		ra_gaia,dec_gaia=ra,dec
+		
+		pos_corrected=CorrectPM([2016,0],[2007,0],ra,dec,pmra,pmdec)
+		ra_galex,dec_galex=pos_corrected[0],pos_corrected[1]
+		
+		pos_corrected=CorrectPM([2016,0],[1991,0],ra,dec,pmra,pmdec)
+		ra_rosat,dec_rosat=pos_corrected[0],pos_corrected[1]
+		
+		pos_corrected=CorrectPM([2016,0],[2017,0],ra,dec,pmra,pmdec)
+		ra_sdss,dec_sdss=pos_corrected[0],pos_corrected[1]
+		
+		pos_corrected=CorrectPM([2016,0],[2010,5],ra,dec,pmra,pmdec)
+		ra_wise,dec_wise=pos_corrected[0],pos_corrected[1]
+		
+		pos_corrected=CorrectPM([2016,0],[1999,0],ra,dec,pmra,pmdec)
+		ra_twomass,dec_twomass=pos_corrected[0],pos_corrected[1]
+		
+	elif pos!=None:
+		ra_panstarrs,dec_panstarrs=pos[0],pos[1]
+		ra_skymapper,dec_skymapper=pos[0],pos[1]
+		ra_gaia,dec_gaia=pos[0],pos[1]
+		ra_galex,dec_galex=pos[0],pos[1]
+		ra_rosat,dec_rosat=pos[0],pos[1]
+		ra_sdss,dec_sdss=pos[0],pos[1]
+		ra_wise,dec_wise=pos[0],pos[1]
+		ra_twomass,dec_twomass=pos[0],pos[1]
+		
+	else:
+		raise Exception('either source or pos input required')
+		return None
+	
+	filter_arr=[]
+	
+	error_arr=[]
+	wavelength_arr=[]
+	flux_arr=[]
+	mag_arr=[]
+	
+	wavelengths=[]
+	flux=[]
+	errors=[]
+	mags=[]
+	photometry=GaiaGetPhotometryCoords(ra_gaia,dec_gaia,radius=radius)
+	if isinstance(photometry,pd.DataFrame):
+		filter_list=['phot_g_mean_mag','phot_rp_mean_mag','phot_bp_mean_mag']
+		error_list=['phot_g_mean_mag_error','phot_rp_mean_mag_error','phot_bp_mean_mag_error']
+		zero_points=[2.5e-9,1.24e-9,4.11e-9]
+
+		#g,rp,bp
+		filter_ref=[5850.88,7690.74,5041.61]
+		
+		for filter, error, zp in zip(filter_list,error_list,zero_points):
+			try:
+				mag=photometry._get_value(0,filter)
+				mags.append(mag)
+
+				errors.append((photometry._get_value(0,error)*u.ABmag))
+				wavelength=(filter_ref[filter_list.index(filter)])
+				wavelengths.append(wavelength)
+
+				fl=zeropoint(mag=mag,zp=zp,wl=wavelength)
+				flux.append(fl)
+
+				filter_arr.append(filter)
+			except:
+				pass
+	wavelength_arr.append(wavelengths)
+	flux_arr.append(flux)
+	error_arr.append(errors)
+	mag_arr.append(mags)
+
+	wavelengths=[]
+	flux=[]
+	errors=[]
+	mags=[]
+	photometry=GALEXGetPhotometryCoords(ra_galex,dec_galex,radius=radius)
+	if isinstance(photometry,pd.DataFrame):
+		filter_list=['NUVmag','FUVmag']
+		error_list=['e_NUVmag','e_FUVmag']
+
+		#nuv,fuv
+		filter_ref=[2303.37,1548.85]
+		
+		for filter,error in zip(filter_list,error_list):
+			try:
+				mag=photometry._get_value(0,filter)
+				mags.append(mag)
+				
+				errors.append((photometry._get_value(0,error)*u.ABmag))
+				wavelength=(filter_ref[filter_list.index(filter)])
+				wavelengths.append(wavelength)
+					
+				zp=10**((5*np.log10(filter_ref[filter_list.index(filter)])+2.406)/-2.5)
+				fl=zeropoint(mag=mag,zp=zp,wl=wavelength)
+				flux.append(fl)
+					
+				filter_arr.append(filter)
+			except:
+				pass
+	wavelength_arr.append(wavelengths)
+	flux_arr.append(flux)
+	error_arr.append(errors)
+	mag_arr.append(mags)
+
+	wavelengths=[]
+	flux=[]
+	errors=[]
+	mags=[]
+	photometry=PanSTARRSGetPhotometryCoords(ra_panstarrs,dec_panstarrs,radius=radius)
+	if isinstance(photometry,pd.DataFrame):
+		filter_list=['gMeanPSFMag','rMeanPSFMag','iMeanPSFMag','zMeanPSFMag','yMeanPSFMag']
+		error_list=['gMeanPSFMagErr','rMeanPSFMagErr','iMeanPSFMagErr','zMeanPSFMagErr','yMeanPSFMagErr']
+
+		#g,r,i,z,y
+		filter_ref=[4810.16,6155.47,7503.03,8668.36,9613.60]
+		
+		for filter,error in zip(filter_list,error_list):
+			try:
+				mag=photometry._get_value(0,filter)
+				mags.append(mag)
+				
+				errors.append((photometry._get_value(0,error)*u.ABmag))
+				wavelength=(filter_ref[filter_list.index(filter)])
+				wavelengths.append(wavelength)
+					
+				zp=10**((5*np.log10(filter_ref[filter_list.index(filter)])+2.406)/-2.5)
+				fl=zeropoint(mag=mag,zp=zp,wl=wavelength)
+				flux.append(fl)
+					
+				filter_arr.append(filter)
+			except:
+				pass
+	wavelength_arr.append(wavelengths)
+	flux_arr.append(flux)
+	error_arr.append(errors)
+	mag_arr.append(mags)
+
+	wavelengths=[]
+	flux=[]
+	errors=[]
+	mags=[]
+	photometry=SkyMapperGetPhotometryCoords(ra_skymapper,dec_skymapper,radius=radius)
+	if isinstance(photometry,pd.DataFrame):
+		filter_list=['g_psf','r_psf','i_psf','z_psf','u_psf','v_psf']
+		error_list=['e_g_psf','e_r_psf','e_i_psf','e_z_psf','e_u_psf','e_v_psf']		
+
+		#g,r,i,z,u,v
+		filter_ref=[5016.05,6076.85,6076.85,9120.25,3500.22,3878.68]
+		
+		for filter,error in zip(filter_list,error_list):
+			try:
+				mag=photometry._get_value(0,filter)
+				mags.append(mag)
+				
+				errors.append((photometry._get_value(0,error)*u.ABmag))
+				wavelength=(filter_ref[filter_list.index(filter)])
+				wavelengths.append(wavelength)
+					
+				zp=10**((5*np.log10(filter_ref[filter_list.index(filter)])+2.406)/-2.5)
+				fl=zeropoint(mag=mag,zp=zp,wl=wavelength)
+				flux.append(fl)
+					
+				filter_arr.append(filter)
+			except:
+				pass
+	wavelength_arr.append(wavelengths)
+	flux_arr.append(flux)
+	error_arr.append(errors)
+	mag_arr.append(mags)
+
+	wavelengths=[]
+	flux=[]
+	errors=[]
+	mags=[]
+	photometry=get_phot_sdss(ra_sdss,dec_sdss,radius=radius)
+	if isinstance(photometry,pd.DataFrame):
+		filter_list=['uPmag','gPmag','rPmag','iPmag','zPmag']
+		error_list=['e_uPmag','e_gPmag','e_rPmag','e_iPmag','e_zPmag']		
+	
+		#u,g,r,i,z
+		filter_ref=[3608.04,4671.78,6141.12,7457.89,8922.78]
+		
+		for filter,error in zip(filter_list,error_list):
+			try:
+				mag=photometry._get_value(0,filter)
+				mags.append(mag)
+				
+				errors.append((photometry._get_value(0,error)*u.ABmag))
+				wavelength=(filter_ref[filter_list.index(filter)])
+				wavelengths.append(wavelength)
+					
+				zp=10**((5*np.log10(filter_ref[filter_list.index(filter)])+2.406)/-2.5)
+				fl=zeropoint(mag=mag,zp=zp,wl=wavelength)
+				flux.append(fl)
+					
+				filter_arr.append(filter)
+			except:
+				pass
+	wavelength_arr.append(wavelengths)
+	flux_arr.append(flux)
+	error_arr.append(errors)
+	mag_arr.append(mags)
+
+	wavelengths=[]
+	flux=[]
+	errors=[]
+	mags=[]
+	photometry=get_phot_wise(ra_wise,dec_wise,radius=radius)
+	if isinstance(photometry,pd.DataFrame):
+		filter_list=['W1mag','W2mag','W3mag','W4mag']
+		error_list=['e_W1mag','e_W2mag','e_W3mag','e_W4mag']
+
+		#w1,w2,w3,w4
+		filter_ref=[33526.00,46028.00,115608.00,220883.00]
+		
+		for filter,error in zip(filter_list,error_list):
+			try:
+				mag=photometry._get_value(0,filter)
+				mags.append(mag)
+				
+				errors.append((photometry._get_value(0,error)*u.ABmag))
+				wavelength=(filter_ref[filter_list.index(filter)])
+				wavelengths.append(wavelength)
+					
+				zp=10**((5*np.log10(filter_ref[filter_list.index(filter)])+2.406)/-2.5)
+				fl=zeropoint(mag=mag,zp=zp,wl=wavelength)
+				flux.append(fl)
+					
+				filter_arr.append(filter)
+			except:
+				pass
+	wavelength_arr.append(wavelengths)
+	flux_arr.append(flux)
+	error_arr.append(errors)
+	mag_arr.append(mags)
+
+	wavelengths=[]
+	flux=[]
+	errors=[]
+	mags=[]
+	photometry=get_phot_twomass(ra_twomass,dec_twomass,radius=radius)
+	if isinstance(photometry,pd.DataFrame):
+		filter_list=['Jmag','Hmag','Kmag']
+		error_list=['e_Jmag','e_Hmag','e_Kmag']
+		
+		#j,h,k
+		filter_ref=[12350.00,16620.00,21590.00]
+
+		for filter,error in zip(filter_list,error_list):
+			try:
+				mag=photometry._get_value(0,filter)
+				mags.append(mag)
+				
+				errors.append((photometry._get_value(0,error)*u.ABmag))
+				wavelength=(filter_ref[filter_list.index(filter)])
+				wavelengths.append(wavelength)
+					
+				zp=10**((5*np.log10(filter_ref[filter_list.index(filter)])+2.406)/-2.5)
+				fl=zeropoint(mag=mag,zp=zp,wl=wavelength)
+				flux.append(fl)
+					
+				filter_arr.append(filter)
+			except:
+				pass
+	wavelength_arr.append(wavelengths)
+	flux_arr.append(flux)
+	error_arr.append(errors)
+	mag_arr.append(mags)
+
+	for i in range(0,len(flux_arr)):
+		for j in range(0,len(flux_arr[i])):
+			try:
+				flux_arr[i][j]=flux_arr[i][j].value
+			except:
+				pass
+	for i in range(0,len(error_arr)):
+		for j in range(0,len(error_arr[i])):
+			error_arr[i][j]=error_arr[i][j].value
+
+	rel_error_arr=[]
+	for i in range(0,len(mag_arr)):	
+		flux_arr[i]=np.asarray(flux_arr[i])
+		error_arr[i]=np.asarray(error_arr[i])
+		mag_arr[i]=np.asarray(mag_arr[i])		
+
+		rel_error_arr.append(flux_arr[i]*error_arr[i]/mag_arr[i])
+
+	if len(flux)==0 and len(wavelengths)==0:
+		return None
+	
+	#Order on plot: GALEX, SDSS, SkyMapper, PanSTARRS, Gaia, 2MASS, WISE
+	#Order in script: Gaia, GALEX, PanSTARRS, SkyMapper, SDSS, WISE, 2MASS
+	colour_arr=['springgreen','royalblue','gold','aquamarine','deepskyblue','orangered','orange']
+	
+	plot=figure(width=400,height=400,title='SED',x_axis_label=r'\[\lambda_{\text{eff}}\text{ }[\text{AA}]\]',y_axis_label=r'\[\text{flux [mJy]}\]',x_axis_type='log',y_axis_type="log")	
+
+	# Bokeh only understands np nan values
+	for i in range(0,len(wavelength_arr)):
+		wavelength_arr[i]=[np.nan if math.isnan(x) else x for x in wavelength_arr[i]]
+		flux_arr[i]=[np.nan if math.isnan(x) else x for x in flux_arr[i]]
+		rel_error_arr[i]=[np.nan if math.isnan(x) else x for x in rel_error_arr[i]]
+
+	sources=[]
+	for i in range(0,len(wavelength_arr)):
+		upper=[x+e for x,e in zip(flux_arr[i],rel_error_arr[i])]
+		lower=[x-e for x,e in zip(flux_arr[i],rel_error_arr[i])]		
+
+		source=ColumnDataSource(data=dict(wavelength=wavelength_arr[i],flux=flux_arr[i],upper=upper,lower=lower))
+		sources.append(source)
+
+	surveys=['Gaia','GALEX','PanSTARRS','SkyMapper','SDSS','WISE','2MASS']
+
+	for i in range(0,len(wavelength_arr)):
+		if len(flux_arr[i])>0:
+			plot.circle(x='wavelength',y='flux',source=sources[i],color=colour_arr[i],legend_label=f'{surveys[i]}')
+			errors=Whisker(source=sources[i],base='wavelength',upper='upper',lower='lower',line_width=0.5,line_color=colour_arr[i])
+			errors.upper_head.line_width,errors.lower_head.line_width=0.5,0.5
+			errors.upper_head.size,errors.lower_head.size=3,3
+			errors.upper_head.line_color,errors.lower_head.line_color=colour_arr[i],colour_arr[i]
+
+	return plot
+
+if __name__=='__main__':
+	img=get_plot(source=6050296829033196032)
+	show(img)

AstroToolkit/Miscellaneous/ProperMotionCorrection.py

@@ -0,0 +1,50 @@
+import math
+
+def PMCorrection(input,target,inputObjRA,inputObjDEC,pmRA,pmDEC,radius=None):
+	#Filter inputs
+	if not isinstance(input[0],int):
+		raise Exception('input year must be an integer')
+	if not isinstance(input[1],int):
+		raise Exception('input month must be an integer')
+	if not isinstance(target[0],int):
+		raise Exception('target year must be an integer')
+	if not isinstance(target[1],int):
+		raise Exception('target month must be an integer')
+	
+	if not isinstance(inputObjRA,float):
+		try:
+			inputObjRA=float(inputObjRA)
+		except:
+			raise Exception('object RA must be a float or an integer')
+	if not isinstance(inputObjDEC,float):
+		try:
+			inputObjDEC=float(inputObjDEC)
+		except:
+			raise Exception('object DEC must be a float or an integer')
+	if not isinstance(pmRA,float):
+		try:
+			pmRA=float(pmRA)
+		except:
+			raise Exception('object PM_RA must be a float or an integer')
+	if not isinstance(pmDEC,float):
+		try:
+			pmDEC=float(pmDEC)
+		except:
+			raise Exception('object PM_DEC must be a float or an integer')
+	
+	inputYear,inputMonth=input[0],input[1]
+	targetYear,targetMonth=target[0],target[1]
+	
+	yearDelta=targetYear-inputYear
+	monthDelta=targetMonth-inputMonth
+
+	inputObjRA+=(yearDelta*pmRA/3600000+monthDelta*pmRA/43200000)*1/math.cos(inputObjDEC/360*2*math.pi)
+	inputObjDEC+=yearDelta*pmDEC/3600000+monthDelta*pmDEC/43200000
+	
+	timeDelta=abs(yearDelta+monthDelta/12)
+	
+	if radius!=None:
+		radius=radius+math.sqrt((pmRA/1000)**2+(pmDEC/1000)**2)*timeDelta
+		return inputObjRA,inputObjDEC,radius
+	else:
+		return inputObjRA,inputObjDEC

AstroToolkit/Miscellaneous/ReadFits.py

@@ -0,0 +1,35 @@
+from astropy.table import Table
+
+def get_source_list(file_name):
+    if not file_name.endswith('.fits'):
+        file_name+='.fits'
+
+    try:
+        data=Table.read(file_name).to_pandas()
+    except:
+        raise Exception(f'File: {file_name} not found, or invalid format.')
+    
+    try:
+        sources=data.loc[:,'source_id'].tolist()
+    except:
+        raise Exception('No source_id column found.')
+
+    return sources
+
+def get_pos_list(file_name):
+    if not file_name.endswith('.fits'):
+        file_name+='.fits'
+
+    try:
+        data=Table.read(file_name).to_pandas()
+    except:
+        raise Exception(f'File: {file_name} not found, or invalid format.')
+    
+    try:
+        ra=data.loc[:,'ra'].tolist()
+        dec=data.loc[:,'dec'].tolist()
+    except:
+        raise Exception('ra/dec columns not found.')
+    
+    pos_list=[list(x) for x in zip(ra,dec)]
+    return pos_list

AstroToolkit/Miscellaneous/TestFile.py

@@ -0,0 +1,138 @@
+from WDPlanetsToolkit.Tools import *
+
+import time
+from bokeh.plotting import show
+
+start_time=time.time()
+
+figs=[]
+
+'''
+ AR Sco (6050296829033196032,245.4470133276900,-22.8862182469700) is used for all imaging functions as it returns images in all surveys, HU Leo (587316166180416640,141.1853304445800,08.0308343220600)
+ is used for everything else.
+'''
+ 
+print('testing imaging functions (source):\n')
+fig1=getpanstarrsimage(source=6050296829033196032)
+fig2=getskymapperimage(source=6050296829033196032)
+fig3=getdssimage(source=6050296829033196032)
+print('success \n\n\n')
+
+figs.append(fig1)
+figs.append(fig2)
+figs.append(fig3)
+
+print('testing imaging functions (pos):\n')
+fig4=getpanstarrsimage(pos=[245.4470133276900,-22.8862182469700])
+fig5=getskymapperimage(pos=[245.4470133276900,-22.8862182469700])
+fig6=getdssimage(pos=[245.4470133276900,-22.8862182469700])
+print('success \n\n\n')
+
+figs.append(fig4)
+figs.append(fig5)
+figs.append(fig6)
+
+print('panstarrsquery (source):\n',panstarrsquery(source=587316166180416640),'\n\n\n')
+print('skymapperquery (source):\n',skymapperquery(source=587316166180416640),'\n\n\n')
+print('gaiaquery: (source):\n',gaiaquery(source=587316166180416640),'\n\n\n')
+print('galexquery (source):\n',galexquery(source=587316166180416640),'\n\n\n')
+print('rosatquery (source):\n',rosatquery(source=587316166180416640),'\n\n\n')
+print('ztfquery (source):\n',ztfquery(source=587316166180416640),'\n\n\n')
+print('sdssquery (source):\n',sdssquery(source=587316166180416640),'\n\n\n')
+print('wisequery (source):\n',wisequery(source=587316166180416640),'\n\n\n')
+print('twomassquery (source):\n',twomassquery(source=587316166180416640),'\n\n\n')
+
+print('panstarrsquery (pos):\n',panstarrsquery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('skymapperquery (pos):\n',skymapperquery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('gaiaquery (pos):\n',gaiaquery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('galexquery (pos):\n',galexquery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('rosatquery (pos):\n',rosatquery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('ztfquery (pos):\n',ztfquery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('sdssquery (pos):\n',sdssquery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('wisequery (pos):\n',wisequery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('twomassquery (pos):\n',twomassquery(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+
+print('getpanstarrsphot (source):\n',getpanstarrsphot(source=587316166180416640),'\n\n\n')
+print('getskymapperphot (source):\n',getskymapperphot(source=587316166180416640),'\n\n\n')
+print('getgaiaphot (source):\n',getgaiaphot(source=587316166180416640),'\n\n\n')
+print('getgalexphot (source):\n',getgalexphot(source=587316166180416640),'\n\n\n')
+print('getsdssphot (source):\n',getsdssphot(source=587316166180416640),'\n\n\n')
+print('getwisephot (source):\n',getwisephot(source=587316166180416640),'\n\n\n')
+print('gettwomassphot (source):\n',gettwomassphot(source=587316166180416640),'\n\n\n')
+
+print('getpanstarrsphot (pos):\n',getpanstarrsphot(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('getskymapperphot (pos):\n',getskymapperphot(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('getgaiaphot (pos):\n',getgaiaphot(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('getgalexphot (pos):\n',getgalexphot(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('getsdssphot (pos):\n',getsdssphot(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('getwisephot (pos):\n',getwisephot(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+print('gettwomassphot (pos):\n',gettwomassphot(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+
+print('getbulkphot (source):\n',getbulkphot(source=587316166180416640),'\n\n\n')
+
+print('getbulkphot (pos):\n',getbulkphot(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+
+print('testing ztf lightcurve functions (source):\n')
+fig7=getztflc(source=587316166180416640)
+print('success \n\n\n')
+
+figs.append(fig7)
+
+print('testing lightcurve functions (pos):\n')
+fig8=getztflc(pos=[141.1853304445800,08.0308343220600])
+print('success \n\n\n')
+
+figs.append(fig8)
+
+print('testing SED functions (source):\n')
+fig9=getsed(source=587316166180416640)
+print('success \n\n\n')
+
+figs.append(fig9)
+
+print('testing SED functions (pos):\n')
+fig10=getsed(pos=[141.1853304445800,08.0308343220600])
+print('success \n\n\n')
+
+figs.append(fig10)
+
+print('getgaiacoords (source):\n',getgaiacoords(source=587316166180416640),'\n\n\n')
+print('getgaiasource (pos):\n',getgaiasource(pos=[141.1853304445800,08.0308343220600]),'\n\n\n')
+
+print('testing spectrum functions (source):\n')
+fig11=getsdssspectrum(source=587316166180416640)
+print('success \n\n\n')
+
+figs.append(fig11)
+
+print('testing spectrum functions (pos):\n')
+fig12=getsdssspectrum(pos=[141.1853304445800,08.0308343220600])
+print('success \n\n\n')
+
+figs.append(fig12)
+
+print('testing HR diagram function (source):\n')
+fig13=gethrd(source=587316166180416640)
+print('success \n\n\n')
+
+figs.append(fig13)
+
+print('testing datapage function (source):\n')
+fig14=getdatapage(source=587316166180416640)
+print('success \n\n\n')
+
+figs.append(fig14)
+
+print('testing datapage function (pos):\n')
+fig15=getdatapage(pos=[141.1853304445800,08.0308343220600])
+print('success \n\n\n')
+
+figs.append(fig15)
+
+delta_time=round((time.time()-start_time),2)
+
+newline='\n'
+print(f'{newline}all tests passed in: {delta_time}s')
+
+for figure in figs:
+    show(figure)