source: trunk/ToBeReviewed/LECTURE/GRIB/read_gribtable.pro @ 67

;+
; NAME:
;       read_gribtable
;
; PURPOSE:
;       Read contents of a gribtable. Gribtables are located
;       in the gribtables subdirectory of HIPHOP
;   
; CATEGORY:
;       HIPHOP, GRIB, ECMWF
;
; CALLING SEQUENCE:
;       read_gribtable,tablename,parmtabl=parmtabl
;
; EXAMPLE:
;       tablename='ectab_128'
;
; INPUTS:
;       tablename : the full path name of a gribtable file
;
; OPTIONAL INPUT PARAMETERS:
;
; KEYWORD INPUT PARAMETERS:
;
; OUTPUTS:
;       parmtable : the parameter table
;
; COMMON BLOCKS:
;
; SIDE EFFECTS:
;
; RESTRICTIONS:
;
; PROCEDURE:
;       
; MODIFICATION HISTORY:
;
; Dominik Brunner, Apr 2000
;-
PRO read_gribtable,tablename,parmtabl=parmtabl,center=center,$
                   subcenter=subcenter,tablnum=tablnum
  
ON_ERROR,2

parmtabl=StrArr(3,256)
center=-1
subcenter=-1
tablnum=-1
; First Subscript (3) is [name,description,units]
; Second  (256) is defined size of a parameter table

IF n_elements(tablename) EQ 0 THEN return

openr,lun,tablename,/get

line=''
; read first line which (eventually) contains information
; about center, subcenter and table number
readf,lun,line
parts=STR_SEP(line,':')
IF n_elements(parts) GT 3 THEN BEGIN
  center=fix(parts[1])
  subcenter=fix(parts[2])
  tablnum=fix(parts[3])
ENDIF ELSE BEGIN
  IF n_elements(parts) GE 3 THEN parmtabl[0:1,fix(parts[0])]=parts[1:2] $
  ELSE IF n_elements(parts) EQ 2 THEN parmtabl[0,fix(parts[0])]=parts[1]
ENDELSE

; loop over remaining lines
REPEAT BEGIN
  readf,lun,line
  parts=STR_SEP(line,':')
  IF n_elements(parts) GE 3 THEN parmtabl[0:1,fix(parts[0])]=parts[1:2] $
  ELSE IF n_elements(parts) EQ 2 THEN parmtabl[0,fix(parts[0])]=parts[1]
END UNTIL EOF(lun)

free_lun,lun

; fill up missing varible names
index=WHERE(parmtabl[0,*] EQ '',count)
IF count GT 0 THEN parmtabl[0,index]='var'+strcompress(index,/rem)

END