Changes from trunk/src/erf.f90 at r4 to branches/FLEXPART_9.1.3/src/erf.f90 at r13

File:

• branches/FLEXPART_9.1.3/src/erf.f90 (modified) (8 diffs, 1 prop)

branches/FLEXPART_9.1.3/src/erf.f90

@@ -42 +42 @@
 
   integer :: j
-  real :: x,tmp,ser,xx,gammln
-  real :: cof(6) = (/ &
+  real(kind=dp) :: x,tmp,ser,xx,gammln
+  real(KIND=dp) :: cof(6) = (/ &
        76.18009173_dp, -86.50532033_dp, 24.01409822_dp, &
        -1.231739516_dp, .120858003e-2_dp, -.536382e-5_dp    /)
-  real :: stp = 2.50662827465_dp
-  real :: half = 0.5_dp, one = 1.0_dp, fpf = 5.5_dp
+  real(KIND=dp) :: stp = 2.50662827465_dp
+  real(KIND=dp) :: half = 0.5_dp, one = 1.0_dp, fpf = 5.5_dp
 
   x=xx-one
@@ -62 +62 @@
 function gammp(a,x)
 
-  implicit none
-
-  real :: a, x, gln, gamser, gammp, gammcf
+  use par_mod, only: dp
+
+  implicit none
+
+  real(KIND=dp) :: a, x, gln, gamser, gammp, gammcf
 
   if(x .lt. 0. .or. a .le. 0.) then
@@ -81 +83 @@
 function gammq(a,x)
 
-  implicit none
-
-  real :: a, x, gln, gamser, gammq, gammcf
+  use par_mod, only: dp
+
+  implicit none
+
+  real(KIND=dp) :: a, x, gln, gamser, gammq, gammcf
 
   if(x.lt.0..or.a.le.0.) then
@@ -100 +104 @@
 subroutine gser(gamser,a,x,gln)
 
+  use par_mod, only: dp
+
   implicit none
 
   integer :: n
-  real :: gamser, a, x, gln, ap, summ, del
-  real, external :: gammln
+  real(KIND=dp) :: gamser, a, x, gln, ap, summ, del
+  real(KIND=dp), external :: gammln
 
   integer,parameter :: itmax=100
@@ -134 +140 @@
 subroutine gcf(gammcf,a,x,gln)
 
+  use par_mod, only: dp
+
   implicit none
 
   integer :: n
-  real :: gammcf, a, x, gln, gold, a0, a1, b0, b1, fac, an, anf, ana, g
-  real, external :: gammln
+  real(KIND=4) :: gammcf, x, a, gln, gold, a0, a1, b0, b1, fac, an, anf, ana, g
+  real(KIND=dp), external :: gammln
 
   integer,parameter :: itmax=100
@@ -172 +180 @@
 function erf(x)
 
-  implicit none
-
-  real :: x, erf
-  real, external :: gammp
+  use par_mod, only: dp
+
+  implicit none
+
+  real(KIND=dp) :: x, erf
+  real(KIND=dp), external :: gammp
 
   if(x.lt.0.)then
@@ -186 +196 @@
 function erfc(x)
 
-  implicit none
-
-  real :: x, erfc
-  real, external :: gammp, gammq
+  use par_mod, only: dp
+
+  implicit none
+
+  real(KIND=dp) :: x, erfc
+  real(KIND=dp), external :: gammp, gammq
 
   if(x.lt.0.)then
@@ -200 +212 @@
 function erfcc(x)
 
-  implicit none
-
-  real :: x, z, t, erfcc
+  use par_mod, only: dp
+
+  implicit none
+
+  real(KIND=dp) :: x, z, t, erfcc
 
   z=abs(x)