关于NCL的逐日距平函数

麦田_smile

正常的形式应该是这样的：

    TIME     = cd_calendar(timec, 0)
    year      = toint( TIME(:,0) )
    month   = toint( TIME(:,1) )
    day       = toint( TIME(:,2) )
    ddd       = day_of_year(year, month, day)
    yyyyddd = year*1000 + ddd     
    yClmDay    = clmDayTLL(hgt, yyyyddd)
    yAnom      = calcDayAnomTLL (y, yyyyddd, yClmDay)

疑惑点：比如我求10年的距平，得到的yClmDay的时间维是366，因为10年里肯定包含闰年。但在计算距平的时候，对闰年来说，应该说没问题，因为时间维都是366。如果是平年，2月28号扣掉2月28号的平均值，那3月1号扣掉的是2月29号的呢还是3月1号的呢？从yyyyddd来看，似乎扣掉的是3月1号的，那后面的时间岂不都错位了？

请问各位怎么理解？

麦田_smile

小沧沧 发表于 2019-10-31 14:06
请问楼主，这个久远的问题你有答案了么？平年3月2日扣掉3月1日均值。。。岂不是相当不准确

http://www.ncl.ucar.edu/Document ... alcDayAnomTLL.shtml 现在版本的函数已经没有这个问题了 时间是366天 具体请看主页

zongqing123

打算计算季节内尺度的，楼主解决问题没有{:eb303:}{:eb303:}

风往北吹

zongqing123 发表于 2016-10-10 18:16
打算计算季节内尺度的，楼主解决问题没有

链接

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Apr 26 2010 - 11:36:26 MDT
The function is hard wired to return (366,:,:)
The input time array is ' yyyyddd '. If the yyyyddd
is correct for the 365 or 360 day calendars, then
    x366 = clmDayTLL(yyyyddd, x)
    xNew = x366(0:364,:,:) ; 365 day year
    xNew = x366(0:359,:,:) ; 360 day year
Good luck
On 4/26/10 11:18 AM, Bridget Thrasher wrote:
> I have 5 files, each with 10 years worth of daily data on a "noleap"
> calendar (3650 days total). I am using clmDayTLL to create the daily
> climatology from these files. The function is creating an array of 366
> days even though I have no leap days in my data. Is this the way it's
> supposed to work? Would it do the same with data on a 360-day calendar?
>
> -Bridget
>
> --
> Bridget Thrasher, PhD
> Postdoctoral Researcher
> Climate Central
> www.climatecentral.org <http://www.climatecentral.org>
>
>
>
>

zongqing123

非常感谢楼主的分享{:eb502:}{:eb502:}

麦田_smile

风往北吹 发表于 2016-10-11 00:41
链接

From: Dennis Shea  

谢谢回复。如果是这样的话，我的理解是，对于平年来讲，3月1日扣掉的多年2月29日的均值，而3月2日扣除的是3月1日的均值，以此类推，这样看似乎还是不太准确的哦？

小沧沧

请问楼主，这个久远的问题你有答案了么？平年3月2日扣掉3月1日均值。。。岂不是相当不准确

小沧沧

麦田_smile 发表于 2019-11-4 17:24
http://www.ncl.ucar.edu/Document/Functions/Contributed/calcDayAnomTLL.shtml 现在版本的函数已经没有 ...

看了官网说明及各种栗子，明白了！谢谢楼主！