load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

begin

  mres                    = True                    ; plot mods desired
  mres@mpMinLatF           = 54.5 ;55.5
  mres@mpMaxLatF           = 55.8 ;57.0           ; choose a different subregion
  mres@mpMinLonF           = 13.8 ;8.0
  mres@mpMaxLonF           = 16. ;9.0
  
  mres@mpDataSetName         = "Earth..4"   ; This new database contains
                                           ; divisions for other countries.
  resltn = (/"Coarsest","Coarse","Medium","Fine","Finest"/)

  mres@mpDataBaseVersion     = "HighRes"    
  mres@mpDataResolution = resltn(3)
  wks = gsn_open_wks("pdf","HovsoreMap")
  gsn_define_colormap(wks,"amwg")

  mres@gsnDraw   = False                             ; don't draw yet
  mres@gsnFrame  = False                             ; don't advance yet

  base_map = gsn_csm_map(wks,mres)

  txres                       = True
  txres@txPerimOn             = True
  txres@txBackgroundFillColor = "Salmon"
  txres@txFontHeightF         = 0.03

  txid = gsn_create_text(wks, "default text settings", txres)

  amres = True
  x_out = new(1,float)
  y_out = new(1,float)
  datatondc(base_map,15.,55.,x_out,y_out)
  print(x_out+"  "+y_out)
  amres@amParallelPosF    = x_out
  amres@amOrthogonalPosF  = y_out
  amres@amJust           = "CenterCenter"
  annoid = gsn_add_annotation(base_map, txid, amres)  ; Attach string to plot
                                                  ; using default values.
  draw(base_map)
  frame(wks)

end