Dear Rupak,<br>you're wellcome.<br>Another suggestion to you is to open the files independently, even the timescales is not the same. Then, if you whish to analyze the difference between 2070-2098 and 1961-90 period (e.g. annual mean), you can proceed in the following way:<br>
1) open the files indenpendently, and find the correspondent T range that correspond two wich period;<br>2) open both files in the script, and from now on, you need to especify the correspondent time when you call the varible. For example:<br>
<br>'open file1.ctl' (file from 1961 to 1990, with 360 time steps)<br>
'open file2.ctl' (file from 2011 to 2098, with 1056 time steps)<br><br>'define PRES=var.1*0'<br>'define FUT=var.2*0'<br><br>count=1<br>while(count<=360)<br>'define PRES=PRES+var.1(t='count')/360'<br>
count=count+1<br>endwhile<br><br>count=697 (the value that you have found in the step 1)<br>
while(count<=1056)<br>
'define FUT=FUT+var.2(t='count')/360'<br>
count=count+1<br>
endwhile<br><br>'define Change=FUT-PRES'<br><br>I hope this can help you.<br><br>Regards,<br>Roger Torres<br><br><br><div class="gmail_quote">2011/9/20 Rupak Rajbhandari <span dir="ltr"><<a href="mailto:rupak.rajbhandari@gmail.com">rupak.rajbhandari@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">Dear Roger,<div><br></div><div>Thanks. You were right. The first file was for base period run from 1961 to 1990 and the second file was for future run from 2011 to 2098. I was trying to see the shifts in rainfall pattern, if there are any. The files are one file for one month and controlled using templates. Since there were some gap on the data (1991-2010), I kept the files in separate folders - hence two ctl files. When I try to view the second file, instead of simply assigning the T dimension value from 1 to 360 (or whatever may be), it extends T dimension backward for second file from 2011 to 1961, which is -500+, and the data obviously not available. Now that I have kept all the files in single folder is working (should work).</div>
<div><br></div><div>Once again thanks.</div><div><br></div><div>-rupak</div><div><br></div><font color="#888888"><div>Rupak Rajbhandari</div><div>Department of Meteorology</div><div>Trichandra Campus, TC/TU</div><div>Kathmandu</div>
</font><div><div></div><div class="h5"><div><br><br>
<div class="gmail_quote">On Mon, Sep 19, 2011 at 5:16 PM, Roger Rodrigues Torres <span dir="ltr"><<a href="mailto:roger.torres@cptec.inpe.br" target="_blank">roger.torres@cptec.inpe.br</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Rupak,<br>one suggestion is to check the declaration of the time in both ctl's. They need to be exactly the same to work in the way that you wish.<br>Regards,<br>Roger Torres<br><br><div class="gmail_quote"><div>
<div></div><div>2011/9/19 Rupak Rajbhandari <span dir="ltr"><<a href="mailto:rupak.rajbhandari@gmail.com" target="_blank">rupak.rajbhandari@gmail.com</a>></span><br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div></div><div>Hi,<div><br></div><div>The following simple script is working well. But when I remove the line 'set dfile 2', then it says data request completely outside file request. why?</div>
<div><br></div><div>Thanks.</div><div>
<br></div><div>-rupak</div><div><br></div><div>****************************************</div><div><div>'reinit'</div><div><br></div><div>'open d:\indus\gridded\precis\rfQ0.6190\wholeindus\rfQ0.6190.WI.ctl'</div>
<div>'open d:\indus\gridded\precis\rfQ0.1198\wholeindus\rfQ0.1198.WI.ctl'</div><div><br></div><div>'exec d:\indus\gs\colors'</div><div>'set gxout grfill'</div><div><br></div><div>'display ave(rf.1,t=1,t=360,12)'</div>
<div>n=1</div><div>while (n<50000) ; n=n+1 ; endwhile</div><div>'set dfile 2'</div><div>'display ave(rf.2,t=1,t=360,12)'</div><div>'cbar'</div></div><div><br></div>
<br></div></div>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><font color="#888888"><br><br clear="all"><br>-- <br>Roger Rodrigues Torres<br>PhD. Student in Meteorology<br>National Institute for Space Research (INPE)<br>Phone: <a href="tel:55%2012%203186-9541" value="+551231869541" target="_blank">55 12 3186-9541</a><br>
Personal Webpage: <a href="http://rtorres.webnode.com.br" target="_blank">http://rtorres.webnode.com.br</a><br>
</font><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org" target="_blank">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br></div>
</div></div><br>_______________________________________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" target="_blank">http://gradsusr.org/mailman/listinfo/gradsusr</a><br>
<br></blockquote></div><br><br clear="all"><br>-- <br>Roger Rodrigues Torres<br>PhD. Student in Meteorology<br>National Institute for Space Research (INPE)<br>Phone: 55 12 3186-9541<br>Personal Webpage: <a href="http://rtorres.webnode.com.br" target="_blank">http://rtorres.webnode.com.br</a><br>