From 00bc0afc20ab8e74f06a64e9e34c38cd8ba86649 Mon Sep 17 00:00:00 2001 From: harriet-godwin Date: Wed, 12 Jun 2019 15:23:06 +0100 Subject: [PATCH] fix strain tensor field rotations in rotations.f90 Change from using longitude to phi --- src/readfields.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/readfields.f90 b/src/readfields.f90 index 5ebb025..20e8504 100644 --- a/src/readfields.f90 +++ b/src/readfields.f90 @@ -1293,7 +1293,7 @@ function load_fw_points(this, coordinates, source_params, model) load_fw_points(:,:,ipoint) = rotate_symm_tensor_voigt_src_to_xyz( & load_fw_points(:,:,ipoint), & - source_params%lon, this%ndumps ) + rotmesh_phi(ipoint), this%ndumps ) load_fw_points(:,:,ipoint) = rotate_symm_tensor_voigt_xyz_src_to_xyz_earth( & load_fw_points(:,:,ipoint), & @@ -1454,7 +1454,7 @@ function load_bw_points(this, coordinates, receiver) ! only need to rotate in case of vs if (this%strain_type.eq.'straintensor_full') then load_bw_points(:,:,ipoint) = rotate_symm_tensor_voigt_src_to_xyz(load_bw_points(:,:,ipoint), & - receiver%lon, this%ndumps) + rotmesh_phi(ipoint), this%ndumps) load_bw_points(:,:,ipoint) = rotate_symm_tensor_voigt_xyz_src_to_xyz_earth(load_bw_points(:,:,ipoint), & receiver%lon, receiver%colat, this%ndumps) end if