Changeset 298


Ignore:
Timestamp:
10/05/07 16:56:26 (17 years ago)
Author:
smasson
Message:

bugfix related to the use of double precision calendar

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro

    r295 r298  
    171171      currentfile = (where(filelist EQ filename))[0] 
    172172      time = (*(extractatt(top_uvalue, 'fileparameters'))[currentfile]).time_counter 
    173       date1 = date2jul(long(beginning[0])) 
    174       if n_elements(ending) NE 0 then date2 = date2jul(long(ending[0])) ELSE date2 = date1 
     173      date1 = date2jul(beginning[0]) 
     174      if n_elements(ending) NE 0 then date2 = date2jul(ending[0]) ELSE date2 = date1 
    175175      firsttps = (where(abs(time - date1) LT 0.9d/86400.d))[0]   ; beware of rounding errors... 
    176176      lasttps = (where(abs(time - date2) LT 0.9d/86400.d))[0] 
     
    185185      ENDIF 
    186186; date1 
    187       date1 = date2jul(long(beginning[0])) 
     187      date1 = date2jul(beginning[0]) 
    188188; date2 
    189       if n_elements(ending) NE 0 then date2 = date2jul(long(ending[0])) ELSE date2 = date1 
     189      if n_elements(ending) NE 0 then date2 = date2jul(ending[0]) ELSE date2 = date1 
    190190; firsttps 
    191191      firsttps = where(time GE (date1 - 0.9d/86400.d)) & firsttps = firsttps[0] 
Note: See TracChangeset for help on using the changeset viewer.