remove dble casts

MFlowCode · Apr 21, 2023 · 8fa3412 · 8fa3412
1 parent 06bcf95
commit 8fa3412
Show file tree

Hide file tree

Showing 5 changed files with 8 additions and 9 deletions.
diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -100,7 +100,6 @@ elseif (CMAKE_Fortran_COMPILER_ID STREQUAL "Flang")
     add_compile_options(
         $<$<COMPILE_LANGUAGE:Fortran>:-Mfreeform>
         $<$<COMPILE_LANGUAGE:Fortran>:-Mpreprocess>
-        $<$<COMPILE_LANGUAGE:Fortran>:-fdefault-real-8>
     )
 elseif (CMAKE_Fortran_COMPILER_ID STREQUAL "Intel")
     add_compile_options($<$<COMPILE_LANGUAGE:Fortran>:-free>)

diff --git a/src/common/m_eigen_solver.f90 b/src/common/m_eigen_solver.f90
@@ -846,8 +846,8 @@ subroutine cbabk2(nm,nl,low,igh,scale,ml,zr,zi)
 !     ------------------------------------------------------------------
 !
       integer i,j,k,ml,nl,ii,nm,igh,low
-      double precision scale(nl),zr(nm,ml),zi(nm,ml)
-      double precision s
+      real(wp) scale(nl),zr(nm,ml),zi(nm,ml)
+      real(wp) s
 
       if (ml .eq. 0) go to 200
       if (igh .eq. low) go to 120

diff --git a/src/post_process/m_global_parameters.f90 b/src/post_process/m_global_parameters.f90
@@ -755,7 +755,7 @@ subroutine s_transcoeff(omega, peclet, Re_trans, Im_trans)
         c3 = (CEXP(c2) - CEXP(-c2))/(CEXP(c2) + CEXP(-c2)) ! TANH(c2)
         trans = ((c2/c3 - 1._wp)**(-1) - 3._wp/c1)**(-1) ! transfer function
 
-        Re_trans = dble(trans)
+        Re_trans = (trans)
         Im_trans = aimag(trans)
 
     end subroutine s_transcoeff
@@ -856,7 +856,7 @@ subroutine s_simpson(Npt)
         ! phi = ln( R0 ) & return R0
         do ir = 1, Npt
             phi(ir) = log(R0mn) &
-                      + dble(ir - 1)*log(R0mx/R0mn)/dble(Npt - 1)
+                      + (ir - 1)*log(R0mx/R0mn)/(Npt - 1)
             R0(ir) = exp(phi(ir))
         end do
         dphi = phi(2) - phi(1)

diff --git a/src/pre_process/m_global_parameters.fpp b/src/pre_process/m_global_parameters.fpp
@@ -759,7 +759,7 @@ contains
         c3 = (CEXP(c2) - CEXP(-c2))/(CEXP(c2) + CEXP(-c2)) ! TANH(c2)
         trans = ((c2/c3 - 1._wp)**(-1) - 3._wp/c1)**(-1) ! transfer function
 
-        Re_trans = dble(trans)
+        Re_trans = (trans)
         Im_trans = aimag(trans)
 
     end subroutine s_transcoeff
@@ -854,7 +854,7 @@ contains
         ! phi = ln( R0 ) & return R0
         do ir = 1, nb
             phi(ir) = log(R0mn) &
-                      + dble(ir - 1)*log(R0mx/R0mn)/dble(nb - 1)
+                      + (ir - 1)*log(R0mx/R0mn)/(nb - 1)
             R0(ir) = exp(phi(ir))
         end do
         dphi = phi(2) - phi(1)

diff --git a/src/simulation/m_global_parameters.fpp b/src/simulation/m_global_parameters.fpp
@@ -951,7 +951,7 @@ contains
         c3 = (CEXP(c2) - CEXP(-c2))/(CEXP(c2) + CEXP(-c2)) ! TANH(c2)
         trans = ((c2/c3 - 1._wp)**(-1) - 3._wp/c1)**(-1) ! transfer function
 
-        Re_trans = dble(trans)
+        Re_trans = (trans)
         Im_trans = aimag(trans)
 
     end subroutine s_transcoeff
@@ -1052,7 +1052,7 @@ contains
         ! phi = ln( R0 ) & return R0
         do ir = 1, nb
             phi(ir) = log(R0mn) &
-                      + dble(ir - 1)*log(R0mx/R0mn)/dble(nb - 1)
+                      + (ir - 1)*log(R0mx/R0mn)/(nb - 1)
             R0(ir) = exp(phi(ir))
         end do
         dphi = phi(2) - phi(1)