Keywords:Phase velocity spectral analysis, Airglow, Gravity wave
Airglow imaging observation, among other ground-based techniques, has been proven to be very effective to study energy and propagation characteristics of atmospheric gravity waves (AGWs). However, the lack of sophisticated analysis methods prevented quantitative studies using huge amount data that have been collected over long period of observation. Matsuda et al.  developed a new phase velocity spectrum analysis method based on 3D FFT algorithm (M-transform) to address this issue. This algorithm can efficiently deal with extensive amounts of imaging data obtained on different years and at various observation sites without bias and treat dynamical/physical effect of AGWs by precisely reflecting amplitude, area and lifetime of each AGW event. Based on Matsuda et al.  method, we have developed a simple and user-friendly function to be used in IDL. The input of this program is a three dimensional array of a time series of 2-D image array where the wave parameters (e.g. horizontal wavelength (λh), wave period (τ), phase speed (c), image resolution in space (dx, dy) and time (dt)) can be customized by the user in one-line command to execute the program. Various simulations by using artificial wave images as test data have been done in order to exhibit the new program’s performance and the characteristics/interpretation of the spectral analysis. The simulation was done by using test data with different time resolution (dt=1 min, 3 min, 5 min and 7 min), image size (256x256, 512x512), duration (30 min, 60 min, 90 min, 120 min) and by changing the horizontal wavelength and wave periods with a fixed phase speed (40 m/s). In addition, by using observational data from Syowa station we conducted an investigation how waves with different horizontal wavelengths behave independently. We divided the horizontal wavelength into three categories: 5< λh <20 km, 20< λh <100 km and 5< λh <100 km. We found that propagation direction for small scale, close to ripple scales (5< λh <20 km), showed omni-directional feature, as previously suggested by Nakamura et al. . The presentation also aims to guide how to use the program within a practical computation time.