# Re: Problem with shaec

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Sun Oct 27 2013 - 08:20:15 MDT

?? Not sure what is meant by 'but "a" is ​​greater than I expected'

Certainly 'a' having an average value of 15931.8 seems reasonable
for an overall average for "global 500hpa Geopotential height"

On 10/27/13 6:49 AM, Wood wrote:
> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Oct 27 08:20:18 2013

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