User: zetienne Date: 2013/10/07 10:24 AM
Removed: /trunk/src/ driver_bcs.F90
Log: Remove dead code.
File Changes:
Directory: /trunk/src/ ======================
File [removed]: driver_bcs.F90 Delta lines: +0 -74 =================================================================== --- trunk/src/driver_bcs.F90 2013-10-01 23:46:41 UTC (rev 6) +++ trunk/src/driver_bcs.F90 2013-10-07 15:24:04 UTC (rev 7) @@ -1,74 +0,0 @@ -!--------------------------- -! Boundary conditions driver -!--------------------------- - -#include "cctk.h" -#include "cctk_Arguments.h" -#include "cctk_Functions.h" -#include "cctk_Parameters.h" - -subroutine IllinoisGRMHD_driver_bcs(CCTK_ARGUMENTS) - - implicit none - - DECLARE_CCTK_ARGUMENTS - DECLARE_CCTK_PARAMETERS - DECLARE_CCTK_FUNCTIONS - - integer, dimension(3) :: ext,fake_ext,thisproc_have_global_bdry_min,thisproc_have_global_bdry_max,syms - real*8 :: dx,dy,dz - integer :: i - integer :: ierr - integer :: levelnumber - - ext = cctk_lsh - - levelnumber = int( log(real(cctk_levfac(1)))/log(2.0) + 1.e-5 ) + 1 - - dx = CCTK_DELTA_SPACE(1) - dy = CCTK_DELTA_SPACE(2) - dz = CCTK_DELTA_SPACE(3) - - ! FIRST SET SYMMETRY GHOSTZONES ON ALL EVOLVED VARIABLES! - call CartSymGN(ierr,cctkGH,'IllinoisGRMHD::grmhd_conservatives') - syms = (/ 1,0, Sym_Bz /) - call IllinoisGRMHD_set_symmetry_gzs_staggered(cctkGH,cctk_lsh,dx,dy,dz,Symmetry, Ax, syms, 0,1,1) - syms = (/ 0,1, Sym_Bz /) - call IllinoisGRMHD_set_symmetry_gzs_staggered(cctkGH,cctk_lsh,dx,dy,dz,Symmetry, Ay, syms, 1,0,1) - syms = (/ 0,0,-Sym_Bz /) - call IllinoisGRMHD_set_symmetry_gzs_staggered(cctkGH,cctk_lsh,dx,dy,dz,Symmetry, Az, syms, 1,1,0) - - if(levelnumber==1 .and. cctk_iteration.gt.0) then - - ! This is a quick and dirty fix: set thisproc_have_global_bdry_min/max so that - ! ALL boundaries are updated, EXCEPT in symmetry ghostzones! - ! We must apply this fix, since have_global_bdry_min/max does not save all the - ! necessary information when being saved to checkpoint. - ! In other words, have_global_bdry_min/max is BROKEN. - thisproc_have_global_bdry_min=1 - if(Symmetry==1) thisproc_have_global_bdry_min(3)=0 - thisproc_have_global_bdry_max=1 - - !---------------------------------------- - ! Impose frozen BC on the "conservative" variables - ! before primitive recovery - if (Matter_BC==2) then - do i=1,cctk_nghostzones(1) - fake_ext = cctk_lsh - cctk_nghostzones + i - - call apply_matter_frozen_bc(ext,fake_ext, & - rho_star,tau,mhd_st_x,mhd_st_y,mhd_st_z, & - Ax,Ay,Az, & - rho_star_p,tau_p,mhd_st_x_p,mhd_st_y_p,mhd_st_z_p, & - Ax_p,Ay_p,Az_p,thisproc_have_global_bdry_min, & - thisproc_have_global_bdry_max) - end do - else - write(*,*) "ERROR: Matter_BC==",Matter_BC,". Only Matter_BC==2 is supported!" - stop - end if - !---------------------------------------- - - end if - -end subroutine IllinoisGRMHD_driver_bcs