;----------------------------------------------------------------------
; maponly_27.ncl
;----------------------------------------------------------------------
; Concepts illustrated:
;   - Drawing six maps with six different land colors
;   - Changing the color of land
;   - Paneling four plots on a page.
;   - Turning off all tickmarks and labels
;   - Turning off borders of a map plot
;----------------------------------------------------------------------
;
; These files are loaded by default in NCL V6.2.0 and newer
; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
; load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   

begin
  wks  = gsn_open_wks("png","maponly1")  ; send graphics to PNG file

 gsn_define_colormap(wks,"OceanLakeLandSnow")
;---Set up common resources for all four maps
  res                   = True

  res@gsnMaximize       = True 
  res@gsnDraw           = False   ; Will panel later
  res@gsnFrame          = False

  res@mpDataBaseVersion  = "MediumRes"
  res@tmBorderThicknessF = 1.       ; 4x thicker
 
;---Turn off all tickmarks and labels
  res@tmXBOn            = True
  res@tmYLOn            = True
  res@tmXTOn            = True
  res@tmYROn            = True
;

  res@mpOceanFillColor            = 3          ; array index in color map
  res@mpLandFillColor             = 30          ; for land, ocean, and inland h20
  res@mpInlandWaterFillColor      = 3

; Start setting resources for individual maps
;
; --Mediterranean Sea
;
  res@mpMinLatF         = 30
  res@mpMaxLatF         = 50
  res@mpMinLonF         = 0
  res@mpMaxLonF         = 60

  res@tmXBBorderOn      = True
  res@tmXTBorderOn      = True
  res@tmYLBorderOn      = True
  res@tmYRBorderOn      = True

  res@vpWidthF          = 0.305214
  res@vpHeightF         = 0.152607

  map4 = gsn_csm_map_ce(wks,res)

;---Central Eurasia

  res@mpMinLatF         = 30
  res@mpMaxLatF         = 50

  res@mpMinLonF         = 60
  res@mpMaxLonF         = 120


  res@tmXBBorderOn      = True
  res@tmXTBorderOn      = True
  res@tmYLBorderOn      = True
  res@tmYRBorderOn      = True

  res@vpWidthF          = 0.305214
  res@vpHeightF         = 0.152607


  map5 = gsn_csm_map_ce(wks,res)

;---East China

  res@mpMinLatF         = 30
  res@mpMaxLatF         = 50
  res@mpMinLonF         = 120
  res@mpMaxLonF         = 180

  res@tmXBBorderOn      = True
  res@tmXTBorderOn      = True
  res@tmYLBorderOn      = True
  res@tmYRBorderOn      = True

  res@vpWidthF          = 0.305214
  res@vpHeightF         = 0.152607


  map6 = gsn_csm_map_ce(wks,res)

;---Sea of okhotsk
  res@mpMinLatF         = 50
  res@mpMaxLatF         = 60
  res@mpMinLonF         = 130
  res@mpMaxLonF         = 150

  res@tmXBBorderOn      = True
  res@tmXTBorderOn      = True
  res@tmYLBorderOn      = True
  res@tmYRBorderOn      = True

  res@vpWidthF          = 0.305214
  res@vpHeightF         = 0.152607

  map3 = gsn_csm_map_ce(wks,res)


;---Arctic Eurasia

  res@mpMinLatF         = 60
  res@mpMaxLatF         = 90
  res@mpMinLonF         = 40
  res@mpMaxLonF         = 180
 
  res@tmXBBorderOn      = True
  res@tmXTBorderOn      = True
  res@tmYLBorderOn      = True
  res@tmYRBorderOn      = True

  res@vpXF              = 0.2                  ; position and size
  res@vpYF              = 0.3

  res@vpWidthF          = 0.305214
  res@vpHeightF         = 0.152607

 
  map2 = gsn_csm_map_ce(wks,res)

;---Scandinavia

  res@mpMinLatF         = 60
  res@mpMaxLatF         = 80
  res@mpMinLonF         = 0
  res@mpMaxLonF         = 40
  
  res@tmXBBorderOn      = True
  res@tmXTBorderOn      = True
  res@tmYLBorderOn      = True
  res@tmYRBorderOn      = True

  res@vpWidthF          = 0.305214
  res@vpHeightF         = 0.152607

 
  map1 = gsn_csm_map_ce(wks,res)


;---Panel all four plots
  pres             = True
  pres@gsnMaximize = True
  pres@gsnPanelRowSpec = True  

  pres@gsnPanelDebug = True 

  gsn_panel(wks,(/map1,map2,map3,map4,map5,map6/),(/3,3/),pres)

end