Hi David
I tried to convert time axis to hours in the ff program, but all my
efforts are not working. I need some help and suggestion on this
convection.I have tried may ways but they are not converting in a
right way.
!p.multi = [0,1,2]
time =[24782555, 24782675, 24782795, 24782915, 24783035, 24783155,
24783275, 24783395, 24783515, 24783635, $
24783755, 24783875, 24783995, 24784115, 24784235, 24784355,
24784475, 24784595, 24784715, 24784835, $
24784955, 24785075, 24785195, 24785315, 24785435, 24785555,
24785675, 24785795, 24785915, 24786035, $
24786155, 24786275, 24786515, 24786635, 24786755, 24786875,
24786995, 24787115, 24787235, 24787355, $
24787475, 24787595, 24787715, 24787835, 24787955, 24788075,
24788195, 24788315, 24788435, 24788555, $
24788675, 24788795, 24788915, 24789035, 24789155, 24789275,
24789395, 24789515, 24789635]
vel = [0.200000, -12.0000, -6.70000, 2.40000, -11.4000, 3.40000,
1.30000, -2.80000, -8.60000, 1.30000, $
-0.600000, -2.70000, -7.40000, -8.40000, 0.100000, 0.000,
23.9000, -5.60000, -2.10000, -10.2000, $
-8.20000, 2.30000, -15.3000, 5.70000, 0.000, -9.80000, 0.00000,
-8.60000, 0.000, 29.6000, $
-32.5000, 8.50000, 0.000, 0.000, 0.0000, 0.0000, 0.0000, 0.000,
0.000, 0.000, $
0.000, 43.4000, 0.000, 7.10000, -2.90000, 22.0000, 13.1000,
-2.60000, 31.5000, 0.0000, $
41.8000, 27.8000, 31.4000, 0.400000, 24.6000, 16.3000,
-3.10000, 4.00000, 4.20000]
U = [24782520, 24782640, 24782760, 24782880,24783000, 24783120,
24783240, 24783360, 24783480, 24783600, $
24783720, 24783840, 24783960, 24784080, 24784200, 24784320, 24784440,
24784560, 24784680, 24784800, $
24784920, 24785040, 24785160, 24785280, 24785400, 24785520, 24785640,
24785760, 24785880, 24786000, $
24786120, 24786240, 24786360, 24786480, 24786600, 24786720, 24786840,
24786960, 24787080, 24787200, $
24787320, 24787440, 24787560, 24787680, 24787800, 24787920, 24788040,
24788160, 24788280, 24788400, $
24788520, 24788640, 24788760, 24788880, 24789000, 24789120, 24789240,
24789360, 24789480, 24789600]
Result = interpol(vel,time,U)
plot, time, vel, linestyle = 2
oplot, u, result
N = n_elements(result)
print, N
T = 120
Spec = fft(result)
freq_axis =findgen(N)/(N*T)
plot, freq_axis[0:N/2],abs(spec[0:N/2]), xrange = [0.001, 0.005]
end
Any help will be highly appreciated
Zolile
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
