load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/home/arhoades/plots/latlonplots/shapefile/shapefile_mask_data.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/cd_string.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"

begin
;******************************************************
;  Open a netCDF file containing the grid and data from the HOMME
;  (High-Order Multiscale Modeling Environment) model, formerly
;  called SEAM (Spectral Element Atmosphere Model).
;******************************************************

  setvalues NhlGetWorkspaceObjectId()
    "wsMaximumSize" : 300000000
  end setvalues

;------------------------------------------------------
;  Type of Plot
;------------------------------------------------------

; type = "x11"
 type = "pdf"
; type = "ps"
; type = "eps"
; type = "ncgm"

;------------------------------------------------------
;  Grab Files and Variables of Interest
;------------------------------------------------------

;UNIGRID 28 - rough topo (HIST)
; MASKED DATA
  fils_HIST_UNI28rough  = systemfunc("ls /home/arhoades/monthlydata/MASKED_SN/RCPvHIST/lnd/SN_mask_HISTandRCP85.nc")
  f_HIST_UNI28rough     = addfile(fils_HIST_UNI28rough,"r")

; UNMASKED DATA
;  fils_HIST_UNI28rough  = systemfunc("ls /home/arhoades/monthlydata/UNIGRID/25km_roughtopo/DJFavg/REGRID/*.avgDJF.box.nc")
;  f_HIST_UNI28rough     = addfiles(fils_HIST_UNI28rough,"r")

  lat_HIST_UNIGRID28_rough = fspan(20,50,240)
  lon_HIST_UNIGRID28_rough = fspan(-130,-100,240)
  HIST_SNOWFALL_UNIGRID28_rough = f_HIST_UNI28rough->DJF_SNOWFALL_HIST      ; 1 kg/m2 = 1 mm
  HIST_SNOWFALL_UNIGRID28_rough@_FillValue = 1e+36
;  HIST_SNOWFALL_UNIGRID28_rough = f_HIST_UNI28rough[:]->H2OSNO      ; 1 kg/m2 = 1 mm
;  HIST_SNOWFALL_UNIGRID28_rough = HIST_SNOWFALL_UNIGRID28_rough-0.1
  HIST_SNOWFALL_UNIGRID28_rough!0   = "time"
  HIST_SNOWFALL_UNIGRID28_rough!1   = "lat"
  HIST_SNOWFALL_UNIGRID28_rough!2   = "lon"
  HIST_SNOWFALL_UNIGRID28_rough&lat = lat_HIST_UNIGRID28_rough
  HIST_SNOWFALL_UNIGRID28_rough&lon = lon_HIST_UNIGRID28_rough
  HIST_SNOWFALL_UNIGRID28_rough&lat@units = "degrees_north"
  HIST_SNOWFALL_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_HIST_SNOWFALL_UNIGRID28_rough   = dimsizes(HIST_SNOWFALL_UNIGRID28_rough)
  nt_HIST_SNOWFALL_UNIGRID28_rough           = dsizes_HIST_SNOWFALL_UNIGRID28_rough(0)
  ny_HIST_SNOWFALL_UNIGRID28_rough           = dsizes_HIST_SNOWFALL_UNIGRID28_rough(1)
  nx_HIST_SNOWFALL_UNIGRID28_rough           = dsizes_HIST_SNOWFALL_UNIGRID28_rough(2)
  printVarSummary(HIST_SNOWFALL_UNIGRID28_rough)
  print ("Minimum="+min(HIST_SNOWFALL_UNIGRID28_rough) +"    Maximum="+max(HIST_SNOWFALL_UNIGRID28_rough))

  lat_HIST_UNIGRID28_rough = fspan(20,50,240)
  lon_HIST_UNIGRID28_rough = fspan(-130,-100,240)
  HIST_SNOWC_UNIGRID28_rough = f_HIST_UNI28rough->DJF_SNOWC_HIST      ; 1 kg/m2 = 1 mm
  HIST_SNOWC_UNIGRID28_rough@_FillValue = 1e+36
;  HIST_SNOWC_UNIGRID28_rough = f_HIST_UNI28rough[:]->FSNO      ; 1 kg/m2 = 1 mm
;  HIST_SNOWC_UNIGRID28_rough = HIST_SNOWC_UNIGRID28_rough-0.001
;  HIST_SNOWC_UNIGRID28_rough = HIST_SNOWC_UNIGRID28_rough*100
  HIST_SNOWC_UNIGRID28_rough!0   = "time"
  HIST_SNOWC_UNIGRID28_rough!1   = "lat"
  HIST_SNOWC_UNIGRID28_rough!2   = "lon"
  HIST_SNOWC_UNIGRID28_rough&lat = lat_HIST_UNIGRID28_rough
  HIST_SNOWC_UNIGRID28_rough&lon = lon_HIST_UNIGRID28_rough
  HIST_SNOWC_UNIGRID28_rough&lat@units = "degrees_north"
  HIST_SNOWC_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_HIST_SNOWC_UNIGRID28_rough   = dimsizes(HIST_SNOWC_UNIGRID28_rough)
  nt_HIST_SNOWC_UNIGRID28_rough           = dsizes_HIST_SNOWC_UNIGRID28_rough(0)
  ny_HIST_SNOWC_UNIGRID28_rough           = dsizes_HIST_SNOWC_UNIGRID28_rough(1)
  nx_HIST_SNOWC_UNIGRID28_rough           = dsizes_HIST_SNOWC_UNIGRID28_rough(2)
  printVarSummary(HIST_SNOWC_UNIGRID28_rough)
  print ("Minimum="+min(HIST_SNOWC_UNIGRID28_rough) +"    Maximum="+max(HIST_SNOWC_UNIGRID28_rough))

  lat_HIST_UNIGRID28_rough = fspan(20,50,240)
  lon_HIST_UNIGRID28_rough = fspan(-130,-100,240)
  HIST_SWE_UNIGRID28_rough = f_HIST_UNI28rough->DJF_SWE_HIST      ; 1 kg/m2 = 1 mm
  HIST_SWE_UNIGRID28_rough@_FillValue = 1e+36
;  HIST_SWE_UNIGRID28_rough = f_HIST_UNI28rough[:]->H2OSNO      ; 1 kg/m2 = 1 mm
  HIST_SWE_UNIGRID28_rough!0   = "time"
  HIST_SWE_UNIGRID28_rough!1   = "lat"
  HIST_SWE_UNIGRID28_rough!2   = "lon"
  HIST_SWE_UNIGRID28_rough&lat = lat_HIST_UNIGRID28_rough
  HIST_SWE_UNIGRID28_rough&lon = lon_HIST_UNIGRID28_rough
  HIST_SWE_UNIGRID28_rough&lat@units = "degrees_north"
  HIST_SWE_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_HIST_SWE_UNIGRID28_rough   = dimsizes(HIST_SWE_UNIGRID28_rough)
  nt_HIST_SWE_UNIGRID28_rough           = dsizes_HIST_SWE_UNIGRID28_rough(0)
  ny_HIST_SWE_UNIGRID28_rough           = dsizes_HIST_SWE_UNIGRID28_rough(1)
  nx_HIST_SWE_UNIGRID28_rough           = dsizes_HIST_SWE_UNIGRID28_rough(2)
  printVarSummary(HIST_SWE_UNIGRID28_rough)
  print ("Minimum="+min(HIST_SWE_UNIGRID28_rough) +"    Maximum="+max(HIST_SWE_UNIGRID28_rough))

  lat_HIST_UNIGRID28_rough = fspan(20,50,240)
  lon_HIST_UNIGRID28_rough = fspan(-130,-100,240)
  HIST_2mST_UNIGRID28_rough = f_HIST_UNI28rough->DJF_2mST_HIST      ; 1 kg/m2 = 1 mm
  HIST_2mST_UNIGRID28_rough@_FillValue = 1e+36
;  HIST_2mST_UNIGRID28_rough = f_HIST_UNI28rough[:]->TSA      ; 1 kg/m2 = 1 mm
  HIST_2mST_UNIGRID28_rough!0   = "time"
  HIST_2mST_UNIGRID28_rough!1   = "lat"
  HIST_2mST_UNIGRID28_rough!2   = "lon"
  HIST_2mST_UNIGRID28_rough&lat = lat_HIST_UNIGRID28_rough
  HIST_2mST_UNIGRID28_rough&lon = lon_HIST_UNIGRID28_rough
  HIST_2mST_UNIGRID28_rough&lat@units = "degrees_north"
  HIST_2mST_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_HIST_2mST_UNIGRID28_rough   = dimsizes(HIST_2mST_UNIGRID28_rough)
  nt_HIST_2mST_UNIGRID28_rough           = dsizes_HIST_2mST_UNIGRID28_rough(0)
  ny_HIST_2mST_UNIGRID28_rough           = dsizes_HIST_2mST_UNIGRID28_rough(1)
  nx_HIST_2mST_UNIGRID28_rough           = dsizes_HIST_2mST_UNIGRID28_rough(2)
  printVarSummary(HIST_2mST_UNIGRID28_rough)
  print ("Minimum="+min(HIST_2mST_UNIGRID28_rough) +"    Maximum="+max(HIST_2mST_UNIGRID28_rough))

;UNIGRID 28 - rough topo (RCP_MIDCENT8.5 - MIDCENT)
; MASKED DATA
  fils_RCP85_MIDCENT_UNI28rough  = systemfunc("ls /home/arhoades/monthlydata/MASKED_SN/RCPvHIST/lnd/SN_mask_HISTandRCP85.nc")
  f_RCP85_MIDCENT_UNI28rough     = addfile(fils_RCP85_MIDCENT_UNI28rough,"r")

; UNMASKED DATA
;  fils_RCP85_MIDCENT_UNI28rough  = systemfunc("ls /group/paullricgrp2/VARRES28km_RCP8.5_MIDCENT/pre-process/lnd/DJFavg/REGRID/*.avgDJF.box.nc")
;  f_RCP85_MIDCENT_UNI28rough     = addfiles(fils_RCP85_MIDCENT_UNI28rough,"r")

  lat_RCP85_MIDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_MIDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough = f_RCP85_MIDCENT_UNI28rough->DJF_SNOWFALL_RCP85_MIDCENT      ; 1 kg/m2 = 1 mm
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough@_FillValue = 1e+36  
;  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough = f_RCP85_MIDCENT_UNI28rough[:]->H2OSNO      ; 1 kg/m2 = 1 mm
;  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough = RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough-0.1
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough!0   = "time"
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough!1   = "lat"
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough!2   = "lon"
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough&lat = lat_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough&lon = lon_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough   = dimsizes(RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough)
  nt_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough(0)
  ny_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough(1)
  nx_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough(2)
  printVarSummary(RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough) +"    Maximum="+max(RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough))

  lat_RCP85_MIDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_MIDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough = f_RCP85_MIDCENT_UNI28rough->DJF_SNOWC_RCP85_MIDCENT      ; 1 kg/m2 = 1 mm
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough@_FillValue = 1e+36
;  RCP85_MIDCENT_SNOWC_UNIGRID28_rough = f_RCP85_MIDCENT_UNI28rough[:]->FSNO      ; 1 kg/m2 = 1 mm
;  RCP85_MIDCENT_SNOWC_UNIGRID28_rough = RCP85_MIDCENT_SNOWC_UNIGRID28_rough-0.001
;  RCP85_MIDCENT_SNOWC_UNIGRID28_rough = RCP85_MIDCENT_SNOWC_UNIGRID28_rough*100
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough!0   = "time"
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough!1   = "lat"
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough!2   = "lon"
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough&lat = lat_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough&lon = lon_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_RCP85_MIDCENT_SNOWC_UNIGRID28_rough   = dimsizes(RCP85_MIDCENT_SNOWC_UNIGRID28_rough)
  nt_RCP85_MIDCENT_SNOWC_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SNOWC_UNIGRID28_rough(0)
  ny_RCP85_MIDCENT_SNOWC_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SNOWC_UNIGRID28_rough(1)
  nx_RCP85_MIDCENT_SNOWC_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SNOWC_UNIGRID28_rough(2)
  printVarSummary(RCP85_MIDCENT_SNOWC_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_MIDCENT_SNOWC_UNIGRID28_rough) +"    Maximum="+max(RCP85_MIDCENT_SNOWC_UNIGRID28_rough))

  lat_RCP85_MIDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_MIDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_MIDCENT_SWE_UNIGRID28_rough = f_RCP85_MIDCENT_UNI28rough->DJF_SWE_RCP85_MIDCENT      ; 1 kg/m2 = 1 mm
  RCP85_MIDCENT_SWE_UNIGRID28_rough@_FillValue = 1e+36
;  RCP85_MIDCENT_SWE_UNIGRID28_rough = f_RCP85_MIDCENT_UNI28rough[:]->H2OSNO      ; 1 kg/m2 = 1 mm
  RCP85_MIDCENT_SWE_UNIGRID28_rough!0   = "time"
  RCP85_MIDCENT_SWE_UNIGRID28_rough!1   = "lat"
  RCP85_MIDCENT_SWE_UNIGRID28_rough!2   = "lon"
  RCP85_MIDCENT_SWE_UNIGRID28_rough&lat = lat_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENT_SWE_UNIGRID28_rough&lon = lon_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENT_SWE_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_MIDCENT_SWE_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_RCP85_MIDCENT_SWE_UNIGRID28_rough   = dimsizes(RCP85_MIDCENT_SWE_UNIGRID28_rough)
  nt_RCP85_MIDCENT_SWE_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SWE_UNIGRID28_rough(0)
  ny_RCP85_MIDCENT_SWE_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SWE_UNIGRID28_rough(1)
  nx_RCP85_MIDCENT_SWE_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_SWE_UNIGRID28_rough(2)
  printVarSummary(RCP85_MIDCENT_SWE_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_MIDCENT_SWE_UNIGRID28_rough) +"    Maximum="+max(RCP85_MIDCENT_SWE_UNIGRID28_rough))

  lat_RCP85_MIDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_MIDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_MIDCENT_2mST_UNIGRID28_rough = f_RCP85_MIDCENT_UNI28rough->DJF_2mST_RCP85_MIDCENT      ; 1 kg/m2 = 1 mm
  RCP85_MIDCENT_2mST_UNIGRID28_rough@_FillValue = 1e+36
;  RCP85_MIDCENT_2mST_UNIGRID28_rough = f_RCP85_MIDCENT_UNI28rough[:]->TSA      ; 1 kg/m2 = 1 mm
  RCP85_MIDCENT_2mST_UNIGRID28_rough!0   = "time"
  RCP85_MIDCENT_2mST_UNIGRID28_rough!1   = "lat"
  RCP85_MIDCENT_2mST_UNIGRID28_rough!2   = "lon"
  RCP85_MIDCENT_2mST_UNIGRID28_rough&lat = lat_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENT_2mST_UNIGRID28_rough&lon = lon_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENT_2mST_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_MIDCENT_2mST_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_RCP85_MIDCENT_2mST_UNIGRID28_rough   = dimsizes(RCP85_MIDCENT_2mST_UNIGRID28_rough)
  nt_RCP85_MIDCENT_2mST_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_2mST_UNIGRID28_rough(0)
  ny_RCP85_MIDCENT_2mST_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_2mST_UNIGRID28_rough(1)
  nx_RCP85_MIDCENT_2mST_UNIGRID28_rough           = dsizes_RCP85_MIDCENT_2mST_UNIGRID28_rough(2)
  printVarSummary(RCP85_MIDCENT_2mST_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_MIDCENT_2mST_UNIGRID28_rough) +"    Maximum="+max(RCP85_MIDCENT_2mST_UNIGRID28_rough))

;UNIGRID 28 - rough topo (RCP_ENDCENT8.5 - ENDCENT)
; MASKED DATA
  fils_RCP85_ENDCENT_UNI28rough  = systemfunc("ls /home/arhoades/monthlydata/MASKED_SN/RCPvHIST/lnd/SN_mask_HISTandRCP85.nc")
  f_RCP85_ENDCENT_UNI28rough     = addfile(fils_RCP85_ENDCENT_UNI28rough,"r")

; UNMASKED DATA
;  fils_RCP85_ENDCENT_UNI28rough  = systemfunc("ls /group/paullricgrp2/VR-CESM_0.25_RCP85_ENDCENT/pre-process/lnd/DJFavg/REGRID/*.avgDJF.box.nc")
;  f_RCP85_ENDCENT_UNI28rough     = addfiles(fils_RCP85_ENDCENT_UNI28rough,"r")

  lat_RCP85_ENDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_ENDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough = f_RCP85_ENDCENT_UNI28rough->DJF_SNOWFALL_RCP85_ENDCENT      ; 1 kg/m2 = 1 mm
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough@_FillValue = 1e+36
;  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough = f_RCP85_ENDCENT_UNI28rough[:]->H2OSNO      ; 1 kg/m2 = 1 mm
;  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough = RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough-0.1
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough!0   = "time"
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough!1   = "lat"
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough!2   = "lon"
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough&lat = lat_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough&lon = lon_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough   = dimsizes(RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough)
  nt_RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough(0)
  ny_RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough(1)
  nx_RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough(2)
  printVarSummary(RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough) +"    Maximum="+max(RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough))

  lat_RCP85_ENDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_ENDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough = f_RCP85_ENDCENT_UNI28rough->DJF_SNOWC_RCP85_ENDCENT      ; 1 kg/m2 = 1 mm
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough@_FillValue = 1e+36
;  RCP85_ENDCENT_SNOWC_UNIGRID28_rough = f_RCP85_ENDCENT_UNI28rough[:]->FSNO      ; 1 kg/m2 = 1 mm
;  RCP85_ENDCENT_SNOWC_UNIGRID28_rough = RCP85_ENDCENT_SNOWC_UNIGRID28_rough-0.001
;  RCP85_ENDCENT_SNOWC_UNIGRID28_rough = RCP85_ENDCENT_SNOWC_UNIGRID28_rough*100
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough!0   = "time"
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough!1   = "lat"
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough!2   = "lon"
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough&lat = lat_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough&lon = lon_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_RCP85_ENDCENT_SNOWC_UNIGRID28_rough   = dimsizes(RCP85_ENDCENT_SNOWC_UNIGRID28_rough)
  nt_RCP85_ENDCENT_SNOWC_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SNOWC_UNIGRID28_rough(0)
  ny_RCP85_ENDCENT_SNOWC_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SNOWC_UNIGRID28_rough(1)
  nx_RCP85_ENDCENT_SNOWC_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SNOWC_UNIGRID28_rough(2)
  printVarSummary(RCP85_ENDCENT_SNOWC_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_ENDCENT_SNOWC_UNIGRID28_rough) +"    Maximum="+max(RCP85_ENDCENT_SNOWC_UNIGRID28_rough))

  lat_RCP85_ENDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_ENDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_ENDCENT_SWE_UNIGRID28_rough = f_RCP85_ENDCENT_UNI28rough->DJF_SWE_RCP85_ENDCENT      ; 1 kg/m2 = 1 mm
  RCP85_ENDCENT_SWE_UNIGRID28_rough@_FillValue = 1e+36
;  RCP85_ENDCENT_SWE_UNIGRID28_rough = f_RCP85_ENDCENT_UNI28rough[:]->H2OSNO      ; 1 kg/m2 = 1 mm
  RCP85_ENDCENT_SWE_UNIGRID28_rough!0   = "time"
  RCP85_ENDCENT_SWE_UNIGRID28_rough!1   = "lat"
  RCP85_ENDCENT_SWE_UNIGRID28_rough!2   = "lon"
  RCP85_ENDCENT_SWE_UNIGRID28_rough&lat = lat_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENT_SWE_UNIGRID28_rough&lon = lon_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENT_SWE_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_ENDCENT_SWE_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_RCP85_ENDCENT_SWE_UNIGRID28_rough   = dimsizes(RCP85_ENDCENT_SWE_UNIGRID28_rough)
  nt_RCP85_ENDCENT_SWE_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SWE_UNIGRID28_rough(0)
  ny_RCP85_ENDCENT_SWE_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SWE_UNIGRID28_rough(1)
  nx_RCP85_ENDCENT_SWE_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_SWE_UNIGRID28_rough(2)
  printVarSummary(RCP85_ENDCENT_SWE_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_ENDCENT_SWE_UNIGRID28_rough) +"    Maximum="+max(RCP85_ENDCENT_SWE_UNIGRID28_rough))

  lat_RCP85_ENDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_ENDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_ENDCENT_2mST_UNIGRID28_rough = f_RCP85_ENDCENT_UNI28rough->DJF_2mST_RCP85_ENDCENT      ; 1 kg/m2 = 1 mm
  RCP85_ENDCENT_2mST_UNIGRID28_rough@_FillValue = 1e+36
;  RCP85_ENDCENT_2mST_UNIGRID28_rough = f_RCP85_ENDCENT_UNI28rough[:]->TSA      ; 1 kg/m2 = 1 mm
  RCP85_ENDCENT_2mST_UNIGRID28_rough!0   = "time"
  RCP85_ENDCENT_2mST_UNIGRID28_rough!1   = "lat"
  RCP85_ENDCENT_2mST_UNIGRID28_rough!2   = "lon"
  RCP85_ENDCENT_2mST_UNIGRID28_rough&lat = lat_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENT_2mST_UNIGRID28_rough&lon = lon_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENT_2mST_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_ENDCENT_2mST_UNIGRID28_rough&lon@units = "degrees_east"
  dsizes_RCP85_ENDCENT_2mST_UNIGRID28_rough   = dimsizes(RCP85_ENDCENT_2mST_UNIGRID28_rough)
  nt_RCP85_ENDCENT_2mST_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_2mST_UNIGRID28_rough(0)
  ny_RCP85_ENDCENT_2mST_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_2mST_UNIGRID28_rough(1)
  nx_RCP85_ENDCENT_2mST_UNIGRID28_rough           = dsizes_RCP85_ENDCENT_2mST_UNIGRID28_rough(2)
  printVarSummary(RCP85_ENDCENT_2mST_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_ENDCENT_2mST_UNIGRID28_rough) +"    Maximum="+max(RCP85_ENDCENT_2mST_UNIGRID28_rough))

;----------------------
; BASE CASE
;----------------------
  HIST_SNOWFALL_UNIGRID28_rough_time = HIST_SNOWFALL_UNIGRID28_rough(lat|:,lon|:,time|:)
  HIST_SNOWFALL_UNIGRID28_rough_base = dim_avg(HIST_SNOWFALL_UNIGRID28_rough_time)
  printVarSummary(HIST_SNOWFALL_UNIGRID28_rough_base)
  print ("Minimum="+min(HIST_SNOWFALL_UNIGRID28_rough_base) +"    Maximum="+max(HIST_SNOWFALL_UNIGRID28_rough_base))

  HIST_SNOWC_UNIGRID28_rough_time = HIST_SNOWC_UNIGRID28_rough(lat|:,lon|:,time|:)
  HIST_SNOWC_UNIGRID28_rough_base = dim_avg(HIST_SNOWC_UNIGRID28_rough_time)
  printVarSummary(HIST_SNOWC_UNIGRID28_rough_base)
  print ("Minimum="+min(HIST_SNOWC_UNIGRID28_rough_base) +"    Maximum="+max(HIST_SNOWC_UNIGRID28_rough_base))

  HIST_SWE_UNIGRID28_rough_time = HIST_SWE_UNIGRID28_rough(lat|:,lon|:,time|:)
  HIST_SWE_UNIGRID28_rough_base = dim_avg(HIST_SWE_UNIGRID28_rough_time)
  printVarSummary(HIST_SWE_UNIGRID28_rough_base)
  print ("Minimum="+min(HIST_SWE_UNIGRID28_rough_base) +"    Maximum="+max(HIST_SWE_UNIGRID28_rough_base))

  HIST_2mST_UNIGRID28_rough_time = HIST_2mST_UNIGRID28_rough(lat|:,lon|:,time|:)
  HIST_2mST_UNIGRID28_rough_base = dim_avg(HIST_2mST_UNIGRID28_rough_time)
  printVarSummary(HIST_2mST_UNIGRID28_rough_base)
  print ("Minimum="+min(HIST_2mST_UNIGRID28_rough_base) +"    Maximum="+max(HIST_2mST_UNIGRID28_rough_base))

;--------------------------------
; DIFFERENCE - RCP85 - BASECASE
;--------------------------------
  HISTvHIST_SNOWFALL_UNIGRID28_rough = new((/nt_HIST_SNOWFALL_UNIGRID28_rough,ny_HIST_SNOWFALL_UNIGRID28_rough,nx_HIST_SNOWFALL_UNIGRID28_rough/),float)

  do i=0,nt_HIST_SNOWFALL_UNIGRID28_rough-1
    HISTvHIST_SNOWFALL_UNIGRID28_rough(i,:,:) = HIST_SNOWFALL_UNIGRID28_rough(i,:,:) - HIST_SNOWFALL_UNIGRID28_rough_base
    printVarSummary(HISTvHIST_SNOWFALL_UNIGRID28_rough)
    printMinMax(HISTvHIST_SNOWFALL_UNIGRID28_rough,True)
  end do

  lat_HIST_UNIGRID28_rough = fspan(20,50,240)
  lon_HIST_UNIGRID28_rough = fspan(-130,-100,240)
  HISTvHIST_SNOWFALL_UNIGRID28_rough!0   = "time"
  HISTvHIST_SNOWFALL_UNIGRID28_rough!1   = "lat"
  HISTvHIST_SNOWFALL_UNIGRID28_rough!2   = "lon"
  HISTvHIST_SNOWFALL_UNIGRID28_rough&lat = lat_HIST_UNIGRID28_rough
  HISTvHIST_SNOWFALL_UNIGRID28_rough&lon = lon_HIST_UNIGRID28_rough 
  HISTvHIST_SNOWFALL_UNIGRID28_rough&lat@units = "degrees_north"
  HISTvHIST_SNOWFALL_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(HISTvHIST_SNOWFALL_UNIGRID28_rough)
  print ("Minimum="+min(HISTvHIST_SNOWFALL_UNIGRID28_rough) +"    Maximum="+max(HISTvHIST_SNOWFALL_UNIGRID28_rough))

  HISTvHIST_SNOWC_UNIGRID28_rough = new((/nt_HIST_SNOWC_UNIGRID28_rough,ny_HIST_SNOWC_UNIGRID28_rough,nx_HIST_SNOWC_UNIGRID28_rough/),float)

  do i=0,nt_HIST_SNOWC_UNIGRID28_rough-1
    HISTvHIST_SNOWC_UNIGRID28_rough(i,:,:) = HIST_SNOWC_UNIGRID28_rough(i,:,:) - HIST_SNOWC_UNIGRID28_rough_base
    printVarSummary(HISTvHIST_SNOWC_UNIGRID28_rough)
    printMinMax(HISTvHIST_SNOWC_UNIGRID28_rough,True)
  end do

  lat_HIST_UNIGRID28_rough = fspan(20,50,240)
  lon_HIST_UNIGRID28_rough = fspan(-130,-100,240)
  HISTvHIST_SNOWC_UNIGRID28_rough!0   = "time"
  HISTvHIST_SNOWC_UNIGRID28_rough!1   = "lat"
  HISTvHIST_SNOWC_UNIGRID28_rough!2   = "lon"
  HISTvHIST_SNOWC_UNIGRID28_rough&lat = lat_HIST_UNIGRID28_rough
  HISTvHIST_SNOWC_UNIGRID28_rough&lon = lon_HIST_UNIGRID28_rough
  HISTvHIST_SNOWC_UNIGRID28_rough&lat@units = "degrees_north"
  HISTvHIST_SNOWC_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(HISTvHIST_SNOWC_UNIGRID28_rough)
  print ("Minimum="+min(HISTvHIST_SNOWC_UNIGRID28_rough) +"    Maximum="+max(HISTvHIST_SNOWC_UNIGRID28_rough))

  HISTvHIST_SWE_UNIGRID28_rough = new((/nt_HIST_SWE_UNIGRID28_rough,ny_HIST_SWE_UNIGRID28_rough,nx_HIST_SWE_UNIGRID28_rough/),float)

  do i=0,nt_HIST_SWE_UNIGRID28_rough-1
    HISTvHIST_SWE_UNIGRID28_rough(i,:,:) = HIST_SWE_UNIGRID28_rough(i,:,:) - HIST_SWE_UNIGRID28_rough_base
    printVarSummary(HISTvHIST_SWE_UNIGRID28_rough)
    printMinMax(HISTvHIST_SWE_UNIGRID28_rough,True)
  end do

  lat_HIST_UNIGRID28_rough = fspan(20,50,240)
  lon_HIST_UNIGRID28_rough = fspan(-130,-100,240)
  HISTvHIST_SWE_UNIGRID28_rough!0   = "time"
  HISTvHIST_SWE_UNIGRID28_rough!1   = "lat"
  HISTvHIST_SWE_UNIGRID28_rough!2   = "lon"
  HISTvHIST_SWE_UNIGRID28_rough&lat = lat_HIST_UNIGRID28_rough
  HISTvHIST_SWE_UNIGRID28_rough&lon = lon_HIST_UNIGRID28_rough
  HISTvHIST_SWE_UNIGRID28_rough&lat@units = "degrees_north"
  HISTvHIST_SWE_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(HISTvHIST_SWE_UNIGRID28_rough)
  print ("Minimum="+min(HISTvHIST_SWE_UNIGRID28_rough) +"    Maximum="+max(HISTvHIST_SWE_UNIGRID28_rough))

  HISTvHIST_2mST_UNIGRID28_rough = new((/nt_HIST_2mST_UNIGRID28_rough,ny_HIST_2mST_UNIGRID28_rough,nx_HIST_2mST_UNIGRID28_rough/),float)

  do i=0,nt_HIST_2mST_UNIGRID28_rough-1
    HISTvHIST_2mST_UNIGRID28_rough(i,:,:) = HIST_2mST_UNIGRID28_rough(i,:,:) - HIST_2mST_UNIGRID28_rough_base
    printVarSummary(HISTvHIST_2mST_UNIGRID28_rough)
    printMinMax(HISTvHIST_2mST_UNIGRID28_rough,True)
  end do

  lat_HIST_UNIGRID28_rough = fspan(20,50,240)
  lon_HIST_UNIGRID28_rough = fspan(-130,-100,240)
  HISTvHIST_2mST_UNIGRID28_rough!0   = "time"
  HISTvHIST_2mST_UNIGRID28_rough!1   = "lat"
  HISTvHIST_2mST_UNIGRID28_rough!2   = "lon"
  HISTvHIST_2mST_UNIGRID28_rough&lat = lat_HIST_UNIGRID28_rough
  HISTvHIST_2mST_UNIGRID28_rough&lon = lon_HIST_UNIGRID28_rough
  HISTvHIST_2mST_UNIGRID28_rough&lat@units = "degrees_north"
  HISTvHIST_2mST_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(HISTvHIST_2mST_UNIGRID28_rough)
  print ("Minimum="+min(HISTvHIST_2mST_UNIGRID28_rough) +"    Maximum="+max(HISTvHIST_2mST_UNIGRID28_rough))

; RCP85 MIDCENT

  RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough = new((/nt_HIST_SNOWFALL_UNIGRID28_rough,ny_HIST_SNOWFALL_UNIGRID28_rough,nx_HIST_SNOWFALL_UNIGRID28_rough/),float)

  do i=0,nt_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough-1
    RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough(i,:,:) = RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough(i,:,:) - HIST_SNOWFALL_UNIGRID28_rough_base
    printVarSummary(RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough)
    printMinMax(RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough,True)
  end do
     
  lat_RCP85_MIDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_MIDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough!0   = "time"
  RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough!1   = "lat"
  RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough!2   = "lon"
  RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough&lat = lat_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough&lon = lon_RCP85_MIDCENT_UNIGRID28_rough    
  RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough) +"    Maximum="+max(RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough))

  RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough = new((/nt_HIST_SNOWC_UNIGRID28_rough,ny_HIST_SNOWC_UNIGRID28_rough,nx_HIST_SNOWC_UNIGRID28_rough/),float)

  do i=0,nt_RCP85_MIDCENT_SNOWC_UNIGRID28_rough-1
    RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough(i,:,:) = RCP85_MIDCENT_SNOWC_UNIGRID28_rough(i,:,:) - HIST_SNOWC_UNIGRID28_rough_base
    printVarSummary(RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough)
    printMinMax(RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough,True)
  end do

  lat_RCP85_MIDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_MIDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough!0   = "time"
  RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough!1   = "lat"
  RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough!2   = "lon"
  RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough&lat = lat_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough&lon = lon_RCP85_MIDCENT_UNIGRID28_rough 
  RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough) +"    Maximum="+max(RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough))

  RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough = new((/nt_HIST_SWE_UNIGRID28_rough,ny_HIST_SWE_UNIGRID28_rough,nx_HIST_SWE_UNIGRID28_rough/),float)

  do i=0,nt_RCP85_MIDCENT_SWE_UNIGRID28_rough-1
    RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough(i,:,:) = RCP85_MIDCENT_SWE_UNIGRID28_rough(i,:,:) - HIST_SWE_UNIGRID28_rough_base
    printVarSummary(RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough)
    printMinMax(RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough,True)
  end do

  lat_RCP85_MIDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_MIDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough!0   = "time"
  RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough!1   = "lat"
  RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough!2   = "lon"
  RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough&lat = lat_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough&lon = lon_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough) +"    Maximum="+max(RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough))

  RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough = new((/nt_HIST_2mST_UNIGRID28_rough,ny_HIST_2mST_UNIGRID28_rough,nx_HIST_2mST_UNIGRID28_rough/),float)

  do i=0,nt_RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough-1
    RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough(i,:,:) = RCP85_MIDCENT_2mST_UNIGRID28_rough(i,:,:) - HIST_2mST_UNIGRID28_rough_base
    printVarSummary(RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough)
    printMinMax(RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough,True)
  end do

  lat_RCP85_MIDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_MIDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough!0   = "time"
  RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough!1   = "lat"
  RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough!2   = "lon"
  RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough&lat = lat_RCP85_MIDCENT_UNIGRID28_rough
  RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough&lon = lon_RCP85_MIDCENT_UNIGRID28_rough 
  RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough) +"    Maximum="+max(RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough))

; RCP85 END CENT

  RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough = new((/nt_HIST_SNOWFALL_UNIGRID28_rough,ny_HIST_SNOWFALL_UNIGRID28_rough,nx_HIST_SNOWFALL_UNIGRID28_rough/),float)

  do i=0,nt_RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough-1
    RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough(i,:,:) = RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough(i,:,:) - HIST_SNOWFALL_UNIGRID28_rough_base
    printVarSummary(RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough)
    printMinMax(RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough,True)
  end do

  lat_RCP85_ENDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_ENDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough!0   = "time"
  RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough!1   = "lat"
  RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough!2   = "lon"
  RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough&lat = lat_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough&lon = lon_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough) +"    Maximum="+max(RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough))

  RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough = new((/nt_HIST_SNOWC_UNIGRID28_rough,ny_HIST_SNOWC_UNIGRID28_rough,nx_HIST_SNOWC_UNIGRID28_rough/),float)

  do i=0,nt_RCP85_ENDCENT_SNOWC_UNIGRID28_rough-1
    RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough(i,:,:) = RCP85_ENDCENT_SNOWC_UNIGRID28_rough(i,:,:) - HIST_SNOWC_UNIGRID28_rough_base
    printVarSummary(RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough)
    printMinMax(RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough,True)
  end do

  lat_RCP85_ENDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_ENDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough!0   = "time"
  RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough!1   = "lat"
  RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough!2   = "lon"
  RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough&lat = lat_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough&lon = lon_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough&lat@units = "degrees_north";
  RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough) +"    Maximum="+max(RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough))

  RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough = new((/nt_HIST_SWE_UNIGRID28_rough,ny_HIST_SWE_UNIGRID28_rough,nx_HIST_SWE_UNIGRID28_rough/),float)

  do i=0,nt_RCP85_ENDCENT_SWE_UNIGRID28_rough-1
    RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough(i,:,:) = RCP85_ENDCENT_SWE_UNIGRID28_rough(i,:,:) - HIST_SWE_UNIGRID28_rough_base
    printVarSummary(RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough)
    printMinMax(RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough,True)
  end do

  lat_RCP85_ENDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_ENDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough!0   = "time"
  RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough!1   = "lat"
  RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough!2   = "lon"
  RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough&lat = lat_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough&lon = lon_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough&lat@units = "degrees_north";
  RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough) +"    Maximum="+max(RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough))

  RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough = new((/nt_HIST_2mST_UNIGRID28_rough,ny_HIST_2mST_UNIGRID28_rough,nx_HIST_2mST_UNIGRID28_rough/),float)

  do i=0,nt_RCP85_ENDCENT_2mST_UNIGRID28_rough-1
    RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough(i,:,:) = RCP85_ENDCENT_2mST_UNIGRID28_rough(i,:,:) - HIST_2mST_UNIGRID28_rough_base
    printVarSummary(RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough)
    printMinMax(RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough,True)
  end do

  lat_RCP85_ENDCENT_UNIGRID28_rough = fspan(20,50,240)
  lon_RCP85_ENDCENT_UNIGRID28_rough = fspan(-130,-100,240)
  RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough!0   = "time"
  RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough!1   = "lat"
  RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough!2   = "lon"
  RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough&lat = lat_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough&lon = lon_RCP85_ENDCENT_UNIGRID28_rough
  RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough&lat@units = "degrees_north"
  RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough&lon@units = "degrees_east"
  printVarSummary(RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough)
  print ("Minimum="+min(RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough) +"    Maximum="+max(RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough)) 

;------------------------------------------------------
;  Titles of plot, text strings, etc.
;------------------------------------------------------
 
; plottitle = "UNIGRID28rough SNOWFALL DJF Climate Average - Units (mm)"
; plottitle2 = "SNOWFALL DJF Climate Average - Diff (2025-2050 - 1980-2005) - Units (mm)"
 filestring = "SN_HOV_1980-2100_DJF_ALL_VRGCM"

;----------------------------------------
;   Open Workstation and Define Colormap
;----------------------------------------

  wks = gsn_open_wks(type,filestring)               ; open workstation

;------------------------------------------------------
; Set Common Mapping Resources for Plots
;------------------------------------------------------

  res                     = True         ; plot modifications desired
;  res@gsnMaximize         = True         ; Maximize size of plot in frame
  res@tiMainOn = False
  res@tiXAxisOn = False
  res@cnInfoLabelOn = False 
  res@gsnSpreadColors     = True         ; Use full colormap, but start
  res@cnFillOn            = True         ; Turn on contour fill  
  res@cnFillMode          = "AreaFill"   ; Style of fill. You can also
;  res@cnFillMode          = "RasterFill" ; use "CellFill" and "RasterFill"
  res@cnLinesOn           = False        ; Turn off contour lines
  res@cnLineLabelsOn      = False        ; Turn off contour line labels
  res@lbLabelAutoStride   = True         ; Clean up labelbar labels.

  res@tmXTOn                = False            ; turn off top labels
  res@tmYROn                = False            ; turn off right labels

;MOVE TEXT STRINGS IN PLOTS AND SET FONT HEIGHTS
  Font = 0.05
  Position = 0.05

  res@tmYLLabelFontHeightF = Font
  res@tmXBLabelFontHeightF = Font
  res@gsnLeftStringFontHeightF = Font
  res@gsnLeftStringOrthogonalPosF = Position       ; Move text up or down from plot boundary
  res@gsnRightStringFontHeightF = Font
  res@gsnRightStringOrthogonalPosF = Position       ; Move text up or down from plot boundary
  res@gsnCenterStringFontHeightF = Font
  res@gsnCenterStringOrthogonalPosF = Position       ; Move text up or down f

; Give 45 degree turn on time tickmark labels
  res@tmXBLabelAngleF = 90.0      ; tilt the XB labels 45 degrees
  res@tmXBLabelDeltaF = 1.5     ; push the XB labels down from plot

; TURN INDIVIDUAL LABEL BARS ON/OFF
  res@tmXBLabelsOn         = True
  res@tmXBOn               = True

;PANEL PLOT SPECIFIC RESOURCES
  res@gsnFrame              = False            ; don't advance frame
  res@gsnDraw               = False           ; don't draw (Panel Plot)
  res@lbLabelBarOn          = False            ; turn off ndvdual cb's

;ASPECT RATIO (B/C OF WRF!!!)
;  res@mpShapeMode   = "FreeAspect"
;  res@vpWidthF      = 0.6
;  res@vpHeightF     = 0.6

; All data is regional
  res@gsnAddCyclic = False

; Reverse Y-Axis for Plot Attachments
;  res@trYReverse = True ; Reverses the y-axis.

;-----------------------------------
;  Set Lat/Lon for CA Domain
;------------------------------------
;  res@gsnMajorLonSpacing = 5.
;  res@gsnMajorLatSpacing = 5.

;  res@mpLimitMode       = "Corners"
;  res@mpLeftCornerLatF  = 35.0
;  res@mpLeftCornerLonF  = -125.0
;  res@mpRightCornerLatF = 45.0
;  res@mpRightCornerLonF = -115.0

;  res@mpLimitMode       = "Corners"
;  res@mpLeftCornerLatF  = 30.0
;  res@mpLeftCornerLonF  = -130.0
;  res@mpRightCornerLatF = 50.0
;  res@mpRightCornerLonF = -110.0

;  res@mpLimitMode       = "Corners"
;  res@mpLeftCornerLatF  = 20.0
;  res@mpLeftCornerLonF  = -130.0
;  res@mpRightCornerLatF = 50.0
;  res@mpRightCornerLonF = -100.0

;-------------------------------------------
; Plot Resources for Each Plot
;-------------------------------------------

  Minvalue_DIFF_SNOWFALL = -5
  Maxvalue_DIFF_SNOWFALL = 5
  Spacing_DIFF_SNOWFALL  = 1

  Minvalue_DIFF_SNOWC = -50
  Maxvalue_DIFF_SNOWC = 0
  Spacing_DIFF_SNOWC  = 5

  Minvalue_DIFF_SWE = -200
  Maxvalue_DIFF_SWE = 0
  Spacing_DIFF_SWE  = 10

  Minvalue_DIFF_2mST = 0
  Maxvalue_DIFF_2mST = 8
  Spacing_DIFF_2mST  = 1

;-------------------------------------------
; DIFFERENCE PLOTS
;-------------------------------------------
; 1980-2005
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon = res
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SNOWFALL,Maxvalue_DIFF_SNOWFALL,Spacing_DIFF_SNOWFALL)
; Set Plot Titles (UNIGRID28)
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@gsnLeftString    = ""
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@gsnCenterString    = ""
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@gsnRightString   = ""
;  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@tmYLMode = "Explicit"
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@tmYLLabels = ispan(1980,2005,1)
  res_SNOWFALL_HIST_DIFF_UNI28rough_lon@tmXBLabelsOn     = False

; 2025-2050
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon = res
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SNOWFALL,Maxvalue_DIFF_SNOWFALL,Spacing_DIFF_SNOWFALL)
; Set Plot Titles (UNIGRID28)
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnLeftString    = ""
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnCenterString    = ""
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnRightString   = ""
;  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLMode = "Explicit"
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabels = ispan(2025,2050,1)
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmXBLabelsOn     = False

; 2075-2100
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon = res
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SNOWFALL,Maxvalue_DIFF_SNOWFALL,Spacing_DIFF_SNOWFALL)
; Set Plot Titles (UNIGRID28)
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnLeftString    = ""
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnCenterString    = ""
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnRightString   = ""
;  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLMode = "Explicit"
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabels = ispan(2075,2100,1)

; 1980-2005
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat = res
; Set Plot Labels
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SNOWFALL,Maxvalue_DIFF_SNOWFALL,Spacing_DIFF_SNOWFALL)
; Set Plot Titles (UNIGRID28)
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@gsnLeftString    = ""
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@gsnCenterString    = ""
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@gsnRightString   = ""
;  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@tmXBMode = "Explicit"
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SNOWFALL_HIST_DIFF_UNI28rough_lat@tmXBLabels = ispan(1980,2005,1)

; 2025-2050
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat = res
; Set Plot Labels
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SNOWFALL,Maxvalue_DIFF_SNOWFALL,Spacing_DIFF_SNOWFALL)
; Set Plot Titles (UNIGRID28)
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnLeftString    = ""
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnCenterString    = ""
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnRightString   = ""
;  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBMode = "Explicit"
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabels = ispan(2025,2050,1)
  res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmYLLabelsOn     = False

; 2075-2100
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat = res
; Set Plot Labels
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SNOWFALL,Maxvalue_DIFF_SNOWFALL,Spacing_DIFF_SNOWFALL)
; Set Plot Titles (UNIGRID28)
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnLeftString    = ""
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnCenterString    = ""
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnRightString   = ""
;  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBMode = "Explicit"
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabels = ispan(2075,2100,1)
  res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmYLLabelsOn     = False

; 1980-2005
  res_SNOWC_HIST_DIFF_UNI28rough_lon = res
  res_SNOWC_HIST_DIFF_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SNOWC_HIST_DIFF_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWC_HIST_DIFF_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SNOWC,Maxvalue_DIFF_SNOWC,Spacing_DIFF_SNOWC)
; Set Plot Titles (UNIGRID28)
  res_SNOWC_HIST_DIFF_UNI28rough_lon@gsnLeftString    = ""
  res_SNOWC_HIST_DIFF_UNI28rough_lon@gsnCenterString    = ""
  res_SNOWC_HIST_DIFF_UNI28rough_lon@gsnRightString   = ""
;  res_SNOWC_HIST_DIFF_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SNOWC_HIST_DIFF_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SNOWC_HIST_DIFF_UNI28rough_lon@tmYLMode = "Explicit"
  res_SNOWC_HIST_DIFF_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SNOWC_HIST_DIFF_UNI28rough_lon@tmYLLabels = ispan(1980,2005,1)
  res_SNOWC_HIST_DIFF_UNI28rough_lon@tmXBLabelsOn     = False

; 2025-2050
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon = res
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SNOWC,Maxvalue_DIFF_SNOWC,Spacing_DIFF_SNOWC)
; Set Plot Titles (UNIGRID28)
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnLeftString    = ""
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnCenterString    = ""
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnRightString   = ""
;  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLMode = "Explicit"
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabels = ispan(2025,2050,1)
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmXBLabelsOn     = False

; 2075-2100
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon = res
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SNOWC,Maxvalue_DIFF_SNOWC,Spacing_DIFF_SNOWC)
; Set Plot Titles (UNIGRID28)
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnLeftString    = ""
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnCenterString    = ""
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnRightString   = ""
;  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLMode = "Explicit"
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabels = ispan(2075,2100,1)

; 1980-2005
  res_SNOWC_HIST_DIFF_UNI28rough_lat = res
; Set Plot Labels
  res_SNOWC_HIST_DIFF_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWC_HIST_DIFF_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SNOWC,Maxvalue_DIFF_SNOWC,Spacing_DIFF_SNOWC)
; Set Plot Titles (UNIGRID28)
  res_SNOWC_HIST_DIFF_UNI28rough_lat@gsnLeftString    = ""
  res_SNOWC_HIST_DIFF_UNI28rough_lat@gsnCenterString    = ""
  res_SNOWC_HIST_DIFF_UNI28rough_lat@gsnRightString   = ""
;  res_SNOWC_HIST_DIFF_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SNOWC_HIST_DIFF_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SNOWC_HIST_DIFF_UNI28rough_lat@tmXBMode = "Explicit"
  res_SNOWC_HIST_DIFF_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SNOWC_HIST_DIFF_UNI28rough_lat@tmXBLabels = ispan(1980,2005,1)

; 2025-2050
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat = res
; Set Plot Labels
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SNOWC,Maxvalue_DIFF_SNOWC,Spacing_DIFF_SNOWC)
; Set Plot Titles (UNIGRID28)
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnLeftString    = ""
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnCenterString    = ""
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnRightString   = ""
;  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBMode = "Explicit"
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabels = ispan(2025,2050,1)
  res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmYLLabelsOn     = False

; 2075-2100
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat = res
; Set Plot Labels
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SNOWC,Maxvalue_DIFF_SNOWC,Spacing_DIFF_SNOWC)
; Set Plot Titles (UNIGRID28)
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnLeftString    = ""
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnCenterString    = ""
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnRightString   = ""
;  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBMode = "Explicit"
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabels = ispan(2075,2100,1)
  res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmYLLabelsOn     = False

; 1980-2005
  res_SWE_HIST_DIFF_UNI28rough_lon = res
  res_SWE_HIST_DIFF_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SWE_HIST_DIFF_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SWE_HIST_DIFF_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SWE,Maxvalue_DIFF_SWE,Spacing_DIFF_SWE)
; Set Plot Titles (UNIGRID28)
  res_SWE_HIST_DIFF_UNI28rough_lon@gsnLeftString    = ""
  res_SWE_HIST_DIFF_UNI28rough_lon@gsnCenterString    = ""
  res_SWE_HIST_DIFF_UNI28rough_lon@gsnRightString   = ""
;  res_SWE_HIST_DIFF_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SWE_HIST_DIFF_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SWE_HIST_DIFF_UNI28rough_lon@tmYLMode = "Explicit"
  res_SWE_HIST_DIFF_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SWE_HIST_DIFF_UNI28rough_lon@tmYLLabels = ispan(1980,2005,1)
  res_SWE_HIST_DIFF_UNI28rough_lon@tmXBLabelsOn     = False

; 2025-2050
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon = res
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SWE,Maxvalue_DIFF_SWE,Spacing_DIFF_SWE)
; Set Plot Titles (UNIGRID28)
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnLeftString    = ""
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnCenterString    = ""
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnRightString   = ""
;  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLMode = "Explicit"
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabels = ispan(2025,2050,1)
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmXBLabelsOn     = False

; 2075-2100
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon = res
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_SWE,Maxvalue_DIFF_SWE,Spacing_DIFF_SWE)
; Set Plot Titles (UNIGRID28)
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnLeftString    = ""
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnCenterString    = ""
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnRightString   = ""
;  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabelStride = 5.0
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLMode = "Explicit"
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabels = ispan(2075,2100,1)

; 1980-2005
  res_SWE_HIST_DIFF_UNI28rough_lat = res
; Set Plot Labels
  res_SWE_HIST_DIFF_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SWE_HIST_DIFF_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SWE,Maxvalue_DIFF_SWE,Spacing_DIFF_SWE)
; Set Plot Titles (UNIGRID28)
  res_SWE_HIST_DIFF_UNI28rough_lat@gsnLeftString    = ""
  res_SWE_HIST_DIFF_UNI28rough_lat@gsnCenterString    = ""
  res_SWE_HIST_DIFF_UNI28rough_lat@gsnRightString   = ""
;  res_SWE_HIST_DIFF_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SWE_HIST_DIFF_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SWE_HIST_DIFF_UNI28rough_lat@tmXBMode = "Explicit"
  res_SWE_HIST_DIFF_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SWE_HIST_DIFF_UNI28rough_lat@tmXBLabels = ispan(1980,2005,1)

; 2025-2050
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat = res
; Set Plot Labels
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SWE,Maxvalue_DIFF_SWE,Spacing_DIFF_SWE)
; Set Plot Titles (UNIGRID28)
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnLeftString    = ""
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnCenterString    = ""
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnRightString   = ""
;  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBMode = "Explicit"
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabels = ispan(2025,2050,1)
  res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmYLLabelsOn     = False

; 2075-2100
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat = res
; Set Plot Labels
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_SWE,Maxvalue_DIFF_SWE,Spacing_DIFF_SWE)
; Set Plot Titles (UNIGRID28)
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnLeftString    = ""
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnCenterString    = ""
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnRightString   = ""
;  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabelStride = 5.0
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBMode = "Explicit"
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabels = ispan(2075,2100,1)
  res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmYLLabelsOn     = False

; 1980-2005
  res_2mST_HIST_DIFF_UNI28rough_lon = res
  res_2mST_HIST_DIFF_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_2mST_HIST_DIFF_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_2mST_HIST_DIFF_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_2mST,Maxvalue_DIFF_2mST,Spacing_DIFF_2mST)
; Set Plot Titles (UNIGRID28)
  res_2mST_HIST_DIFF_UNI28rough_lon@gsnLeftString    = ""
  res_2mST_HIST_DIFF_UNI28rough_lon@gsnCenterString    = ""
  res_2mST_HIST_DIFF_UNI28rough_lon@gsnRightString   = ""
;  res_2mST_HIST_DIFF_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_2mST_HIST_DIFF_UNI28rough_lon@tmYLLabelStride = 5.0
  res_2mST_HIST_DIFF_UNI28rough_lon@tmYLMode = "Explicit"
  res_2mST_HIST_DIFF_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_2mST_HIST_DIFF_UNI28rough_lon@tmYLLabels = ispan(1980,2005,1)
  res_2mST_HIST_DIFF_UNI28rough_lon@tmXBLabelsOn     = False

; 2025-2050
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon = res
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_2mST,Maxvalue_DIFF_2mST,Spacing_DIFF_2mST)
; Set Plot Titles (UNIGRID28)
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnLeftString    = ""
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnCenterString    = ""
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@gsnRightString   = ""
;  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabelStride = 5.0
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLMode = "Explicit"
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmYLLabels = ispan(2025,2050,1)
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon@tmXBLabelsOn     = False

; 2075-2100
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon = res
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@trYReverse = True ; Reverses the y-axis.
; Set Plot Labels
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@cnLevelSelectionMode = "ExplicitLevels"
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@cnLevels = ispan(Minvalue_DIFF_2mST,Maxvalue_DIFF_2mST,Spacing_DIFF_2mST)
; Set Plot Titles (UNIGRID28)
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnLeftString    = ""
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnCenterString    = ""
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@gsnRightString   = ""
;  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabelJust   = "BottomCenter"
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabelStride = 5.0
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLMode = "Explicit"
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLValues = ispan(0,24,1)
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon@tmYLLabels = ispan(2075,2100,1)

; 1980-2005
  res_2mST_HIST_DIFF_UNI28rough_lat = res
; Set Plot Labels
  res_2mST_HIST_DIFF_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_2mST_HIST_DIFF_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_2mST,Maxvalue_DIFF_2mST,Spacing_DIFF_2mST)
; Set Plot Titles (UNIGRID28)
  res_2mST_HIST_DIFF_UNI28rough_lat@gsnLeftString    = ""
  res_2mST_HIST_DIFF_UNI28rough_lat@gsnCenterString    = ""
  res_2mST_HIST_DIFF_UNI28rough_lat@gsnRightString   = ""
;  res_2mST_HIST_DIFF_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_2mST_HIST_DIFF_UNI28rough_lat@tmXBLabelStride = 5.0
  res_2mST_HIST_DIFF_UNI28rough_lat@tmXBMode = "Explicit"
  res_2mST_HIST_DIFF_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_2mST_HIST_DIFF_UNI28rough_lat@tmXBLabels = ispan(1980,2005,1)

; 2025-2050
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat = res
; Set Plot Labels
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_2mST,Maxvalue_DIFF_2mST,Spacing_DIFF_2mST)
; Set Plot Titles (UNIGRID28)
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnLeftString    = ""
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnCenterString    = ""
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@gsnRightString   = ""
;  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabelStride = 5.0
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBMode = "Explicit"
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmXBLabels = ispan(2025,2050,1)
  res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat@tmYLLabelsOn     = False

; 2075-2100
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat = res
; Set Plot Labels
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@cnLevelSelectionMode = "ExplicitLevels"
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@cnLevels = ispan(Minvalue_DIFF_2mST,Maxvalue_DIFF_2mST,Spacing_DIFF_2mST)
; Set Plot Titles (UNIGRID28)
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnLeftString    = ""
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnCenterString    = ""
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@gsnRightString   = ""
;  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabelJust   = "BottomCenter"
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabelStride = 5.0
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBMode = "Explicit"
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBValues = ispan(0,24,1)
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmXBLabels = ispan(2075,2100,1)
  res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat@tmYLLabelsOn     = False

;-------------------------------------------
; Average Across Latitude/Longitude to make 2D array and Reorder for Plot Function 
;-------------------------------------------
;DIFFERENCE PLOTS
; EXAMPLE ARRAY - RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough

; 1980-2005
  HIST_SNOWFALL_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(HISTvHIST_SNOWFALL_UNIGRID28_rough,(/1/))
  HIST_SNOWFALL_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(HISTvHIST_SNOWFALL_UNIGRID28_rough,(/2/))
  HIST_SNOWFALL_UNIGRID28_rough_DIFF_lat = HIST_SNOWFALL_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)  

; 2025-2050
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough,(/1/))
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(RCP85_MIDCENTvHIST_SNOWFALL_UNIGRID28_rough,(/2/))
  RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lat = RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 2075-2100
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough,(/1/))
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(RCP85_ENDCENTvHIST_SNOWFALL_UNIGRID28_rough,(/2/))
  RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lat = RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 1980-2005
  HIST_SNOWC_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(HISTvHIST_SNOWC_UNIGRID28_rough,(/1/))
  HIST_SNOWC_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(HISTvHIST_SNOWC_UNIGRID28_rough,(/2/))
  HIST_SNOWC_UNIGRID28_rough_DIFF_lat = HIST_SNOWC_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 2025-2050
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough,(/1/))
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(RCP85_MIDCENTvHIST_SNOWC_UNIGRID28_rough,(/2/))
  RCP85_MIDCENT_SNOWC_UNIGRID28_rough_DIFF_lat = RCP85_MIDCENT_SNOWC_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 2075-2100
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough,(/1/))
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(RCP85_ENDCENTvHIST_SNOWC_UNIGRID28_rough,(/2/))
  RCP85_ENDCENT_SNOWC_UNIGRID28_rough_DIFF_lat = RCP85_ENDCENT_SNOWC_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 1980-2005
  HIST_SWE_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(HISTvHIST_SWE_UNIGRID28_rough,(/1/))
  HIST_SWE_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(HISTvHIST_SWE_UNIGRID28_rough,(/2/))
  HIST_SWE_UNIGRID28_rough_DIFF_lat = HIST_SWE_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 2025-2050
  RCP85_MIDCENT_SWE_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough,(/1/))
  RCP85_MIDCENT_SWE_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(RCP85_MIDCENTvHIST_SWE_UNIGRID28_rough,(/2/))
  RCP85_MIDCENT_SWE_UNIGRID28_rough_DIFF_lat = RCP85_MIDCENT_SWE_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 2075-2100
  RCP85_ENDCENT_SWE_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough,(/1/))
  RCP85_ENDCENT_SWE_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(RCP85_ENDCENTvHIST_SWE_UNIGRID28_rough,(/2/))
  RCP85_ENDCENT_SWE_UNIGRID28_rough_DIFF_lat = RCP85_ENDCENT_SWE_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 1980-2005
  HIST_2mST_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(HISTvHIST_2mST_UNIGRID28_rough,(/1/))
  HIST_2mST_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(HISTvHIST_2mST_UNIGRID28_rough,(/2/))
  HIST_2mST_UNIGRID28_rough_DIFF_lat = HIST_2mST_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 2025-2050
  RCP85_MIDCENT_2mST_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough,(/1/))
  RCP85_MIDCENT_2mST_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(RCP85_MIDCENTvHIST_2mST_UNIGRID28_rough,(/2/))
  RCP85_MIDCENT_2mST_UNIGRID28_rough_DIFF_lat = RCP85_MIDCENT_2mST_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

; 2075-2100
  RCP85_ENDCENT_2mST_UNIGRID28_rough_DIFF_lon = dim_avg_n_Wrap(RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough,(/1/))
  RCP85_ENDCENT_2mST_UNIGRID28_rough_DIFF_lat_avg = dim_avg_n_Wrap(RCP85_ENDCENTvHIST_2mST_UNIGRID28_rough,(/2/))
  RCP85_ENDCENT_2mST_UNIGRID28_rough_DIFF_lat = RCP85_ENDCENT_2mST_UNIGRID28_rough_DIFF_lat_avg(lat|:,time|:)

;------------------------------------------
; PLOT
;------------------------------------------
; SNOWFALL
; DIFF

;  gsn_define_colormap(wks,"WhiteYellowOrangeRed")
  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)


  plot_SNOWFALL_HIST_DIFF_lon = new(3,graphic)
  plot_SNOWFALL_HIST_DIFF_lon(0) = gsn_csm_hov(wks,HIST_SNOWFALL_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SNOWFALL_HIST_DIFF_UNI28rough_lon)
  plot_SNOWFALL_HIST_DIFF_lon(1) = gsn_csm_hov(wks,RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lon)
  plot_SNOWFALL_HIST_DIFF_lon(2) = gsn_csm_hov(wks,RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lon)

  plot_SNOWFALL_HIST_DIFF_lat = new(3,graphic)
  plot_SNOWFALL_HIST_DIFF_lat(0) = gsn_csm_hov(wks,HIST_SNOWFALL_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SNOWFALL_HIST_DIFF_UNI28rough_lat)
  plot_SNOWFALL_HIST_DIFF_lat(1) = gsn_csm_lat_time(wks,RCP85_MIDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SNOWFALL_RCP85_DIFF_MIDCENT_UNI28rough_lat)
  plot_SNOWFALL_HIST_DIFF_lat(2) = gsn_csm_lat_time(wks,RCP85_ENDCENT_SNOWFALL_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SNOWFALL_RCP85_DIFF_ENDCENT_UNI28rough_lat)

; SNOWC
; DIFF

  gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  plot_SNOWC_HIST_DIFF_lon = new(3,graphic)
  plot_SNOWC_HIST_DIFF_lon(0) = gsn_csm_hov(wks,HIST_SNOWC_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SNOWC_HIST_DIFF_UNI28rough_lon)
  plot_SNOWC_HIST_DIFF_lon(1) = gsn_csm_hov(wks,RCP85_MIDCENT_SNOWC_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lon)
  plot_SNOWC_HIST_DIFF_lon(2)  = gsn_csm_hov(wks,RCP85_ENDCENT_SNOWC_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lon)

  plot_SNOWC_HIST_DIFF_lat = new(3,graphic)
  plot_SNOWC_HIST_DIFF_lat(0) = gsn_csm_hov(wks,HIST_SNOWC_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SNOWC_HIST_DIFF_UNI28rough_lat)
  plot_SNOWC_HIST_DIFF_lat(1) = gsn_csm_lat_time(wks,RCP85_MIDCENT_SNOWC_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SNOWC_RCP85_DIFF_MIDCENT_UNI28rough_lat)
  plot_SNOWC_HIST_DIFF_lat(2) = gsn_csm_lat_time(wks,RCP85_ENDCENT_SNOWC_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SNOWC_RCP85_DIFF_ENDCENT_UNI28rough_lat)

; ATTACH PLOTS
;  attachres1 = True
;  attachres1@gsnAttachPlotsXAxis = True ; attaches along x-axis 
;  attachres2 = True
;  attachid1  = 
;  gsn_attach_plots(plot_SNOWFALL_HIST,(/plot_SNOWFALL_RCP85_MIDCENT,plot_SNOWFALL_RCP85_ENDCENT/),attachres1,attachres2)

; SWE
; DIFF

  gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  plot_SWE_HIST_DIFF_lon = new(3,graphic)
  plot_SWE_HIST_DIFF_lon(0) = gsn_csm_hov(wks,HIST_SWE_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SWE_HIST_DIFF_UNI28rough_lon)
  plot_SWE_HIST_DIFF_lon(1) = gsn_csm_hov(wks,RCP85_MIDCENT_SWE_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lon)
  plot_SWE_HIST_DIFF_lon(2) = gsn_csm_hov(wks,RCP85_ENDCENT_SWE_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lon)

  plot_SWE_HIST_DIFF_lat = new(3,graphic)
  plot_SWE_HIST_DIFF_lat(0) = gsn_csm_hov(wks,HIST_SWE_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SWE_HIST_DIFF_UNI28rough_lat)
  plot_SWE_HIST_DIFF_lat(1) = gsn_csm_lat_time(wks,RCP85_MIDCENT_SWE_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SWE_RCP85_DIFF_MIDCENT_UNI28rough_lat)
  plot_SWE_HIST_DIFF_lat(2) = gsn_csm_lat_time(wks,RCP85_ENDCENT_SWE_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_SWE_RCP85_DIFF_ENDCENT_UNI28rough_lat)

; 2mST
; DIFF

  gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  plot_2mST_HIST_DIFF_lon = new(3,graphic)
  plot_2mST_HIST_DIFF_lon(0) = gsn_csm_hov(wks,HIST_2mST_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_2mST_HIST_DIFF_UNI28rough_lon)
  plot_2mST_HIST_DIFF_lon(1) = gsn_csm_hov(wks,RCP85_MIDCENT_2mST_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lon)
  plot_2mST_HIST_DIFF_lon(2) = gsn_csm_hov(wks,RCP85_ENDCENT_2mST_UNIGRID28_rough_DIFF_lon(time|:,{lon|-121.8:-118}),res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lon)

  plot_2mST_HIST_DIFF_lat = new(3,graphic)
  plot_2mST_HIST_DIFF_lat(0) = gsn_csm_hov(wks,HIST_2mST_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_2mST_HIST_DIFF_UNI28rough_lat)
  plot_2mST_HIST_DIFF_lat(1) = gsn_csm_lat_time(wks,RCP85_MIDCENT_2mST_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_2mST_RCP85_DIFF_MIDCENT_UNI28rough_lat)
  plot_2mST_HIST_DIFF_lat(2) = gsn_csm_lat_time(wks,RCP85_ENDCENT_2mST_UNIGRID28_rough_DIFF_lat({lat|34:40.5},time|:),res_2mST_RCP85_DIFF_ENDCENT_UNI28rough_lat)

; ATTACH PLOTS
;  attachres1 = True
;  attachres1@gsnAttachPlotsXAxis = True ; attaches along x-axis
;  attachres2 = True
;  attachid1  =
;  gsn_attach_plots(plot_SNOWFALL_HIST,(/plot_SNOWFALL_RCP85_MIDCENT,plot_SNOWFALL_RCP85_ENDCENT/),attachres1,attachres2)

;---------------------------
; Panel Plots Along X-Axis
;---------------------------
 resP_DIFF = True
 resP_DIFF@gsnMaximize = True

; CREATE WHITE SPACE BETWEEN PLOTS IN X/Y

 resP_DIFF@gsnPanelYWhiteSpacePercent = 0.00
 resP_DIFF@gsnPanelXWhiteSpacePercent = 0.00

; LABEL BAR LOCATION AND TEXT LOCATION
 resP_DIFF@gsnPanelLabelBar       = True                ; add common colorbar
 resP_DIFF@gsnPanelRowSpec = True
 resP_DIFF@lbOrientation =   "vertical"          ; vertical label bar
 resP_DIFF@lbTitlePosition      = "Right"                           ; title location
 resP_DIFF@lbTitleDirection     = "Across"                          ; letter angle
 resP_DIFF@lbTitleAngleF        = 90.                               ; title angle

; LABEL BAR FONT
  resP_DIFF@gsnPanelFigureStringsFontHeightF = 0.0075
  resP_DIFF@lbTitleFontHeightF = 0.0075       ; tweak the size of the labelbar title
  resP_DIFF@lbLabelFontHeightF               = 0.010                       ; label bar tickmark font

; LABEL BAR AND TITLE OFFSET
  resP_DIFF@lbTitleOffsetF = 0.15           ; positive values are up when the labelbar is under the plot
  resP_DIFF@pmLabelBarOrthogonalPosF = 0.15 ; move the labelbar down a bit so we have room for the overall x-axis

; PANEL PLOT SPECIFIC RESOURCES
  resP_DIFF@gsnFrame              = False            ; don't advance frame
  resP_DIFF@gsnDraw               = True           ; draw (Panel Plot)

; SNOWFALL 
 resP_SNOWFALL = resP_DIFF 
 resP_SNOWFALL@cnLevelSelectionMode = "ExplicitLevels"
 resP_SNOWFALL@cnLevels = ispan(Minvalue_DIFF_SNOWFALL,Maxvalue_DIFF_SNOWFALL,Spacing_DIFF_SNOWFALL)
 resP_SNOWFALL@gsnPanelFigureStrings = (/"a)","b)","c)"/) ; add strings to panel
 resP_SNOWFALL@lbTitleString = "Snowfall Anomaly (mm)" ; give the labelbar a title
 resP_SNOWFALL@gsnPanelTop = 0.950
 resP_SNOWFALL@gsnPanelBottom = 0.725

; SNOWC
 resP_SNOWC = resP_DIFF
 resP_SNOWC@cnLevelSelectionMode = "ExplicitLevels"
 resP_SNOWC@cnLevels = ispan(Minvalue_DIFF_SNOWC,Maxvalue_DIFF_SNOWC,Spacing_DIFF_SNOWC)
 resP_SNOWC@gsnPanelFigureStrings = (/"d)","e)","f)"/) ; add strings to panel
 resP_SNOWC@lbTitleString = "Snow Cover Anomaly (%)" ; give the labelbar a title
 resP_SNOWC@gsnPanelTop = 0.725
 resP_SNOWC@gsnPanelBottom = 0.500

; SWE
 resP_SWE = resP_DIFF
 resP_SWE@cnLevelSelectionMode = "ExplicitLevels"
 resP_SWE@cnLevels = ispan(Minvalue_DIFF_SWE,Maxvalue_DIFF_SWE,Spacing_DIFF_SWE)
 resP_SWE@gsnPanelFigureStrings = (/"g)","h)","i)"/) ; add strings to panel
 resP_SWE@lbTitleString = "Snow Water Equivalent Anomaly (mm)" ; give the labelbar a title
 resP_SWE@gsnPanelTop = 0.500
 resP_SWE@gsnPanelBottom = 0.275

; 2mST
 resP_2mST = resP_DIFF
 resP_2mST@cnLevelSelectionMode = "ExplicitLevels"
 resP_2mST@cnLevels = ispan(Minvalue_DIFF_2mST,Maxvalue_DIFF_2mST,Spacing_DIFF_2mST)
 resP_2mST@gsnPanelFigureStrings = (/"j)","k)","l)"/) ; add strings to panel
 resP_2mST@lbTitleString = "Surface Temperature Anomaly (K)" ; give the labelbar a title
 resP_2mST@gsnPanelTop = 0.275
 resP_2mST@gsnPanelBottom = 0.05

; LAT PLOTS
; SNOWFALL

;  gsn_define_colormap(wks,"WhiteYellowOrangeRed")
  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  gsn_panel(wks,(/plot_SNOWFALL_HIST_DIFF_lat/),(/3,1/),resP_SNOWFALL)

 gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  gsn_panel(wks,(/plot_SNOWC_HIST_DIFF_lat/),(/3,1/),resP_SNOWC)

 gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  gsn_panel(wks,(/plot_SWE_HIST_DIFF_lat/),(/3,1/),resP_SWE)

 gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
;  gsn_reverse_colormap(wks)

  gsn_panel(wks,(/plot_2mST_HIST_DIFF_lat/),(/3,1/),resP_2mST)

 frame(wks)

; LON PLOTS

;  gsn_define_colormap(wks,"WhiteYellowOrangeRed")
  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  gsn_panel(wks,(/plot_SNOWFALL_HIST_DIFF_lon/),(/1,1,1/),resP_SNOWFALL)

 gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  gsn_panel(wks,(/plot_SNOWC_HIST_DIFF_lon/),(/1,1,1/),resP_SNOWC)

 gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
  gsn_reverse_colormap(wks)

  gsn_panel(wks,(/plot_SWE_HIST_DIFF_lon/),(/1,1,1/),resP_SWE)

 gsn_define_colormap(wks,"WhiteYellowOrangeRed")
;  gsn_define_colormap(wks,"GMT_polar")   ; choose colormap
;  gsn_define_colormap(wks,"topo_15lev")   ; choose colormap
;  gsn_reverse_colormap(wks)

  gsn_panel(wks,(/plot_2mST_HIST_DIFF_lon/),(/1,1,1/),resP_2mST)

  print("Done Creating Climate Plot")

end