*** JAGUAR-T639L340 2018/19 SSW hindcast Data *** The coordinate values of the JAGUAR-T639L340 are described in 'zahyo_high.f90'. The potential temperature coordinate is described in 'zahyo_theta.f90'. All variables are formatted as used to draw the figures in the paper to reduce data size: zonal- and/or 4-day mean and extracted for the latitude regions in the figures. Missing values are filed by -999. File format : binary The descriptions of the data in the files are as follows. ./ 1) ZM_T_tz_70N-80N.bin : (Fig 1) Zonal-mean T (unit:K) averaged over 70°N-80°N 2) ZM_u_4d_NH_50N-70N.bin : (Fig 1) Zonal-mean u (unit:m/s) averaged over 50°N-70°N 3) EPF_(all,PW,MW,GW)_tz_50N-70N.bin : (Fig 2) EPFD (unit:m/s/d) of (all wave,PW,MW,GW) components averaged over 50°N-70°N 4) ZM_T_4d_NH.bin : (Fig 3&5) Zonal- & 4-day mean T (unit:K) at 20°S-90°N 5) ZM_u_4d_NH.bin : (Fig 3&5) Zonal- & 4-day mean u (unit:m/s) at 20°S-90°N 6) ZM_MPV_4d_NH.bin : (Fig 3&5) Zonal- & 4-day mean MPV (unit:PVU) in the theta coordinate at 20°S-90°N 7) EPF_(all,PW,GW)_4d_NH.bin : (Fig 3&5) Zonal- & 4-day mean Eliassen-Palm flux (unit:m^2/s^2) and its divergence (unit:m/s/d) of (all wave,PW,GW) components at 20°S-90°N. var1:EPFD, var2:Fy/(ρ0*a), var3:Fz/(ρ0*a) 8) (GPH,u,MPV,N2,duwdz,dwdz,w,3D-flux-M)_NH-map_*(hPa,K)_*.bin : (Fig 4,7,8&9) 4-day mean [GPH(unit:m),u(m/s),MPV(PVU),N2(/s^2),-du'w'/dz(m/s/d) of GWs,dw*/dz(/s),w*(cm/s),-(∇・F_1)/ρ0(m/s/d)] in the Northern Hemisphere (90°N-EQ) at each pressure/potential temperature 9) (u,3D-flux-W)_4d_xz_60N-70N_20190106.bin : (Fig 6) 4-day mean [u(unit:m/s),3D-flux-W(var1:Fx/ρ0,var2:Fz/ρ0,m^2/s^2)] averaged over 60°N-70°N during 6-9 Jan 2019 10) T_tz_70N.bin : (Fig 8) Vertical profiles of T (unit:K) at 8 stations located at a latitude of 70°N 11) GPH_xt_60N-70N_80km_20181217-20190117.bin : (Fig 10a) Hourly data of GPH (unit:m) averaged over 60°N-70°N at 80km from 17 Dec 2018-17 Jan 2019 12) (GPH,3D-flux-W)_xz_(65N,70N)_*.bin : (Fig 10b-d) Daily mean [GPH(unit:m),3D-flux-W(var1:Fx/ρ0,var2:Fz/ρ0,m^2/s^2)] at a latitude of (65°N,70°N) 13) uw_cz_*_*.bin : (Fig 11) u'w' (unit:m^2/s^2) of GWs at each latitude & ground-based phase velocity (c) ***-----------------------------*** *** SAMPLE OF DATA READING(F90) *** ***-----------------------------*** Files 1-3: do itime = 1, 24*44 !-- hourly data for 44 days read(ifile) (var(ilev),ilev=1,nlev) !-- nlev=340 end do Files 4-7: do itime = 1, 11 !-- 4 days * 11 periods do ilev = 1, nlev(or ntheta in File 6) !-- nlev=340, ntheta=138 read(ifile) (var1(ilat),ilat=nlat1,nlat2) !-- nlat1=1(90°N), nlat2=587(20°S) ( read(ifile) (var2(ilat),ilat=nlat1,nlat2) ) ( ... ) end do end do Files 8: do ilon = 1, nlon !-- nlon=1920 read(ifile) (var(ilat),ilat=1,nlat/2) !-- nlat=960 end do Files 9: do ilev = 1, nlev !-- nlev=340 read(ifile) (var1(ilon),ilon=1,nlon) !-- nlon=1920 ( read(ifile) (var2(ilon),ilon=1,nlon) ) ( ... ) end do File 10: do istation = 1, 8 read(ifile) longitude do itime = 1, 24*44 !-- hourly data for 44 days read(ifile) (var(ilev),ilev=1,nlev) end do end do File 11: do itime = 1, 24*32 !-- hourly data for 32 days read(ifile) (var(ilon),ilon=1,nlon) end do Files 12: do itime = 1, 4 !-- daily data for 4 days do ilev = 1, nlev !-- nlev=340 read(ifile) (var1(ilon),ilon=1,nlon) !-- nlon=1920 ( read(ifile) (var2(ilon),ilon=1,nlon) ) ( ... ) end do end do Files 13: read(ifile) (c(ic),ic=-nc,nc) !-- c: ground-based phase velocity, nc=100 do ilev = 1, nlev read(ifile) (var(ic,ilev),ic=-nc,nc) end do