;;;
begin
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  print("Part1: Reading .csv files (2015_2017.csv)") 
  fname_R24     = "2015_2017.csv"
; asciiread (filepath, dimensions, datatype) ; If -1 is given for the dimensions parameter, 
;                      all values in the file will be read into a one-dimensional variable. 
  lines     = asciiread(fname_R24,-1,"string")
  print(lines(0))
  print(dimsizes(lines))

  delim = ","

;---Read fields 1-6 from fname_Report
  st_name         =           str_get_field(lines(1:),1,delim)
  station_id      = tointeger(str_get_field(lines(1:),2,delim))
  lon             =   tofloat(str_get_field(lines(1:),3,delim))
  lat             =   tofloat(str_get_field(lines(1:),4,delim))
  date            =           str_get_field(lines(1:),6,delim)
  Vart             =   tofloat(str_get_field(lines(1:),7,delim))
  
;---Print the information
  ;print("dimsizes(st_name) = "+dimsizes(st_name))
  ;print("dimsizes(station_id) = "+dimsizes(station_id))
  ;print("sum_R24 = "+Vart)
  ;print("station_id is '" + station_id + "', st_name is " + st_name)
  station_id_u   = get_unique_values(station_id)
  station_id_uc = count_unique_values(station_id)
  
  printVarSummary(station_id_u)
  printVarSummary(station_id_uc)
  print("station_id_uc = "+station_id_uc)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

end