Esempio n. 1
0
    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)
Esempio n. 2
0
    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')