Change "1.d0" to "1d0"

Author: cphyc

amr/adaptive_loop.f90

@@ -205,7 +205,7 @@ subroutine adaptive_loop
               output_now=.true.
               if(myid==1) write(*,*) 'Dumping snapshot before walltime runs out'
               ! Now set walltime to a negative number so we don't keep printing outputs
-              walltime_hrs = -1.d0
+              walltime_hrs = -1d0
            endif
         endif
      endif

amr/amr_parameters.f90

@@ -260,9 +260,9 @@ module amr_parameters
   real(dp),dimension(1:MAXREGION)   ::x_center=0
   real(dp),dimension(1:MAXREGION)   ::y_center=0
   real(dp),dimension(1:MAXREGION)   ::z_center=0
-  real(dp),dimension(1:MAXREGION)   ::length_x=1.d10
-  real(dp),dimension(1:MAXREGION)   ::length_y=1.d10
-  real(dp),dimension(1:MAXREGION)   ::length_z=1.d10
+  real(dp),dimension(1:MAXREGION)   ::length_x=1d10
+  real(dp),dimension(1:MAXREGION)   ::length_y=1d10
+  real(dp),dimension(1:MAXREGION)   ::length_z=1d10
   real(dp),dimension(1:MAXREGION)   ::exp_region=2
 
   ! Boundary conditions parameters

amr/constants.f90

@@ -3,7 +3,7 @@ module constants
 
   ! Numerical constants
   real(dp),parameter ::twopi        = 6.2831853d0
-  real(dp),parameter ::pi           = twopi/2.d0
+  real(dp),parameter ::pi           = twopi/2d0
 
   real(dp),parameter ::mu_mol       = 1.2195d0
 

amr/init_time.f90

@@ -45,7 +45,7 @@ subroutine init_time
      ! Compute Friedman model look up table
      if(myid==1)write(*,*)'Computing Friedman model'
      call friedman(dble(omega_m),dble(omega_l),dble(omega_k), &
-          & 1.d-6,dble(aexp_ini), &
+          & 1d-6,dble(aexp_ini), &
           & aexp_frw,hexp_frw,tau_frw,t_frw,n_frw)
 
      ! Compute initial conformal time
@@ -152,7 +152,7 @@ subroutine init_time
      if(cosmo) then
         my_grackle_units%comoving_coordinates = 1
         ! Reonization redshift has to be later than starting redshift
-        z_reion=min(1.d0/(1.1d0*aexp_ini)-1.d0,z_reion)
+        z_reion=min(1d0/(1.1d0*aexp_ini)-1d0,z_reion)
         ! Approximate initial temperature
         T2_start=1.356d-2/aexp_ini**2
         if(nrestart==0)then
@@ -232,7 +232,7 @@ subroutine init_time
         endif
         if(cosmo)then
            ! Reonization redshift has to be later than starting redshift
-           z_reion=min(1.d0/(1.1d0*aexp_ini)-1.d0,z_reion)
+           z_reion=min(1d0/(1.1d0*aexp_ini)-1d0,z_reion)
            call set_model(Nmodel,dble(J21*1d-21),-1.0d0,dble(a_spec),-1.0d0,dble(z_reion), &
                 & -1,2, &
                 & dble(h0/100),dble(omega_b),dble(omega_m),dble(omega_l), &
@@ -244,7 +244,7 @@ subroutine init_time
         else
            call set_model(Nmodel,dble(J21*1d-21),-1.0d0,dble(a_spec),-1.0d0,dble(z_reion), &
                 & -1,2, &
-                & dble(70.d0/100.d0),0.04d0,0.3d0,0.7d0, &
+                & dble(70d0/100d0),0.04d0,0.3d0,0.7d0, &
                 & dble(aexp_ini),T2_sim)
         endif
      end if
@@ -258,7 +258,7 @@ subroutine init_time
      endif
      if(cosmo)then
         ! Reonization redshift has to be later than starting redshift
-        z_reion=min(1.d0/(1.1d0*aexp_ini)-1.d0,z_reion)
+        z_reion=min(1d0/(1.1d0*aexp_ini)-1d0,z_reion)
         call set_model(Nmodel,dble(J21*1d-21),-1.0d0,dble(a_spec),-1.0d0,dble(z_reion), &
              & -1,2, &
              & dble(h0/100.),dble(omega_b),dble(omega_m),dble(omega_l), &
@@ -284,7 +284,7 @@ subroutine init_time
      endif
      if(cosmo)then
         ! Reonization redshift has to be later than starting redshift
-        z_reion=min(1.d0/(1.1d0*aexp_ini)-1.d0,z_reion)
+        z_reion=min(1d0/(1.1d0*aexp_ini)-1d0,z_reion)
         call rt_set_model(dble(h0/100.),dble(omega_b),dble(omega_m),dble(omega_l), &
              & dble(aexp_ini),T2_sim)
         T2_start=T2_sim
@@ -553,7 +553,7 @@ subroutine init_cosmo
      endif
      omega_m = gadgetheader%omega0
      omega_l = gadgetheader%omegalambda
-     h0 = gadgetheader%hubbleparam * 100.d0
+     h0 = gadgetheader%hubbleparam * 100d0
      boxlen_ini = gadgetheader%boxsize
      aexp = gadgetheader%time
      aexp_ini = aexp
@@ -578,7 +578,7 @@ subroutine init_cosmo
      write(*,'(" omega_m=",F7.3," omega_l=",F7.3," omega_b=",F7.3)')omega_m,omega_l,omega_b
      write(*,'(" box size=",1pe10.3," h-1 Mpc")')boxlen_ini
   end if
-  omega_k=1.d0-omega_l-omega_m
+  omega_k=1d0-omega_l-omega_m
 
   ! Compute linear scaling factor between aexp and astart(ilevel)
   do ilevel=levelmin,nlevelmax_part
@@ -631,8 +631,8 @@ function fy(a)
     real(dp)::fy
     real(dp)::y,a
 
-    y=omega_m*(1.d0/a-1.d0) + omega_l*(a*a-1.d0) + 1.d0
-    fy=1.d0/y**1.5d0
+    y=omega_m*(1d0/a-1d0) + omega_l*(a*a-1d0) + 1d0
+    fy=1d0/y**1.5d0
 
     return
   end function fy
@@ -649,7 +649,7 @@ function d1a(a)
        write(*,*)'a=',a
        call clean_stop
     end if
-    y=omega_m*(1.d0/a-1.d0) + omega_l*(a*a-1.d0) + 1.d0
+    y=omega_m*(1d0/a-1d0) + omega_l*(a*a-1d0) + 1d0
     if(y .lt. 0.0D0)then
        write(*,*)'y=',y
        call clean_stop
@@ -672,7 +672,7 @@ end function d1a
 !!$    a=1.e-7
 !!$    niter=0
 !!$10  niter=niter+1
-!!$    da=(d1/d1a(a)-1.d0)/fpeebl(a)*a
+!!$    da=(d1/d1a(a)-1d0)/fpeebl(a)*a
 !!$    a=a+da
 !!$    if (abs(da).gt.1.0e-8.and.niter.lt.10) go to 10
 !!$    ad1=a
@@ -686,9 +686,9 @@ function fpeebl(a)
     real(dp) :: fact,y,eps
 
     eps=1.0d-6
-    y=omega_m*(1.d0/a-1.d0) + omega_l*(a*a-1.d0) + 1.d0
+    y=omega_m*(1d0/a-1d0) + omega_l*(a*a-1d0) + 1d0
     fact=rombint(eps,a,eps)
-    fpeebl=(omega_l*a*a-0.5d0*omega_m/a)/y - 1.d0 + a*fy(a)/fact
+    fpeebl=(omega_l*a*a-0.5d0*omega_m/a)/y - 1d0 + a*fy(a)/fact
     return
   end function fpeebl
   !cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
@@ -796,12 +796,12 @@ subroutine friedman(O_mat_0,O_vac_0,O_k_0,alpha,axp_min, &
 
      nstep = nstep + 1
      dtau = alpha * axp_tau / dadtau(axp_tau,O_mat_0,O_vac_0,O_k_0)
-     axp_tau_pre = axp_tau - dadtau(axp_tau,O_mat_0,O_vac_0,O_k_0)*dtau/2.d0
+     axp_tau_pre = axp_tau - dadtau(axp_tau,O_mat_0,O_vac_0,O_k_0)*dtau/2d0
      axp_tau = axp_tau - dadtau(axp_tau_pre,O_mat_0,O_vac_0,O_k_0)*dtau
      tau = tau - dtau
 
      dt = alpha * axp_t / dadt(axp_t,O_mat_0,O_vac_0,O_k_0)
-     axp_t_pre = axp_t - dadt(axp_t,O_mat_0,O_vac_0,O_k_0)*dt/2.d0
+     axp_t_pre = axp_t - dadt(axp_t,O_mat_0,O_vac_0,O_k_0)*dt/2d0
      axp_t = axp_t - dadt(axp_t_pre,O_mat_0,O_vac_0,O_k_0)*dt
      t = t - dt
 
@@ -814,10 +814,10 @@ subroutine friedman(O_mat_0,O_vac_0,O_k_0,alpha,axp_min, &
 
   nskip=nstep/ntable
 
-  axp_t = 1.d0
-  t = 0.d0
-  axp_tau = 1.d0
-  tau = 0.d0
+  axp_t = 1d0
+  t = 0d0
+  axp_tau = 1d0
+  tau = 0d0
   nstep = 0
   nout=0
   t_out(nout)=t
@@ -829,12 +829,12 @@ subroutine friedman(O_mat_0,O_vac_0,O_k_0,alpha,axp_min, &
 
      nstep = nstep + 1
      dtau = alpha * axp_tau / dadtau(axp_tau,O_mat_0,O_vac_0,O_k_0)
-     axp_tau_pre = axp_tau - dadtau(axp_tau,O_mat_0,O_vac_0,O_k_0)*dtau/2.d0
+     axp_tau_pre = axp_tau - dadtau(axp_tau,O_mat_0,O_vac_0,O_k_0)*dtau/2d0
      axp_tau = axp_tau - dadtau(axp_tau_pre,O_mat_0,O_vac_0,O_k_0)*dtau
      tau = tau - dtau
 
      dt = alpha * axp_t / dadt(axp_t,O_mat_0,O_vac_0,O_k_0)
-     axp_t_pre = axp_t - dadt(axp_t,O_mat_0,O_vac_0,O_k_0)*dt/2.d0
+     axp_t_pre = axp_t - dadt(axp_t,O_mat_0,O_vac_0,O_k_0)*dt/2d0
      axp_t = axp_t - dadt(axp_t_pre,O_mat_0,O_vac_0,O_k_0)*dt
      t = t - dt
 

amr/light_cone.f90

@@ -38,8 +38,8 @@ subroutine output_cone()
 
   if(nstep_coarse.lt.2)return
 
-  z2=1/aexp_old-1.d0
-  z1=1/aexp-1.d0
+  z2=1/aexp_old-1d0
+  z1=1/aexp-1d0
 
   if(z2.gt.zmax_cone)return
   if(abs(z2-z1)<1d-6)return
@@ -738,7 +738,7 @@ function coord_distance(zz,Omega0,OmegaL,OmegaR,coverH0)
   real(kind=8) :: z,res,coord_distance,zz
   real(kind=8) :: Omega0,OmegaL,OmegaR,coverH0
   z=abs(zz)
-  call qromb(0.d0,z,res,omega0,omegaL,OmegaR)
+  call qromb(0d0,z,res,omega0,omegaL,OmegaR)
   coord_distance=coverH0*res
   if (zz.lt.0) coord_distance=-coord_distance
 end function coord_distance
@@ -750,15 +750,15 @@ function funcE(z,Omega0,OmegaL,OmegaR)
   real(kind=8) :: funcE,z,HsurH0
   real(kind=8) :: omega0,omegaL,OmegaR
 
-  funcE=1.d0/HsurH0(z,Omega0,OmegaL,OmegaR)
+  funcE=1d0/HsurH0(z,Omega0,OmegaL,OmegaR)
 end function funcE
 
 !===========================================================================
 function HsurH0(z,omega0,omegaL,OmegaR)
   !===========================================================================
   implicit none
   real(kind=8) :: z,omega0,omegaL,OmegaR,HsurH0
-  HsurH0=sqrt(Omega0*(1.d0+z)**3+OmegaR*(1.d0+z)**2+OmegaL)
+  HsurH0=sqrt(Omega0*(1d0+z)**3+OmegaR*(1d0+z)**2+OmegaL)
 end function HsurH0
 
 
@@ -768,15 +768,15 @@ SUBROUTINE qromb(a,b,ss,omega0,omegaL,OmegaR)
   implicit none
   INTEGER :: JMAX,JMAXP,K,KM
   REAL(kind=8) :: a,b,ss,EPS,omega0,omegaL,OmegaR
-  PARAMETER (EPS=1.d-6, JMAX=20, JMAXP=JMAX+1, K=5, KM=K-1)
+  PARAMETER (EPS=1d-6, JMAX=20, JMAXP=JMAX+1, K=5, KM=K-1)
   !  USES polint,trapzd
   INTEGER :: j
   REAL(kind=8) :: dss,h(JMAXP),s(JMAXP)
   h(1)=1.
   do j=1,JMAX
      call trapzd(a,b,s(j),j,omega0,omegaL,OmegaR)
      if (j.ge.K) then
-        call polint(h(j-KM),s(j-KM),K,0.d0,ss,dss)
+        call polint(h(j-KM),s(j-KM),K,0d0,ss,dss)
         if (abs(dss).le.EPS*abs(ss)) return
      endif
      s(j+1)=s(j)

amr/load_balance.f90

@@ -742,18 +742,18 @@ subroutine cmp_ordering(x,order,nn)
      zc=boxlen/2
 #endif
      do i=1,nn
-        xx=x(i,1)-xc+1.d-10
+        xx=x(i,1)-xc+1d-10
         yy=x(i,2)-yc
 #if NDIM>2
         zz=x(i,3)
 #endif
         if(xx>0.)then
-           order(i)=atan(yy/xx)+acos(-1.d0)/2.d0
+           order(i)=atan(yy/xx)+acos(-1d0)/2d0
         else
-           order(i)=atan(yy/xx)+acos(-1.d0)*3.d0/2.d0
+           order(i)=atan(yy/xx)+acos(-1d0)*3d0/2d0
         endif
 #if NDIM>2
-        if(zz.gt.zc)order(i)=order(i)+2.d0*acos(-1.d0)
+        if(zz.gt.zc)order(i)=order(i)+2d0*acos(-1d0)
 #endif
      end do
   end if
@@ -861,53 +861,53 @@ subroutine cmp_minmaxorder(x,order_min,order_max,dx,nn)
      do i=1,nn
         if(dx==boxlen)then
            order_min(i)=0
-           order_max(i)=4.d0*acos(-1.d0)
+           order_max(i)=4d0*acos(-1d0)
         else
            ! x- y-
            yy=x(i,2)-yc-dxx
            xx=x(i,1)-xc-dxx
            if(xx.ge.0.)then
               xx=xx+1d-10
-              theta1=atan(yy/xx)+acos(-1.d0)/2
+              theta1=atan(yy/xx)+acos(-1d0)/2
            else
               xx=xx-1d-10
-              theta1=atan(yy/xx)+acos(-1.d0)*3.d0/2.d0
+              theta1=atan(yy/xx)+acos(-1d0)*3d0/2d0
            endif
            ! x+ y-
            xx=x(i,1)-xc+dxx
            if(xx.gt.0.)then
               xx=xx+1d-10
-              theta2=atan(yy/xx)+acos(-1.d0)/2
+              theta2=atan(yy/xx)+acos(-1d0)/2
            else
               xx=xx-1d-10
-              theta2=atan(yy/xx)+acos(-1.d0)*3.d0/2.d0
+              theta2=atan(yy/xx)+acos(-1d0)*3d0/2d0
            endif
 
            ! x+ y+
            yy=x(i,2)-yc+dxx
            if(xx.gt.0.)then
               xx=xx+1d-10
-              theta3=atan(yy/xx)+acos(-1.d0)/2.d0
+              theta3=atan(yy/xx)+acos(-1d0)/2d0
            else
               xx=xx-1d-10
-              theta3=atan(yy/xx)+acos(-1.d0)*3.d0/2.d0
+              theta3=atan(yy/xx)+acos(-1d0)*3d0/2d0
            endif
            ! x- y+
            xx=x(i,1)-xc-dxx
            if(xx.ge.0.)then
               xx=xx+1d-10
-              theta4=atan(yy/xx)+acos(-1.d0)/2.d0
+              theta4=atan(yy/xx)+acos(-1d0)/2d0
            else
               xx=xx-1d-10
-              theta4=atan(yy/xx)+acos(-1.d0)*3.d0/2.d0
+              theta4=atan(yy/xx)+acos(-1d0)*3d0/2d0
            endif
            order_min(i)=min(theta1,theta2,theta3,theta4)
            order_max(i)=max(theta1,theta2,theta3,theta4)
 #if NDIM>2
            zz=x(i,3)
            if(zz.gt.zc)then
-              order_min(i)=order_min(i)+2.d0*acos(-1.d0)
-              order_max(i)=order_max(i)+2.d0*acos(-1.d0)
+              order_min(i)=order_min(i)+2d0*acos(-1d0)
+              order_max(i)=order_max(i)+2d0*acos(-1d0)
            endif
 #endif
         endif

amr/movie.f90

@@ -281,9 +281,9 @@ subroutine output_frame()
   if(cosmo) then
      if(tend_theta_camera(proj_ind).le.0d0) tend_theta_camera(proj_ind) = aendmov
      if(tend_phi_camera(proj_ind).le.0d0) tend_phi_camera(proj_ind) = aendmov
-     theta_cam  = theta_camera(proj_ind)*pi/180.d0                                                                                 &
+     theta_cam  = theta_camera(proj_ind)*pi/180d0                                                                                 &
                 +min(max(aexp-tstart_theta_camera(proj_ind),0d0),tend_theta_camera(proj_ind))*dtheta_camera(proj_ind)*pi/180./(aendmov-astartmov)
-     phi_cam    = phi_camera(proj_ind)*pi/180.d0                                                                                   &
+     phi_cam    = phi_camera(proj_ind)*pi/180d0                                                                                   &
                 +min(max(aexp-tstart_theta_camera(proj_ind),0d0),tend_phi_camera(proj_ind))*dphi_camera(proj_ind)*pi/180./(aendmov-astartmov)
      dist_cam   = dist_camera(proj_ind)+min(max(aexp-tstart_theta_camera(proj_ind),0d0),tend_theta_camera(proj_ind))*ddist_camera(proj_ind)/(aendmov-astartmov)
   else

amr/output_amr.f90

@@ -595,9 +595,9 @@ subroutine savegadget(filename)
   header%npart = 0
   header%npart(2) = npart
   header%mass = 0
-  header%mass(2) = omega_m*RHOcrit*(boxlen_ini)**3/npart_tot/1.d10
+  header%mass(2) = omega_m*RHOcrit*(boxlen_ini)**3/npart_tot/1d10
   header%time = aexp
-  header%redshift = 1.d0/aexp-1.d0
+  header%redshift = 1d0/aexp-1d0
   header%flag_sfr = 0
   header%nparttotal = 0
 #ifndef LONGINT

amr/update_time.f90

@@ -560,7 +560,7 @@ SUBROUTINE getProperTime(tau,tproper)
   implicit none
   real(dp)::tau, tproper
   integer::i
-  if(.not. cosmo .or. tau .eq. 0.d0) then ! this might happen quite often
+  if(.not. cosmo .or. tau .eq. 0d0) then ! this might happen quite often
      tproper = tau
      return
   endif