From 138e9c43eb8d19ea3abfe23e00a7be59107be20d Mon Sep 17 00:00:00 2001 From: Martin Diehl Date: Fri, 25 Feb 2011 16:24:23 +0000 Subject: [PATCH] some comments (talking to Ricardo) --- code/mpie_spectral.f90 | 6 +++--- code/prec_single.f90 | 2 +- 2 files changed, 4 insertions(+), 4 deletions(-) diff --git a/code/mpie_spectral.f90 b/code/mpie_spectral.f90 index 6e3017bef..54669ce01 100644 --- a/code/mpie_spectral.f90 +++ b/code/mpie_spectral.f90 @@ -283,7 +283,7 @@ program mpie_spectral xi(1) = real(k_s(1), pReal)/geomdimension(1) if (any(xi /= 0.0_pReal)) then do l = 1,3; do m = 1,3 - xinormdyad(l,m) = xi(l)*xi(m)/sum(xi**2) + xinormdyad(l,m) = xi(l)*xi(m)/sum(xi**2) ! Unit sphere, unit vectors in Fourier space enddo; enddo temp33_Real = math_inv3x3(math_mul3333xx33(c0, xinormdyad)) else @@ -511,7 +511,7 @@ program mpie_spectral temp33_Complex(m,n) = sum(gamma_hat(1,1,1,m,n,:,:) *(workfft(i*2-1,j,k,:,:)& +workfft(i*2 ,j,k,:,:)*img)) enddo; enddo - workfft(i*2-1,j,k,:,:) = real (temp33_Complex) + workfft(i*2-1,j,k,:,:) = real (temp33_Complex) ! change of strain workfft(i*2 ,j,k,:,:) = aimag(temp33_Complex) enddo; enddo; enddo else !use precalculated gamma-operator @@ -520,7 +520,7 @@ program mpie_spectral temp33_Complex(m,n) = sum(gamma_hat(i,j,k, m,n,:,:) *(workfft(i*2-1,j,k,:,:)& + workfft(i*2 ,j,k,:,:)*img)) enddo; enddo - workfft(i*2-1,j,k,:,:) = real (temp33_Complex) + workfft(i*2-1,j,k,:,:) = real (temp33_Complex) ! change of strain workfft(i*2 ,j,k,:,:) = aimag(temp33_Complex) enddo; enddo; enddo endif diff --git a/code/prec_single.f90 b/code/prec_single.f90 index ae86adf11..428fe59a9 100644 --- a/code/prec_single.f90 +++ b/code/prec_single.f90 @@ -9,7 +9,7 @@ integer, parameter :: pReal = selected_real_kind(6,37) ! 6 significant digits, up to 1e+-37 integer, parameter :: pInt = selected_int_kind(9) ! up to +- 1e9 integer, parameter :: pLongInt = 8 ! should be 64bit - real(pReal), parameter :: tol_math_check = 1.0e-6_pReal + real(pReal), parameter :: tol_math_check = 1.0e-5_pReal real(pReal), parameter :: tol_gravityNodePos = 1.0e-36_pReal type :: p_vec