;*******************************************************************************
;  Pintar solo los datos de terreno sobre el mapa, con la máxima resolución de
;  costas y en blanco y negro
;  Modificado para BL Meteor
;*******************************************************************************

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

begin

  a = addfile("./geo_em.d03.nc","r")                              ; Open a file


;   type = "x11"
 type = "pdf"
; type = "ps"
; type = "eps"
; type = "png"
 wks = gsn_open_wks(type,"Domain03_bn")            ; Create a plot workstation
 gsn_define_colormap(wks,"gsdtol")          ; Paleta de colores para ter en BN
 gsn_reverse_colormap(wks)                     ; Invertir la escala de colores
;  gsn_draw_colormap(wks)

  opts = True                                    ; Set some Basic Plot options
  opts@MainTitle = "Island of Fuerteventura, Domain 03"
  opts@InitTime = False                          ; Do not plot time or footers
  opts@Footer = False    

  ter = wrf_user_getvar(a,"ter",0)               ; Read the variable to memory 
;***************** opciones para terrain ******************************************
  res = opts                                   ; Use basic options for this field
  res@cnFillOn = True                          ; Create a color fill plot
  res@cnLevelSelectionMode = "ExplicitLevels"
  res@cnLevels = (/ 25., 50., 100., 150., 200., 250., 300., 350., 400., \
                   450., 500., 550., 600., 650., 700., 750., 800. /) ; Set the levels
  res@gsnSpreadColorStart = 2
  res@gsnSpreadColorEnd   = 32
  contour = wrf_contour(a,wks,ter,res)
;**********************************************************************************
;****************** opciones del mapa *********************************************
  pltres = True                                                  ; Set plot options
  mpres = True                                                    ; Set map options
  mpres@gsnDraw            = False                                     ; don't draw
  mpres@gsnFrame           = False                            ; don't advance frame
  mpres@mpProjection = "Mercator"                                       ;Projection
  mpres@mpGeophysicalLineColor      = "Black"  
  mpres@mpGridLineColor             = "Black"
  mpres@mpPerimLineColor            = "Black"
  mpres@mpGeophysicalLineThicknessF = 3.0
  mpres@mpGridLineThicknessF        = 0.5
  mpres@mpGridSpacingF              = 0.5
  mpres@mpDataBaseVersion           = "HighRes"
  mpres@mpDataResolution             = "Fine" ; La linea de costa con mucho detalle
;**********************************************************************************
  plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres) ; Plot field over map

; Add marks where the stations are located
  resp = True                                                    ; Poly res desired
  resp@tfPolyDrawOrder = "PostDraw" 
  resp@gsMarkerColor = "Black"
  resp@gsMarkerIndex = 16
  resp@gsMarkerSizeF = 0.01
  dummy = gsn_add_polymarker(wks,plot,-14.012,28.432,resp)        ; Molino de Antigua, GC08
  
; Add the name of the station
  txres               = True                                    ; text mods desired
  txres@txFontHeightF = 0.015                                     ; text font height  
  gsn_text(wks,plot,"Antigua",-14.012,28.472,txres)
  
 draw(plot)
 frame(wks)
end