''' COTD WOOOOOO!!!!!!! !!!!!!600+ UPVOTES!!!! THANK YOU SO MUCH!!!!!!!! vvv Check Below for Helpful Hints! vvv Scroll to the bottom to Explore the Multiverse! Thanks for viewing my code! It took many hours to gather data on the specifics of the universe to make this simulation as realistic as possible, but I loved every minute. I hope you enjoy playing with the simulator! Hit run to see how it works and then scroll to the bottom to explore the multiverse! Let me know if you want me to add anything or have any concerns/bugs. More to come! -nl ''' ''' PLEASE ENTER YOUR SYSTEM INTO OUR DATABASE HERE!!! https://goo.gl/forms/Mxfz4D4afuVQcOjR2 I WILL FEATURE YOUR NAME AND SYSTEM FOR A DAY WHEN THE CODE IS RUN! VIEW CURRENT DATABASE HERE: https://docs.google.com/spreadsheets/d/1CMLx9CX9pWvqZYBzXCieHIiTZ8OprwJ8kz3aC8P7tlM ''' ''' Helpful Hints: @ == Galaxy + == Empty Sector Dimension names can be ANYTHING! Don't restrict yourself to numbers! Each Universe Dimension is infinite. Travel across the Sectors of the dimension to expand your universal view! Each Galaxy is restrained by size. Look to see the Radius of the Galaxy, it will tell you how many Sub-Sectors the Galaxy reaches out in positive and negative coordinate directions! Galaxies are Circular, you will be more likely find stars towards the center! Stars are measured based on their temperature from hottest to coolest (O,B,A,F,G,K,M) and are plotted on the galaxy sub-sector in this way. 1 AU is the average distance from the Sun to Earth. Earth's gravity is about 9.8m/s2. To find your weight on a different planet, do this>> (Planet Gravity/9.8)*(Your Weight) In output, the Mass, Radius and Luminosity of stars are measured compared to our Sun (M,R,L) respectively. The Mass and Radius of planets are measured compared to Earth (M+,R+) respectively. Sectors cover a distance of 3,250,000ly across and down Sub-Sectors cover a distance of 65ly across and down ''' import sys, codecs import random as rand #sys.stdout = codecs.getwriter('utf_16')(sys.stdout.buffer, 'strict') seed = '' def universe(draw=0,dim=0,X=0,Y=0): global seed seed = str(dim)+'/'+str(X)+','+str(Y) rand.seed(seed) spread = rand.randint(1,50) if draw == 1: sector = create('U',spread) num_gal = 0 for i in sector: for j in i: if j != '+': num_gal += 1 print('Dimension: {}'.format(dim)) print('Sector: {:,} , {:,}'.format(X,Y)) print('Spread: {}%'.format(spread)) print('Number of Galaxies in Sector: {}'.format(num_gal)) print() print(' 0123456789abc') x=0 for i in sector: if x == 10: i = 'a '+i elif x == 11: i = 'b '+i elif x == 12: i = 'c '+i else: i = str(x)+' '+i print(i) x+=1 print('--------------------') else: return create('U',spread) def create(ctype,spread,gx=0,gy=0): global seed sd = seed if ctype == 'G': seed = '/'.join(seed.split('/')[:3]) rand.seed(seed) stars = 100000000 low = .00033 h = low n = .0033 m = rand.uniform(low,1-(h+n)-low*6) k = rand.uniform(low,1-(n+h+m)-low*5) g = rand.uniform(low,1-(n+h+k+m)-low*4) f = rand.uniform(low,1-(n+h+g+k+m)-low*3) a = rand.uniform(low,1-(n+h+f+g+k+m)-low*2) b = rand.uniform(low,1-(n+h+a+f+g+k+m)-low) o = 1-(n+h+b+a+f+g+k+m) star_stats = {m:'M',k:'K',g:'G',f:'F',a:'A',b:'B',o:'O',h:'H',n:'N'} types = [m,k,g,f,a,b,o,h,n] m*=stars k*=stars g*=stars f*=stars a*=stars b*=stars o*=stars h*=stars n*=stars sector = [] for i in range(13): line = '' for j in range(13): check = rand.randint(1,100) if check <= spread: if ctype == 'U': line = line + '@' elif ctype == 'G': star = rand.randint(1,stars) if gx == 0 and gy == 0 and i==6 and j==6: line = line + 'H' else: if star<=m: line = line + 'M' elif m 10000000: highrad -= 1 else: super = 0 rad = round(rand.uniform(30,highrad),2) lum = round((rad**2)*(temp/5776)**4,2) else: rad = round(rand.uniform(lowr,highr),4) lum = (rad**2)*(temp/5776)**4 if lum<1: lum = round(lum,5) else: lum = round(lum,2) elif s_type == 'H': mass = round(rand.uniform(lowm,highm),2) rad = 2.95*mass if mass < 100: classify = 'Stellar' rad = round(rad,2) elif mass < 10**5: classify = 'Intermediate' else: classify = 'Supermassive' p_chance = rand.randint(1,10) if p_chance <= 3: planets = rand.randint(6,8) else: planets = rand.randint(0,5) print('Dimension: {}'.format(dim)) print('Sector: {:,} , {:,}'.format(X,Y)) print('Galaxy: {} , {}'.format(x,y)) print('Sub-Sector: {:,} , {:,}'.format(gx,gy)) print('System: {} , {}'.format(a,b)) user_database = {'902/234,141/0,10/-4584,-1740/12,1':'Erik Rasolov','286/-111/837/7/4/264/268/3/9':'Toodles','90/90,90/3,5/4,4/4,4':'Madhur','305/-271,883/1,11/-5649,-3041/7,12':'Julia M'} try: print('\n[[{}\'s System]]'.format(user_database[seed])) except: pass print() if s_type not in 'HN': print('Stellar Classification: {}{}'.format(s_type,classify)) print('Radius: {:,}R'.format(rad)) elif s_type == 'H': print('Stellar Classification: {} Black Hole'.format(classify)) if rad >= 347850: rad = round(rad/695700,2) lowh,highh = 0,0 print('Radius: {:,}R'.format(rad)) else: print('Radius: {:,}km'.format(rad)) print('Mass: {:,}M'.format(mass)) elif s_type == 'N': print('Radius: {:,}km'.format(rad)) if temp != 0: print('Mass: {:,}M'.format(mass)) print('Surface Temperature: {:,}K'.format(temp)) print('Luminosity: {:,}L'.format(lum)) lowh = round(.95*(lum**.5),3) highh = round(1.37*(lum**.5),3) print('Habitable Zone: {:,}AU--{:,}AU'.format(lowh,highh)) else: lowh,highh = 0,0 print('Planets: {}'.format(planets)) print('--------------------') smass = mass*(1.989*10**30) mass = mass * 333000 highd = 0 if planets != 0: for i in range(1,planets+1): ptype = rand.randint(1,8) if ptype <= 4: ptype = 'Terrestrial' elif 4 < ptype <= 6: ptype = 'Ice Giant' else: ptype = 'Gas Giant' if ptype == 'Terrestrial': pmass = rand.randint(1,4) if pmass <= 3: pmass = round(rand.uniform(.01,4),2) else: pmass = round(rand.uniform(4,16),2) den = round(rand.uniform(3,12),2) sat = rand.randint(0,3) elif ptype == 'Ice Giant': pmass = round(rand.uniform(10,50),2) den = round(rand.uniform(1,3),2) sat = rand.randint(0,7) elif ptype == 'Gas Giant': pmass = round(rand.uniform(50,10000)) den = round(rand.uniform(.01,2),2) sat = rand.randint(0,70) vol = round((pmass*(5.972*10**27)/den)/10**15,2) prad = round((((vol*3)/12.56)**(1/3))/6371,2) roche = (prad*6371*(2*(mass/pmass))**(1/3))/149597871 if i == 0: far = 1.5 else: far = 5 lowd = round(rand.uniform(roche+highd,roche+highd+far),3) highd = rand.randint(1,4) if highd <= 3: long = 1 else: long = 3 highd = round(rand.uniform(lowd,lowd+long),3) vel = (((6.76*10**-11)*smass)/((lowd+highd)*1.496*10**11/2))**.5 vel = round(vel,2) grav = (6.67408*10**-11)*((pmass*5.979*(10**24))/(prad*(6376000**2))) grav = round(grav,2) a = ((1.496*10**11)*(highd+lowd)/2) time = 6.28*((a**3/(smass*(6.674*10**-11))))**.5 per = 'Earth Days' time = round(time/86400,2) rot = rand.randint(1,6) if rot <= 4: rot = rand.uniform(.1,3) elif 4 < rot <= 5: rot = rand.uniform(.1,50) else: rot = rand.uniform(.1,400) rot = round(rot,2) rotper = per if rot>=365: rot = round(rot/365,2) rotper = 'Earth Years' if time>=365: time = round(time/365,2) per = 'Earth Years' print('Planet from Star: {}'.format(i)) print('Planet Type: {}'.format(ptype)) print('Radius: {:,}R+'.format(prad)) print('Mass: {:,}M+'.format(pmass)) print('Volume: {:0.3e}km3'.format(vol)) print('Density: {:,}g/cm3'.format(den)) print('Gravity: {:,}m/s2'.format(grav)) print('Min Distance from Star: {:,}AU'.format(lowd)) print('Max Distance from Star: {:,}AU'.format(highd)) print('Average Orbital Velocity: {:,}m/s'.format(vel)) print('Orbital Period: {:,} {}'.format(time,per)) print('Rotational Period: {:,} {}'.format(rot,rotper)) if lowd >= lowh and highd <= highh: hab = 'Orbits' elif lowd >= lowh and lowd <= highh or highd >= lowh and highd <= highh or lowd <= lowh and highd >= highh: hab = 'Enters' else: hab = 'Outside' print('Enters/Orbits in Habitable Zone: {}'.format(hab)) if ptype in ('Gas Giant','Ice Giant'): ring = rand.randint(1,2) else: ring = rand.randint(1,10) if ring: ring = 'Yes' else: ring = 'No' print('Rings: {}'.format(ring)) print('Number of Satellites: {}'.format(sat)) if i != planets: print() mass = pmass print('--------------------\n') def random_system(): print('--Random System--\n') dim = rand.randint(-1000,1000) X = rand.randint(-1000,1000) Y = rand.randint(-1000,1000) gx,gy= 0,0 x,y = 0,0 a,b = 0,0 y,x = rand.choice(galaxy(x,y,gx,gy,X,Y,dim,0,1)) universe(1,dim,X,Y) size = int(galaxy(x,y,gx,gy,X,Y,dim,0,0,1)*.5) gx,gy = rand.randint(-size,size),rand.randint(-size,size) galaxy(x,y,gx,gy,X,Y,dim,1) a,b = rand.choice(solar(a,b,x,y,gx,gy,X,Y,dim,1)) solar(a,b,x,y,gx,gy,X,Y,dim) '''PLEASE ENTER YOUR SYSTEM INTO OUR DATABASE HERE!!! https://goo.gl/forms/Mxfz4D4afuVQcOjR2 ''' #----Alter These---- #a=10,b=11,c=12 for (sub-)sector coordinates run_type = 0 #0 = Random System, 1 = Configured (here), 2 = Seed (Below) dim = 902 #dimension (Any str or num) X = 234 #Sector X Y = 141 #Sector Y x = 0 #Galaxy X +Numbers >>> y = 10 #Galaxy Y +Numbers vvv gx = -4584 #Sub-Sector X gy = -1740 #Sub-Sector Y a = 12 #Solar System X +Numbers >>> b = 1 #Solar System Y +Numbers vvv #------------------- #X,Y and gx,gy can be positive or negative #x,y and a,b can only be 0-12 #----Leave These---- if run_type == 0: try: random_system() except: print('No Systems Exist: Run Again') elif run_type == 1: universe(1,dim,X,Y) galaxy(x,y,gx,gy,X,Y,dim,1) solar(a,b,x,y,gx,gy,X,Y,dim) elif run_type == 2: #Place your known seed here seed = '305/-271,883/1,11/-5649,-3041/7,12' seed = seed.split('/') dim = int(seed[0]) X,Y = map(int,seed[1].split(',')) x,y = map(int,seed[2].split(',')) gx,gy = map(int,seed[3].split(',')) a,b = map(int,seed[4].split(',')) universe(1,dim,X,Y) galaxy(x,y,gx,gy,X,Y,dim,1) solar(a,b,x,y,gx,gy,X,Y,dim) #------------------- '''PLEASE ENTER YOUR SYSTEM INTO OUR DATABASE HERE!!! https://goo.gl/forms/Mxfz4D4afuVQcOjR2 '''

Procedurally Generated Universe