| Contour with Map_Set after reading ASCII column data [message #66520] |
Fri, 15 May 2009 14:09 |
jacques.frastain
Messages: 1
Registered: May 2009
|
Junior Member |
|
|
Hi guys,
I am new to IDL using it to process data for my master project. I have
the following problem when it comes to using CONTOUR together with
MAP_SET after reading in an ASCII data set.
My data file looks in principle as follows:
lon lat temp ; this header actually doesn't exists in the .dat
file
10 10 5
20 10 10
30 10 15
10 15 6
20 15 11
30 15 16
....
I use the following code to read the .dat file:
OPENR,lun,'lon_lat_temp.dat',/GET_LUN
rows = 1000
data = fltarr(3,rows)
READF,lun,data
lon = data(0,*)
lat = data(1,*)
temp = data(2,*)
lon2 = REFORM(lon)
lat2 = REFORM(lat)
temp2 = REFORM(temp)
This way I create a data array 'DATA FLOAT = Array[3, 1000]' and by
using REFORM I am able to create for each of the columns within
my .dat file a row vector. However, from this point onwards I am in
trouble - say I would like to plot a world map showing contours of the
temperature (= temp) depending on my latitudes (lat) and longitudes
(lon).
Can anyone help me solving this issue?
Many thanks!
Jacques
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
