#! /usr/bin/env python2.7
import matplotlib
matplotlib.use('Agg')
import denudationRateAnalysis as dra
from matplotlib import pyplot as plt
import numpy as np
portenga_data = dra.read_csv('portenga_filtered_steepness.csv')
dr = [r[3] for r in portenga_data]
dr_sig = [r[4] for r in portenga_data]
ks = dict()
r2 = dict()
ks['0.4'] = [
np.float64(r[10].replace('[', '').replace(']', '')) for r in portenga_data
]
r2['0.4'] = [
np.float64(r[11].replace('[', '').replace(']', '')) for r in portenga_data
]
ks['0.5'] = [
np.float64(r[12].replace('[', '').replace(']', '')) for r in portenga_data
]
r2['0.5'] = [
np.float64(r[13].replace('[', '').replace(']', '')) for r in portenga_data
]
ks['0.6'] = [
#! /usr/bin/env python2.6
import denudationRateAnalysis as dra
import numpy as np
data = dra.read_csv('portengadata.csv')
del (data[0])
ksn_vec, area_vec = dra.calculate_ksn_for_data(data, 1000000, 0.4)
np.savez_compressed('ksn_area_data_0_4.npz',
ksn_vec=ksn_vec,
area_vec=area_vec)