diff options
author | jthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-11-17 14:55:50 +0000 |
---|---|---|
committer | jthorn <jthorn@e296648e-0e4f-0410-bd07-d597d9acff87> | 2002-11-17 14:55:50 +0000 |
commit | ac099068da91d9528088af849caa7aa2e6e1758a (patch) | |
tree | f8837bb58d475eb0492f4a6a604a0d51bba94423 /src/boundary.F | |
parent | d72eb352c6182603bde73125f71c4313f22cea28 (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.F | 125 |
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 |