# Problem with shaec

From: Wood <564104718_at_nyahnyahspammersnyahnyah>
Date: Sun Oct 27 2013 - 06:49:59 MDT

Hello ,
I Computes spherical harmonic analysis of a global 500hpa Geopotential height on a fixed grid(73x144),but "a" is ​​greater than I expected,I wonder if someone can help me
nlat = 73
nlon =144
N = nlat
M = nlon
if (M%2 .eq.0) then ; note % is NCL's modulus operator
N = min((/ nlat, (M+2)/2 /))
else ; nlon must be odd
N = min((/ nlat, (M+1)/2 /))
end if
a = new ( (/nlat,N/), float)
b = new ( (/nlat,N/), float)

shaec(g,a,b)
print(a)
Variable: a (subsection)
Type: float
Total Size: 292 bytes
73 values
Number of Dimensions: 1
Dimensions and sizes: [73]
Coordinates:
Number Of Attributes: 1
_FillValue : 9.96921e+36
(0) 15931.8

(1) -115.5847
(2) -737.8795
(3) 70.40801
(4) -119.614
(5) 33.65253
(6) 45.77473
(7) -82.48395
(8) -9.278745
(9) -22.75922
(10) 0.8107558
(11) 18.12147
(12) -4.767975
(13) 5.310674
(14) 2.703082

.....

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Oct 27 06:50:12 2013

This archive was generated by hypermail 2.1.8 : Fri Nov 01 2013 - 08:58:14 MDT