;************************************************* ; ; Simple example of how to draw selected geometry from a shapefile, ; based upon properties of an associated non-spatial variable. ; ; This example color-fills the states based upon "percent unemployment", ; which is calculated from several of the non-spatial variables in the ; file. ; ; "states.shp" is from the National Atlas (http://www.nationalatlas.gov/) ; ; You must also have the files "states.dbf" and "states.shx" for this ; example to run. ;************************************************* 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" load "/home/arhoades/plots/latlonplots/shapefile/shapefile_mask_data.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 setvalues NhlGetWorkspaceObjectId() "wsMaximumSize" : 300000000 end setvalues ;------------------------------------------------------ ; Type of Plot ;------------------------------------------------------ ; type = "x11" type = "pdf" ; type = "ps" ; type = "eps" ; type = "ncgm" ;--------------------------------------------------- ; Open Workstation and Name File ;--------------------------------------------------- filename = "SN_TRENDLINE_ALLVAR_HISTvRCP85" wks = gsn_open_wks (type,filename) ; open workstation ;------------------------------------------------------ ; ADD CLIMATE AVERAGE MODEL DATA ;------------------------------------------------------ ; VR_CESM 0.25 - 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_VR_CESM_28_rough = fspan(20,50,240) lon_HIST_VR_CESM_28_rough = fspan(-130,-110,160) HIST_SNOWFALL_VR_CESM_28_rough = f_HIST_UNI28rough->DJF_SNOWFALL_HIST ; 1 kg/m2 = 1 mm HIST_SNOWFALL_VR_CESM_28_rough@_FillValue = 1e+36 ; HIST_SNOWFALL_VR_CESM_28_rough = f_HIST_UNI28rough[:]->H2OSNO ; 1 kg/m2 = 1 mm ; HIST_SNOWFALL_VR_CESM_28_rough = HIST_SNOWFALL_VR_CESM_28_rough-0.1 HIST_SNOWFALL_VR_CESM_28_rough!0 = "time" HIST_SNOWFALL_VR_CESM_28_rough!1 = "lat" HIST_SNOWFALL_VR_CESM_28_rough!2 = "lon" HIST_SNOWFALL_VR_CESM_28_rough&lat = lat_HIST_VR_CESM_28_rough HIST_SNOWFALL_VR_CESM_28_rough&lon = lon_HIST_VR_CESM_28_rough HIST_SNOWFALL_VR_CESM_28_rough&lat@units = "degrees_north" HIST_SNOWFALL_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_HIST_SNOWFALL_VR_CESM_28_rough = dimsizes(HIST_SNOWFALL_VR_CESM_28_rough) nt_HIST_SNOWFALL_VR_CESM_28_rough = dsizes_HIST_SNOWFALL_VR_CESM_28_rough(0) ny_HIST_SNOWFALL_VR_CESM_28_rough = dsizes_HIST_SNOWFALL_VR_CESM_28_rough(1) nx_HIST_SNOWFALL_VR_CESM_28_rough = dsizes_HIST_SNOWFALL_VR_CESM_28_rough(2) printVarSummary(HIST_SNOWFALL_VR_CESM_28_rough) print ("Minimum="+min(HIST_SNOWFALL_VR_CESM_28_rough) +" Maximum="+max(HIST_SNOWFALL_VR_CESM_28_rough)) lat_HIST_VR_CESM_28_rough = fspan(20,50,240) lon_HIST_VR_CESM_28_rough = fspan(-130,-110,160) HIST_SNOWC_VR_CESM_28_rough = f_HIST_UNI28rough->DJF_SNOWC_HIST ; 1 kg/m2 = 1 mm HIST_SNOWC_VR_CESM_28_rough@_FillValue = 1e+36 ; HIST_SNOWC_VR_CESM_28_rough = f_HIST_UNI28rough[:]->FSNO ; 1 kg/m2 = 1 mm ; HIST_SNOWC_VR_CESM_28_rough = HIST_SNOWC_VR_CESM_28_rough-0.001 ; HIST_SNOWC_VR_CESM_28_rough = HIST_SNOWC_VR_CESM_28_rough*100 HIST_SNOWC_VR_CESM_28_rough!0 = "time" HIST_SNOWC_VR_CESM_28_rough!1 = "lat" HIST_SNOWC_VR_CESM_28_rough!2 = "lon" HIST_SNOWC_VR_CESM_28_rough&lat = lat_HIST_VR_CESM_28_rough HIST_SNOWC_VR_CESM_28_rough&lon = lon_HIST_VR_CESM_28_rough HIST_SNOWC_VR_CESM_28_rough&lat@units = "degrees_north" HIST_SNOWC_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_HIST_SNOWC_VR_CESM_28_rough = dimsizes(HIST_SNOWC_VR_CESM_28_rough) nt_HIST_SNOWC_VR_CESM_28_rough = dsizes_HIST_SNOWC_VR_CESM_28_rough(0) ny_HIST_SNOWC_VR_CESM_28_rough = dsizes_HIST_SNOWC_VR_CESM_28_rough(1) nx_HIST_SNOWC_VR_CESM_28_rough = dsizes_HIST_SNOWC_VR_CESM_28_rough(2) printVarSummary(HIST_SNOWC_VR_CESM_28_rough) print ("Minimum="+min(HIST_SNOWC_VR_CESM_28_rough) +" Maximum="+max(HIST_SNOWC_VR_CESM_28_rough)) lat_HIST_VR_CESM_28_rough = fspan(20,50,240) lon_HIST_VR_CESM_28_rough = fspan(-130,-110,160) HIST_SWE_VR_CESM_28_rough = f_HIST_UNI28rough->DJF_SWE_HIST ; 1 kg/m2 = 1 mm HIST_SWE_VR_CESM_28_rough@_FillValue = 1e+36 ; HIST_SWE_VR_CESM_28_rough = f_HIST_UNI28rough[:]->SWE ; 1 kg/m2 = 1 mm HIST_SWE_VR_CESM_28_rough!0 = "time" HIST_SWE_VR_CESM_28_rough!1 = "lat" HIST_SWE_VR_CESM_28_rough!2 = "lon" HIST_SWE_VR_CESM_28_rough&lat = lat_HIST_VR_CESM_28_rough HIST_SWE_VR_CESM_28_rough&lon = lon_HIST_VR_CESM_28_rough HIST_SWE_VR_CESM_28_rough&lat@units = "degrees_north" HIST_SWE_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_HIST_SWE_VR_CESM_28_rough = dimsizes(HIST_SWE_VR_CESM_28_rough) nt_HIST_SWE_VR_CESM_28_rough = dsizes_HIST_SWE_VR_CESM_28_rough(0) ny_HIST_SWE_VR_CESM_28_rough = dsizes_HIST_SWE_VR_CESM_28_rough(1) nx_HIST_SWE_VR_CESM_28_rough = dsizes_HIST_SWE_VR_CESM_28_rough(2) printVarSummary(HIST_SWE_VR_CESM_28_rough) print ("Minimum="+min(HIST_SWE_VR_CESM_28_rough) +" Maximum="+max(HIST_SWE_VR_CESM_28_rough)) lat_HIST_VR_CESM_28_rough = fspan(20,50,240) lon_HIST_VR_CESM_28_rough = fspan(-130,-110,160) HIST_TSA_VR_CESM_28_rough = f_HIST_UNI28rough->DJF_2mST_HIST ; 1 kg/m2 = 1 mm HIST_TSA_VR_CESM_28_rough@_FillValue = 1e+36 ; HIST_TSA_VR_CESM_28_rough = f_HIST_UNI28rough[:]->TSA ; 1 kg/m2 = 1 mm HIST_TSA_VR_CESM_28_rough!0 = "time" HIST_TSA_VR_CESM_28_rough!1 = "lat" HIST_TSA_VR_CESM_28_rough!2 = "lon" HIST_TSA_VR_CESM_28_rough&lat = lat_HIST_VR_CESM_28_rough HIST_TSA_VR_CESM_28_rough&lon = lon_HIST_VR_CESM_28_rough HIST_TSA_VR_CESM_28_rough&lat@units = "degrees_north" HIST_TSA_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_HIST_TSA_VR_CESM_28_rough = dimsizes(HIST_TSA_VR_CESM_28_rough) nt_HIST_TSA_VR_CESM_28_rough = dsizes_HIST_TSA_VR_CESM_28_rough(0) ny_HIST_TSA_VR_CESM_28_rough = dsizes_HIST_TSA_VR_CESM_28_rough(1) nx_HIST_TSA_VR_CESM_28_rough = dsizes_HIST_TSA_VR_CESM_28_rough(2) printVarSummary(HIST_TSA_VR_CESM_28_rough) print ("Minimum="+min(HIST_TSA_VR_CESM_28_rough) +" Maximum="+max(HIST_TSA_VR_CESM_28_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_VR_CESM_28_rough = fspan(20,50,240) lon_RCP85_MIDCENT_VR_CESM_28_rough = fspan(-130,-110,160) RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough = f_RCP85_MIDCENT_UNI28rough->DJF_SNOWFALL_RCP85_MIDCENT ; 1 kg/m2 = 1 mm RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough@_FillValue = 1e+36 ; RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough = f_RCP85_MIDCENT_UNI28rough[:]->H2OSNO ; 1 kg/m2 = 1 mm ; RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough = RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough-0.1 RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough!0 = "time" RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough!1 = "lat" RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough!2 = "lon" RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough&lat = lat_RCP85_MIDCENT_VR_CESM_28_rough RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough&lon = lon_RCP85_MIDCENT_VR_CESM_28_rough RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough&lat@units = "degrees_north" RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough = dimsizes(RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough) nt_RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough(0) ny_RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough(1) nx_RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough(2) printVarSummary(RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough) print ("Minimum="+min(RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough) +" Maximum="+max(RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough)) lat_RCP85_MIDCENT_VR_CESM_28_rough = fspan(20,50,240) lon_RCP85_MIDCENT_VR_CESM_28_rough = fspan(-130,-110,160) RCP85_MIDCENT_SNOWC_VR_CESM_28_rough = f_RCP85_MIDCENT_UNI28rough->DJF_SNOWC_RCP85_MIDCENT ; 1 kg/m2 = 1 mm RCP85_MIDCENT_SNOWC_VR_CESM_28_rough@_FillValue = 1e+36 ; RCP85_MIDCENT_SNOWC_VR_CESM_28_rough = f_RCP85_MIDCENT_UNI28rough[:]->FSNO ; 1 kg/m2 = 1 mm ; RCP85_MIDCENT_SNOWC_VR_CESM_28_rough = RCP85_MIDCENT_SNOWC_VR_CESM_28_rough-0.001 ; RCP85_MIDCENT_SNOWC_VR_CESM_28_rough = RCP85_MIDCENT_SNOWC_VR_CESM_28_rough*100 RCP85_MIDCENT_SNOWC_VR_CESM_28_rough!0 = "time" RCP85_MIDCENT_SNOWC_VR_CESM_28_rough!1 = "lat" RCP85_MIDCENT_SNOWC_VR_CESM_28_rough!2 = "lon" RCP85_MIDCENT_SNOWC_VR_CESM_28_rough&lat = lat_RCP85_MIDCENT_VR_CESM_28_rough RCP85_MIDCENT_SNOWC_VR_CESM_28_rough&lon = lon_RCP85_MIDCENT_VR_CESM_28_rough RCP85_MIDCENT_SNOWC_VR_CESM_28_rough&lat@units = "degrees_north" RCP85_MIDCENT_SNOWC_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_RCP85_MIDCENT_SNOWC_VR_CESM_28_rough = dimsizes(RCP85_MIDCENT_SNOWC_VR_CESM_28_rough) nt_RCP85_MIDCENT_SNOWC_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SNOWC_VR_CESM_28_rough(0) ny_RCP85_MIDCENT_SNOWC_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SNOWC_VR_CESM_28_rough(1) nx_RCP85_MIDCENT_SNOWC_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SNOWC_VR_CESM_28_rough(2) printVarSummary(RCP85_MIDCENT_SNOWC_VR_CESM_28_rough) print ("Minimum="+min(RCP85_MIDCENT_SNOWC_VR_CESM_28_rough) +" Maximum="+max(RCP85_MIDCENT_SNOWC_VR_CESM_28_rough)) lat_RCP85_MIDCENT_VR_CESM_28_rough = fspan(20,50,240) lon_RCP85_MIDCENT_VR_CESM_28_rough = fspan(-130,-110,160) RCP85_MIDCENT_SWE_VR_CESM_28_rough = f_RCP85_MIDCENT_UNI28rough->DJF_SWE_RCP85_MIDCENT ; 1 kg/m2 = 1 mm RCP85_MIDCENT_SWE_VR_CESM_28_rough@_FillValue = 1e+36 ; RCP85_MIDCENT_SWE_VR_CESM_28_rough = f_RCP85_MIDCENT_UNI28rough[:]->SWE ; 1 kg/m2 = 1 mm RCP85_MIDCENT_SWE_VR_CESM_28_rough!0 = "time" RCP85_MIDCENT_SWE_VR_CESM_28_rough!1 = "lat" RCP85_MIDCENT_SWE_VR_CESM_28_rough!2 = "lon" RCP85_MIDCENT_SWE_VR_CESM_28_rough&lat = lat_RCP85_MIDCENT_VR_CESM_28_rough RCP85_MIDCENT_SWE_VR_CESM_28_rough&lon = lon_RCP85_MIDCENT_VR_CESM_28_rough RCP85_MIDCENT_SWE_VR_CESM_28_rough&lat@units = "degrees_north" RCP85_MIDCENT_SWE_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_RCP85_MIDCENT_SWE_VR_CESM_28_rough = dimsizes(RCP85_MIDCENT_SWE_VR_CESM_28_rough) nt_RCP85_MIDCENT_SWE_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SWE_VR_CESM_28_rough(0) ny_RCP85_MIDCENT_SWE_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SWE_VR_CESM_28_rough(1) nx_RCP85_MIDCENT_SWE_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_SWE_VR_CESM_28_rough(2) printVarSummary(RCP85_MIDCENT_SWE_VR_CESM_28_rough) print ("Minimum="+min(RCP85_MIDCENT_SWE_VR_CESM_28_rough) +" Maximum="+max(RCP85_MIDCENT_SWE_VR_CESM_28_rough)) lat_RCP85_MIDCENT_VR_CESM_28_rough = fspan(20,50,240) lon_RCP85_MIDCENT_VR_CESM_28_rough = fspan(-130,-110,160) RCP85_MIDCENT_TSA_VR_CESM_28_rough = f_RCP85_MIDCENT_UNI28rough->DJF_2mST_RCP85_MIDCENT ; 1 kg/m2 = 1 mm RCP85_MIDCENT_TSA_VR_CESM_28_rough@_FillValue = 1e+36 ; RCP85_MIDCENT_TSA_VR_CESM_28_rough = f_RCP85_MIDCENT_UNI28rough[:]->TSA ; 1 kg/m2 = 1 mm RCP85_MIDCENT_TSA_VR_CESM_28_rough!0 = "time" RCP85_MIDCENT_TSA_VR_CESM_28_rough!1 = "lat" RCP85_MIDCENT_TSA_VR_CESM_28_rough!2 = "lon" RCP85_MIDCENT_TSA_VR_CESM_28_rough&lat = lat_RCP85_MIDCENT_VR_CESM_28_rough RCP85_MIDCENT_TSA_VR_CESM_28_rough&lon = lon_RCP85_MIDCENT_VR_CESM_28_rough RCP85_MIDCENT_TSA_VR_CESM_28_rough&lat@units = "degrees_north" RCP85_MIDCENT_TSA_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_RCP85_MIDCENT_TSA_VR_CESM_28_rough = dimsizes(RCP85_MIDCENT_TSA_VR_CESM_28_rough) nt_RCP85_MIDCENT_TSA_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_TSA_VR_CESM_28_rough(0) ny_RCP85_MIDCENT_TSA_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_TSA_VR_CESM_28_rough(1) nx_RCP85_MIDCENT_TSA_VR_CESM_28_rough = dsizes_RCP85_MIDCENT_TSA_VR_CESM_28_rough(2) printVarSummary(RCP85_MIDCENT_TSA_VR_CESM_28_rough) print ("Minimum="+min(RCP85_MIDCENT_TSA_VR_CESM_28_rough) +" Maximum="+max(RCP85_MIDCENT_TSA_VR_CESM_28_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_28_RCP85_ENDCENT/pre-process/lnd/DJFavg/REGRID/*.avgDJF.box.nc") ; f_RCP85_ENDCENT_UNI28rough = addfiles(fils_RCP85_ENDCENT_UNI28rough,"r") lat_RCP85_ENDCENT_VR_CESM_28_rough = fspan(20,50,240) lon_RCP85_ENDCENT_VR_CESM_28_rough = fspan(-130,-110,160) RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough = f_RCP85_ENDCENT_UNI28rough->DJF_SNOWFALL_RCP85_ENDCENT ; 1 kg/m2 = 1 mm RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough@_FillValue = 1e+36 ; RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough = f_RCP85_ENDCENT_UNI28rough[:]->H2OSNO ; 1 kg/m2 = 1 mm ; RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough = RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough-0.1 RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough!0 = "time" RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough!1 = "lat" RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough!2 = "lon" RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough&lat = lat_RCP85_ENDCENT_VR_CESM_28_rough RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough&lon = lon_RCP85_ENDCENT_VR_CESM_28_rough RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough&lat@units = "degrees_north" RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough = dimsizes(RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough) nt_RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough(0) ny_RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough(1) nx_RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough(2) printVarSummary(RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough) print ("Minimum="+min(RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough) +" Maximum="+max(RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough)) lat_RCP85_ENDCENT_VR_CESM_28_rough = fspan(20,50,240) lon_RCP85_ENDCENT_VR_CESM_28_rough = fspan(-130,-110,160) RCP85_ENDCENT_SNOWC_VR_CESM_28_rough = f_RCP85_ENDCENT_UNI28rough->DJF_SNOWC_RCP85_ENDCENT ; 1 kg/m2 = 1 mm RCP85_ENDCENT_SNOWC_VR_CESM_28_rough@_FillValue = 1e+36 ; RCP85_ENDCENT_SNOWC_VR_CESM_28_rough = f_RCP85_ENDCENT_UNI28rough[:]->FSNO ; 1 kg/m2 = 1 mm ; RCP85_ENDCENT_SNOWC_VR_CESM_28_rough = RCP85_ENDCENT_SNOWC_VR_CESM_28_rough-0.001 ; RCP85_ENDCENT_SNOWC_VR_CESM_28_rough = RCP85_ENDCENT_SNOWC_VR_CESM_28_rough*100 RCP85_ENDCENT_SNOWC_VR_CESM_28_rough!0 = "time" RCP85_ENDCENT_SNOWC_VR_CESM_28_rough!1 = "lat" RCP85_ENDCENT_SNOWC_VR_CESM_28_rough!2 = "lon" RCP85_ENDCENT_SNOWC_VR_CESM_28_rough&lat = lat_RCP85_ENDCENT_VR_CESM_28_rough RCP85_ENDCENT_SNOWC_VR_CESM_28_rough&lon = lon_RCP85_ENDCENT_VR_CESM_28_rough RCP85_ENDCENT_SNOWC_VR_CESM_28_rough&lat@units = "degrees_north" RCP85_ENDCENT_SNOWC_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_RCP85_ENDCENT_SNOWC_VR_CESM_28_rough = dimsizes(RCP85_ENDCENT_SNOWC_VR_CESM_28_rough) nt_RCP85_ENDCENT_SNOWC_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SNOWC_VR_CESM_28_rough(0) ny_RCP85_ENDCENT_SNOWC_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SNOWC_VR_CESM_28_rough(1) nx_RCP85_ENDCENT_SNOWC_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SNOWC_VR_CESM_28_rough(2) printVarSummary(RCP85_ENDCENT_SNOWC_VR_CESM_28_rough) print ("Minimum="+min(RCP85_ENDCENT_SNOWC_VR_CESM_28_rough) +" Maximum="+max(RCP85_ENDCENT_SNOWC_VR_CESM_28_rough)) lat_RCP85_ENDCENT_VR_CESM_28_rough = fspan(20,50,240) lon_RCP85_ENDCENT_VR_CESM_28_rough = fspan(-130,-110,160) RCP85_ENDCENT_SWE_VR_CESM_28_rough = f_RCP85_ENDCENT_UNI28rough->DJF_SWE_RCP85_ENDCENT ; 1 kg/m2 = 1 mm RCP85_ENDCENT_SWE_VR_CESM_28_rough@_FillValue = 1e+36 ; RCP85_ENDCENT_SWE_VR_CESM_28_rough = f_RCP85_ENDCENT_UNI28rough[:]->SWE ; 1 kg/m2 = 1 mm RCP85_ENDCENT_SWE_VR_CESM_28_rough!0 = "time" RCP85_ENDCENT_SWE_VR_CESM_28_rough!1 = "lat" RCP85_ENDCENT_SWE_VR_CESM_28_rough!2 = "lon" RCP85_ENDCENT_SWE_VR_CESM_28_rough&lat = lat_RCP85_ENDCENT_VR_CESM_28_rough RCP85_ENDCENT_SWE_VR_CESM_28_rough&lon = lon_RCP85_ENDCENT_VR_CESM_28_rough RCP85_ENDCENT_SWE_VR_CESM_28_rough&lat@units = "degrees_north" RCP85_ENDCENT_SWE_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_RCP85_ENDCENT_SWE_VR_CESM_28_rough = dimsizes(RCP85_ENDCENT_SWE_VR_CESM_28_rough) nt_RCP85_ENDCENT_SWE_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SWE_VR_CESM_28_rough(0) ny_RCP85_ENDCENT_SWE_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SWE_VR_CESM_28_rough(1) nx_RCP85_ENDCENT_SWE_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_SWE_VR_CESM_28_rough(2) printVarSummary(RCP85_ENDCENT_SWE_VR_CESM_28_rough) print ("Minimum="+min(RCP85_ENDCENT_SWE_VR_CESM_28_rough) +" Maximum="+max(RCP85_ENDCENT_SWE_VR_CESM_28_rough)) lat_RCP85_ENDCENT_VR_CESM_28_rough = fspan(20,50,240) lon_RCP85_ENDCENT_VR_CESM_28_rough = fspan(-130,-110,160) RCP85_ENDCENT_TSA_VR_CESM_28_rough = f_RCP85_ENDCENT_UNI28rough->DJF_2mST_RCP85_ENDCENT ; 1 kg/m2 = 1 mm RCP85_ENDCENT_TSA_VR_CESM_28_rough@_FillValue = 1e+36 ; RCP85_ENDCENT_TSA_VR_CESM_28_rough = f_RCP85_ENDCENT_UNI28rough[:]->TSA ; 1 kg/m2 = 1 mm RCP85_ENDCENT_TSA_VR_CESM_28_rough!0 = "time" RCP85_ENDCENT_TSA_VR_CESM_28_rough!1 = "lat" RCP85_ENDCENT_TSA_VR_CESM_28_rough!2 = "lon" RCP85_ENDCENT_TSA_VR_CESM_28_rough&lat = lat_RCP85_ENDCENT_VR_CESM_28_rough RCP85_ENDCENT_TSA_VR_CESM_28_rough&lon = lon_RCP85_ENDCENT_VR_CESM_28_rough RCP85_ENDCENT_TSA_VR_CESM_28_rough&lat@units = "degrees_north" RCP85_ENDCENT_TSA_VR_CESM_28_rough&lon@units = "degrees_east" dsizes_RCP85_ENDCENT_TSA_VR_CESM_28_rough = dimsizes(RCP85_ENDCENT_TSA_VR_CESM_28_rough) nt_RCP85_ENDCENT_TSA_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_TSA_VR_CESM_28_rough(0) ny_RCP85_ENDCENT_TSA_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_TSA_VR_CESM_28_rough(1) nx_RCP85_ENDCENT_TSA_VR_CESM_28_rough = dsizes_RCP85_ENDCENT_TSA_VR_CESM_28_rough(2) printVarSummary(RCP85_ENDCENT_TSA_VR_CESM_28_rough) print ("Minimum="+min(RCP85_ENDCENT_TSA_VR_CESM_28_rough) +" Maximum="+max(RCP85_ENDCENT_TSA_VR_CESM_28_rough)) ;------------------------------------------------------ ; CHANGE ARRAY TO 1D FOR STAT DISPERSION ;------------------------------------------------------ oned_SNOWFALL_VR_CESM_28_HIST = ndtooned(HIST_SNOWFALL_VR_CESM_28_rough) oned_SNOWFALL_VR_CESM_28_RCP85_MIDCENT = ndtooned(RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough) oned_SNOWFALL_VR_CESM_28_RCP85_ENDCENT = ndtooned(RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough) oned_SNOWC_VR_CESM_28_HIST = ndtooned(HIST_SNOWC_VR_CESM_28_rough) oned_SNOWC_VR_CESM_28_RCP85_MIDCENT = ndtooned(RCP85_MIDCENT_SNOWC_VR_CESM_28_rough) oned_SNOWC_VR_CESM_28_RCP85_ENDCENT = ndtooned(RCP85_ENDCENT_SNOWC_VR_CESM_28_rough) oned_SWE_VR_CESM_28_HIST = ndtooned(HIST_SWE_VR_CESM_28_rough) oned_SWE_VR_CESM_28_RCP85_MIDCENT = ndtooned(RCP85_MIDCENT_SWE_VR_CESM_28_rough) oned_SWE_VR_CESM_28_RCP85_ENDCENT = ndtooned(RCP85_ENDCENT_SWE_VR_CESM_28_rough) oned_TSA_VR_CESM_28_HIST = ndtooned(HIST_SNOWC_VR_CESM_28_rough) oned_TSA_VR_CESM_28_RCP85_MIDCENT = ndtooned(RCP85_MIDCENT_SNOWC_VR_CESM_28_rough) oned_TSA_VR_CESM_28_RCP85_ENDCENT = ndtooned(RCP85_ENDCENT_SNOWC_VR_CESM_28_rough) ;------------------------------------------------------ ; STAT DISPERSION ;------------------------------------------------------ opt = True opt@PrintStat = True print("VR_CESM 0.25 - HISTORICAL - AVG DJF Sierra SNOWFALL (mm)") Sierra_STAT_VR_CESM_28_HIST = stat_dispersion(oned_SNOWFALL_VR_CESM_28_HIST,opt) print("VR_CESM 0.25 - MIDCENT - AVG DJF Sierra SNOWFALL (mm)") Sierra_STAT_VR_CESM_28_RCP85_MIDCENT = stat_dispersion(oned_SNOWFALL_VR_CESM_28_RCP85_MIDCENT,opt) print("VR_CESM 0.25 - ENDCENT - AVG DJF Sierra SNOWFALL (mm)") Sierra_STAT_VR_CESM_28_RCP85_ENDCENT = stat_dispersion(oned_SNOWFALL_VR_CESM_28_RCP85_ENDCENT,opt) print("VR_CESM 0.25 - HISTORICAL - AVG DJF Sierra SNOWC (%)") Sierra_STAT_VR_CESM_28_HIST = stat_dispersion(oned_SNOWC_VR_CESM_28_HIST,opt) print("VR_CESM 0.25 - MIDCENT - AVG DJF Sierra SNOWC (%)") Sierra_STAT_VR_CESM_28_RCP85_MIDCENT = stat_dispersion(oned_SNOWC_VR_CESM_28_RCP85_MIDCENT,opt) print("VR_CESM 0.25 - ENDCENT - AVG DJF Sierra SNOWC (%)") Sierra_STAT_VR_CESM_28_RCP85_ENDCENT = stat_dispersion(oned_SNOWC_VR_CESM_28_RCP85_ENDCENT,opt) print("VR_CESM 0.25 - HISTORICAL - AVG DJF Sierra SWE (mm)") Sierra_STAT_VR_CESM_28_HIST = stat_dispersion(oned_SWE_VR_CESM_28_HIST,opt) print("VR_CESM 0.25 - MIDCENT - AVG DJF Sierra SWE (mm)") Sierra_STAT_VR_CESM_28_RCP85_MIDCENT = stat_dispersion(oned_SWE_VR_CESM_28_RCP85_MIDCENT,opt) print("VR_CESM 0.25 - ENDCENT - AVG DJF Sierra SWE (mm)") Sierra_STAT_VR_CESM_28_RCP85_ENDCENT = stat_dispersion(oned_SWE_VR_CESM_28_RCP85_ENDCENT,opt) print("VR_CESM 0.25 - HISTORICAL - AVG DJF Sierra TSA (K)") Sierra_STAT_VR_CESM_28_HIST = stat_dispersion(oned_TSA_VR_CESM_28_HIST,opt) print("VR_CESM 0.25 - MIDCENT - AVG DJF Sierra TSA (K)") Sierra_STAT_VR_CESM_28_RCP85_MIDCENT = stat_dispersion(oned_TSA_VR_CESM_28_RCP85_MIDCENT,opt) print("VR_CESM 0.25 - ENDCENT - AVG DJF Sierra TSA (K)") Sierra_STAT_VR_CESM_28_RCP85_ENDCENT = stat_dispersion(oned_TSA_VR_CESM_28_RCP85_ENDCENT,opt) ;------------------------------------------------------ ; CREATE SEASONAL CLIMATE AVERAGE FOR Sierra Nevada ;------------------------------------------------------ SNOWFALL_MASK_LATLON_VR_CESM_28_HIST = dim_avg_n_Wrap(HIST_SNOWFALL_VR_CESM_28_rough, (/1,2/)) printVarSummary(SNOWFALL_MASK_LATLON_VR_CESM_28_HIST) SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT = dim_avg_n_Wrap(RCP85_MIDCENT_SNOWFALL_VR_CESM_28_rough, (/1,2/)) printVarSummary(SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT = dim_avg_n_Wrap(RCP85_ENDCENT_SNOWFALL_VR_CESM_28_rough, (/1,2/)) printVarSummary(SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) SNOWC_MASK_LATLON_VR_CESM_28_HIST = dim_avg_n_Wrap(HIST_SNOWC_VR_CESM_28_rough, (/1,2/)) printVarSummary(SNOWC_MASK_LATLON_VR_CESM_28_HIST) SNOWC_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT = dim_avg_n_Wrap(RCP85_MIDCENT_SNOWC_VR_CESM_28_rough, (/1,2/)) printVarSummary(SNOWC_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) SNOWC_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT = dim_avg_n_Wrap(RCP85_ENDCENT_SNOWC_VR_CESM_28_rough, (/1,2/)) printVarSummary(SNOWC_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) SWE_MASK_LATLON_VR_CESM_28_HIST = dim_avg_n_Wrap(HIST_SWE_VR_CESM_28_rough, (/1,2/)) printVarSummary(SWE_MASK_LATLON_VR_CESM_28_HIST) SWE_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT = dim_avg_n_Wrap(RCP85_MIDCENT_SWE_VR_CESM_28_rough, (/1,2/)) printVarSummary(SWE_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) SWE_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT = dim_avg_n_Wrap(RCP85_ENDCENT_SWE_VR_CESM_28_rough, (/1,2/)) printVarSummary(SNOWC_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) TSA_MASK_LATLON_VR_CESM_28_HIST = dim_avg_n_Wrap(HIST_TSA_VR_CESM_28_rough, (/1,2/)) printVarSummary(TSA_MASK_LATLON_VR_CESM_28_HIST) TSA_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT = dim_avg_n_Wrap(RCP85_MIDCENT_TSA_VR_CESM_28_rough, (/1,2/)) printVarSummary(TSA_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) TSA_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT = dim_avg_n_Wrap(RCP85_ENDCENT_TSA_VR_CESM_28_rough, (/1,2/)) printVarSummary(TSA_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) ;------------------------------------------------------ ; LINEAR REGRESSION - Sierra Nevada ;------------------------------------------------------ yyyy_HIST = fspan(1980,2005,25) yyyy_MIDCENT = fspan(2025,2050,25) yyyy_ENDCENT = fspan(2075,2100,25) SEASREG_SNOWFALL_VR_CESM_28_HIST_SNOWFALLAVG = regline(yyyy_HIST,SNOWFALL_MASK_LATLON_VR_CESM_28_HIST) print(SEASREG_SNOWFALL_VR_CESM_28_HIST_SNOWFALLAVG) SEASREG_SNOWFALL_VR_CESM_28_RCP85_MIDCENT_SNOWFALLAVG = regline(yyyy_MIDCENT,SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) print(SEASREG_SNOWFALL_VR_CESM_28_RCP85_MIDCENT_SNOWFALLAVG) SEASREG_SNOWFALL_VR_CESM_28_RCP85_ENDCENT_SNOWFALLAVG = regline(yyyy_ENDCENT,SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) print(SEASREG_SNOWFALL_VR_CESM_28_RCP85_ENDCENT_SNOWFALLAVG) SEASREG_SNOWC_VR_CESM_28_HIST_SNOWCAVG = regline(yyyy_HIST,SNOWC_MASK_LATLON_VR_CESM_28_HIST) print(SEASREG_SNOWC_VR_CESM_28_HIST_SNOWCAVG) SEASREG_SNOWC_VR_CESM_28_RCP85_MIDCENT_SNOWCAVG = regline(yyyy_MIDCENT,SNOWC_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) print(SEASREG_SNOWC_VR_CESM_28_RCP85_MIDCENT_SNOWCAVG) SEASREG_SNOWC_VR_CESM_28_RCP85_ENDCENT_SNOWCAVG = regline(yyyy_ENDCENT,SNOWC_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) print(SEASREG_SNOWC_VR_CESM_28_RCP85_ENDCENT_SNOWCAVG) SEASREG_SWE_VR_CESM_28_HIST_SWEAVG = regline(yyyy_HIST,SWE_MASK_LATLON_VR_CESM_28_HIST) print(SEASREG_SWE_VR_CESM_28_HIST_SWEAVG) SEASREG_SWE_VR_CESM_28_RCP85_MIDCENT_SWEAVG = regline(yyyy_MIDCENT,SWE_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) print(SEASREG_SWE_VR_CESM_28_RCP85_MIDCENT_SWEAVG) SEASREG_SWE_VR_CESM_28_RCP85_ENDCENT_SWEAVG = regline(yyyy_ENDCENT,SWE_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) print(SEASREG_SWE_VR_CESM_28_RCP85_ENDCENT_SWEAVG) SEASREG_TSA_VR_CESM_28_HIST_TSAAVG = regline(yyyy_HIST,TSA_MASK_LATLON_VR_CESM_28_HIST) print(SEASREG_TSA_VR_CESM_28_HIST_TSAAVG) SEASREG_TSA_VR_CESM_28_RCP85_MIDCENT_TSAAVG = regline(yyyy_MIDCENT,TSA_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) print(SEASREG_TSA_VR_CESM_28_RCP85_MIDCENT_TSAAVG) SEASREG_TSA_VR_CESM_28_RCP85_ENDCENT_TSAAVG = regline(yyyy_ENDCENT,TSA_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) print(SEASREG_TSA_VR_CESM_28_RCP85_ENDCENT_TSAAVG) ;------------------------------------------------------ ; NORMALIZE VARIABLES AND COMPARE CORRELATION COEFFICIENTS ;------------------------------------------------------ SNOWFALL_NORMAL_VR_CESM_28_HIST = dim_standardize_n_Wrap(SNOWFALL_MASK_LATLON_VR_CESM_28_HIST,0,0) print(SNOWFALL_NORMAL_VR_CESM_28_HIST) SNOWFALL_NORMAL_VR_CESM_28_RCP85_MIDCENT = dim_standardize_n_Wrap(SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT,0,0) print(SNOWFALL_NORMAL_VR_CESM_28_RCP85_MIDCENT) SNOWFALL_NORMAL_VR_CESM_28_RCP85_ENDCENT = dim_standardize_n_Wrap(SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT,0,0) print(SNOWFALL_NORMAL_VR_CESM_28_RCP85_ENDCENT) SNOWC_NORMAL_VR_CESM_28_HIST = dim_standardize_n_Wrap(SNOWC_MASK_LATLON_VR_CESM_28_HIST,0,0) print(SNOWC_NORMAL_VR_CESM_28_HIST) SNOWC_NORMAL_VR_CESM_28_RCP85_MIDCENT = dim_standardize_n_Wrap(SNOWC_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT,0,0) print(SNOWC_NORMAL_VR_CESM_28_RCP85_MIDCENT) SNOWC_NORMAL_VR_CESM_28_RCP85_ENDCENT = dim_standardize_n_Wrap(SNOWC_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT,0,0) print(SNOWC_NORMAL_VR_CESM_28_RCP85_ENDCENT) SWE_NORMAL_VR_CESM_28_HIST = dim_standardize_n_Wrap(SWE_MASK_LATLON_VR_CESM_28_HIST,0,0) print(SWE_NORMAL_VR_CESM_28_HIST) SWE_NORMAL_VR_CESM_28_RCP85_MIDCENT = dim_standardize_n_Wrap(SWE_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT,0,0) print(SWE_NORMAL_VR_CESM_28_RCP85_MIDCENT) SWE_NORMAL_VR_CESM_28_RCP85_ENDCENT = dim_standardize_n_Wrap(SWE_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT,0,0) print(SWE_NORMAL_VR_CESM_28_RCP85_ENDCENT) TSA_NORMAL_VR_CESM_28_HIST = dim_standardize_n_Wrap(TSA_MASK_LATLON_VR_CESM_28_HIST,0,0) print(TSA_NORMAL_VR_CESM_28_HIST) TSA_NORMAL_VR_CESM_28_RCP85_MIDCENT = dim_standardize_n_Wrap(TSA_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT,0,0) print(TSA_NORMAL_VR_CESM_28_RCP85_MIDCENT) TSA_NORMAL_VR_CESM_28_RCP85_ENDCENT = dim_standardize_n_Wrap(TSA_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT,0,0) print(TSA_NORMAL_VR_CESM_28_RCP85_ENDCENT) ;------------------------------------------------------ ; CALCULATE REGRESSION COEFFICIENTS ;------------------------------------------------------ SNOWFALL_REGCOEFF_VR_CESM_28_HIST = regline_stats(yyyy_HIST,SNOWFALL_MASK_LATLON_VR_CESM_28_HIST) print(SNOWFALL_REGCOEFF_VR_CESM_28_HIST) SNOWFALL_REGCOEFF_VR_CESM_28_RCP85_MIDCENT = regline_stats(yyyy_MIDCENT,SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) print(SNOWFALL_REGCOEFF_VR_CESM_28_RCP85_MIDCENT) SNOWFALL_REGCOEFF_VR_CESM_28_RCP85_ENDCENT = regline_stats(yyyy_ENDCENT,SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) print(SNOWFALL_REGCOEFF_VR_CESM_28_RCP85_ENDCENT) SNOWC_REGCOEFF_VR_CESM_28_HIST = regline_stats(yyyy_HIST,SNOWC_MASK_LATLON_VR_CESM_28_HIST) print(SNOWC_REGCOEFF_VR_CESM_28_HIST) SNOWC_REGCOEFF_VR_CESM_28_RCP85_MIDCENT = regline_stats(yyyy_MIDCENT,SNOWC_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) print(SNOWC_REGCOEFF_VR_CESM_28_RCP85_MIDCENT) SNOWC_REGCOEFF_VR_CESM_28_RCP85_ENDCENT = regline_stats(yyyy_ENDCENT,SNOWC_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) print(SNOWC_REGCOEFF_VR_CESM_28_RCP85_ENDCENT) SWE_REGCOEFF_VR_CESM_28_HIST = regline_stats(yyyy_HIST,SWE_MASK_LATLON_VR_CESM_28_HIST) print(SWE_REGCOEFF_VR_CESM_28_HIST) SWE_REGCOEFF_VR_CESM_28_RCP85_MIDCENT = regline_stats(yyyy_MIDCENT,SWE_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) print(SWE_REGCOEFF_VR_CESM_28_RCP85_MIDCENT) SWE_REGCOEFF_VR_CESM_28_RCP85_ENDCENT = regline_stats(yyyy_ENDCENT,SWE_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) print(SWE_REGCOEFF_VR_CESM_28_RCP85_ENDCENT) TSA_REGCOEFF_VR_CESM_28_HIST = regline_stats(yyyy_HIST,TSA_MASK_LATLON_VR_CESM_28_HIST) print(TSA_REGCOEFF_VR_CESM_28_HIST) TSA_REGCOEFF_VR_CESM_28_RCP85_MIDCENT = regline_stats(yyyy_MIDCENT,TSA_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT) print(TSA_REGCOEFF_VR_CESM_28_RCP85_MIDCENT) TSA_REGCOEFF_VR_CESM_28_RCP85_ENDCENT = regline_stats(yyyy_ENDCENT,TSA_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT) print(TSA_REGCOEFF_VR_CESM_28_RCP85_ENDCENT) ;----------------------------------------------------- ; STORE STANDARDIZED REGRESSION COEFFICIENTS ;------------------------------------------------------ SNOWFALL_STANDREG_VR_CESM_28_HIST = SNOWFALL_REGCOEFF_VR_CESM_28_HIST@bstd(1) SNOWFALL_STANDREG_VR_CESM_28_RCP85_MIDCENT = SNOWFALL_REGCOEFF_VR_CESM_28_RCP85_MIDCENT@bstd(1) SNOWFALL_STANDREG_VR_CESM_28_RCP85_ENDCENT = SNOWFALL_REGCOEFF_VR_CESM_28_RCP85_ENDCENT@bstd(1) SNOWC_STANDREG_VR_CESM_28_HIST = SNOWC_REGCOEFF_VR_CESM_28_HIST@bstd(1) SNOWC_STANDREG_VR_CESM_28_RCP85_MIDCENT = SNOWC_REGCOEFF_VR_CESM_28_RCP85_MIDCENT@bstd(1) SNOWC_STANDREG_VR_CESM_28_RCP85_ENDCENT = SNOWC_REGCOEFF_VR_CESM_28_RCP85_ENDCENT@bstd(1) SWE_STANDREG_VR_CESM_28_HIST = SWE_REGCOEFF_VR_CESM_28_HIST@bstd(1) SWE_STANDREG_VR_CESM_28_RCP85_MIDCENT = SWE_REGCOEFF_VR_CESM_28_RCP85_MIDCENT@bstd(1) SWE_STANDREG_VR_CESM_28_RCP85_ENDCENT = SWE_REGCOEFF_VR_CESM_28_RCP85_ENDCENT@bstd(1) TSA_STANDREG_VR_CESM_28_HIST = TSA_REGCOEFF_VR_CESM_28_HIST@bstd(1) TSA_STANDREG_VR_CESM_28_RCP85_MIDCENT = TSA_REGCOEFF_VR_CESM_28_RCP85_MIDCENT@bstd(1) TSA_STANDREG_VR_CESM_28_RCP85_ENDCENT = TSA_REGCOEFF_VR_CESM_28_RCP85_ENDCENT@bstd(1) ;--------------------------------------------------- ;STORE ARRAY VALUES FOR SIERRA NEVADA MASK REGION ;PLOT ARRAY MUST BE (/6,time/) ; 6 is for the MODEL/OBS SNOWFALL (sum or avg), and trendline coefficients*time (Pearson's and/or Spearman's) ;--------------------------------------------------- ; AVERAGE SNOWFALL nyrs = (/25/) ydata_SNOWFALL_AVG = new((/6,nyrs/),float,1e+36) print(ydata_SNOWFALL_AVG) ydata_SNOWFALL_AVG(0,:) = SNOWFALL_MASK_LATLON_VR_CESM_28_HIST ydata_SNOWFALL_AVG(1,:) = SEASREG_SNOWFALL_VR_CESM_28_HIST_SNOWFALLAVG*(yyyy_HIST-SEASREG_SNOWFALL_VR_CESM_28_HIST_SNOWFALLAVG@xave) + SEASREG_SNOWFALL_VR_CESM_28_HIST_SNOWFALLAVG@yave ydata_SNOWFALL_AVG(2,:) = SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT ydata_SNOWFALL_AVG(3,:) = SEASREG_SNOWFALL_VR_CESM_28_RCP85_MIDCENT_SNOWFALLAVG*(yyyy_MIDCENT-SEASREG_SNOWFALL_VR_CESM_28_RCP85_MIDCENT_SNOWFALLAVG@xave) + SEASREG_SNOWFALL_VR_CESM_28_RCP85_MIDCENT_SNOWFALLAVG@yave ydata_SNOWFALL_AVG(4,:) = SNOWFALL_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT ydata_SNOWFALL_AVG(5,:) = SEASREG_SNOWFALL_VR_CESM_28_RCP85_ENDCENT_SNOWFALLAVG*(yyyy_ENDCENT-SEASREG_SNOWFALL_VR_CESM_28_RCP85_ENDCENT_SNOWFALLAVG@xave) + SEASREG_SNOWFALL_VR_CESM_28_RCP85_ENDCENT_SNOWFALLAVG@yave print(ydata_SNOWFALL_AVG) ydata_SNOWC_AVG = new((/6,nyrs/),float,1e+36) print(ydata_SNOWC_AVG) ydata_SNOWC_AVG(0,:) = SNOWC_MASK_LATLON_VR_CESM_28_HIST ydata_SNOWC_AVG(1,:) = SEASREG_SNOWC_VR_CESM_28_HIST_SNOWCAVG*(yyyy_HIST-SEASREG_SNOWC_VR_CESM_28_HIST_SNOWCAVG@xave) + SEASREG_SNOWC_VR_CESM_28_HIST_SNOWCAVG@yave ydata_SNOWC_AVG(2,:) = SNOWC_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT ydata_SNOWC_AVG(3,:) = SEASREG_SNOWC_VR_CESM_28_RCP85_MIDCENT_SNOWCAVG*(yyyy_MIDCENT-SEASREG_SNOWC_VR_CESM_28_RCP85_MIDCENT_SNOWCAVG@xave) + SEASREG_SNOWC_VR_CESM_28_RCP85_MIDCENT_SNOWCAVG@yave ydata_SNOWC_AVG(4,:) = SNOWC_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT ydata_SNOWC_AVG(5,:) = SEASREG_SNOWC_VR_CESM_28_RCP85_ENDCENT_SNOWCAVG*(yyyy_ENDCENT-SEASREG_SNOWC_VR_CESM_28_RCP85_ENDCENT_SNOWCAVG@xave) + SEASREG_SNOWC_VR_CESM_28_RCP85_ENDCENT_SNOWCAVG@yave print(ydata_SNOWC_AVG) ydata_SWE_AVG = new((/6,nyrs/),float,1e+36) print(ydata_SWE_AVG) ydata_SWE_AVG(0,:) = SWE_MASK_LATLON_VR_CESM_28_HIST ydata_SWE_AVG(1,:) = SEASREG_SWE_VR_CESM_28_HIST_SWEAVG*(yyyy_HIST-SEASREG_SWE_VR_CESM_28_HIST_SWEAVG@xave) + SEASREG_SWE_VR_CESM_28_HIST_SWEAVG@yave ydata_SWE_AVG(2,:) = SWE_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT ydata_SWE_AVG(3,:) = SEASREG_SWE_VR_CESM_28_RCP85_MIDCENT_SWEAVG*(yyyy_MIDCENT-SEASREG_SWE_VR_CESM_28_RCP85_MIDCENT_SWEAVG@xave) + SEASREG_SWE_VR_CESM_28_RCP85_MIDCENT_SWEAVG@yave ydata_SWE_AVG(4,:) = SWE_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT ydata_SWE_AVG(5,:) = SEASREG_SWE_VR_CESM_28_RCP85_ENDCENT_SWEAVG*(yyyy_ENDCENT-SEASREG_SWE_VR_CESM_28_RCP85_ENDCENT_SWEAVG@xave) + SEASREG_SWE_VR_CESM_28_RCP85_ENDCENT_SWEAVG@yave print(ydata_SWE_AVG) ydata_TSA_AVG = new((/6,nyrs/),float,1e+36) print(ydata_TSA_AVG) ydata_TSA_AVG(0,:) = TSA_MASK_LATLON_VR_CESM_28_HIST ydata_TSA_AVG(1,:) = SEASREG_TSA_VR_CESM_28_HIST_TSAAVG*(yyyy_HIST-SEASREG_TSA_VR_CESM_28_HIST_TSAAVG@xave) + SEASREG_TSA_VR_CESM_28_HIST_TSAAVG@yave ydata_TSA_AVG(2,:) = TSA_MASK_LATLON_VR_CESM_28_RCP85_MIDCENT ydata_TSA_AVG(3,:) = SEASREG_TSA_VR_CESM_28_RCP85_MIDCENT_TSAAVG*(yyyy_MIDCENT-SEASREG_TSA_VR_CESM_28_RCP85_MIDCENT_TSAAVG@xave) + SEASREG_TSA_VR_CESM_28_RCP85_MIDCENT_TSAAVG@yave ydata_TSA_AVG(4,:) = TSA_MASK_LATLON_VR_CESM_28_RCP85_ENDCENT ydata_TSA_AVG(5,:) = SEASREG_TSA_VR_CESM_28_RCP85_ENDCENT_TSAAVG*(yyyy_ENDCENT-SEASREG_TSA_VR_CESM_28_RCP85_ENDCENT_TSAAVG@xave) + SEASREG_TSA_VR_CESM_28_RCP85_ENDCENT_TSAAVG@yave print(ydata_TSA_AVG) ;--------------------------------------------------- ; Set Plot Resources ;--------------------------------------------------- res = True ; plot mods desired ; res@gsnMaximize = True res@gsnPaperOrientation = "portrait" ; X AND Y LABELS res@tiYAxisFontHeightF = 0.025 ; Y-axis string size (X/0.025 for scaling) res@tiXAxisFontHeightF = 0.025 ; Y-axis string size (X/0.025 for scaling) res@tiMainFontHeightF = 0.025 ; Title string size (X/0.025 for scaling) res@tmXBLabelFontHeightF = 0.025 ; Label font height (X/0.020 for scaling) res@tmYLLabelFontHeightF = 0.025 ; Label font height (X/0.020 for scaling) res@xyDashPatterns = 0 ; solid line res@xyLineColors = (/"black", "red"/) res@xyLineThicknesses = (/1,3/) ; res@vpHeightF = 0.4 ; change aspect ratio of plot ; res@vpWidthF = 0.8 ; res@vpXF = 0.1 ; start plot at x ndc coord ; ADD GRID TO PLOT ; res@tmXMajorGrid = True ; implement x grid ; res@tmXMajorGridThicknessF = 0.5 ; 2.0 is default ; res@tmXMajorGridLineDashPattern = 2 ; select short dash lines res@tmYMajorGrid = True ; implement y grid res@tmYMajorGridThicknessF = 0.5 ; 2.0 is default res@tmYMajorGridLineDashPattern = 2 ; select short dash lines ; SET MAX AND MINS FOR X AND Y AXIS ; res@trXMinF = 1980 ; res@trXMaxF = 2005 ; res@trYMinF = 0 ; res@trYMaxF = 10 ;----------------------------------------------------- ; INDIVIDUAL PLOTTING RESOURCES ;----------------------------------------------------- res@gsnCenterString = "Sierra Nevada" plot_SNOWFALL_avg = new (3,"graphic") plot_SNOWC_avg = new (3,"graphic") plot_SWE_avg = new (3,"graphic") plot_TSA_avg = new (3,"graphic") res@gsnDraw = False res@gsnFrame = False res@tmXBLabelAngleF = 45 ; tilt the XB labels 45 degrees res@tmXBLabelDeltaF = 1.5 ; push the XB labels down from plot res_HIST = res res_HIST@trXMinF = 1980 res_HIST@trXMaxF = 2005 res_HIST@tmXBMode = "Explicit" res_HIST@tmXBValues = (/1980,1985,1990,1995,2000,2005/) res_HIST@tmXBLabels = "" + res_HIST@tmXBValues res_HIST@tmXBMinorValues = ispan(1980,2005,1) res_HIST_SNOWFALL = res_HIST res_HIST_SNOWFALL@tiMainString = "VR-CESM 0.25~S~o~N~ - Historical" res_HIST_SNOWFALL@tiYAxisString = "Snowfall (mm/day)" ; add an Y axis title res_HIST_SNOWFALL@trYMinF = 0 res_HIST_SNOWFALL@trYMaxF = 10 plot_SNOWFALL_avg(0) = gsn_csm_xy(wks,yyyy_HIST,ydata_SNOWFALL_AVG(0:1,:),res_HIST_SNOWFALL) ; create plot_SNOWFALL_avg res_HIST_SNOWC = res_HIST res_HIST_SNOWC@tiMainString = "" res_HIST_SNOWC@tiYAxisString = "Snow Cover (%)" ; add an Y axis title res_HIST_SNOWC@trYMinF = 0 res_HIST_SNOWC@trYMaxF = 100 plot_SNOWC_avg(0) = gsn_csm_xy(wks,yyyy_HIST,ydata_SNOWC_AVG(0:1,:),res_HIST_SNOWC) ; create plot_SNOWC_avg res_HIST_SWE = res_HIST res_HIST_SWE@tiMainString = "" res_HIST_SWE@trYMinF = 0 res_HIST_SWE@trYMaxF = 400 res_HIST_SWE@tiYAxisString = "Snow Water Equivalent (mm)" ; add an Y axis title plot_SWE_avg(0) = gsn_csm_xy(wks,yyyy_HIST,ydata_SWE_AVG(0:1,:),res_HIST_SWE) ; create plot_SWE_avg res_HIST_TSA = res_HIST res_HIST_TSA@tiMainString = "" res_HIST_TSA@tiYAxisString = "Surface Temperature (K)" ; add an Y axis title res_HIST_TSA@trYMinF = 268 res_HIST_TSA@trYMaxF = 285 plot_TSA_avg(0) = gsn_csm_xy(wks,yyyy_HIST,ydata_TSA_AVG(0:1,:),res_HIST_TSA) ; create plot_TSA_avg res_MIDCENT = res res_MIDCENT@tmYLLabelsOn = False ; add an Y axis title res_MIDCENT@trXMinF = 2025 res_MIDCENT@trXMaxF = 2050 res_MIDCENT@tmXBMode = "Explicit" res_MIDCENT@tmXBValues = (/2025,2030,2035,2040,2045,2050/) res_MIDCENT@tmXBLabels = "" + res_MIDCENT@tmXBValues res_MIDCENT@tmXBMinorValues = ispan(2025,2050,1) res_MIDCENT_SNOWFALL = res_MIDCENT res_MIDCENT_SNOWFALL@tiMainString = "VR-CESM 0.25~S~o~N~ - RCP85 Mid-century" res_MIDCENT_SNOWFALL@tiYAxisString = "" ; add an Y axis title res_MIDCENT_SNOWFALL@trYMinF = 0 res_MIDCENT_SNOWFALL@trYMaxF = 10 plot_SNOWFALL_avg(1) = gsn_csm_xy(wks,yyyy_MIDCENT,ydata_SNOWFALL_AVG(2:3,:),res_MIDCENT_SNOWFALL) ; create plot_SNOWFALL_avg res_MIDCENT_SNOWC = res_MIDCENT res_MIDCENT_SNOWC@tiMainString = "" res_MIDCENT_SNOWC@tiYAxisString = "" ; add an Y axis title res_MIDCENT_SNOWC@trYMinF = 0 res_MIDCENT_SNOWC@trYMaxF = 100 plot_SNOWC_avg(1) = gsn_csm_xy(wks,yyyy_MIDCENT,ydata_SNOWC_AVG(2:3,:),res_MIDCENT_SNOWC) ; create plot_SNOWC_avg res_MIDCENT_SWE = res_MIDCENT res_MIDCENT_SWE@tiMainString = "" res_MIDCENT_SWE@tiYAxisString = "" ; add an Y axis title res_MIDCENT_SWE@trYMinF = 0 res_MIDCENT_SWE@trYMaxF = 400 plot_SWE_avg(1) = gsn_csm_xy(wks,yyyy_MIDCENT,ydata_SWE_AVG(2:3,:),res_MIDCENT_SWE) ; create plot_SWE_avg res_MIDCENT_TSA = res_MIDCENT res_MIDCENT_TSA@tiMainString = "" res_MIDCENT_TSA@tiYAxisString = "" ; add an Y axis title res_MIDCENT_TSA@trYMinF = 268 res_MIDCENT_TSA@trYMaxF = 285 plot_TSA_avg(1) = gsn_csm_xy(wks,yyyy_MIDCENT,ydata_TSA_AVG(2:3,:),res_MIDCENT_TSA) ; create plot_TSA_avg res_ENDCENT = res res_ENDCENT@tmYLLabelsOn = False ; add an Y axis title res_ENDCENT@trXMinF = 2075 res_ENDCENT@trXMaxF = 2100 res_ENDCENT@tmXBMode = "Explicit" res_ENDCENT@tmXBValues = (/2075,2080,2085,2090,2095,2100/) res_ENDCENT@tmXBLabels = "" + res_ENDCENT@tmXBValues res_ENDCENT@tmXBMinorValues = ispan(2075,2100,1) res_ENDCENT_SNOWFALL = res_ENDCENT res_ENDCENT_SNOWFALL@tiMainString = "VR-CESM 0.25~S~o~N~ - RCP85 End-century" res_ENDCENT_SNOWFALL@tiYAxisString = "" ; add an Y axis title res_ENDCENT_SNOWFALL@trYMinF = 0 res_ENDCENT_SNOWFALL@trYMaxF = 10 plot_SNOWFALL_avg(2) = gsn_csm_xy(wks,yyyy_ENDCENT,ydata_SNOWFALL_AVG(4:5,:),res_ENDCENT_SNOWFALL) ; create plot_SNOWFALL_avg res_ENDCENT_SNOWC = res_ENDCENT res_ENDCENT_SNOWC@tiMainString = "" res_ENDCENT_SNOWC@tiYAxisString = "" ; add an Y axis title res_ENDCENT_SNOWC@trYMinF = 0 res_ENDCENT_SNOWC@trYMaxF = 100 plot_SNOWC_avg(2) = gsn_csm_xy(wks,yyyy_ENDCENT,ydata_SNOWC_AVG(4:5,:),res_ENDCENT_SNOWC) ; create plot_SNOWC_avg res_ENDCENT_SWE = res_ENDCENT res_ENDCENT_SWE@tiMainString = "" res_ENDCENT_SWE@tiYAxisString = "" ; add an Y axis title res_ENDCENT_SWE@trYMinF = 0 res_ENDCENT_SWE@trYMaxF = 400 plot_SWE_avg(2) = gsn_csm_xy(wks,yyyy_ENDCENT,ydata_SWE_AVG(4:5,:),res_ENDCENT_SWE) ; create plot_SWE_avg res_ENDCENT_TSA = res_ENDCENT res_ENDCENT_TSA@tiMainString = "" res_ENDCENT_TSA@tiYAxisString = "" ; add an Y axis title res_ENDCENT_TSA@trYMinF = 268 res_ENDCENT_TSA@trYMaxF = 285 plot_TSA_avg(2) = gsn_csm_xy(wks,yyyy_ENDCENT,ydata_TSA_AVG(4:5,:),res_ENDCENT_TSA) ; create plot_TSA_avg ;--------------------------------------------------- ; PANEL PLOT RESOURCES ;--------------------------------------------------- resP = True resP@gsnPaperOrientation = "portrait" resP@gsnMaximize = True resP@gsnPanelRowSpec = True resP@gsnPanelYWhiteSpacePercent = 0.0 resP@gsnPanelXWhiteSpacePercent = 0.0 ; ADD LABELS TO PLOTS resP@gsnPanelRowSpec = True ; resP@gsnPanelFigureStringsFontHeightF = 0.0085 resP@amJust = "TopLeft" resP_SNOWFALL = resP resP_SNOWFALL@gsnPanelFigureStrings = (/"a) Standardized Coefficient = "+ SNOWFALL_STANDREG_VR_CESM_28_HIST,"b) Standardized Coefficient = "+ SNOWFALL_STANDREG_VR_CESM_28_RCP85_MIDCENT,"c) Standardized Coefficient = "+ SNOWFALL_STANDREG_VR_CESM_28_RCP85_ENDCENT/) ; add strings to panel resP_SNOWFALL@gsnPanelFigureStringsFontHeightF = 0.0060 resP_SNOWFALL@gsnPanelLeft = 0.05 resP_SNOWFALL@gsnPanelRight = 0.95 resP_SNOWFALL@gsnPanelTop = 0.90 resP_SNOWFALL@gsnPanelBottom = 0.70 resP_SNOWFALL@gsnFrame = False ; resP_SNOWFALL@gsnDraw = False gsn_panel(wks,plot_SNOWFALL_avg,(/3,1/),resP_SNOWFALL) ; now draw as one plot_SNOWFALL resP_SNOWC = resP resP_SNOWC@gsnPanelFigureStrings = (/"a) Standardized Coefficient = "+ SNOWC_STANDREG_VR_CESM_28_HIST,"b) Standardized Coefficient = "+ SNOWC_STANDREG_VR_CESM_28_RCP85_MIDCENT,"c) Standardized Coefficient = "+ SNOWC_STANDREG_VR_CESM_28_RCP85_ENDCENT/) ; add strings to panel resP_SNOWC@gsnPanelFigureStringsFontHeightF = 0.0060 resP_SNOWC@gsnPanelLeft = 0.05 resP_SNOWC@gsnPanelRight = 0.95 resP_SNOWC@gsnPanelTop = 0.70 resP_SNOWC@gsnPanelBottom = 0.50 resP_SNOWC@gsnFrame = False ; resP_SNOWC@gsnDraw = False gsn_panel(wks,plot_SNOWC_avg,(/3,1/),resP_SNOWC) ; now draw as one plot_SNOWC resP_SWE = resP resP_SWE@gsnPanelFigureStrings = (/"a) Standardized Coefficient = "+ SWE_STANDREG_VR_CESM_28_HIST,"b) Standardized Coefficient = "+ SWE_STANDREG_VR_CESM_28_RCP85_MIDCENT,"c) Standardized Coefficient = "+ SWE_STANDREG_VR_CESM_28_RCP85_ENDCENT/) ; add strings to panel resP_SWE@gsnPanelFigureStringsFontHeightF = 0.0060 resP_SWE@gsnPanelLeft = 0.05 resP_SWE@gsnPanelRight = 0.95 resP_SWE@gsnPanelTop = 0.50 resP_SWE@gsnPanelBottom = 0.30 resP_SWE@gsnFrame = False ; resP_SWE@gsnDraw = False gsn_panel(wks,plot_SWE_avg,(/3,1/),resP_SWE) ; now draw as one plot_SWE resP_TSA = resP resP_TSA@gsnPanelFigureStrings = (/"a) Standardized Coefficient = "+ TSA_STANDREG_VR_CESM_28_HIST,"b) Standardized Coefficient = "+ TSA_STANDREG_VR_CESM_28_RCP85_MIDCENT,"c) Standardized Coefficient = "+ TSA_STANDREG_VR_CESM_28_RCP85_ENDCENT/) ; add strings to panel resP_TSA@gsnPanelFigureStringsFontHeightF = 0.0060 resP_TSA@gsnPanelLeft = 0.05 resP_TSA@gsnPanelRight = 0.95 resP_TSA@gsnPanelTop = 0.30 resP_TSA@gsnPanelBottom = 0.10 resP_TSA@gsnFrame = False ; resP_TSA@gsnDraw = False gsn_panel(wks,plot_TSA_avg,(/3,1/),resP_TSA) ; now draw as one plot_TSA ;---This will maximize the size of all the paneled stuff. ; maximize_output(wks,True) ;------------------------------------------------------ print("Done Calculating Statistics") end