aboutsummaryrefslogtreecommitdiff
path: root/src/boundary.F
diff options
context:
space:
mode:
authorjthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87>2002-11-17 14:55:50 +0000
committerjthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87>2002-11-17 14:55:50 +0000
commitac099068da91d9528088af849caa7aa2e6e1758a (patch)
treef8837bb58d475eb0492f4a6a604a0d51bba94423 /src/boundary.F
parentd72eb352c6182603bde73125f71c4313f22cea28 (diff)
Modified Files:
make.code.defn Added Files: Bona_Masso_data.F77 blended_boundary.F77 boundary.F77 decode_pars.F77 gauge.F77 initialize.F77 metric.F77 xyz_blended_boundary.F77 Removed Files: Bona_Masso_data.F blended_boundary.F boundary.F decode_pars.F gauge.F initialize.F metric.F xyz_blended_boundary.F Rename files which really are Fortran 77 from foo.F (= Cactus Fortran 90 fixed form) to foo.F77 (= Fortran 77) This means they're now compiled with a Fortran 77 compiler. This should make no difference to the semantics (they were already Fortran 77 code), but makes it easier to compile this thorn on platforms which have a Fortran 77 compiler but no Fortran 90 compiler. Also small bugfix: blended_boundary.F --> .F77 xyz_blended_boundary.F --> .F77 * change declaration of ierr from CCTK_REAL --> integer git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/Exact/trunk@135 e296648e-0e4f-0410-bd07-d597d9acff87
Diffstat (limited to 'src/boundary.F')
-rw-r--r--src/boundary.F125
1 files changed, 0 insertions, 125 deletions
diff --git a/src/boundary.F b/src/boundary.F
deleted file mode 100644
index 1862e43..0000000
--- a/src/boundary.F
+++ /dev/null
@@ -1,125 +0,0 @@
-C $Header$
-
-#include "cctk.h"
-#include "cctk_Parameters.h"
-#include "cctk_Arguments.h"
-
- subroutine Exact__boundary(CCTK_ARGUMENTS)
-
- implicit none
-
- DECLARE_CCTK_ARGUMENTS
-
- integer i,j,k
- integer nx,ny,nz
-
- CCTK_REAL tplusone
- CCTK_REAL
- $ dxgxxjunk, dxgyyjunk, dxgzzjunk,
- $ dxgxyjunk, dxgyzjunk, dxgxzjunk,
- $ dygxxjunk, dygyyjunk, dygzzjunk,
- $ dygxyjunk, dygyzjunk, dygxzjunk,
- $ dzgxxjunk, dzgyyjunk, dzgzzjunk,
- $ dzgxyjunk, dzgyzjunk, dzgxzjunk,
- $ axjunk, ayjunk, azjunk,
- $ bxxjunk, bxyjunk, bxzjunk,
- $ byxjunk, byyjunk, byzjunk,
- $ bzxjunk, bzyjunk, bzzjunk
-
-C Grid parameters.
-
- nx = cctk_lsh(1)
- ny = cctk_lsh(2)
- nz = cctk_lsh(3)
-
-C Set all initial data including dijk and vi on all points which
-C are on the boundary of the domain if it really is the boundary
-C of the complete grid. Treat all six sides of the grid cube this way.
-
-c Set t = time + dt. This is necessary here because by the time
-c we reach this point the geometry has been evolved one time step
-c but the variable `time' still hasn't been updated.
-
- tplusone = cctk_time + cctk_delta_time
-
-C Note we also always set the lapse and shift at the boundaries at
-C time t+1. This is to provide boundary conditions for testing
-C elliptic gauge conditions. If they are not used, they will be
-C overwritten by Exact__gauge.
-
-#define EXACTDATAPOINT \
- call Exact__Bona_Masso_data( \
- decoded_exact_model, \
- x(i,j,k), y(i,j,k), z(i,j,k), tplusone, \
- gxx(i,j,k), gyy(i,j,k), gzz(i,j,k), \
- gxy(i,j,k), gyz(i,j,k), gxz(i,j,k), \
- kxx(i,j,k), kyy(i,j,k), kzz(i,j,k), \
- kxy(i,j,k), kyz(i,j,k), kxz(i,j,k), \
- dxgxxjunk, dxgyyjunk, dxgzzjunk, \
- dxgxyjunk, dxgyzjunk, dxgxzjunk, \
- dygxxjunk, dygyyjunk, dygzzjunk, \
- dygxyjunk, dygyzjunk, dygxzjunk, \
- dzgxxjunk, dzgyyjunk, dzgzzjunk, \
- dzgxyjunk, dzgyzjunk, dzgxzjunk, \
- alp(i,j,k), axjunk, ayjunk, azjunk, \
- betax(i,j,k), betay(i,j,k), betaz(i,j,k), \
- bxxjunk, bxyjunk, bxzjunk, \
- byxjunk, byyjunk, byzjunk, \
- bzxjunk, bzyjunk, bzzjunk)
-
- if (cctk_bbox(1) .eq. 1) then
- i=1
- do j=1,ny
- do k=1,nz
- EXACTDATAPOINT
- end do
- end do
- end if
-
- if (cctk_bbox(2) .eq. 1) then
- i=nx
- do j=1,ny
- do k=1,nz
- EXACTDATAPOINT
- end do
- end do
- end if
-
- if (cctk_bbox(3) .eq. 1) then
- j=1
- do i=1,nx
- do k=1,nz
- EXACTDATAPOINT
- end do
- end do
- end if
-
- if (cctk_bbox(4) .eq. 1) then
- j=ny
- do i=1,nx
- do k=1,nz
- EXACTDATAPOINT
- end do
- end do
- end if
-
- if (cctk_bbox(5) .eq. 1) then
- k=1
- do j=1,ny
- do i=1,nx
- EXACTDATAPOINT
- end do
- end do
- end if
-
- if (cctk_bbox(6) .eq. 1) then
- k=nz
- do j=1,ny
- do i=1,nx
- EXACTDATAPOINT
- end do
- end do
- end if
-
- return
- end