begin
; file path & name
  f = addfile("../3Dyn_vitt/dist_solver/build/testout.nc","r")

; specifc which time index to use (if only ne timestep the it is 0  
  itime2get = 0 
  
; get the variables     
  pot  = f->ELECPOTEN(:,:,:)
  potHl= f->HILAT_POT(:,:,:)
  pot  = pot*1.e-3    ; V to kV
  potHl= potHl*1.e-3  ; V to kV
  facHl= f->HILAT_FAC(:,:,:)
  
  date = f->date    ; YYYYMMDD
  ut   = tofloat(f->datesec) ; current seconds of current day
  printVarSummary(ut)
  ut   = ut/60.  ; sec to hours

;------------------------------------------------------------------------  
  wks = gsn_open_wks("png","ce")
  gsn_define_colormap(wks,"BlueDarkOrange18")
  plot = new(5,graphic)                          ; create a plot array
  
  res                       = True            ; plot mods desired
  res@gsnDraw               = False           ; don't draw
  res@gsnFrame              = False           ; don't advance frame

  res@cnFillOn              = True            ; turn on color
  res@cnLinesOn             = False           ; no contour lines
  res@cnFillMode            = "RasterFill"    ; Raster Mode
  res@cnLineLabelsOn        = False           ; no line labels
  res@lbOrientation         = "vertical"
  res@gsnYAxisIrregular2Linear = True
  res@tiXAxisString   = "mag. longitude"  
  res@tiYAxisString   = "mag. apex latitude"
  res@gsnLeftString = " " 
  res@gsnRightString= " " 
  res@vpWidthF      = 0.8
  res@vpHeightF     = 0.4
  res@gsnLeftString   = " "
; control the color bar
;  res@cnLevelSelectionMode = "ManualLevels"	; manually set the contour levels with the following 3 resources
;  res@cnMinLevelValF  = -2.			; set the minimum contour level
;  res@cnMaxLevelValF  =  2.			; set the maximum contour level
;  res@cnLevelSpacingF = 0.25			; set the interval between contours

;  res@trYMinF =  40.0   ; Set your minimum Y-axis value
;  res@trYMaxF =  90.0   ; Set your maximum Y-axis value
;  res@trXMinF = -50.0   ; Set your minimum Y-axis value
;  res@trXMaxF =  50.0   ; Set your maximum Y-axis value
  
   it =itime2get
   ; define what is plotted	
    res@tiMainString          = "potential [kV]: date="+date(it)+" ut="+ut(it)
    plot(0) = gsn_csm_contour(wks,pot(it,:,:),res)
    res@tiMainString          = "potential [kV] high lat: date="+date(it)+" ut="+ut(it)
    plot(1) = gsn_csm_contour(wks,potHl(it,:,:),res)

    res@tiMainString          = "FAC [A/m2?] high latitude date="+date(it)+" ut="+ut(it)
    plot(2) = gsn_csm_contour(wks,facHl(it,:,:),res)
 
    gsn_panel(wks,plot(0:2),(/3,1/),False)             ; now draw as one plot 
          
end