slopes = ImgData(os.path.join(folder,'big_spider_slopes.tif'))
aspects = ImgData(os.path.join(folder,'big_spider_aspects.tif'))
spider1 = ImgData(os.path.join(folder,'ESP_022607_0985_cropped_big_spider.cub'))
# spider2 = ImgData(os.path.join(folder,'ESP_022699_0985_cropped_big_spider.cub'))
# read in required data
dem.read_all()
slopes.read_all()
aspects.read_all()
spider1.read_all()
correct_azimuth(dem, aspects)
# create MarsSpicer object
mspice = MarsSpicer()
mspice.goto_ls_0()
mspice.advance_time_by(24*3600*356)
utc = mspice.utc
# seconds per step
timestep = 600
# how many steps for the time series
nsteps = 500
mspice.goto('inca')
# getting l_s times
times, _ = mspice.time_series('F_aspect',timestep, nsteps, provide_times='l_s')
insol = np.zeros_like(dem.data)
to_append = mspicer.solar_constant * dt.seconds * math.cos(diff_angle)
energies.append(to_append)
mspicer.time += dt
return np.array((ets, energies))
def save_data_to_file(data, fname, start, dt):
with open(fname,'w') as f:
f.write('# Start: {1}, dt: {0} s, Intensities in J/(dt*m**2)\n'.format(dt.seconds,start.isoformat()))
np.savetxt(f, data.T)
# location setup
mspicer = MarsSpicer()
mspicer.set_spoint_by(lat=85,lon=0)
# timing setup
mspicer.goto_ls_0()
mspicer.time -= dt.timedelta(days=10,hours=12)
# dt, the time step
my_dt = dt.timedelta(hours=1)
# saving for restart later
start_time = mspicer.time
out_flat = create_arrays(mspicer, 50, my_dt)
save_data_to_file(out_flat[1], 'insolation_flat.txt', start_time, my_dt)
mspicer.get_tilted_normal(30)
mspicer.rotate_tnormal(180)
mspicer.time = start_time
