Re: Smoothing a line

Thank you all!! Karin, I wasn't very sure about using "runave_n_Wrap" for the red line when I was planning to use a polynomial interpolation (http://goo.gl/tVBWR5), so I would avoid the borders effects I have in the running average. But I must admit that your red line (running average) really fits better than my red line (polynomial interpolation - http://goo.gl/tVBWR5). After some consideration on what to do, I realized that, because the 365 values in the array refers to days in a year, I can duplicate January and put at the end of December and duplicate December and put at the before January. This would avoid losing information from the end and beginning of the year.   Furthermore, I don't have the same (correct) results of Dennis Shea even using his own code. I must assume that Mary is right and the problem is the NCL version I'm running on (6.1.2). I will upload my NCL version.  Again, thank you very, very much. You are the best!  Best wishes, Rodrigo. Em Quarta-feira, 16 de Abril de 2014 13:29, Dennis Shea <shea_at_ucar.edu> escreveu: The attached matches your Excel [1] 13 pt running average [2] polynomial fit to [1] [3] plot On 4/16/14, 7:57 AM, Mary Haley wrote: > Hi Rodrigo, > > Please see the attached PNG image.  Is this what you’re expecting? > > NCL V6.1.2 gives me a plot like the one you are seeing, and V6.2.0 gives me the > attached (which I converted to PNG so I could trim it). > > “lspoly” is the culprit, and if you look at the “what’s new” page for NCL > V6.2.0, you’ll see that “lspoly” was replaced with a new routine: > > http://www.ncl.ucar.edu/current_release.shtml > > You may want to consider upgrading to NCL V6.2.0. > > —Mary > > > On Apr 15, 2014, at 9:21 PM, Rodrigo Colpo <rodrigocolpo2_at_yahoo.com.br > <mailto:rodrigocolpo2_at_yahoo.com.br>> wrote: > >> Dear Alan (and other users able to help), >> >> My plot still looks very different from the red line I had in Excel (a 4th >> order polynomial interpolation). >> What I'm doing wrong? >> >> The full code and data I'm using is bellow and the graph is attached. >> >> Thank you very, very much! >> Rodrigo. >> >> >> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" >> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" >> >> y = (/24.9264,-8.361866, -32.16906, 20.60294, 7.384959, 40.00099, -92.79871, >> 50.42744, -5.070408, 40.27705, -27.88171, -15.03915, 11.47537, -29.03874, >> -1.756587, -20.16744, 21.60254, -22.96009, 20.59936, -28.42203, 27.0571, >> 12.97648, -88.9543, 85.17473, -55.08803, 45.47628, 24.15122, -58.92344, >> 43.25131, -51.89326, 20.54425, 33.07741, -39.74909, -11.37656, 51.00746, >> -15.13431, 36.12533, -1.521067, -52.99825,12.57256, 28.42521, 3.828768, >> -53.74265, 26.17332, -9.199653, -2.309943, -29.50837, 41.74754, -10.51091, >> 31.03505, 4.09079, -43.97025, 31.43508, 10.02173, 41.16813, -2.80632, >> -79.53146, -20.74533, 59.44746, -49.75961, 39.98883, 10.05335, 0.8409821, >> -35.33345, 18.40772, -12.2347, 1.105038, 31.84602, -13.8735, -7.891514, >> 0.3682705, -5.585191, -28.66333, 52.90786, -5.409936, 7.680257, 0.7249982, >> -40.92494, -25.33436, 23.74692, 57.71513, -32.00339, -16.75868, 23.42031, >> -4.792629, 37.20881, -4.456624, -23.72248, 7.767448, -21.8565, 26.44745, >> -29.24916, -6.314954, 30.94873, -13.91431, 21.38537, -37.93886, 26.56374, >> -32.10296, 47.78776, -0.4801184, -25.8139, 58.57794, -0.7972997, -40.70943, >> -19.4354, 6.80695, -15.61182, 16.48778, 16.05497, 16.69331, 5.19874, >> -23.85107, 27.78862, 40.81667, -84.96756, -12.97966, -9.382995, 7.511517, >> 23.06323,-14.97894, 15.4442, 6.898167, -17.49207, 14.06445, 11.29572, >> 37.50677, -26.72132, 35.59035, -4.198732, -54.52291, -16.66267, 41.58054, >> 32.10245, -10.67811, 22.81898, 3.371699, -20.39278, -2.959632, -12.20442, >> 7.349309, 12.05357, 19.43616, -4.506427, -38.90607, 2.41556, >> 56.17836,-14.21334, 16.37525, 26.71167, 9.09359, 39.39703, 21.15053, >> -71.65178, -8.757985, -22.7589, 14.46975, 33.85746, -17.84287, 23.09438, >> 17.78081, 25.97645, -7.781834, -26.64459, -6.69821, -20.50336, 37.03909, >> 9.835811, 2.60719, 5.255515, -12.98845, -7.794468, -24.1282, 9.487386, >> 6.582564, 25.11583, 15.52696, 37.84516, 1.754072, -16.12908, -16.88242, >> 36.66131, -1.120736, -26.04581, 23.76035, 32.78828, -31.48173, 36.0989, >> 9.240784, -14.55615, 16.70981, 20.86206, -17.27525, 39.56828, 13.07125, >> -12.99599, -17.7401, -46.12418, 19.58386, 56.07545, -0.5685795, -0.9008529, >> -25.9442, -1.752784, 10.51259, 9.666623, 47.25177, -17.35054, 9.756303, >> -30.42878, 39.36865, 3.802191, 10.75947, -32.46304, 4.70647, -0.3541092, >> 16.2095, -25.91894, 13.74797, 26.89765, -6.783369, -3.999881, -4.279597, >> -0.8352942, -29.25169, -13.72089, 20.788, -0.8843041, -9.306261, -1.17004, >> 8.377829, 40.41809, 18.00184, 23.30745, -21.75604, -51.99948, -8.58614, >> -16.91891, 24.7299, 3.172685, -19.49131, -14.85031, 10.22103, -51.69467, >> -0.7927365, 38.09153, 1.860253, 30.74613, -3.257262, 4.002953, 0.2218886, >> 13.85386, -31.56291, -12.19289, 32.99896, -32.52431, -19.88405, -54.25632, >> 20.47882, -8.080286, 35.64525, -15.34531, 13.49685, -13.06714, -3.831783, >> 19.64373, -68.76797, 12.5983, 3.35871, -18.91072, 10.27991, -27.42312, >> -15.79099, 29.58066, -9.914533, -0.6830215, -9.077754, -9.185384, -25.31908, >> 67.3791, -30.95538, -22.8384, -25.2431, 17.0081, 66.54207, -16.88133, >> -16.24662, -36.7114, -25.94501, 23.17612, -45.73708, 30.64136, -18.41788, >> 40.28629, -54.11975, 44.79348, -35.77917, -7.19503, -13.19775, -12.92189, >> 24.29116, -16.98229, 49.63627, -6.452826, 11.50777, -64.53977, -5.951408, >> 37.638, -15.28019, -26.81924, 63.38461, -23.64748, 8.656753, -21.12774, >> 60.26141, -81.49531, 34.91807, -16.63806, 13.47313, -27.20161, -5.661833, >> -5.45722, 12.96912, 14.15574, -13.23332, -45.49319, 28.04377, 0.7724994, >> -17.6555, -20.47201, 28.1443, -14.44097, 23.40522, -0.1788197, -68.70046, >> 16.8583, 21.12821, 19.91207, -17.77668, -9.989548, -49.73676, 54.32707, >> -2.146452, -32.35265, 32.95036, 2.574606, -9.437281, -15.56939, -0.3845669, >> 12.12393, 14.93732, -28.40279, -14.69085, 49.71231, -1.143925, -6.234774, >> -55.59104, 61.16213, -68.4799, 16.58861, 14.23054, -47.54554, 49.72101, >> -48.24809, 32.33456/) >> x  = ispan(1,365,1) >> c  = lspoly(x, y, 1, 4) >> yo  = c(0) + c(1)*x + c(2)*x^2 + c(3)*x^3; >> xo_new = fspan(1.,365.,365*5) >> yo_new = ftcurv(x, yo, xo_new) >> >> ;---Start the graphics. >> wks = gsn_open_wks ("ps","plot")            ; open workstation >> res                  = True                    ; plot mods desired >> res_at_xyMarkLineModes  = "Lines" >> plot = gsn_csm_xy (wks, xo_new, yo_new, res) >>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Apr 17 11:29:38 2014