formatting and adding missings d0s in eos (trivial)

Author: pmocz

eos/private/eos_helm_eval.f90

@@ -98,12 +98,12 @@ subroutine get_helm_for_eosdt( &
          d_dlnd(i_phase:i_latent_ddlnRho) = 0d0
 
          ! zero all components
-         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0
+         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
 
          ! mark this one
-         res(i_frac_HELM) = 1.0
+         res(i_frac_HELM) = 1.0d0
 
       end subroutine get_helm_for_eosdt
 

eos/private/eoscms_eval.f90

@@ -295,12 +295,12 @@ subroutine get_CMS_for_eosdt( &
       d_dlnd(i_phase:i_latent_ddlnRho) = 0d0
 
       ! zero all components
-      res(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-      d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-      d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0
+      res(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+      d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+      d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
 
       ! mark this one
-      res(i_frac_CMS) = 1.0
+      res(i_frac_CMS) = 1.0d0
 
    end subroutine get_CMS_for_eosdt
 

eos/private/eosdt_eval.f90

@@ -311,9 +311,9 @@ subroutine get_other_for_eosdt( &
             ierr)
 
          ! zero all frac components
-         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0
+         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
 
          skip = .false.
 
@@ -1013,12 +1013,12 @@ subroutine get_opal_scvh_for_eosdt( &
          d_dlnd(i_phase:i_latent_ddlnRho) = 0d0
 
          ! zero all components
-         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0
+         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
 
          ! mark this one
-         res(i_frac_OPAL_SCVH) = 1.0
+         res(i_frac_OPAL_SCVH) = 1.0d0
 
       end subroutine get_opal_scvh_for_eosdt
 
@@ -1055,12 +1055,12 @@ subroutine get_FreeEOS_for_eosdt( &
          d_dlnd(i_phase:i_latent_ddlnRho) = 0d0
 
          ! zero all components
-         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0
+         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
 
          ! mark this one
-         res(i_frac_FreeEOS) = 1.0
+         res(i_frac_FreeEOS) = 1.0d0
 
       end subroutine get_FreeEOS_for_eosdt
 

eos/private/eospc_eval.f90

@@ -139,12 +139,12 @@ subroutine get_PC_for_eosdt( &
          d_dlnd(i_latent_ddlnT:i_latent_ddlnRho) = 0d0
 
          ! zero all components
-         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0
+         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
 
          ! mark this one
-         res(i_frac_PC) = 1.0
+         res(i_frac_PC) = 1.0d0
 
       end subroutine get_PC_for_eosdt
 

eos/private/ideal.f90

@@ -51,12 +51,12 @@ subroutine get_ideal_for_eosdt(handle, dbg, Z, X, abar, zbar, species, chem_id,
    skip = .false.
 
    ! zero all components
-   res(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-   d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-   d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0
+   res(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+   d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+   d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
 
    ! mark this one
-   res(i_frac_ideal) = 1.0
+   res(i_frac_ideal) = 1.0d0
 
    end subroutine get_ideal_for_eosdt
 

eos/private/skye.f90

@@ -177,12 +177,12 @@ subroutine get_Skye_for_eosdt(handle, dbg, Z, X, abar, zbar, species, chem_id, n
          skip = .false.
 
          ! zero all components
-         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0
-         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0
+         res(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnd(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
+         d_dlnT(i_frac:i_frac+num_eos_frac_results-1) = 0.0d0
 
          ! mark this one
-         res(i_frac_Skye) = 1.0
+         res(i_frac_Skye) = 1.0d0
 
       end subroutine get_Skye_for_eosdt
 

eos/private/skye_thermodynamics.f90

@@ -77,7 +77,7 @@ subroutine compute_derived_quantities(temp, dens, s, e, p, cv, cp, chit, chid, g
       cp = cv * gam1 / chid
 
       ! Sound speed
-      cs = clight * sqrt(gam1 / (1d0 + (dens / p) * (e + clight**2)))
+      cs = clight * sqrt(gam1 / (1d0 + (dens / p) * (e + clight*clight)))
    end subroutine compute_derived_quantities
 
    !> Computes thermodynamic quantities from Skye and packs them into the EOS return vectors.

eos/public/eos_def.f90

@@ -646,6 +646,4 @@ end subroutine free_EosDT_XZ_Info
 
       end subroutine eos_def_shutdown
 
-
       end module eos_def
-

eos/public/eos_lib.f90

@@ -32,9 +32,7 @@ module eos_lib
 
       implicit none
 
-      contains  ! the procedure interface for the library
-      ! client programs should only call these routines.
-
+      contains
 
       subroutine eos_init( &
             eosDT_cache_dir, use_cache, info)
@@ -108,7 +106,6 @@ elemental real(dp) function Radiation_Energy(T)
       end function Radiation_Energy
 
 
-
       ! eos evaluation
 
       ! you can call these routines after you've allocated a handle.
@@ -171,9 +168,7 @@ subroutine eosDT_get_component( &
          ! INPUT
 
          integer, intent(in) :: handle  ! eos handle; from star, pass s% eos_handle
-
          integer, intent(in) :: which_eos  ! see eos_def: i_eos_<component>
-
          integer, intent(in) :: species  ! number of species
          integer, pointer :: chem_id(:)  ! maps species to chem id
             ! index from 1 to species
@@ -904,7 +899,6 @@ subroutine eos_get_control_namelist(handle, name, val, ierr)
          call get_eos_ptr(handle,rq,ierr)
          if(ierr/=0) return
          call get_eos_controls(rq, name, val, ierr)
-
       end subroutine eos_get_control_namelist
 
       subroutine eos_set_control_namelist(handle, name, val, ierr)
@@ -919,7 +913,6 @@ subroutine eos_set_control_namelist(handle, name, val, ierr)
          call get_eos_ptr(handle,rq,ierr)
          if(ierr/=0) return
          call set_eos_controls(rq, name, val, ierr)
-
       end subroutine eos_set_control_namelist
 
 

kap/private/op_def.f90

@@ -73,7 +73,7 @@ module op_def
    integer, dimension(17, 91, 25), save :: ne1p, ne2p, np, kp1, kp2, kp3, npp
    real, dimension(-1:28, 28, 91, 25), save :: fionp
    real, allocatable, dimension(:), save :: yy2, yx
-   INTEGER, allocatable, dimension(:), save :: nx
+   integer, allocatable, dimension(:), save :: nx
 
    integer, parameter :: op_cache_version = 1