<html><head></head><body><div style="font-family: Verdana;font-size: 12.0px;"><div>Hi all!</div>

<div>&nbsp;</div>

<div>I am trying to make as scatter plot. I have monthly data, the dimensions of the variables x and y are time,lat,lon.</div>

<div>I would like to have scatter plots for different regions and to have monthly data&nbsp;de-seasonalized.</div>

<div>The script is working, but it seems like I have annual means.</div>

<div>Here is my script:</div>

<div>&nbsp;</div>

<div>
<div>in_x = addfile(&quot;~/DATA/x.nc&quot;,&quot;r&quot;)<br/>
x = in_x-&gt;x(:,{-10:-20},{180:190})</div>

<div>in_y = addfile(&quot;~/DATA/y.nc&quot;,&quot;r&quot;)<br/>
y = in_y-&gt;y(:,{-10:-20},{180:190})</div>

<div>&nbsp;</div>
;*********************************************************************************

<div>; remove annual cycle<br/>
;*********************************************************************************</div>

<div>x_ave = rmMonAnnCycTLL(x)<br/>
y = rmMonAnnCycTLL(y)</div>

<div>&nbsp;</div>

<div>;********************************************************************************<br/>
; calculate average LON LAT<br/>
;********************************************************************************</div>

<div>x_ave_new = dim_avg_n(x_ave,(/1,2/))</div>

<div>y_ave_new = dim_avg_n(y_ave,(/1,2/))</div>

<div>&nbsp;</div>

<div>;********************************************************************************<br/>
; regression line<br/>
;********************************************************************************</div>

<div>rc = regline(x_ave_new,y_ave_new)</div>

<div>print(rc)</div>

<div>&nbsp;</div>

<div>;********************************************************************************<br/>
; create array<br/>
;********************************************************************************</div>

<div>x = x_ave_new<br/>
y = y_ave_new</div>

<div>data1 = new((/2,dimsizes(y)/), typeof(y))</div>

<div><br/>
data1(0,:) = y</div>

<div>data1(1,:) = rc*(x-rc@xave) + rc@yave</div>

<div>&nbsp;</div>

<div>&nbsp;</div>

<div><br/>
;*********************************************************************************<br/>
; plotting parameters<br/>
;*********************************************************************************</div>

<div>wks = gsn_open_wks(&quot;pdf&quot;,&quot;scatter&quot;)</div>

<div>res &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= True<br/>
;res@gsnMinimize&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;= True</div>

<div>...</div>

<div>&nbsp;</div>

<div>plot = gsn_csm_xy (wks,x,data1,res)</div>

<div>&nbsp;</div>

<div>end</div>

<div>;********************************************************************</div>

<div>&nbsp;</div>

<div>&nbsp;</div>

<div>Thanks.</div>

<div>Best regards,</div>

<div>Lena</div>
</div></div></body></html>