def __init__(self,num_nodes,dx,bslip_start, bslip_stop, b): #initializes the model's fields self.dx = dx self.generate_bslip_array(bslip_start, bslip_stop) self.num_nodes = num_nodes + self.nodes_past_divide self.time = 0 self.x= np.array(range(0,((self.num_nodes)*self.dx),self.dx)) self.ice_thickness= np.zeros(self.num_nodes) self.bed_elev = b for i in range(1210, 410, -40): self.bed_elev.append(i) self.surface_elev= self.bed_elev #start with no ice self.mass_balance = tools.load_mbal(index=1)
def __init__(self,num_nodes,dx,bslip_1, bslip_2, bslip_3, b): #initializes the model's fields self.dx = dx self.bslip_1 = bslip_1 self.bslip_2 = bslip_2 self.bslip_3 = bslip_3 self.num_nodes = num_nodes + 20 self.time = 0 self.x= np.array(range(0,((self.num_nodes)*self.dx),self.dx)) self.ice_thickness= np.zeros(self.num_nodes) self.bed_elev = b for i in range(1210, 410, -40): self.bed_elev.append(i) self.surface_elev= self.bed_elev #start with no ice self.mass_balance = tools.load_mbal('reduced_smoothed_mbal.csv')