Changeset 474 for trunk/SRC/ReadWrite/ncdf_gettime.pro

Timestamp:

03/14/12 09:45:50 (12 years ago)

Author:

smasson

Message:

take into accound hh:mm:ss when reading time in netcdf files

File:

• trunk/SRC/ReadWrite/ncdf_gettime.pro (modified) (3 diffs)

trunk/SRC/ReadWrite/ncdf_gettime.pro

@@ -72 +72 @@
     tobeclosed = 1
   ENDIF
+  IF n_elements(caller) EQ 0 THEN caller = 'read_ncdf'
   inq = ncdf_inquire(cdfid)
 ;----------------------------------------------------
@@ -204 +205 @@
     ENDIF
     start = str_sep(words[2], '-')
+    IF stregex(value, '[^ ]* since ([0-9]){1,4}-([0-9]){1,2}-([0-9]){1,2} ([0-9]){1,2}:([0-9]){1,2}:([0-9]){1,2}', /boolean) EQ 0 THEN BEGIN
+      print, 'hh:mm:ss of the calendar origin not found... -> start calendar at 00:00:00'
+      start = [start, '0', '0', '0']
+    ENDIF ELSE start = [start, str_sep(words[3], ':')]
   ENDIF ELSE units = value
 ;----------------------------------------------------
@@ -212 +217 @@
   case units OF
     'true julian day':
-    'second':time = julday(start[1], start[2], start[0], 0, 0, 0)+time/86400.d
-    'minute':time = julday(start[1], start[2], start[0], 0, 0, 0)+time/1440.d
-    'hour':time = julday(start[1], start[2], start[0], 0, 0, 0)+time/24.d
-    'day':time = julday(start[1], start[2], start[0], 0, 0, 0)+time
+    'second':time = julday(start[1], start[2], start[0], start[3], start[4], start[5])+time/86400.d
+    'minute':time = julday(start[1], start[2], start[0], start[3], start[4], start[5])+time/1440.d
+    'hour':time = julday(start[1], start[2], start[0], start[3], start[4], start[5])+time/24.d
+    'day':time = julday(start[1], start[2], start[0], start[3], start[4], start[5])+time
     'month':BEGIN
       if total(fix(time) NE time) NE 0 then $ ; we switch to days with 30d/m
-         time = julday(start[1], start[2], start[0], 0, 0, 0) + round(time*30) $
-      ELSE time = julday(start[1]+fix(time), replicate(start[2], jpt), replicate(start[0], jpt), 0, 0, 0)
+         time = julday(start[1], start[2], start[0], start[3], start[4], start[5]) + round(time*30) $
+      ELSE time = julday(start[1]+fix(time), replicate(start[2], jpt), replicate(start[0], jpt), start[3], start[4], start[5])
     END
     'year':BEGIN
       if total(fix(time) NE time) NE 0 then $ ; we switch to days with 365d/y
-         time = julday(start[1], start[2], start[0], 0, 0, 0) + round(time*365) $
-      ELSE time = julday(replicate(start[1], jpt), replicate(start[2], jpt), start[0]+fix(time), 0, 0, 0)
+         time = julday(start[1], start[2], start[0], start[3], start[4], start[5]) + round(time*365) $
+      ELSE time = julday(replicate(start[1], jpt), replicate(start[2], jpt), start[0]+fix(time), start[3], start[4], start[5])
     END
   ENDCASE