load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
;-- read the data and define

  diri = "/home/model-user/Desktop/"
  fili = "IAGOS_timeseries_2013101819404151.nc" 
  pthi = diri+fili
  f = addfile(pthi, "r")
  
  
  var = f->O3_PM(0)  
  lat1d = f->lat
  lon1d = f->lon
  printVarSummary(var)
  print(var)
  print ("min/max var = " + min(var) + "/" + max(var)) 
  

wks = gsn_open_wks("png","plot_contour_filled_map")
;-- set resources
res = True
res@gsnMaximize = True
res@cnFillOn = True ;-- turn on contour fill
res@cnFillPalette  = "BlueWhiteOrangeRed" ;-- choose color map
res@cnLinesOn= False ;-- turn off contour lines
res@cnLineLabelsOn = False ;-- turn off line labels
res@cnLevelSelectionMode = "ManualLevels"
res@cnMinLevelValF = 0
res@cnMaxLevelValF = 250
res@gsnAddCyclic = False
res@lbBoxMinorExtentF = 0.15 ;-- decrease the height of the labelbar
res@tiMainString = "IAGOS data" ;-- title string
res@tiMainFontHeightF = 0.02
;---Lat/lon arrays
res@sfXArray = lon1d
res@sfYArray = lat1d

cnConstFEnableFill = True
   step = 5
   res@sfXArray = lon1d(:)
   res@sfYArray = lat1d(:) 
;-- draw the contour map
plot = gsn_csm_contour_map(wks, var, res)
end