Commit 8c155b73 authored by Jorn Bruggeman's avatar Jorn Bruggeman
Browse files

propagate missing_value of conserved quantities

parent bb892252
...@@ -2258,6 +2258,7 @@ contains ...@@ -2258,6 +2258,7 @@ contains
consvar%target_hz => get_variable_by_standard_variable(self, standard_variable%at_interfaces()) consvar%target_hz => get_variable_by_standard_variable(self, standard_variable%at_interfaces())
_ASSERT_(associated(consvar%target_hz), 'classify_variables', 'Conserved quantity ' // trim(standard_variable%name) // ' not found at interfaces.') _ASSERT_(associated(consvar%target_hz), 'classify_variables', 'Conserved quantity ' // trim(standard_variable%name) // ' not found at interfaces.')
end select end select
consvar%missing_value = consvar%target%missing_value
standard_variable_node => standard_variable_node%next standard_variable_node => standard_variable_node%next
end do end do
call standard_variable_set%finalize() call standard_variable_set%finalize()
......
...@@ -1257,7 +1257,7 @@ contains ...@@ -1257,7 +1257,7 @@ contains
#endif #endif
call model%get_interior_conserved_quantities(_PREARG_INTERIOR_IN_ total_int _INTERIOR_SLICE_RANGE_PLUS_1_) call model%get_interior_conserved_quantities(_PREARG_INTERIOR_IN_ total_int _INTERIOR_SLICE_RANGE_PLUS_1_)
do ivar = 1, size(model%conserved_quantities) do ivar = 1, size(model%conserved_quantities)
call check_interior_slice_plus_1(total_int _INTERIOR_SLICE_RANGE_PLUS_1_, ivar, -2.e20_rke, & call check_interior_slice_plus_1(total_int _INTERIOR_SLICE_RANGE_PLUS_1_, ivar, model%conserved_quantities(ivar)%missing_value, &
(interior_state_offset + 0.5_rke * (test_model%nstate + 1)) * test_model%nstate _POSTARG_INTERIOR_IN_) (interior_state_offset + 0.5_rke * (test_model%nstate + 1)) * test_model%nstate _POSTARG_INTERIOR_IN_)
end do end do
_END_OUTER_INTERIOR_LOOP_ _END_OUTER_INTERIOR_LOOP_
...@@ -1275,7 +1275,7 @@ contains ...@@ -1275,7 +1275,7 @@ contains
#endif #endif
call model%get_horizontal_conserved_quantities(_PREARG_HORIZONTAL_IN_ total_hz _HORIZONTAL_SLICE_RANGE_PLUS_1_) call model%get_horizontal_conserved_quantities(_PREARG_HORIZONTAL_IN_ total_hz _HORIZONTAL_SLICE_RANGE_PLUS_1_)
do ivar = 1, size(model%conserved_quantities) do ivar = 1, size(model%conserved_quantities)
call check_horizontal_slice_plus_1(total_hz _HORIZONTAL_SLICE_RANGE_PLUS_1_, ivar, -2.e20_rke, & call check_horizontal_slice_plus_1(total_hz _HORIZONTAL_SLICE_RANGE_PLUS_1_, ivar, model%conserved_quantities(ivar)%missing_value, &
(surface_state_offset + 0.5_rke * (test_model%nsurface_state + 1)) * test_model%nsurface_state + & (surface_state_offset + 0.5_rke * (test_model%nsurface_state + 1)) * test_model%nsurface_state + &
(bottom_state_offset + 0.5_rke * (test_model%nbottom_state + 1)) * test_model%nbottom_state _POSTARG_HORIZONTAL_IN_) (bottom_state_offset + 0.5_rke * (test_model%nbottom_state + 1)) * test_model%nbottom_state _POSTARG_HORIZONTAL_IN_)
end do end do
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment