[cig-commits] [commit] : added comment block for outputting rho for old socal model (bb18476)

[cig_noreply at geodynamics.org]

Repository : ssh://geoshell/specfem3d

On branch  : 
Link       : https://github.com/geodynamics/specfem2d/compare/1e201257d91c794056b990a43329e05d04f77454...0000000000000000000000000000000000000000

>---------------------------------------------------------------

commit bb184765c02be1a0fd33197d49f4560d4c8cc098
Author: Carl Tape <carltape at gi.alaska.edu>
Date:   Thu Feb 3 23:14:32 2011 +0000

    added comment block for outputting rho for old socal model


>---------------------------------------------------------------

bb184765c02be1a0fd33197d49f4560d4c8cc098
 save_arrays_solver.f90 | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/save_arrays_solver.f90 b/save_arrays_solver.f90
index b248ff6..5190f01 100644
--- a/save_arrays_solver.f90
+++ b/save_arrays_solver.f90
@@ -218,6 +218,11 @@
 !!$  open(unit=27,file=prname(1:len_trim(prname))//'vs.bin',status='unknown',form='unformatted')
 !!$  write(27) mustore / rho_vs
 !!$  close(27)
+!!$
+!!$! rho (for checking the mesh and model)
+!!$  open(unit=27,file=prname(1:len_trim(prname))//'rho.bin',status='unknown',form='unformatted')
+!!$  write(27) rho_vp**2 / (FOUR_THIRDS * mustore + kappastore)
+!!$  close(27)
 
 ! kappa
   open(unit=27,file=prname(1:len_trim(prname))//'kappa.bin',status='unknown',form='unformatted')