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

[becker at geodynamics.org]

Author: becker
Date: 2013-10-14 10:42:47 -0700 (Mon, 14 Oct 2013)
New Revision: 22961

Modified:
   mc/1D/hc/trunk/calc_dyn_topo
   mc/1D/hc/trunk/hc_invert_dtopo.c
Log:


Modified: mc/1D/hc/trunk/calc_dyn_topo
===================================================================
--- mc/1D/hc/trunk/calc_dyn_topo	2013-10-13 23:01:39 UTC (rev 22960)
+++ mc/1D/hc/trunk/calc_dyn_topo	2013-10-14 17:42:47 UTC (rev 22961)
@@ -36,3 +36,6 @@
     gawk -v fac=$fac '{print($1,$2,$3*fac)}' | \
     xyz2grd -fg -Rg -I1 -Nnan  -Gdtopo.grd
 grd2map dtopo.grd
+
+
+hc_extract_sh_layer dscaled.sol.bin 2 4
\ No newline at end of file

Modified: mc/1D/hc/trunk/hc_invert_dtopo.c
===================================================================
--- mc/1D/hc/trunk/hc_invert_dtopo.c	2013-10-13 23:01:39 UTC (rev 22960)
+++ mc/1D/hc/trunk/hc_invert_dtopo.c	2013-10-14 17:42:47 UTC (rev 22961)
@@ -88,16 +88,15 @@
   {
     /* compute solution */
     hc_solve(model,p->free_slip,p->solution_mode,sol_spectral,
-	     (solved)?(FALSE):(TRUE), /* density changed? */
+	     TRUE, /* density changed? */
 	     (solved)?(FALSE):(TRUE), /* plate velocity changed? */
 	     TRUE,			/* viscosity changed */
 	     FALSE,p->compute_geoid,
 	     pvel,model->dens_anom,geoid,
 	     p->verbose);
     /* extract the top tractions */
-    hc_compute_dynamic_topography(model,sol_spectral,&dtopo,FALSE,p->verbose);
-    
-    sh_single_par_and_exp_to_file(dtopo,"dtopo.ab",TRUE,p->verbose);
+    hc_compute_dynamic_topography(model,sol_spectral,&dtopo,TRUE,p->verbose);
+    //sh_single_par_and_exp_to_file(dtopo,"dtopo.ab",TRUE,p->verbose);
 
     /* geoid correlation */
     hc_compute_correlation(geoid,p->ref_geoid,(gcorr),0,p->verbose); /* full correlation */