Page Comparison

...

1. Spherical harmonic transform available in NCL and pyNGL for lat/lon data.

2. Interpolate EAM native grid output to a (N+1)x2N lat/lon cap grid with TempestRemap’s “highorder” (native SE shape function) algorithm.

3. For cubed sphere grid with resolution NE, take N>=NE*6 (degrees of freedom, pole to pole)

4. Need instantaneous output of (U,V) or (vor,div) and PS on the GLL grid (not the PG2) grid. PS is needed if interpolating to pressure levels

   1. in EAM, U and V may be on the PG2 grid in EAMv2. To get GLL output, use DYN_U and DYN_DYN and V_DYN. VOR and DIV output is always on the GLL grid.

   2. fincl3 = ‘VOR:I','DIV:I','DYN_PS:I' or: ‘U‘DYN_DYNU:I’,'DYN_V_DYN:I','DYN_PS:I'

5. For smoother results, need spectra from ~60 snapshots. 2x per day for 1 month will produce nice results in Aqua planet. What about hourly over 3 days? Are there seasonal effects in Real planet?

6. At NE256, interpolating and computing the spectra for each snapshot takes about 5min and 30GB of memory. At NE1024 this will probably require ~500GB of memory and could take 320min per snapshot.

...