-
Notifications
You must be signed in to change notification settings - Fork 0
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
[WIP] Create quarter administration #360
Conversation
@@ -78,7 +78,7 @@ subroutine set_2d_computational_nodes_lines(origin, lgrmin, kmax, mmax, nmax, dx | |||
deallocate(area_mask_padded) | |||
endif | |||
write(*,*) '** INFO: Number of 2D nodes is: ', nod - 1 | |||
write(*,*) '** INFO: Number of 2D lines is: ', l_u + l_v | |||
write(*,*) '** INFO: Number of 2D lines is: ', l_u + (l_v - l_u) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
interesting!
if (ku==kd) then | ||
quarter_line(get_quarter_idx(nodd, 2), 1) = l - 1 ! This needs to be an index in Python again. | ||
quarter_line(get_quarter_idx(nodd, 4), 1) = l - 1 | ||
quarter_line(get_quarter_idx(nodu, 1), 1) = l - 1 | ||
quarter_line(get_quarter_idx(nodu, 3), 1) = l - 1 | ||
quarter_neighbour(get_quarter_idx(nodd, 2), 1) = nodu - 1 | ||
quarter_neighbour(get_quarter_idx(nodd, 4), 1) = nodu - 1 | ||
quarter_neighbour(get_quarter_idx(nodu, 1), 1) = nodd - 1 | ||
quarter_neighbour(get_quarter_idx(nodu, 3), 1) = nodd - 1 | ||
elseif (ku == kd + 1) then | ||
if (nd == 2 * nu - 1) then |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Some comments here explaining the logic would be nice I think.
if (ku==kd) then | |
quarter_line(get_quarter_idx(nodd, 2), 1) = l - 1 ! This needs to be an index in Python again. | |
quarter_line(get_quarter_idx(nodd, 4), 1) = l - 1 | |
quarter_line(get_quarter_idx(nodu, 1), 1) = l - 1 | |
quarter_line(get_quarter_idx(nodu, 3), 1) = l - 1 | |
quarter_neighbour(get_quarter_idx(nodd, 2), 1) = nodu - 1 | |
quarter_neighbour(get_quarter_idx(nodd, 4), 1) = nodu - 1 | |
quarter_neighbour(get_quarter_idx(nodu, 1), 1) = nodd - 1 | |
quarter_neighbour(get_quarter_idx(nodu, 3), 1) = nodd - 1 | |
elseif (ku == kd + 1) then | |
if (nd == 2 * nu - 1) then | |
if (ku==kd) then | |
! ? | |
quarter_line(get_quarter_idx(nodd, 2), 1) = l - 1 ! This needs to be an index in Python again. | |
quarter_line(get_quarter_idx(nodd, 4), 1) = l - 1 | |
quarter_line(get_quarter_idx(nodu, 1), 1) = l - 1 | |
quarter_line(get_quarter_idx(nodu, 3), 1) = l - 1 | |
quarter_neighbour(get_quarter_idx(nodd, 2), 1) = nodu - 1 | |
quarter_neighbour(get_quarter_idx(nodd, 4), 1) = nodu - 1 | |
quarter_neighbour(get_quarter_idx(nodu, 1), 1) = nodd - 1 | |
quarter_neighbour(get_quarter_idx(nodu, 3), 1) = nodd - 1 | |
elseif (ku == kd + 1) then | |
if (nd == 2 * nu - 1) then | |
! ? |
threedigrid_builder/application.py
Outdated
@@ -225,6 +228,7 @@ def make_gridadmin( | |||
progress_callback: Optional[Callable[[float, str], None]] = None, | |||
upgrade: bool = False, | |||
convert_to_geopackage: bool = False, | |||
create_triangulation: bool = False, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This is never used?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
2 small comments/questions
No description provided.