Changeset 284

Timestamp:

09/13/07 17:34:43 (17 years ago)

Author:

smasson

Message:

make sur dates are given as long integer

File:

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

trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

@@ -155 +155 @@
 ;------------------------------------------------------------
   IF n_elements(beginning) EQ 0 THEN BEGIN
-    beginning = 0
-    timestep = 1
+    beginning = 0L
+    timestep = 1L
   ENDIF
 ; define time and jpt
   CASE 1 OF
     keyword_set(timestep):BEGIN 
-      firsttps = beginning[0]
-      IF n_elements(ending) NE 0 THEN lasttps = ending[0] ELSE lasttps = firsttps
+      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)
@@ -172 +172 @@
       currentfile = (where(filelist EQ filename))[0]
       time = (*(extractatt(top_uvalue, 'fileparameters'))[currentfile]).time_counter
-      date1 = date2jul(beginning[0])
-      if n_elements(ending) NE 0 then date2 = date2jul(ending[0]) ELSE date2 = date1
+      date1 = date2jul(long(beginning[0]))
+      if n_elements(ending) NE 0 then date2 = date2jul(long(ending[0])) ELSE date2 = date1
       firsttps = (where(abs(time - date1) LT 0.9d/86400.d))[0]   ; beware of rounding errors...
       lasttps = (where(abs(time - date2) LT 0.9d/86400.d))[0]
@@ -186 +186 @@
       ENDIF
 ; date1
-      date1 = date2jul(beginning[0])
+      date1 = date2jul(long(beginning[0]))
 ; date2
-      if n_elements(ending) NE 0 then date2 = date2jul(ending[0]) ELSE date2 = date1
+      if n_elements(ending) NE 0 then date2 = date2jul(long(ending[0])) ELSE date2 = date1
 ; firsttps
       firsttps = where(time GE (date1 - 0.9d/86400.d)) & firsttps = firsttps[0]