Re: time coordinate in netcdf data

From: Kai Zhang <k.zhang.iap_at_nyahnyahspammersnyahnyah>
Date: Thu, 14 May 2009 17:48:28 +0200

Excellent solution! Thanks!

2009/5/14 Dennis Shea <shea_at_ucar.edu>:
> Hi Kai .... Try a combination of
> http://www.ncl.ucar.edu/Document/Functions/Built-in/ut_calendar.shtml
> http://www.ncl.ucar.edu/Document/Functions/Built-in/ut_inv_calendar.shtml
>
>
> time1 = f->time ; "days since 1998-10-01 00:00"
> time2 = g->time ; "hours since 2000-03-01 00:00
>
> date2 = ut_calendar(time2, 0) ; datex(:,6) [array]
>
> year = *floattointeger*
> <http://www.ncl.ucar.edu/Document/Functions/Built-in/floattointeger.shtml>(utc_date(:,0))
> ; Convert to integer for ut_inv_calendar
> month = *floattointeger*
> <http://www.ncl.ucar.edu/Document/Functions/Built-in/floattointeger.shtml>(date2(:,1))
> day = *floattointeger*
> <http://www.ncl.ucar.edu/Document/Functions/Built-in/floattointeger.shtml>(date2(:,2))
> hour = *floattointeger*
> <http://www.ncl.ucar.edu/Document/Functions/Built-in/floattointeger.shtml>(date2(:,3))
> minute = *floattointeger*
> <http://www.ncl.ucar.edu/Document/Functions/Built-in/floattointeger.shtml>(date2(:,4))
> second = date2(:,5) ; must be float or double
>
> ; note "time1"
> TIME2 = *ut_inv_calendar*(year,month,day,hour,minute,second,time1_at_units,
> 0)
> print(time1+" "+ time2+" "+TIME2)
> or, print(sprintf("%12.2f",time1)+" " \
> +sprintf("%12.2f",time2)+" " \
> +sprintf("%12.2f",TIME2) )
>
>
>
> Good Luck
> D
> =======================
> Kai Zhang wrote:
>>
>> Dear all,
>>
>> I have two netcdf data files which use different time coordinates. One
>> file is
>>
>> time:units = "days since 1998-10-01 00:00" and the other file is
>>
>> time:units = "hours since 2000-03-01 00:00".
>>
>> Is there a simple or well-developed method to convert the time
>> coordinate in second netcdf data to the same time coordinate in the
>> first data?
>>
>> Thanks!
>>
>> Best regards, Kai
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>

-- 
---------------------------------------------------------------
Kai Zhang
Max-Planck-Institut fuer Meteorologie
Bundesstrasse 53
20146 Hamburg
-------------------------------------------
Tel: +49 040 / 41173-130
E-Mail: kai.zhang_at_zmaw.de
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu May 14 2009 - 09:48:28 MDT

This archive was generated by hypermail 2.2.0 : Tue May 26 2009 - 10:17:16 MDT