[cig-commits] r12354 - mc/1D/hc/trunk

[becker at geodynamics.org]

Author: becker
Date: 2008-06-27 16:23:17 -0700 (Fri, 27 Jun 2008)
New Revision: 12354

Modified:
   mc/1D/hc/trunk/hc_init.c
Log:
Modified output and added check for low densities (e.g. in PREM water layer).



Modified: mc/1D/hc/trunk/hc_init.c
===================================================================
--- mc/1D/hc/trunk/hc_init.c	2008-06-27 22:31:06 UTC (rev 12353)
+++ mc/1D/hc/trunk/hc_init.c	2008-06-27 23:23:17 UTC (rev 12354)
@@ -704,7 +704,7 @@
 	  }
 	  reported = TRUE;
 	  if(verbose >= 2)
-	    fprintf(stderr,"hc_assign_density: non_dim radius                 %% factor    PREM \\rho                  layer #            depth[km]\n");
+	    fprintf(stderr,"hc_assign_density: non_dim radius                 %% factor    PREM \\rho/mean_rho          layer #            depth[km]\n");
 	}
 
 	/* 
@@ -729,6 +729,9 @@
 	  */
 	  prem_get_rho(&rho0,hc->rden[hc->inho],hc->prem);
 	  rho0 /= 1000.0;
+	  if(rho0 < 3)
+	    fprintf(stderr,"\nhc_assign_density: WARNING: using small (%g) density from PREM for layer at depth %g\n\n",
+		    rho0*1000,HC_Z_DEPTH(hc->rden[hc->inho]));
 	}else{
 	  /* mean value */
 	  rho0 =  hc->avg_den_mantle;
@@ -742,7 +745,7 @@
 	  
 	  fprintf(stderr,"hc_assign_density: r: %11g anom scales: %11g x %11g = %11g\t%5i out of %i, z: %11g\n",
 		  hc->rden[hc->inho],
-		  HC_DENSITY_SCALING,rho0,dens_scale[0],hc->inho+1,nset,zlabel);
+		  HC_DENSITY_SCALING,rho0/ hc->avg_den_mantle,dens_scale[0],hc->inho+1,nset,zlabel);
 	}
 	if(hc->inho){	
 	  /*