begin

ifil1="correl.rain.with.slp.WEST_of_Goyder_Region.nc"
ifil2="correl.rain.with.slp.CENTRAL_of_Goyder_Region.nc"
ifil3="correl.rain.with.slp.NEAST_of_Goyder_Region.nc"
ifil4="correl.rain.with.slp.SEAST_of_Goyder_Region.nc"

dfh1=addfile(ifil1,"r")
dfh2=addfile(ifil2,"r")
dfh3=addfile(ifil3,"r")
dfh4=addfile(ifil4,"r")

lat=dfh1->lat
lon=dfh1->lon

data1=dfh1->correl
data2=dfh2->correl
data3=dfh3->correl
data4=dfh4->correl

; wks = gsn_open_wks("x11","plot4")
 wks = gsn_open_wks("ps","plot4")               ; open a ps file
  gsn_define_colormap(wks,"precip_diff_12lev")  ; choose a colormap
  
  plots = new(4,graphic)                          ; create a plot array

  res                     = True

;  gsn_reverse_colormap(wks) 

  res@mpDataSetName         = "Earth..4"   ; This new database contains
                                           ; divisions for other countries.
  res@mpDataBaseVersion     = "MediumRes"  ; Medium resolution database
  res@mpOutlineOn           = True         ; Turn on map outlines
  res@mpOutlineSpecifiers   = (/"Australia:states"/)

  res@gsnDraw             = False           ; don't draw
  res@gsnFrame            = False           ; don't advance frame
  res@cnInfoLabelOn       = False           ; turn off cn info label
  res@cnFillOn            = True            ; turn on color
  res@cnLinesOn           = False       ; do not draw contour lines

  res@lbLabelBarOn        = True

  res@mpCenterLonF         = 130.
  res@mpMinLonF            = 100.
  res@mpMaxLonF            = 160.0
  res@mpMinLatF            = -45.0
  res@mpMaxLatF            = -7.5
  res@mpFillOn             = False    ; turn off modern continents
  res@gsnAddCyclic = False
  res@cnFillMode = "CellFill"


  res@gsnLeftStringFontHeightF=.03
  res@gsnLeftString  = "WESTERN SA"

  res@cnLevelSelectionMode = "ExplicitLevels"

  res@cnLevels             = (/-0.3,-0.2,-0.1,0.,0.1,0.2/)

  plots(0)=gsn_csm_contour_map_ce(wks,data1,res)
  res@gsnLeftString  = "CENTRAL SA"

  delete(res@cnLevels)
  delete(res@cnFillColors)
  res@cnLevels             = (/-0.3,-0.2,-0.1,0.,0.1,0.2/)


  plots(1)=gsn_csm_contour_map_ce(wks,data2,res)

  delete(res@cnLevels)

  res@cnLevels             = (/-0.3,-0.2,-0.1,0.,0.1,0.2/)

  res@gsnLeftString  = "N/EAST SA"
  plots(2)=gsn_csm_contour_map_ce(wks,data3,res)

  delete(res@cnLevels)

  res@cnLevels             = (/-0.3,-0.2,-0.1,0.,0.1,0.2/)

  res@gsnLeftString  = "S/EAST SA"
  plots(3)=gsn_csm_contour_map_ce(wks,data4,res)

  resP                  = True                   ; modify the panel plot
  
  resP@lbLabelFontHeightF  = 0.009           ; label bar font height

  resP@txFontHeightF = .02
  resP@txString         = "Correlation between PSL and RAINFALL"
  resP@gsnPanelBottom   = 0.05                   ; add space at bottom
  resP@gsnPanelFigureStrings= (/"a)","b)","c)","d)"/) ; add strings to panel
  res@txFontHeightF     = .24 
  gsn_panel(wks,plots,(/2,2/),resP)               ; now draw as one plot


end