No commit message

geodynamics · Oct 14, 2013 · 5d89a4b · 5d89a4b
1 parent d8da888
commit 5d89a4b
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 4 deletions.
diff --git a/calc_dyn_topo b/calc_dyn_topo
@@ -36,3 +36,6 @@ cat tmp.ab | \
     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
diff --git a/hc_invert_dtopo.c b/hc_invert_dtopo.c
@@ -88,16 +88,15 @@ int main(int argc, char **argv)
   {
     /* 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 */