[cig-commits] [commit] Hiro_latest: Fix initial temperature (ca32efb)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/calypso

On branch  : Hiro_latest
Link       : https://github.com/geodynamics/calypso/compare/93e9f8f974c7a247c8f02e54ec18de063f86c8fb...3c548304673360ddedd7d68c8095b3fb74a2b9ce

>---------------------------------------------------------------

commit ca32efbfa8b2cea3879f6ce2543b58444e597e20
Author: Hiroaki Matsui <h_kemono at mac.com>
Date:   Mon Sep 23 16:09:13 2013 -0700

    Fix initial temperature


>---------------------------------------------------------------

ca32efbfa8b2cea3879f6ce2543b58444e597e20
 .../MHD_src/sph_MHD/set_initial_sph_dynamo.f90               | 12 ++++++------
 1 file changed, 6 insertions(+), 6 deletions(-)

diff --git a/src/Fortran_libraries/MHD_src/sph_MHD/set_initial_sph_dynamo.f90 b/src/Fortran_libraries/MHD_src/sph_MHD/set_initial_sph_dynamo.f90
index 9841c0e..584d9a9 100644
--- a/src/Fortran_libraries/MHD_src/sph_MHD/set_initial_sph_dynamo.f90
+++ b/src/Fortran_libraries/MHD_src/sph_MHD/set_initial_sph_dynamo.f90
@@ -23,7 +23,7 @@
       implicit none
 !
       private :: set_initial_temp_sph, set_initial_magne_sph
-      private :: set_initial_scalar_sph, set_initial_velo_sph
+      private :: set_initial_light_sph, set_initial_velo_sph
       private :: set_ini_reference_temp_sph
 !
 !-----------------------------------------------------------------------
@@ -64,7 +64,7 @@
           call set_initial_temp_sph(isig)
         end if
         if(ipol%i_light .gt. 0) then
-          call set_initial_scalar_sph(isig, ipol%i_light)
+          call set_initial_light_sph(isig, ipol%i_light)
         end if
 !
         if(iflag_restart .eq. i_rst_dbench1) then
@@ -211,7 +211,7 @@
           inod = jj + (k-1)*nidx_rj(2)
 !
           d_rj(inod,ipol%i_temp) = (one-three*xr**2+three*xr**4-xr**6)  &
-     &                            * 0.1d0 * six / (sqrt(pi))
+     &                            * 0.1d0 * three / (two*sqrt(pi))
         end do
       end if
 !
@@ -219,7 +219,7 @@
 !
 !-----------------------------------------------------------------------
 !
-      subroutine set_initial_scalar_sph(isig, is_fld)
+      subroutine set_initial_light_sph(isig, is_fld)
 !
       use m_spheric_parameter
       use m_sph_spectr_data
@@ -257,11 +257,11 @@
           xr = two * radius_1d_rj_r(k) - one * (r_CMB+r_ICB) / shell
           inod = jj + (k-1)*nidx_rj(2)
           d_rj(inod,is_fld) = (one-three*xr**2+three*xr**4-xr**6)      &
-     &                       * 0.1d0 * six / (sqrt(pi))
+     &                            * 0.1d0 * three / (two*sqrt(pi))
         end do
       end if
 !
-      end subroutine set_initial_scalar_sph
+      end subroutine set_initial_light_sph
 !
 !-----------------------------------------------------------------------
 !