diff options
| author | schnetter <schnetter@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2008-08-25 20:07:03 +0000 |
|---|---|---|
| committer | schnetter <schnetter@a678b1cf-93e1-4b43-a69d-d43939e66649> | 2008-08-25 20:07:03 +0000 |
| commit | 35e170b820fe9aac5158d3e0aa3672433f2be9a1 (patch) | |
| tree | 392ca2db81c03e89586b4c1666a58e85503f8bac /src/brillq.F | |
| parent | 9eef156a017128bf2800c455e2238b46130cb235 (diff) | |
Don't use explicit double precision Fortran functions; use type
generic functions instead
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDBrillData/trunk@110 a678b1cf-93e1-4b43-a69d-d43939e66649
Diffstat (limited to 'src/brillq.F')
| -rw-r--r-- | src/brillq.F | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/src/brillq.F b/src/brillq.F index c3e8d6d..98e3ddd 100644 --- a/src/brillq.F +++ b/src/brillq.F @@ -91,11 +91,11 @@ c "exp" brill data if (qtype.eq.0) then - brillq = exp_a*dabs(rho)**(2 + exp_b) - . / dexp((rho - exp_rho0)**2)/(rho**2 + z**2) + brillq = exp_a*abs(rho)**(2 + exp_b) + . / exp((rho - exp_rho0)**2)/(rho**2 + z**2) if (exp_sigmaz.ne.0.0D0) then - brillq = brillq/dexp((z/exp_sigmaz)**2) + brillq = brillq/exp((z/exp_sigmaz)**2) end if else if (qtype.eq.1) then @@ -114,7 +114,7 @@ c smooth "pure quadrupole" q should be, plus the choice of c Holz et al. brillq = gundlach_a*(rho/gundlach_sigmarho)**gundlach_b - & / dexp(((rho**2 + z**2 - gundlach_r0**2)/ + & / exp(((rho**2 + z**2 - gundlach_r0**2)/ & gundlach_sigmar**2)**(gundlach_c/2)) else |