[cig-commits] r15804 - mc/3D/CitcomS/trunk/lib

[becker at geodynamics.org]

Author: becker
Date: 2009-10-13 09:37:49 -0700 (Tue, 13 Oct 2009)
New Revision: 15804

Modified:
   mc/3D/CitcomS/trunk/lib/Ggrd_handling.c
Log:
Fixed typo.



Modified: mc/3D/CitcomS/trunk/lib/Ggrd_handling.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Ggrd_handling.c	2009-10-12 21:40:46 UTC (rev 15803)
+++ mc/3D/CitcomS/trunk/lib/Ggrd_handling.c	2009-10-13 16:37:49 UTC (rev 15804)
@@ -701,13 +701,13 @@
     if(use_codes)
       vscale *=  E->data.radius_km*1e3/1e6*1e2*M_PI/180.;		/* for deg/Myr -> cm/yr conversion */
     if(E->parallel.me == 0)
-      fprintf(stderr,"ggrd_read_vtop_from_file: expecting velocity grids in cm/yr\n");
+      fprintf(stderr,"ggrd_read_vtop_from_file: expecting velocity grids in cm/yr\n, scaling factor: %g",vscale);
   }else{
     /* stress scale, from MPa */
-    vscale =  1e-6/(E->data.ref_viscosity*E->data.therm_diff/(E->data.radius_km*E->data.radius_km*1e6));
+    vscale =  1e6/(E->data.ref_viscosity*E->data.therm_diff/(E->data.radius_km*E->data.radius_km*1e6));
     if((!mode_warned) && (E->parallel.me == 0)){
       fprintf(stderr,"ggrd_read_vtop_from_file: WARNING: make sure traction control is what you want, not free slip\n");
-      fprintf(stderr,"ggrd_read_vtop_from_file: expecting traction grids in MPa\n");
+      fprintf(stderr,"ggrd_read_vtop_from_file: expecting traction grids in MPa, scaling factor: %g\n",vscale);
       mode_warned = TRUE;
     }
   }