Changeset 298
- Timestamp:
- 10/05/07 16:56:26 (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SRC/ToBeReviewed/LECTURE/read_ncdf.pro
r295 r298 171 171 currentfile = (where(filelist EQ filename))[0] 172 172 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 = date1173 date1 = date2jul(beginning[0]) 174 if n_elements(ending) NE 0 then date2 = date2jul(ending[0]) ELSE date2 = date1 175 175 firsttps = (where(abs(time - date1) LT 0.9d/86400.d))[0] ; beware of rounding errors... 176 176 lasttps = (where(abs(time - date2) LT 0.9d/86400.d))[0] … … 185 185 ENDIF 186 186 ; date1 187 date1 = date2jul( long(beginning[0]))187 date1 = date2jul(beginning[0]) 188 188 ; date2 189 if n_elements(ending) NE 0 then date2 = date2jul( long(ending[0])) ELSE date2 = date1189 if n_elements(ending) NE 0 then date2 = date2jul(ending[0]) ELSE date2 = date1 190 190 ; firsttps 191 191 firsttps = where(time GE (date1 - 0.9d/86400.d)) & firsttps = firsttps[0]
Note: See TracChangeset
for help on using the changeset viewer.