aselimov/CACmb
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Working redef box option

Browse Source
This commit is contained in:
Alex Selimov 2020-10-13 11:57:20 -04:00
parent b658202a1e
commit c42db27f57
2 changed files with 5 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
src/call_option.f90
View File

@ -5,6 +5,7 @@ subroutine call_option(option, arg_pos)
use opt_orient use opt_orient
use opt_deform use opt_deform
use opt_delete use opt_delete
use opt_redef_box
use box use box
implicit none implicit none
@ -38,6 +39,8 @@ subroutine call_option(option, arg_pos)
call run_delete(arg_pos) call run_delete(arg_pos)
case('-set_cac') case('-set_cac')
arg_pos=arg_pos +3 arg_pos=arg_pos +3
case('-redef_box')
call redef_box(arg_pos)
case default case default
print *, 'Option ', trim(adjustl(option)), ' is not currently accepted.' print *, 'Option ', trim(adjustl(option)), ' is not currently accepted.'
stop 3 stop 3

4
src/opt_redef_box.f90
View File

@ -35,9 +35,9 @@ module opt_redef_box
delete_list(:) = 0 delete_list(:) = 0
do i = 1, ele_num do i = 1, ele_num
!Determine if all nodes are within the new boundaries !Determine if all nodes are within the new boundaries
node_out(:) = .false.
do inod = 1, ng_node(lat_ele(i)) do inod = 1, ng_node(lat_ele(i))
do ibasis = 1, basisnum(lat_ele(i)) do ibasis = 1, basisnum(lat_ele(i))
node_out(:) = .false.
if(.not.in_block_bd(r_node(:,ibasis,inod,i), new_bd)) then if(.not.in_block_bd(r_node(:,ibasis,inod,i), new_bd)) then
node_out(inod) = .true. node_out(inod) = .true.
end if end if
@ -66,7 +66,7 @@ module opt_redef_box
end if end if
end do end do
call delete_elements(delete_num, delete_list) call delete_elements(delete_num, delete_list(1:delete_num))
box_bd=new_bd box_bd=new_bd
box_bc = new_bc box_bc = new_bc