Square components (issue andwatson#20)

JackDMcGrath · Jul 27, 2023 · 56ce0eb · 56ce0eb
1 parent 44c55d0
commit 56ce0eb
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/util/vel_decomp_vE_vUN.m b/util/vel_decomp_vE_vUN.m
@@ -142,9 +142,9 @@
 
 % same for the uncertainty
 if ~isempty(gnss_sN)
-    UN2U_var = (1 - sind(inc).^2 .* cosd(az).^2) ./ cosd(inc);
+    UN2U_var = (1 - sind(inc).^2 .* cosd(az).^2) ./ cosd(inc).^2;
     N2U_var = ((sind(az).^2).*(sind(inc).^2)) ./ (cosd(inc).^2);
-    var_up = sqrt((var_UN.^2 .* UN2U_var) + (gnss_sN .* N2U_var));
+    var_up = sqrt((var_UN.^2 .* UN2U_var) + ((gnss_sN.^2) .* N2U_var));
 else
     var_up = ones(size(m_up));
 end