[Commits] [svn:einsteintoolkit] IllinoisGRMHD/trunk/src/ (Rev. 7)
zachetie at gmail.com
zachetie at gmail.com
Mon Oct 7 10:24:05 CDT 2013
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
More information about the Commits
mailing list