;   Example script display location of model domains
;   Only works for ARW domains
;   Beta code used
;   November 2008

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

begin
;

; We generate plots, but what kind do we prefer?
  type = "x11"
; type = "pdf"
; type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"wps_show_dom_ll-rot")

  colors = (/"white","black","White","ForestGreen","DeepSkyBlue"/)
  gsn_define_colormap(wks, colors)  

  mpres = True
  mpres@mpFillOn = True
  mpres@mpFillColors  = (/"background","DeepSkyBlue","ForestGreen","DeepSkyBlue", "transparent"/)
  mpres@mpGeophysicalLineColor      = "Black"
  mpres@mpGridLineColor             = "Black"
  mpres@mpLimbLineColor             = "Black"
  mpres@mpNationalLineColor         = "Black"
  mpres@mpPerimLineColor            = "Black"
  mpres@mpUSStateLineColor          = "Black"

  lnres = True 
  lnres@gsLineThicknessF = 2.5

  txres = True
  txres@txFont = "helvetica-bold"
  txres@txJust = "BottomLeft"
  txres@txPerimOn = False
  txres@txFontHeightF = 0.015

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;


  mpres@max_dom           = 2
  mpres@parent_id         = (/   1,   1,   1,   3,   1 /)
  mpres@parent_grid_ratio = (/   1,   3,   3,   5,   5 /)
  mpres@i_parent_start    = (/   1, 191,  25,  25,  10 /)
  mpres@j_parent_start    = (/   1,  52,  52,  25,  10 /)
  mpres@e_we              = (/ 361, 460, 490, 101,  51/)
  mpres@e_sn              = (/ 181, 268, 238, 101,  51 /)
  mpres@dx                = .85
  mpres@dy                = .85
  mpres@map_proj          = "lat-lon"
  mpres@ref_lat           =  0.00
  mpres@ref_lon           =  180.00
  mpres@truelat1          =  0.0
  mpres@truelat2          =  0.0
  mpres@stand_lon         =  100.0
  mpres@pole_lat = 45.0
  mpres@pole_lon = -100.0

  mp = wrf_wps_dom (wks,mpres,lnres,txres)

  frame(wks)


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

end