write vp rho

xumi1993 · Jan 10, 2025 · ba0ff65 · ba0ff65
1 parent e44bbb6
commit ba0ff65
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/model.f90 b/src/model.f90
@@ -324,6 +324,10 @@ subroutine write_model(this, subname)
       call h%add('/lat',this%ygrids)
       call h%add('/dep',this%zgrids)
       call h%add('/vs',transpose_3(this%vs3d))
+      if (ap%inversion%use_alpha_beta_rho) then
+        call h%add('/vp',transpose_3(this%vp3d))
+        call h%add('/rho',transpose_3(this%rho3d))
+      endif
       call h%close()
     endif
     call synchronize_all()