[ncl-talk] delete 31 march for leap year
Marcus Morgan
mmorgan2009 at my.fit.edu
Mon May 11 20:59:12 MDT 2015
Here you would use the modulus function. Every four years a leap year
occurs, so:
leap=ind(mod(Time(:,0), 4) .eq. 0 .and. TIME(:,1).eq.3 .and.
TIME(:,2).eq.31) ; Chooses March 31st From those Leap Years
keep=ind(.not.Time(leap,0))
FOO = foo(keep,:,:)
delete(foo)
Haven't tested it, but it should work.
Good Luck
*Marcus N. Morgan*Graduate of Meteorology
*Department of Marine and Environmental Systems* Florida Institute of
Technology
150 W. University Blvd.
Melbourne, FL 32901
On Fri, Apr 17, 2015 at 1:52 AM, Xi Chang <xi.chang01 at gmail.com> wrote:
> Thanks Dennis,
>
> Yes because i just need to analyze JFM (90 days) for 40 years. I know the
> above trick but i just want delete 31 march only for LEAP YEAR not for
> whole years. Do you know how to do that?
>
> thanks
>
>
> On Friday, April 17, 2015, Dennis Shea <shea at ucar.edu> wrote:
>
>> [1] Yes ... there is a way
>> [2] But ... why would ypu want to do that? It introduces a 'spectral gap'
>> in the time series.
>>
>> ===
>> foo(time,lat,lon)
>>
>> time has units understood by cd_calendar
>> https://www.ncl.ucar.edu/Document/Functions/Built-in/cd_calendar.shtml
>>
>> TIME = cd_calendar(time, 0)
>>
>> ikeep = ind(.not.(TIME(:,1).eq.3 .and. TIME(:,2).eq.31) )
>> FOO = foo(ikeep,:,:)
>> delete(foo)
>>
>>
>>
>> On Thu, Apr 16, 2015 at 4:16 PM, Xi Chang <xi.chang01 at gmail.com> wrote:
>>
>>> Hi NCL,
>>> I have 40 years daily data, i need to delete 31 march only for the leap
>>> year, is there any quick way to do that?
>>>
>>> Regards
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150512/64c9c923/attachment.html
More information about the ncl-talk
mailing list