<div dir="ltr"><div><span style="font-size:12.8px">Dear all,</span><br></div><div><span style="font-size:12.8px">I want to remove ENSO index from all SST data and use the residual SST to plot correlation between residual SST and area averaged (15S- 08N, 95E-130E) rainfall index.</span></div><div><span style="font-size:12.8px">My residual SST has only time dimension so I cannot plot the correlation map.I tried to find the residual SST as follows. Is there any other method or how can I plot the correlation map. </span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Thank you for your help.</span><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px">Priyanka.</span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"><br></span></div><div><span style="font-size:12.8px"><br></span></div><div><br></div><div><br></div><div>begin</div><div>;*******************************************************************************</div><div>    ;Seasonal mean precipitation averaged over the domain</div><div>;*****************************************************************************</div><div>f    = "<a href="http://mon.precip.nc">mon.precip.nc</a>"</div><div>fin  = addfile(f, "r")</div><div>time = fin->time</div><div>ymd = cd_calendar(time,2)</div><div>iTim = ind (ymd.ge.19790101.and.ymd.lt.20130101)</div><div>preci = fin->precip(iTim,{8:-15},{95:130})</div><div>;printVarSummary(precip)</div><div>pano = rmMonAnnCycTLL(preci)</div><div>;*******************************************************************</div><div>               ;area averaged rainfall index</div><div>;********************************************************************</div><div>index =dim_avg_Wrap (dim_avg_Wrap (pano(time|:,lat|:,lon|:)))</div><div>;printVarSummary(index)</div><div>xNew = dim_standardize(index, 1)</div><div>xNew!0="time"</div><div>xNew&time = index&time</div><div>printVarSummary(xNew)</div><div>;------------------------------------------------------------------------</div><div>f1 = "HadISST_ssta.nc"<br></div><div>fin1 = addfile(f1, "r")</div><div>time2 = fin1->time</div><div>ymd2 = cd_calendar(time2,2)</div><div>iTim2 = ind (ymd2.ge.19790101.and.ymd2.lt.20130101)</div><div>sst = fin1->ssta(iTim2,:,:)</div><div><br></div><div>;*******************************************************************</div><div>                  ;ENSO index</div><div>;******************************************************************</div><div>NINO = fin1->ssta(iTim2,{-5:5},{-120:-170})</div><div>Nino_index = dim_avg_Wrap (dim_avg_Wrap(NINO(time|:,latitude|:,longitude|:)))</div><div>;printVarSummary(Nino_index)</div><div><br></div><div>;---------------Remove ENSO index from SST;---------------------------</div><div>SST = dim_avg_Wrap (dim_avg_Wrap(sst(time|:,latitude|:,longitude|:)))</div><div>Residula_SST = SST - Nino_index</div><div>Residula_SST!0 = "time"</div><div>printVarSummary(Residula_SST)</div><div>end</div><div><br></div></div>