Title: Observed AMV index Description: | Observed AMV index, defined as the 10-year low-pass-filtered area-weighted average of residual SST anomalies over the North Atlantic (80°W to 0°, 0°-65°N). (detrend process) The residual SST anomaly at each grid point is computed by removing the local component regressed on the global mean SST anomaly. 'REF: https://doi.org/10.1029/2019RG000644' COMMENT: | Observed AMV index, defined as the 10-year low-pass-filtered area-weighted average of residual SST anomalies over the North Atlantic (80°W to 0°, 0°-65°N). (detrend process) The residual SST anomaly at each grid point is computed by removing the local component regressed on the global mean SST anomaly. 'REF: https://doi.org/10.1029/2019RG000644' #Thumbnail: maxML_var_ratio_average_case_thumb.png COMPONENT: mpiom MONTHS: 1,2,3,4,5,6,7,8,9,10,11,12 VN: sst VARIABLE: VN LATBE: 0.,65. LONBE: -80.,0. TSMAINPROC : | varyzll = read_mpiexm_season_YZLL(case, casedir,component,varname,years,months,modeln) varyzll := varyzll(:,0,:,:) ;; lev varyzll := dim_rmvmean_n_Wrap(varyzll,0) gts = mean_with_latbe_lonbe(varyzll,(/-90.,90./),(/0.,360./)) gts!0 = varyzll!0 gts&$gts!0$ = varyzll&$varyzll!0$ ;; detrend regmap = regCoef_n(gts,varyzll,0,0) resualSST = varyzll dims = dimsizes(varyzll) do t = 0, dims(0)-1 resualSST(t,:,:) = varyzll(t,:,:) - (regmap(:,:)*gts(t)) end do ts = mean_with_latbe_lonbe(resualSST,latbe,lonbe) ts!0 = varyzll!0 ts&$ts!0$ = varyzll&$varyzll!0$ delete(varyzll) TSPOSTPROC: | ;; bandpass wgt = filwgts_lanczos(21,0,1./10.,-999,1.) ts := wgt_runave(ts,wgt,0) x = ispan(min(ybe),max(ybe),1) TSRES: | res@vpHeightF = .3 res@trYMinF = -0.4 res@trYMaxF = 0.4 res@tmXBFormat = "f3.0" ;;ts := ts - avg(ts) REGRES: | ;; REGRES res@mpLimitMode = "LatLon" res@mpMinLatF = -90. res@mpMaxLatF = 90. res@mpMinLonF = -330. res@mpMaxLonF = 30. res@pmTickMarkDisplayMode = "Always" res@mpCenterLonF = -150. res@cnLinesOn = False res@cnFillOn = True res@lbOrientation = "Vertical" res@pmLabelBarWidthF = 0.1 res@pmLabelBarOrthogonalPosF = 0.01 res@cnLevelSelectionMode = "ExplicitLevels" res@cnLevels = fspan(-1.2,1.2,13) RERRES: | res@mpLimitMode = "LatLon" res@mpMinLatF = 10. res@mpMaxLatF = 80. res@mpMinLonF = -100. res@mpMaxLonF = 30. res@pmTickMarkDisplayMode = "Always" res@mpCenterLatF = 50. res@mpCenterLonF = -35. res@mpProjection = "LambertConformal" res@mpGridAndLimbOn = True res@gsnMaskLambertConformal = True res@cnLinesOn = False res@cnFillOn = True res@lbOrientation = "Vertical" res@pmLabelBarWidthF = 0.1 res@pmLabelBarOrthogonalPosF = 0.01 res@cnLevelSelectionMode = "ExplicitLevels" res@cnLevels = fspan(-1.2,1.2,13) VRERES: | ;; REGRES res@mpLimitMode = "LatLon" res@mpMinLatF = -90. res@mpMaxLatF = 90. res@mpMinLonF = -330. res@mpMaxLonF = 30. res@pmTickMarkDisplayMode = "Always" res@mpCenterLonF = -150. res@mpOutlineOn = True res@mpLandFillColor = "grey90" res@mpGeophysicalLineColor = "grey60" ;res@cnLinesOn = False ;res@cnFillOn = True ;res@lbOrientation = "Vertical" ;res@pmLabelBarWidthF = 0.1 ;res@pmLabelBarOrthogonalPosF = 0.01 ;res@cnLevelSelectionMode = "ExplicitLevels" ;res@cnLevels = fspan(-1.2,1.2,13) res@vcMinDistanceF = .02 res@vcRefLengthF = .05 res@vcRefMagnitudeF = 3. res@vcGlyphStyle = "CurlyVector" res@vcMonoLineArrowColor = True res@vcLineArrowColor = "black" res@vcLineArrowThicknessF = 1.5 res@vcRefAnnoString1 = res@vcRefMagnitudeF+" (m/s)/K" res@lbBoxSeparatorLinesOn = False res@vcRefAnnoOrthogonalPosF = -0.143 res@vcRefAnnoParallelPosF = 0.176 res@vcRefAnnoSide = "Left" SIGM: | ## apply to VRERES is needed ;; plot signifinicant vector only siglev = .05 utval = onedtond(abs(ureg@tval),dimsizes(ureg)) vtval = onedtond(abs(vreg@tval),dimsizes(vreg)) windtval = onedtond(abs(windreg@tval),dimsizes(windreg)) ;; no windreg now... uP = student_t(utval,dimsizes(years)) vP = student_t(vtval,dimsizes(years)) windP = student_t(windtval,dimsizes(years)) tmask = windP.le.siglev ureg = where(tmask,ureg,ureg@_FillValue) vreg = where(tmask,vreg,vreg@_FillValue) CASEDIR: /projects/NS9207K/pgchiu/CASE ## CASE define later CTN: f08_1AGCM ## proxy of CCASE 'control case' CCASE: CTN CONTROLDIR: /projects/NS9207K/pgchiu/CCASE YBEC: 1101,1250 CONTROLYEARBE: YBEC NYC: 150 CAA: f08_5AGCMs_NASPG_average CNA: IE-average YBEA: 1101,1250 NYA: 150 CAS: f08_5AGCMs_NASPG_sample CNS: CTL2 YBES: 1101,1250 NYS: 150 Scripts: - plotScript: mpiexm/plot_hadisst_obsamv.ncl TITLE: "HadISST observed AMV index" FIGFILENAME: hadisst_obsamv - plotScript: mpiexm/plot_hadisst_obsamv_sst_reg.ncl TITLE: "HadISST ObsAMV index reg." FIGFILENAME: hadisst_obsamv_regmap TSNETCDFFILE: hadisst_obsamv TSVARNAME: ts PLOTRES: REGRES - plotScript: mpiexm/plot_mpiexm_ts.ncl CASE: CTN YEARBE: YBEC PLOTRES: TSRES TITLE: "ObsAMV index CTL (model year: NYC)" FIGFILENAME: CTN_VN_obsamv - plotScript: mpiexm/plot_mpiexm_ts.ncl CASE: CAA YEARBE: YBEA PLOTRES: TSRES TITLE: "ObsAMV index CNA (model year: NYA)" FIGFILENAME: CAA_VN_obsamv - plotScript: mpiexm/plot_mpiexm_ts.ncl CASE: CAS YEARBE: YBES PLOTRES: TSRES TITLE: "ObsAMV index CNS (model year: NYS)" FIGFILENAME: CAS_VN_obsamv - plotScript: mpiexm/plot_mpiexm_regmap.ncl CASE: CTN YEARBE: YBEC TITLE: "ObsAMV index reg. CTL (model year: NYC)" FIGFILENAME: CASE_VN_regmap TSNETCDFFILE: CASE_VN_obsamv TSVARNAME: ts PLOTRES: REGRES - plotScript: mpiexm/plot_mpiexm_regmap.ncl CASE: CAA YEARBE: YBEA TITLE: "ObsAMV index reg. CNA (model year: NYA)" FIGFILENAME: CASE_VN_regmap TSNETCDFFILE: CASE_VN_obsamv TSVARNAME: ts PLOTRES: REGRES - plotScript: mpiexm/plot_mpiexm_regmap.ncl CASE: CAS YEARBE: YBES TITLE: "ObsAMV index reg. CNS (model year: NYS)" FIGFILENAME: CASE_VN_regmap TSNETCDFFILE: CASE_VN_obsamv TSVARNAME: ts PLOTRES: REGRES - plotScript: mpiexm/plot_mpiexm_regmap.ncl CASE: CTN YEARBE: YBEC TITLE: "ObsAMV index reg. CTL (model year: NYC)" FIGFILENAME: CASE_VN_regmap_natl CACHEFILENAME: CASE_VN_regmap TSNETCDFFILE: CASE_VN_obsamv TSVARNAME: ts PLOTRES: RERRES - plotScript: mpiexm/plot_mpiexm_regmap.ncl CASE: CAA YEARBE: YBEA TITLE: "ObsAMV index reg. CNA (model year: NYA)" FIGFILENAME: CASE_VN_regmap_natl CACHEFILENAME: CASE_VN_regmap TSNETCDFFILE: CASE_VN_obsamv TSVARNAME: ts PLOTRES: RERRES - plotScript: mpiexm/plot_mpiexm_regmap.ncl CASE: CAS YEARBE: YBES TITLE: "ObsAMV index reg. CNS (model year: NYS)" FIGFILENAME: CASE_VN_regmap_natl CACHEFILENAME: CASE_VN_regmap TSNETCDFFILE: CASE_VN_obsamv TSVARNAME: ts PLOTRES: RERRES ##- ## plotScript: mpiexm/plot_mpiexm_regmap_vect.ncl ## CASE: CTN ## YEARBE: YBEC ## TITLE: "ObsAMV index 10m wind reg. CTL (model year: NYC)" ## COMPONENT: echam6 ## UVARIABLE: u10 ## VVARIABLE: v10 ## MODELN: 0 ## FIGFILENAME: CASE_VN_sfc_uv_reg ## TSNETCDFFILE: CASE_VN_obsamv ## TSVARNAME: ts ## PLOTRES: VRERES ##- ## plotScript: mpiexm/plot_mpiexm_regmap_vect.ncl ## CASE: CAA ## YEARBE: YBEA ## TITLE: "ObsAMV index 10m wind reg. CNA (model year: NYA)" ## COMPONENT: echam6 ## UVARIABLE: u10 ## VVARIABLE: v10 ## MODELN: "(/1,2,3,4,5/)" ## FIGFILENAME: CASE_VN_sfc_uv_reg ## TSNETCDFFILE: CASE_VN_obsamv ## TSVARNAME: ts ## PLOTRES: VRERES ##- ## plotScript: mpiexm/plot_mpiexm_regmap_vect.ncl ## CASE: CAA ## YEARBE: YBEA ## TITLE: "ObsAMV index 10m wind reg. CNA (model year: NYA)" ## COMPONENT: echam6 ## UVARIABLE: u10 ## VVARIABLE: v10 ## MODELN: 1 ## FIGFILENAME: CASE_VN_sfc_uv_reg_1agcm ## TSNETCDFFILE: CASE_VN_obsamv ## TSVARNAME: ts ## PLOTRES: VRERES ##- ## plotScript: mpiexm/plot_mpiexm_regmap_vect.ncl ## CASE: CAS ## YEARBE: YBES ## TITLE: "ObsAMV index 10m wind reg. CNS (model year: NYS)" ## COMPONENT: echam6 ## UVARIABLE: u10 ## VVARIABLE: v10 ## MODELN: 1 ## FIGFILENAME: CASE_VN_sfc_uv_reg ## TSNETCDFFILE: CASE_VN_obsamv ## TSVARNAME: ts ## PLOTRES: VRERES - plotScript: mk_html_page.sh FIGFILES: hadisst_obsamv hadisst_obsamv_regmap CTN_VN_obsamv CAA_VN_obsamv CAS_VN_obsamv CTN_VN_regmap CAA_VN_regmap CAS_VN_regmap CTN_VN_regmap_natl CAA_VN_regmap_natl CAS_VN_regmap_natl CTN_VN_sfc_uv_reg CAA_VN_sfc_uv_reg CAS_VN_sfc_uv_reg CAA_VN_sfc_uv_reg_1agcm COLUMN: "'HadISST Observed AMV index' tr 2 'Observed AMV index Model' 'tr' CTL CNA CNS 'tr' 3 'Regression map' 'tr' 3 3 3 '-' 1 "