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_2002021617203349.nc" pthi = diri+fili f = addfile(pthi, "r") var = f->O3_PM printVarSummary(var) print ("min/max var = " + min(var) + "/" + max(var)) time = ispan (62400,64260,5) long_name = "UTC time" wks = gsn_open_wks("png", "O3_PM lagos to Lome Flight_20020216") res = True res@xyExplicitLegendLabels = (/"Ozone at Lagos"/) res@xyLineColor = (/"red"/) res@xyLineThicknessF = 3 res@tiYAxisString = "O3_PM ppb" res@tiYAxisFont = 21 res@tiYAxisFontAspectF = 1.3 res@tiYAxisFontHeightF = 0.012 res@tiXAxisString = "Time" res@tiXAxisFont = 21 res@tiXAxisFontAspectF = 1.3 res@tiXAxisFontHeightF = 0.012 res@tiXAxisOffsetYF = 0.0 res@tiMainOffsetYF = 0.11 res@tiMainString = "Ozone_PM: Lagos to Lome flight" res@tmXBLabelFontAspectF = 1.1 res@tmXBLabelFontHeightF = 0.012 res@trYMinF = 15. res@trYMaxF = 100.0 res@gsnMaximize = True res@gsnDraw = False res@gsnFrame = False ;-- create plot plot = gsn_xy(wks,time,var,res) draw(wks) frame(wks) end