mehwish ramzan mehwish.ramzan at gmail.com
Sat Nov 18 08:09:57 EST 2017

```Dear GrADS users,

I need help in fixing the following problem.

I am trying to create an inter-annual time series for present and future
climate. Now, when I run the script it gives me overlapping time series.
What is want to obtain is the distinction between present and future
precipitation change.

Also, how can I obtain a linear regression within the same graph?Kindly
enlighten me.

Please see the script and an image as attached files.

Thank you

With Best Regards,

Mehwish
-------------- next part --------------
An HTML attachment was scrubbed...
-------------- next part --------------
{\rtf1\ansi\ansicpg1252\cocoartf1265\cocoasubrtf210
{\fonttbl\f0\fswiss\fcharset0 Helvetica;}
{\colortbl;\red255\green255\blue255;}
\margl1440\margr1440\vieww25400\viewh14520\viewkind0
\pard\tx720\tx1440\tx2160\tx2880\tx3600\tx4320\tx5040\tx5760\tx6480\tx7200\tx7920\tx8640\pardirnatural

\f0\fs24 \cf0 'reinit'\
'open rcp85.ctl'\
'open rcp45.ctl'\
'open 20C.ctl'\
\
'set display color white'\
'c'\
\
'set vrange 0 7'\
'set xaxis 1980 2100'\
\
'set dfile 1'\
'set x 1'\
'set y 1'\
'set t 1 972'\
'define p1 = aave(pratesfc*60*60*24, lon=12, lon=123, lat=-15, lat=45)'\
'set cmark 0'\
'set ccolor 1'\
'set cstyle 1'\
'set line 1 1 5'\
'd p1'\
\
'set dfile 2'\
'set x 1'\
'set y 1'\
'set t 1 972'\
'define p2 = aave(pratesfc*60*60*24, lon=12, lon=123, lat=-15, lat=45)'\
'set cmark 0'\
'set ccolor 3'\
'set cstyle 1'\
'set line 3 1 5'\
'd p2'\
\
'set dfile 3'\
'set cmark 0'\
'set ccolor 4'\
'set x 1'\
'set y 1'\
'set t 1 312'\
'define p3= aave(prate*60*60*24, lon=12, lon=123, lat=-15, lat=45)'\
'set style 1'\
'set line 4 1 5'\
'd p\'923\
'cbar_line -x 2.5 -y 2.5 -c 1 3 4 -m 0 -t "RCP8.5" "RCP4.5" "20C" -p'\
'draw title Inter-annual time series\'92\
\
'gxprint inter-annual_prec_21C_20C.png'\
}
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inter-annual_timeseries.png
Type: image/png
Size: 152538 bytes
Desc: not available