14c14 < sign = -1. --- > sign = 1. 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 = "G_JRA55_test1_vdc013l" > ; datadir = "/work/yiwen89419/CCSM_ncl/G_JRA55_test1_vdc013l/" > ; netcdf = "MOCtot.306-366.eof.atl.lz10+.nodt.nc" 20c20,24 < casename = "taiesm" ; --- > ; casename = "b_20ct_ctrl_1" ;EOF1*(-1) > ; datadir = "/work/yiwen89419/CCSM_ncl/" > ; netcdf = "./Monthly/MOCtot.1921-2000.eof.atl.lz10+.nodt.nc" > > casename = "taiesm" ;EOF1*(-1) 22,25c26,29 < ; netcdf = "./f09.B-PI.tn15.cmip6.j01.re/yearly/MOCtot.500-700.eof.atl.lz10+.nodt.nc" < ; netcdf = "MOCtot.500-700.eof.atl.nc" ;detrend & smooth(x) < netcdf = "MOCtot.500-700.eof.atl.lz10+.nodt.nc" ;not detrend & smooth[default] < ; netcdf = "MOCtot.500-700.eof.atl.lz10+.nc" ;detrend & smooth --- > netcdf = "./f09.B-PI.tn15.cmip6.j01.re/yearly/MOCtot.500-700.eof.atl.lz10+.nodt.nc" > > ; sign = -1. > 28c32 < psout = "Fig2_AMOCeofs" --- > psout = sfx@fBase 30a35 > 40,41c45,46 < eof(0,:,:) = eof(0,:,:)*sign < eofts(0,:) = eofts(0,:)*sign --- > eof = eof*sign > eofts = eofts*sign 43,45c48,49 < mocz = mocz/100. < mocz = mocz/1000. < mocz@units = "km" --- > mocz = mocz/100./1000. > mocz@units = "m" 62,63c66 < ;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!) 70,74c73,75 < gsn_define_colormap(wks,"gsdtol") < blk = NhlNewColor(wks,0.,0.,0.) ; add gray to map < grey = NhlNewColor(wks,0.5,0.5,0.5) ; add gray to map < ltgrey = NhlNewColor(wks,0.8,0.8,0.8) ; add gray to map < wt = NhlNewColor(wks,1.,1.,1.) ; add gray to map --- > ; gsn_define_colormap(wks,"amwg_blueyellowred") ; choose colormap > gsn_define_colormap(wks,"BlueWhiteOrangeRed") ; choose colormap > i = NhlNewColor(wks,0.8,0.8,0.8) ; add gray to map 83c84 < res@cnMissingValFillColor = blk --- > res@cnMissingValFillColor = "grey" 86,87c87,88 < res@cnLinesOn = False ; turn off contour lines < res@cnLineLabelsOn = False ; turn the line labels off --- > res@cnLinesOn = True ; turn off contour lines > res@cnLineLabelsOn = True ; turn the line labels off 89,92c90,91 < res@gsnSpreadColorStart = 32 ; start at color 10 < res@gsnSpreadColorEnd = 8 ; end at color 96 < res@gsnLeftStringFontHeightF = 0.025 < res@gsnRightStringFontHeightF = 0.025 --- > res@gsnSpreadColorStart = 2 ; start at color 10 > res@gsnSpreadColorEnd = 255 ; end at color 96 96,97c95,96 < res@cnMinLevelValF = 0.0 < res@cnMaxLevelValF = 1.0 --- > res@cnMinLevelValF = -1.4 > res@cnMaxLevelValF = 1.4 108,133d106 < res@tmXTOn = False < res@tmYROn = False < res@tmXBLabelsOn = True < res@tmXBMode = "Explicit" < res@tmXBValues = (/-30.,0.,30.,60.,90./) < res@tmXBLabels = (/"30~S~o~N~S","0~S~o~N~N","30~S~o~N~N","60~S~o~N~N","90~S~o~N~N"/) < < lineres = True < lineres@gsnDraw = False < lineres@gsnFrame = False < lineres@cnFillOn = False ; turn on color fill < lineres@cnLinesOn = True ; turn off contour lines < lineres@cnLineLabelFontHeightF = 0.018 < lineres@cnLineLabelsOn = True ; turn the line labels off < lineres@cnLevelSelectionMode = "ManualLevels" ; manually set the contour levels < lineres@cnMinLevelValF = -1.0 < lineres@cnMaxLevelValF = 1.0 < lineres@cnLevelSpacingF = 0.1 < lineres@cnMonoLineThickness = False < lineres@cnMonoLineDashPattern = False < lineres@cnLineThicknesses = (/1,1,1,1,1,1,1,1,1,1,2,1,1,1,1,1,1,1,1,1,1/) < lineres@cnLineDashPatterns = (/1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0/) < lineres@cnInfoLabelOn = False < < < 144d116 < resxy@gsnLeftStringFontHeightF = 0.025 155,156d126 < resxy@tmXTOn = False < resxy@tmYROn = False 164c134 < ;spec@tiYAxisString = "Frequency x Power (Sv~S~2~N~)" --- > spec@tiYAxisString = "Frequency x Power (Sv~S~2~N~)" 170,171d139 < spec@gsnLeftStringFontHeightF = 0.025 < spec@gsnRightStringFontHeightF = 0.025 174c142 < spec@xyLineColors = (/blk,grey,grey,grey/) --- > spec@xyLineColors = (/"black","red","red","red"/) 207,215c175,177 < ; res@gsnLeftString = "ci = "+flt2string(res@cnLevelSpacingF)+" Sv" < res@gsnLeftString = "a. AMOC EOF1" < res@gsnRightString = sprintf("%3.1f", eof@pcvar(0))+"%" < tmp = eof(0,:,:) < tmp = abs(eof(0,:,:)) < plot(0) = gsn_csm_contour(wks,tmp,res) < ploto = gsn_csm_contour(wks,eof(0,:,:),lineres) < overlay (plot(0), ploto) < delete(tmp) --- > res@gsnLeftString = "ci = "+flt2string(res@cnLevelSpacingF)+" Sv" > res@gsnRightString = "EOF 1, "+sprintf("%4.2f",eof@pcvar(0) ) + "%" > plot(0) = gsn_csm_contour(wks,eof(0,:,:),res) 217,218c179 < resxy@gsnLeftString = "b. PC1 Time Series" < ; resxy@gsnCenterString = "" --- > resxy@gsnCenterString = "PC 1 time series" 222,224c183,184 < spectrum = specx_ci(dof,0.95,0.) < spec@gsnLeftString = "c. PC1 Power Spectrum" < spec@gsnCenterString = "" --- > spectrum = specx_ci(dof,0.95,0.99) > spec@gsnCenterString = "Power Spectrum" 255,265c215,217 < res@gsnLeftString = "d. AMOC EOF2" < res@gsnRightString = sprintf("%3.1f", eof@pcvar(1))+"%" < res@tiXAxisString = "Latitude" < tmp = eof(1,:,:) < tmp = abs(eof(1,:,:)) < plot(3) = gsn_csm_contour(wks,tmp,res) < ploto3 = gsn_csm_contour(wks,eof(1,:,:),lineres) < delete(tmp) < overlay (plot(3), ploto3) < resxy@gsnLeftString = "e. PC2 Time Series" < resxy@tiXAxisString = "Time (years)" --- > res@gsnRightString = "EOF 2, "+sprintf("%4.2f",eof@pcvar(1) ) + "%" > plot(3) = gsn_csm_contour(wks,eof(1,:,:),res) > resxy@gsnCenterString = "PC 2 time series" 269,270c221,222 < spectrum = specx_ci(dof,0.95,0.) < spec@gsnLeftString = "f. PC2 Power Spectrum" --- > spectrum = specx_ci(dof,0.95,0.99) > spec@gsnCenterString = "" 301a254 > pres@txString = casename 307d259 < pres@gsnPanelRight = 0.95