fixed gaussian weights

martalmeida · Jan 18, 2021 · 6173835 · 6173835
1 parent 7c00ef3
commit 6173835
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 1 deletion.
diff --git a/nmf3d_mat/hvf_baroclinic.m b/nmf3d_mat/hvf_baroclinic.m
@@ -38,6 +38,7 @@
   x   = sind(LAT);
 elseif isequal(latType,'gaussian')
   [x,w] = lgwt(dlat,-1,1);
+  x=x(end:-1:1);
 end
 
 disp('Part I');

diff --git a/nmf3d_mat/lgwt.m b/nmf3d_mat/lgwt.m
@@ -59,7 +59,7 @@
 
 % changed nov 2020:
 % note that w is always symmetric, no need to change
-x=x(end:-1:1);
+%x=x(end:-1:1);
 
 % Compute the weights
 w=(b-a)./((1-y.^2).*Lp.^2)*(N2/N1)^2;