16,18c16,18
< ; casename = "b40.1850.track1.1deg.006"
< ; datadir = "/ccsm/ocn/yeager/amoc_ts/"+casename+"/eofs/"
< ; netcdf = "MOC.700-1299.eof.atl.nc"
---
> ;  casename = "b_20ct_ctrl_1"
> ;  datadir = "/work/yiwen89419/CESM_ncl/"
> ;  netcdf = "./Monthly/MOCtot.1921-2000.eof.atl.lz10+.nodt.nc"
24c24
<   fin= datadir+netcdf
---
>   fin= netcdf
61,62c61,62
< ;jave = 5 ; number of periodograms to average (if < 3, no smoothing!)
< jave = 1 ; number of periodograms to average (if < 3, no smoothing!)
---
> jave = 5 ; number of periodograms to average (if < 3, no smoothing!)
> ;jave = 1 ; number of periodograms to average (if < 3, no smoothing!)
163c163
< ;spec@tiYAxisString = "Frequency x Power (Sv~S~2~N~)"
---
> ;spec@tiYAxisString = "Power (Sv~S~2~N~)"
178c178
< spec@xyXStyle = "Log"
---
> ;spec@xyXStyle = "Log"
185c185
<  spec@trXMaxF = 0.5
---
>  spec@trXMaxF = 0.05
187c187
<  spec@trYMaxF = 1.
---
>  spec@trYMaxF = 30.
224c224,225
<   plot(2)=gsn_csm_xy(wks,dof@frq,spectrum*conform(spectrum,dof@frq,(/1/)),spec)
---
> ; plot(2)=gsn_csm_xy(wks,dof@frq,spectrum*conform(spectrum,dof@frq,(/1/)),spec)
>   plot(2)=gsn_csm_xy(wks,dof@frq,spectrum,spec)
271c272
<   plot(5)=gsn_csm_xy(wks,dof@frq,spectrum*conform(spectrum,dof@frq,(/1/)),spec)
---
>   plot(5)=gsn_csm_xy(wks,dof@frq,spectrum,spec)