;Id = NCDF_OPEN('/aim/d/waccm/smedyear/03/wa319_2x_smedyear.cam2.h1.1951-03-21-00000.nc', /NOWRITE)

Id = NCDF_OPEN('/net/aim/d/waccm/sminyear/03/wa319_2x_sminyear.cam2.h1.1951-03_MeanTZ3.nc', /NOWRITE)
IdUV = NCDF_OPEN('/net/aim/d/waccm/sminyear/03/wa319_2x_sminyear.cam2.h1.1951-03_MeanUV.nc', /NOWRITE)
IdWAM = NCDF_OPEN('/net/aim/d/waccm/wa319_TIPHYS/WAM/wam_t_mar.nc', /NOWRITE)
IdUWAM = NCDF_OPEN('/net/aim/d/waccm/wa319_TIPHYS/WAM/wam_u_mar.nc', /NOWRITE)
IdVWAM = NCDF_OPEN('/net/aim/d/waccm/wa319_TIPHYS/WAM/wam_v_mar.nc', /NOWRITE)
;Id = NCDF_OPEN('/aim/d/waccm/sminyear/03/wa319_2x_sminyear.cam2.h1.1951-03_MeanUV.nc', /NOWRITE)
;Id = NCDF_OPEN('/aim/d/waccm/wa319_TIPHYS/MTM/WAX3548GWSMnQUP50_2x.cam2.h1.1954_MarchMeanTZ3.nc', /NOWRITE)
;Id = NCDF_OPEN('/aim/d/waccm/wa319_TIPHYS/MTM/wa319_2x_sminyear.cam2.h1.1951-06_MeanTZ3.nc', /NOWRITE)
;Id = NCDF_OPEN('/aim/d/waccm/sminyear/03/wa319_2x_sminyear.cam2.h1.1951-03-11-00000.nc', /NOWRITE)
;Id = NCDF_OPEN('/aim/d/waccm/sminyear/03/wa319_2x_sminyear.cam2.h1.1951-03-27-00000.nc', /NOWRITE)
;Id = NCDF_OPEN('/aim/d/waccm/wa319_TIPHYS/MTM/WAX3548GWSMnQUP50_2x.cam2.h1.1954-03-11-00000.nc', /NOWRITE)
;Id = NCDF_OPEN('/aim/d/waccm/wa319_TIPHYS/MTM/WAX3548GWSMnQUP50_2x.cam2.h1.1954-03-27-00000.nc', /NOWRITE)

;; Coordinates

ncdf_varget, Id, 'lon', lon         ; longitude
ncdf_varget, Id, 'lat', lat         ; latitude
ncdf_varget, Id, 'lev', zplev       ; log pressure, -ln(5.e-7/p) (p in hPa)
ncdf_varget, Id, 'time', time_min   ; time in minites

ncdf_varget, IdWAM, 'lon', lonWAM         ; longitude
ncdf_varget, IdWAM, 'lat', latWAM         ; latitude
ncdf_varget, IdWAM, 'lev', zplevWAM       ; log pressure, -ln(5.e-7/p) (p in hPa)
ncdf_varget, IdWAM, 'time', time_minWAM   ; time in minites

;; Field variables

;ncdf_varget, Id, 'Z3', zg     ;geopotential height
;zgm = zg * (1. + zg/6370000.0)

nx = n_elements(lon)
ny = n_elements(lat)
nz = n_elements(zplev)
ntime = n_elements(time_min)

nxWAM = n_elements(lonWAM)
nyWAM = n_elements(latWAM)
nzWAM = n_elements(zplevWAM)
ntimeWAM = n_elements(time_minWAM)

; Dimensions lon,lat,lev,time

ncdf_varget, Id, 'T', tn      ;neutral temperature

ncdf_varget, IdUV, 'U', un      ;neutral zonal wind

ncdf_varget, IdUV, 'V', vn     ;neutral meridional wind

;ncdf_varget, Id, 'OMEGA', om     ;neutral vertical wind

; Dimensions lon,lat,lev,time

ncdf_varget, IdWAM, 'temp', tnWAM      ;neutral temperature
ncdf_varget, IdUWAM, 'u', unWAM      ;neutral zonal wind
ncdf_varget, IdVWAM, 'v', vnWAM      ;neutral zonal wind

NCDF_CLOSE, Id
NCDF_CLOSE, IdUV

NCDF_CLOSE, IdWAM
NCDF_CLOSE, IdUWAM
NCDF_CLOSE, IdVWAM

ntime2=ntime/2+1
nshift=(nx+1)/2-1
frq=findgen(ntime2)
xwvn=-findgen(nx)+nshift

;zonalfreq0,tn,ntime,nx,ny,nz,nshift,ntime2,uspc,uphs
;zonalfreq0,un,ntime,nx,ny,nz,nshift,ntime2,uspc,uphs
;zonalfreq0,vn,ntime,nx,ny,nz,nshift,ntime2,uspc,uphs
;
;  The following are for the reconstruction and don't apply to amplitude plots made here of uspc only ucomp
;
nfrq = -1
nwv = 1 ;migrating diurnal wavenumber 1

;zonalfreq1,tn,ntime,nx,ny,nz,nshift,ntime2,uspc,nwv,nfrq,tcompd
;zonalfreq1,un,ntime,nx,ny,nz,nshift,ntime2,uspc,nwv,nfrq,tcompd
;zonalfreq1,vn,ntime,nx,ny,nz,nshift,ntime2,uspc,nwv,nfrq,tcompd

ntime2WAM=ntimeWAM/2+1
nshiftWAM=(nxWAM+1)/2-1
frqWAM=findgen(ntime2WAM)
xwvnWAM=-findgen(nxWAM)+nshiftWAM
;
;  The following are for the reconstruction and don't apply to amplitude plots made here of uspc only ucomp
;
nfrq = -1
nwv = 1 ;migrating diurnal wavenumber 1

;zonalfreq0,tnWAM,ntimeWAM,nxWAM,nyWAM,nzWAM,nshiftWAM,ntime2WAM,uspcWAM,uphsWAM
;zonalfreq0,unWAM,ntimeWAM,nxWAM,nyWAM,nzWAM,nshiftWAM,ntime2WAM,uspcWAM,uphsWAM
;zonalfreq0,vnWAM,ntimeWAM,nxWAM,nyWAM,nzWAM,nshiftWAM,ntime2WAM,uspcWAM,uphsWAM

;zonalfreq1,tnWAM,ntimeWAM,nxWAM,nyWAM,nzWAM,nshiftWAM,ntime2WAM,uspcWAM,nwv,nfrq,tcompdWAM
;zonalfreq1,unWAM,ntimeWAM,nxWAM,nyWAM,nzWAM,nshiftWAM,ntime2WAM,uspcWAM,nwv,nfrq,tcompdWAM
;zonalfreq1,vnWAM,ntimeWAM,nxWAM,nyWAM,nzWAM,nshiftWAM,ntime2WAM,uspcWAM,nwv,nfrq,tcompdWAM

;plot_waccm_wam_t_prof_eq, zplev, zplevWAM, tn, tnWAM, un, unWAM, vn, vnWAM

tnZM = FLTARR(ny,nz,ntime)
unZM = FLTARR(ny,nz,ntime)
vnZM = FLTARR(ny,nz,ntime)

FOR iTime = 0, ntime-1 DO BEGIN
  FOR iLev = 0, nz-1 DO BEGIN
    FOR iLat = 0, ny-1 DO BEGIN

     tnZM(iLat,iLev,iTime) = MEAN(tn(*,iLat,iLev,iTime))
     unZM(iLat,iLev,iTime) = MEAN(un(*,iLat,iLev,iTime))
     vnZM(iLat,iLev,iTime) = MEAN(vn(*,iLat,iLev,iTime))

    ENDFOR
  ENDFOR
ENDFOR

tnZMDev = FLTARR(nx,ny,nz,ntime)
unZMDev = FLTARR(nx,ny,nz,ntime)
vnZMDev = FLTARR(nx,ny,nz,ntime)

FOR iTime = 0, ntime-1 DO BEGIN
  FOR iLev = 0, nz-1 DO BEGIN
    FOR iLat = 0, ny-1 DO BEGIN

     tnZMDev(*,iLat,iLev,iTime) = tn(*,iLat,iLev,iTime) - tnZM(iLat,iLev,iTime)
     unZMDev(*,iLat,iLev,iTime) = un(*,iLat,iLev,iTime) - unZM(iLat,iLev,iTime)
     vnZMDev(*,iLat,iLev,iTime) = vn(*,iLat,iLev,iTime) - vnZM(iLat,iLev,iTime)

    ENDFOR
  ENDFOR
ENDFOR

tnWAMZM = FLTARR(nyWAM,nzWAM,ntimeWAM)
unWAMZM = FLTARR(nyWAM,nzWAM,ntimeWAM)
vnWAMZM = FLTARR(nyWAM,nzWAM,ntimeWAM)

FOR iTime = 0, ntimeWAM-1 DO BEGIN
  FOR iLev = 0, nzWAM-1 DO BEGIN
    FOR iLat = 0, nyWAM-1 DO BEGIN

     tnWAMZM(iLat,iLev,iTime) = MEAN(tnWAM(*,iLat,iLev,iTime))
     unWAMZM(iLat,iLev,iTime) = MEAN(unWAM(*,iLat,iLev,iTime))
     vnWAMZM(iLat,iLev,iTime) = MEAN(vnWAM(*,iLat,iLev,iTime))

    ENDFOR
  ENDFOR
ENDFOR

tnWAMZMDev = FLTARR(nxWAM,nyWAM,nzWAM,ntimeWAM)
unWAMZMDev = FLTARR(nxWAM,nyWAM,nzWAM,ntimeWAM)
vnWAMZMDev = FLTARR(nxWAM,nyWAM,nzWAM,ntimeWAM)

FOR iTime = 0, ntimeWAM-1 DO BEGIN
  FOR iLev = 0, nzWAM-1 DO BEGIN
    FOR iLat = 0, nyWAM-1 DO BEGIN

     tnWAMZMDev(*,iLat,iLev,iTime) = tnWAM(*,iLat,iLev,iTime) - tnWAMZM(iLat,iLev,iTime)
     unWAMZMDev(*,iLat,iLev,iTime) = unWAM(*,iLat,iLev,iTime) - unWAMZM(iLat,iLev,iTime)
     vnWAMZMDev(*,iLat,iLev,iTime) = vnWAM(*,iLat,iLev,iTime) - vnWAMZM(iLat,iLev,iTime)

    ENDFOR
  ENDFOR
ENDFOR

plot_waccm_wam_t_prof_zmdev, zplev, zplevWAM, tnZMDev, tnWAMZMDev, unZMDev, unWAMZMDev, vnZMDev, vnWAMZMDev

;plot_waccm_wam_t_prof_phase, zplev, zplevWAM, uphs, uphsWAM
;plot_waccm_wam_u_prof_phase, zplev, zplevWAM, uphs, uphsWAM
;plot_waccm_wam_v_prof_phase, zplev, zplevWAM, uphs, uphsWAM

;plot_waccm_wam_t_prof, zplev, zplevWAM, uspc, uspcWAM
;plot_waccm_wam_u_prof, zplev, zplevWAM, uspc, uspcWAM
;plot_waccm_wam_v_prof, zplev, zplevWAM, uspc, uspcWAM

;plotting_waccm, tn, lat, zplev, uspc, tcompd

;plotting_waccm_u, un, lat, zplev, uspc, tcompd

;plotting_waccm_v, vn, lat, zplev, uspc, tcompd
 
;plot_waccm_t_prof, tn, zplev, uspc
;plot_waccm_u_prof, un, zplev, uspc
;plot_waccm_v_prof, un, zplev, uspc
 
;plot_recon_waccm,tcompsd,lat,lon,zplev,time_min,nx,ny,nz,ntime

;plotting_recon, tn, lat, zplev, uspc, tcompd_sd_td_qd, uwvtotal

;plot_recon_lon, tcompd, lon, lev

;zonalwvn,tn,nx,ny,nz,uspc4,uphs4,uwv4

;plotting_wvn, tn, uwv1UT, lon

;plot_utave_lonave_t, uwvtotal_7to12, ntime, nx, ny, nz,lat


;plot_utave_lonave_z_t, uphs, 5, nx, ny, nz, zgm

END