Skip to content
Snippets Groups Projects
Commit c94f31a4 authored by David Seus's avatar David Seus
Browse files

debug calc_gli_term

parent 2dc9b5a4
No related branches found
No related tags found
No related merge requests found
[interf_ind]"""Class definitions for subdomains representing soil patches for LDD-TPR simulation """Class definitions for subdomains representing soil patches for LDD-TPR simulation
This file contains class definitions for the class DomainPatch which defines the This file contains class definitions for the class DomainPatch which defines the
subdomains used for the LDD simulation. Each subdomain of this class holds subdomains used for the LDD simulation. Each subdomain of this class holds
...@@ -581,7 +581,7 @@ class DomainPatch(df.SubDomain): ...@@ -581,7 +581,7 @@ class DomainPatch(df.SubDomain):
phi_a = self.testfunction[phase] phi_a = self.testfunction[phase]
gli_pressure_part = df.assemble(-2*dt*Lambda[phase]*pa_prev*phi_a*ds) gli_pressure_part = df.assemble(-2*dt*Lambda[phase]*pa_prev*phi_a*ds)
gli_p_prev_part = gli_pressure_part.get_local() gli_p_prev_part = gli_pressure_part.get_local()
gli = gli_p_prev_part - gli_prev_neighbour gli = gli_p_prev_part - gli_prev_of_neighbour
# check if there are shared dofs with another interface. # check if there are shared dofs with another interface.
if dofs_in_common_with_other_interfaces[phase].size != 0: if dofs_in_common_with_other_interfaces[phase].size != 0:
for common_dof in dofs_in_common_with_other_interfaces[phase]: for common_dof in dofs_in_common_with_other_interfaces[phase]:
......
...@@ -85,10 +85,10 @@ isRichards = { ...@@ -85,10 +85,10 @@ isRichards = {
} }
##################### MESH ##################################################### ##################### MESH #####################################################
mesh_resolution = 40 mesh_resolution = 6
##################### TIME ##################################################### ##################### TIME #####################################################
timestep_size = 4*0.0001 timestep_size = 4*0.0001
number_of_timesteps = 5000 number_of_timesteps = 200
# decide how many timesteps you want analysed. Analysed means, that we write out # decide how many timesteps you want analysed. Analysed means, that we write out
# subsequent errors of the L-iteration within the timestep. # subsequent errors of the L-iteration within the timestep.
number_of_timesteps_to_analyse = 11 number_of_timesteps_to_analyse = 11
...@@ -115,8 +115,8 @@ L = {# ...@@ -115,8 +115,8 @@ L = {#
lambda_param = {# lambda_param = {#
# subdom_num : lambda parameter for the L-scheme # subdom_num : lambda parameter for the L-scheme
1 : {'wetting' :100},# 1 : {'wetting' :4},#
2 : {'wetting' :100} 2 : {'wetting' :4}
} }
## relative permeabilty functions on subdomain 1 ## relative permeabilty functions on subdomain 1
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment