def main():
b0 = [-53.86014283247338, -96.55148545680116, -209.56378728526175, -329.5285447843742, -276.1078577736731, -99.4682311472077, -74.60838257235906, 103.76931702628877, 199.59759885797695, 325.49387992076936, 294.9387516971606, 190.74077827555544, 95.99117577626988, 252.4851292030383, 410.3736616110993, 500, 500, 444.71357635367707, 326.35184433516275, 223.25791552817049, 98.55471776145947, 134.18167926051117, 30.336550811088184, -88.39136174706994, -121.36948314990242, -18.45822543176338, -66.01158355801323, 52.93628093781889, -92.66509252409524, 15.70587677548626, -22.132177754463765, 6.391541935525751, -68.42339808954196, 39.59640808291448, 53.464027842215515, 136.78502122848863, 254.14520471035019, 192.5665648106059, 65.30439277418417, 14.10079843084735, -114.02424000705798, -289.4634592194706, -303.0305745281302, -266.62019841090347, -198.4732468197231, -210.421471160374, -143.2408448139692, -66.7972463258127, -18.769521382966776, 46.9742033841991, 105.26603075720368, 164.98279637713722, -43.50064538327361, -79.11047511726842, -3.8891318593645394, 12.390579201231496, -59.56957705146683, 0]
base = tools.load_nolan_bedrock()
b0 = map(operator.add, base, b0)
run1 = isothermalISM(58, 1000, 0.0015, .0005, 0.00022, b0) #55 nodes, 1000-meter spacing, basal slip was .0005
run1.openOutput('run1.nc')
for i in range(1500): #5000 years
run1.timestep(1)
if(i%100==0):
print 'on timestep', i
run1.write()
#run1.calculate_velocity()
run1.close()
tools.plot_model_run('run1.nc')
def main():
base = tools.load_nolan_bedrock()
b0 = base#map(operator.add, base, b0)
run1 = isothermalISM(58, 1000, 1e-16, 1e-16, b0)#.001125, .000025, b0) #55 nodes, 1000-meter spacing, basal slip was .0005
run1.openOutput('run1.nc')
for i in range(5000): #5000 years
run1.timestep(31536000)#1 yr in seconds??
if(i%100==0):
print 'on timestep', i
run1.write()
#basal = run1.calculate_basal_velocity()
#deformation = run1.calculate_velocity()
run1.close()
tools.plot_model_run('run1.nc')
print run1.get_ice_thickness()
