[cig-commits] r16621 - in mc/1D/hc/trunk: . prem

[becker at geodynamics.org]

Author: becker
Date: 2010-05-03 22:36:41 -0700 (Mon, 03 May 2010)
New Revision: 16621

Modified:
   mc/1D/hc/trunk/hc_init.c
   mc/1D/hc/trunk/prem.h
   mc/1D/hc/trunk/prem/prem.dat
Log:
Minor radius compatibility fix



Modified: mc/1D/hc/trunk/hc_init.c
===================================================================
--- mc/1D/hc/trunk/hc_init.c	2010-05-04 04:25:20 UTC (rev 16620)
+++ mc/1D/hc/trunk/hc_init.c	2010-05-04 05:36:41 UTC (rev 16621)
@@ -273,8 +273,10 @@
 
   */
   hc->re = hc->prem->r0;
-  if(fabs(hc->re - (HC_RE_KM * 1e3)) > 1e-7)
-    HC_ERROR("hc_init_constants","Earth radius mismatch")
+  if(fabs(hc->re - (HC_RE_KM * 1e3)) > 1e-7){
+    fprintf(stderr,"%.7e %.7e\n",HC_RE_KM * 1e3,hc->re);
+    HC_ERROR("hc_init_constants","Earth radius mismatch");
+  }
 
   hc->secyr = HC_SECYR;	/* seconds/year  */
 

Modified: mc/1D/hc/trunk/prem/prem.dat
===================================================================
--- mc/1D/hc/trunk/prem/prem.dat	2010-05-04 04:25:20 UTC (rev 16620)
+++ mc/1D/hc/trunk/prem/prem.dat	2010-05-04 05:36:41 UTC (rev 16621)
@@ -1,5 +1,5 @@
         13 4
-          6371000.
+          6371008.7714
           1221500.  3480000.  3630000.  5600000.  5701000.  5771000.
             5971000.  6151000.  6291000.  6346600.  6356000.  6368000.
             6371000.

Modified: mc/1D/hc/trunk/prem.h
===================================================================
--- mc/1D/hc/trunk/prem.h	2010-05-04 04:25:20 UTC (rev 16620)
+++ mc/1D/hc/trunk/prem.h	2010-05-04 05:36:41 UTC (rev 16621)
@@ -66,7 +66,7 @@
    radius of earth in km 
 
 */
-#define PREM_RE_KM 6371.0
+#define PREM_RE_KM 6371.0087714
 
 /*