麻烦前辈指点用IDL如何读取探空数据

liuqinglong

我拿到的探空数据是txt格式的，但是不会用IDL读取，我将数据前面的版本信息，测站和仪器信息等全删了，变成纯数据的文本，还是读不了，请前辈指点~

topmad

鄙人也不会IDL  但是根据你的问题来看 应该还是IDL操作的问题 这个东西你还是看看IDL的基本操作  文件读取那部分

chenhongfengye

你要读所有的数据吗？
我读过水汽压，相对湿度之类的

chenhongfengye

ii=10000
  
  for i=0l,ii-1 do begin
;  for i=0,0 do begin
  
  head_line=' '
  time_year=' '
  time_month=' '
  time_day=' '
  time_time=' '
  
  elevation=0.
   
  readf,lun,head_line
;  print,head_line
  
  time_year=strmid(head_line,5,4)
  time_month=strmid(head_line,9,2)
  time_day=strmid(head_line,11,2)
  time_time=strmid(head_line,13,2)
   
  num_data=fix(strmid(head_line,19,4))


  data_arr=fltarr(15,num_data)

  readf,lun,data_arr;,format=format_data
  

  elevation=data_arr[1,0]        ;elevation of station


if (data_arr[1,0] ne -99999) and (num_data ge 10) then begin

;if (data_arr[1,0] ne -99999) and (num_data ge 10) and (time_month eq '06' or $
;    time_month eq '07' or time_month eq '08') then begin
     
  press_arr=data_arr[0,*]/100.       ;observed pressure (mb)=100 Pa= 100*(kg*m/s2)/m2
  
  OBSgph_arr=data_arr[1,*]           ;observed geopotential Height (m)
  CALCgph_arr=data_arr[2,*]          ;calculated geopotential Height (m)
  
  AGL_Height_arr=CALCgph_arr-elevation   ;above ground level
  
  temp_arr=data_arr[3,*]/10.         ;observed temperature (K)
  tempgrad_arr=data_arr[4,*]/10.     ;temperature gradient (K/km)

  ptemp_arr=data_arr[5,*]/10.        ;potential temperature (K)
  ptempgrad_arr=data_arr[6,*]/10.    ;potential temperature gradient (K/km)
  
  VAPpress_arr=data_arr[7,*]/100.     ;vapor pressure (mb)
  SATvap_arr=data_arr[8,*]/1000.       ;saturation vapor pressure (mb)
  RH_arr=data_arr[9,*]/10.          ;relative humidity (%,percent)
  RHgrad_arr=data_arr[10,*]/10.      ;relative humidity gradient (%/km)


ind_99999=where(data_arr[0,*] ne -99999 and data_arr[3,*] ne -99999 and $
          data_arr[7,*] ne -99999 and data_arr[9,*] ne -99999,cc_99999)

VAPpress=VAPpress_arr(ind_99999)
press=press_arr(ind_99999)
temp=temp_arr(ind_99999)
ptemp=ptemp_arr(ind_99999)
RH=RH_arr(ind_99999)
CALCgph=CALCgph_arr(ind_99999)
AGL_Height=AGL_Height_arr(ind_99999)

num_data=n_elements(CALCgph)

liuqinglong

非常感谢您无私的帮助！
我想读的数据包括测站的基本信息，比如经纬度，还有温度，高度，气压等这些数据。

liuqinglong

chenhongfengye 发表于 2013-1-2 11:17

                               
登录/注册后可看大图

ii=10000
  
  for i=0l,ii-1 do begin

非常感谢您无私的帮助！
我想读的数据包括有温度，高度，气压等，还有测站的经纬度这些数据。

liuqinglong

topmad 发表于 2013-1-2 10:19

                               
登录/注册后可看大图

鄙人也不会IDL  但是根据你的问题来看 应该还是IDL操作的问题 这个东西你还是看看IDL的基本操作  文件读取那 ...

谢谢你的热心提示~

chenhongfengye

不客气，加油哈。。。