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

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/calypso

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

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

commit 6a778c42c8d02afaae9c8cd4b55e525317aaf635
Author: Hiroaki Matsui <h_kemono at mac.com>
Date:   Mon Sep 23 16:10:51 2013 -0700

    Fix initial temperature


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

6a778c42c8d02afaae9c8cd4b55e525317aaf635
 .../data_utilities/INITIAL_FIELD/const_sph_initial_spectr.f90         | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/programs/data_utilities/INITIAL_FIELD/const_sph_initial_spectr.f90 b/src/programs/data_utilities/INITIAL_FIELD/const_sph_initial_spectr.f90
index 5857d3b..ac349b4 100644
--- a/src/programs/data_utilities/INITIAL_FIELD/const_sph_initial_spectr.f90
+++ b/src/programs/data_utilities/INITIAL_FIELD/const_sph_initial_spectr.f90
@@ -169,7 +169,7 @@
 !    set initial temperature
           xr = two * rr - one * (r_CMB+r_ICB) / shell
           d_rj(inod,ipol%i_temp) = (one-three*xr**2+three*xr**4-xr**6)  &
-     &                            * A_temp * six / (sqrt(pi))
+     &                            * A_temp * three / (two*sqrt(pi))
         end do
       end if
 !
@@ -217,7 +217,7 @@
           inod = local_sph_data_address(k,jj)
           xr = two * radius_1d_rj_r(k) - one * (r_CMB+r_ICB) / shell
           d_rj(inod,ipol%i_light) = (one-three*xr**2+three*xr**4-xr**6) &
-     &                       * A_light * six / (sqrt(pi))
+     &                            * A_light * three / (two*sqrt(pi))
         end do
       end if
 !