*** Processed Data Used in the Following Paper ***
! Okui, H., Sato, K., & Watanabe, S., Contribution of gravity waves to universal vertical wavenumber (m^(-3)) spectra revealed by a gravity-wave-permitting general circulation model. submitted to Journal of Geophysical Research: Atmosphere.
! Corresponding author: Okui, H. (okui@eps.s.u-tokyo.ac.jp)
*** ------------------------------------------ ***

Coordinates: zahyo_high.f90

*** ------------------------------------------ ***
Figures 1a-f: u', v', & T'g/ToN spectra at Shigaraki, Japan
*** ------------------------------------------ ***
!** Spectral density **!
./{u,v,T}-GW_PSDm_lat35lon136
./{u,v,T}-ALL_PSDm_lat35lon136
!-- How to read --!
read(10) range_str,(PSD_str(i),i=1,26)
read(10) range_mes,(PSD_mes(i),i=1,26)
! real :: range_str !! Height range for calculation in the stratosphere [m]
! real :: range_mes !! Height range for calculation in the mesosphere [m]
! real :: PSD_str(26) !! Spectral density in the stratosphere [(m/s)^2/(cyc/m)]
! real :: PSD_mes(26) !! Spectral density in the mesosphere [(m/s)^2/(cyc/m)]
! m(1:26)=(/(i/range,i=0,25)/)

!** N2 values for theoretical spectra [Smith et al. 1987] **!
./N2_lat35lon136
!-- How to read --!
read(10) rN2_str
read(10) rN2_mes
! real :: rN2_str !! N2 in the stratosphere [s^-2]
! real :: rN2_mes !! N2 in the mesosphere [s^-2]

*** ------------------------------------------ ***
Figure 1g: Zonal-mean v' spectra
*** ------------------------------------------ ***
!** Spectral density **!
./v-GW_PSDm_lat
./v-ALL_PSDm_lat
!-- How to read --!
do ilat = 1,18  ! (85N,75N,...,75S,85S)
  read(10) range,(PSD(i),i=1,26)
end do
! real :: range !! Height range for calculation in the stratosphere [m]
! real :: PSD(26) !! Spectral density in the stratosphere [(m/s)^2/(cyc/m)]
! m(1:26)=(/(i/range,i=0,25)/)

!** N2 values for theoretical spectra [Smith et al. 1987] **!
./N2_lat
!-- How to read --!
read(10) (rN2(ilat),ilat=1,18)

*** ------------------------------------------ ***
Figure 2: Zonal-mean mg*, t, & F0
*** ------------------------------------------ ***
!** mg*, t, & F0 **!
./{u,v,T}_zonal-mean-mgstar
./{u,v,T}_zonal-mean-t
./{u,v,T}_zonal-mean-F0
!-- How to read --!
do ialt = 1, 4 ! z=(18-33,33-45,45-60,60-90)km
  read(10) (data(ialt,ilat),ilat=1,nlat/4) !ilat:(lat(i*4-3)+lat(i*4-2))/2
end do

*** ------------------------------------------ ***
Figures 3, 4: Global mg*, t, & F0
*** ------------------------------------------ ***
!** mg*, t, & F0 **!
./mgstar-t-F0_xy_{18-33,33-45,45-60,60-90}km
!-- How to read --!
do i = 1, nlon/4
  read(10) (data_m(i,j),j=1,nlat/4)
  read(10) (data_t(i,j),j=1,nlat/4)
  read(10) (data_F(i,j),j=1,nlat/4)
end do

!** Zonal wind **!
./U_xy_{18-33,33-45,45-60,60-90}km
!-- How to read --!
do i = 1, nlon/16
  read(10) (data(i,j),j=1,nlat/16)
end do

*** ------------------------------------------ ***
Figure 5: u' spectra & zonal-mean U
*** ------------------------------------------ ***
!** Model-simulated spectra **!
./u-GW_PSDm_50-60{N,S}
!-- How to read --!
do ialt = 1, 4 ! z=(18-33,33-45,45-60,60-90)km
  read(10) range,(PSD(i),i=1,25)
end do
! real :: range !! Height range for calculation in the stratosphere [m]
! real :: PSD(25) !! Spectral density in the stratosphere [(m/s)^2/(cyc/m)]
! m(1:26)=(/(i/range,i=1,25)/)

!** N2 values for theoretical spectra [Smith et al. 1987] **!
./N2_fig5
!-- How to read --!
read(10) (rN2(ialt),ialt=1,4)

!** U profile for shear-effect calculation **!
./U-profile_50-60{N,S}
!-- How to read --!
read(10) (U(ialt),ialt=1,nlev)

*** ------------------------------------------ ***
Figure 6: GW-KE
*** ------------------------------------------ ***
./GW-KE_yz
!-- How to read --!
do j = 1, nlat/3
  read(10) (data(j,k),k=1,nlev)
end do