Changeset 462 for trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

Timestamp:

04/15/11 19:49:48 (13 years ago)

Author:

smasson

Message:

get time axis even when using /timestep keyword in read_ncdfSRC/ToBeReviewed/LECTURE/read_ncdf.pro

File:

• trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro (modified) (2 diffs)

trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

@@ -189 +189 @@
       firsttps = long(beginning[0])
       IF n_elements(ending) NE 0 THEN lasttps = long(ending[0]) ELSE lasttps = firsttps
-      jpt = lasttps-firsttps+1
-      IF NOT keyword_set(callitself) then time = julday(1, 1, 1) + lindgen(jpt)
+      IF NOT keyword_set(callitself) then BEGIN 
+        time = ncdf_gettime(filename, cdfid, caller = 'read_ncdf', err = err, _extra = ex)
+        jpt = lasttps-firsttps+1
+        IF time[0] LT 0 then time = julday(1, 1, 1) + lindgen(jpt) $
+        ELSE time = time[firsttps:lasttps]
+      ENDIF ELSE jpt = lasttps-firsttps+1
     END
     keyword_set(parent):BEGIN
@@ -204 +208 @@
       jpt = lasttps-firsttps+1
     END
-    keyword_set(allrecords):BEGIN
+    keyword_set(allrecords) OR n_elements(beginning) EQ 0:BEGIN
       time = ncdf_gettime(filename, cdfid, caller = 'read_ncdf', err = err, _extra = ex)
       IF time[0] LT 0 then BEGIN