Python Model.init_beta_stats示例

示例#1

    def eval(self, values, x):
        values[0] = self.Tmin + self.St * sqrt((x[0] - x0)**2 + (x[1] - y0)**2)


T_s = SurfaceTemperature(element=model.Q.ufl_element())

model.set_geometry(S, B, deform=False)
model.initialize_variables()

model.init_adot(adot)
model.init_T_surface(T_s)
model.init_H(thklim)
model.init_H_bounds(thklim, 1e4)
model.init_q_geo(model.ghf)
model.init_beta_stats('Ubar')

model.save_pvd(model.S, 'S')
model.save_pvd(model.B, 'B')
model.save_pvd(model.adot, 'adot')
model.save_pvd(model.T_surface, 'T_surface')

model.eps_reg = 1e-10

T = HybridTransientSolver(model, config)
T.solve()

File(out_dir + 'Ts.xml') << model.Ts
File(out_dir + 'Tb.xml') << model.Tb
File(out_dir + 'Mb.xml') << model.Mb
File(out_dir + 'H.xml') << model.H

示例#2

    def eval(self, values, x):
        values[0] = self.Tmin + self.St * sqrt((x[0] - x0)**2 + (x[1] - y0)**2)


T_s = SurfaceTemperature(element=model.Q.ufl_element())

model.set_geometry(S, B, deform=False)
model.initialize_variables()

model.init_adot(adot)
model.init_T_surface(T_s)
model.init_H(thklim)
model.init_H_bounds(thklim, 1e4)
model.init_q_geo(model.ghf)
model.init_beta_stats('U')

model.save_pvd(model.S, 'S')
model.save_pvd(model.B, 'B')
model.save_pvd(model.adot, 'adot')
model.save_pvd(model.T_surface, 'T_surface')

model.eps_reg = 1e-10

T = HybridTransientSolver(model, config)
T.solve()

File(out_dir + 'Ts.xml') << model.Ts
File(out_dir + 'Tb.xml') << model.Tb
File(out_dir + 'Mb.xml') << model.Mb
File(out_dir + 'H.xml') << model.H