;***************************************************************************
; read CXI surface met data. format: NASA-Ames
; PSD provides this code to make it easier to read the data. We do
; not support NCL and advise those with questions to read NCL webpages
; and/or contact their support staff.
;
; This code reads the CXISLAND Sfc met data  into arrays and plots a variable.
; See http://www.esrl.noaa.gov/psd/enso/rapid_response/data_pub/
;*************************************************************************** 
; written by CAS 8/31/2016. V1.0
;*************************************************************************** 

begin

; --- Read Data ----------------------------------------;

  filein="SHOUT_ENRR-UASO3_GH_20160202_RA.ict"
  titledate = filein
  cnlvl=systemfunc("/usr/local/bin/wc "+ filein)
  nlvl=stringtointeger(cnlvl) 
  print(nlvl+" obs+70")
  skiprows=39
  ncol = 2
  nlvl=nlvl-skiprows
  print(nlvl+" lines of obs")

  TestData = readAsciiTable (filein , (/ncol/), "float",skiprows) 
  TestData@_FillValue=-999.0
  TestData@missing_value=-999.0
  
  time   = TestData (:,0)     ; time [UTCday)
  ozone  = TestData (:,1)     ; air pressure     [mb] or [hPa] 


  i=0

;*************************
; create plot
;*************************

  res=True
  wks  = gsn_open_wks ("png", "SHOUT_ozone")
  res@tiMainString     = "Time(UTC day) vs ozone(ppb) plane"          ; add title
  plot  = gsn_csm_xy (wks,time,ozone,res) ; create plot
  max_ozone = max(ozone)
  min_ozone = min(ozone)   
  print(max_ozone+" max ozone")
  print(min_ozone+" min ozone")  

 end