<div dir="ltr">Hi Will,<div><br></div><div>Thanks for your help. However, the span command was not recognized.</div><div><br></div><div>Regards,</div><div>James</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">
On 14 August 2014 09:08, Will Hobbs <span dir="ltr"><<a href="mailto:Will.Hobbs@utas.edu.au" target="_blank">Will.Hobbs@utas.edu.au</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word;color:rgb(0,0,0);font-size:14px;font-family:Calibri,sans-serif">
<div>James</div>
<div><br>
</div>
<div>As the name of function might suggest month_to_season() is designed to convert monthly data to seasonal data; it's not designed to work with 4xtimes daily data. If in doubt about any function, your first step should ALWAYS be to consult the NCL help page
for that function (<a href="http://www.ncl.ucar.edu/Document/Functions/Contributed/month_to_season.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/Contributed/month_to_season.shtml</a>)</div>
<div><br>
</div>
<div>There is a page listing all of NCL's functions: <a href="http://www.ncl.ucar.edu/Document/Functions/climo.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/climo.shtml</a>, at the moment they only work with monthly or daily data, so you will have to use
a different approach.</div>
<div><br>
</div>
<div>My advise would be to use one of NCL's date functions (<a href="http://www.ncl.ucar.edu/Document/Functions/date.shtml" target="_blank">http://www.ncl.ucar.edu/Document/Functions/date.shtml</a>) to find all the data points for a given season. I've given an example below to calculate the seasonal mean for Jan-Mar (JFM), whcih you can adjust
to your own purposes. This assumes that your data has an appropriate time coord system, which I'm calling 'time' here. Lt's say you have an array, x, dimensioned (time,lat,lon)</div>
<div><br>
</div>
<div>> nyr = yren - yrst+1 ;yrst and yren are the forst and last years in x</div>
<div>> xSeas = x(:nyr-1,:,:) ; nyr is the number of years in your data set</div>
<div>>year = span(yrst, yren,1) ; integer array of all years in the range</div>
<div>> tdat = cd_calendar(x&time, -5) ; converts time to a multidimensional array giving year, month, day, time etc for each time step</div>
<div><br>
</div>
<div>> do yy = 0, nyr-1</div>
<div>> ii = ind(tdat(:,0).eq.year(yy).and.tdat(:,1).ge.1.and.tdat(:,1).le.3) ;location of all times in given year, that are between Jan (1) and March(3)</div>
<div>> xSeas(yy,:,:) = dim_avg_n(x(ii,:,:), 0) ; get average of all those data points</div>
<div>> delete(ii)</div>
<div>> end do</div>
<div><br>
</div>
<div><br>
</div>
<div>Hope that helps. Don't forget that there is an enormous amount of help to be gained by on the NCL website; it has a search function which includes archived 'ncl-talk' posts. Nine times out of ten somebody else will have come across the same problem. </div>
<div><br>
</div>
<div>Will </div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<div><br>
</div>
<span>
<div style="font-family:Calibri;font-size:11pt;text-align:left;color:black;BORDER-BOTTOM:medium none;BORDER-LEFT:medium none;PADDING-BOTTOM:0in;PADDING-LEFT:0in;PADDING-RIGHT:0in;BORDER-TOP:#b5c4df 1pt solid;BORDER-RIGHT:medium none;PADDING-TOP:3pt">
<span style="font-weight:bold">From: </span>James Ciarlo` <<a href="mailto:james.ciarlo@physics.org" target="_blank">james.ciarlo@physics.org</a>><br>
<span style="font-weight:bold">Date: </span>Thursday, 14 August 2014 4:43 PM<br>
<span style="font-weight:bold">To: </span>"<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>" <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>><br>
<span style="font-weight:bold">Subject: </span>[ncl-talk] eof script problems<br>
</div><div><div class="h5">
<div><br>
</div>
<div>
<div>
<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Dear NCL users,</span>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">I am new to NCL and am trying to modify an eof script to suit a different dataset. </div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">The data the script is designed for uses a command month_to_season but since my data is not monthly but rather 4 time daily, I need an alternative method. And in fact running the script gives me the following
error (to start with, as I am sure fixing this will result in numerous other errors later but lets start with this)</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">(0) contributed.ncl: month_to_season: dimension must be a multiple of 12</div>
<div style="font-family:arial,sans-serif;font-size:13px"><br>
</div>
<div style="font-family:arial,sans-serif;font-size:13px">Regards,</div>
<div style="font-family:arial,sans-serif;font-size:13px">James</div>
</div>
</div>
</div>
</div></div></span>
</div>
</blockquote></div><br></div>