;********
srcFileName="tasmax_rcp85_land-rcm_uk_12km_time10.nc"

;---Read File
f=addfile(srcFileName,"r")
print(f)

;---Get the source lat/lon grid
tas   = f->tasmax(0,0,:,:)
lat2d = f->grid_latitude
lon2d = f->grid_longitude
printMinMax(lat2d,0)
printMinMax(lon2d,0)

itm   = f->transverse_mercator

print(itm)
print("itm@latitude_of_projection_origin="+itm@latitude_of_projection_origin)
print("itm@longitude_of_central_meridian="+itm@longitude_of_central_meridian)

lat_shift = itm@latitude_of_projection_origin
lon_shift = itm@longitude_of_central_meridian

lat2d = lat2d + lat_shift
lon2d = lon2d + lon_shift
printMinMax(lat2d,0)
printMinMax(lon2d,0)

tas@lat2d = lat2d             ; use 'shifted' coordinates
tas@lon2d = lon2d    

;---Plot data
  wks = gsn_open_wks("png","SBR_map_v1")

  res              = True
  res@gsnMaximize  = True     ; maximize plot in frame
  res@cnFillOn     = True     ; turn on contour fill
  res@cnLinesOn    = False    ; turn off contour fill
  res@gsnAddCyclic = False
  res@mpMinLatF    = min(lat2d)
  res@mpMaxLatF    = max(lat2d)
  res@mpMinLonF    = min(lon2d)
  res@mpMaxLonF    = max(lon2d)
  res@mpDataBaseVersion = "MediumRes"
  res@mpLandFillColor = "white"

print(res)

  plot = gsn_csm_contour_map(wks,tas,res)