;************************************************
; Program: read.xband.nxl
; Purpose: show the variableas and dimenstions in the xband radar files
; Provided as a courtesy. Please see NCL webpages for questions.
;************************************************
; to run
; ncl -Qn <read.xband.ncl
; by Cathy Smith oct 2016.

begin
  infile="csu.rxm25-20160203-224453.qc.netcdf"
  f = addfile(infile,"r")
  vNames = getfilevarnames(f) ; get variable names of all groups on file
; alternately you can uncomment out the next line to get a simple list.
;  print (vNames)               ; print all variable names on file
; here we loop through the variables and print the dimenstions (1 or 2).

  nNames = dimsizes (vNames)
  print("List of varaibles in file "+infile+":")
  print(" ")
  do n=0,nNames-1              ; loop thru each variable
      v = f->$vNames(n)$
      dimv = dimsizes(v)
      rank = dimsizes(dimv)     ; rank [ie: number of dimensions]
      if(rank.eq.1)then
        print(vNames(n)+" ("+dimv+")")
      else
        print(vNames(n)+" ("+dimv(0)+","+dimv(1)+")")
      end if
      delete(dimv)
      delete(v)
   end do
end