From 74a992da87c5a312669bfffd9e6dff41dbe321ee Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Fri, 12 Jan 2007 19:14:00 +0000 Subject: IDScalarWaveFO: Use double precison when calling random_number Create random numbers explicitly in a double precision variable, since gfortran does not support yet higher precisions for the intrinsic random_number. darcs-hash:20070112191430-dae7b-bf8b0a500cbf4770b588f4c4a2efa0de972934d6.gz --- CarpetExtra/IDScalarWaveFO/src/initialdata.F90 | 13 +++++++++---- 1 file changed, 9 insertions(+), 4 deletions(-) (limited to 'CarpetExtra') diff --git a/CarpetExtra/IDScalarWaveFO/src/initialdata.F90 b/CarpetExtra/IDScalarWaveFO/src/initialdata.F90 index 3530fdfbc..bd8edf30b 100644 --- a/CarpetExtra/IDScalarWaveFO/src/initialdata.F90 +++ b/CarpetExtra/IDScalarWaveFO/src/initialdata.F90 @@ -13,6 +13,7 @@ CCTK_REAL pi parameter (pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068d0) CCTK_REAL omega + double precision rphi, rpsix, rpsiy, rpsiz integer i, j, k if (CCTK_EQUALS(initial_data, "plane")) then @@ -49,10 +50,14 @@ do k=1,cctk_lsh(3) do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) - call random_number (phi(i,j,k)) - call random_number (psix(i,j,k)) - call random_number (psiy(i,j,k)) - call random_number (psiz(i,j,k)) + call random_number (rphi) + call random_number (rpsix) + call random_number (rpsiy) + call random_number (rpsiz) + phi(i,j,k) = rphi + psix(i,j,k) = rpsix + psiy(i,j,k) = rpsiy + psiz(i,j,k) = rpsiz end do end do end do -- cgit v1.2.3