;************************************************* ; ; 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" ;------------------------------------------------------ ; ADD CLIMATE AVERAGE MODEL DATA ;------------------------------------------------------ fils = systemfunc("ls /home/arhoades/monthlydata/MASKED_SN/SN_mask_SWE_ALL.nc") fils1 = systemfunc("ls /home/arhoades/monthlydata/MASKED_SN/SN_mask_SWE_GFDL.nc") f_data = addfile(fils,"r") f_GFDL = addfile(fils1,"r") ;------------------------------------------------------ ; GRAB SWE VARIABLES ;------------------------------------------------------ ;UNIGRID14 - SWE SWE_mask_array_UNIGRID14 = f_data->SWE_mask_array_UNIGRID14 SWE_mask_array_UNIGRID14@_FillValue = 1e+36 SWE_mask_array_UNIGRID14!0 = "time" SWE_mask_array_UNIGRID14!1 = "lat" SWE_mask_array_UNIGRID14!2 = "lon" lat1d_mask_array_UNIGRID14 = fspan(33,43,81) lon1d_mask_array_UNIGRID14 = fspan(-125,-115,81) SWE_mask_array_UNIGRID14&lat = lat1d_mask_array_UNIGRID14 SWE_mask_array_UNIGRID14&lon = lon1d_mask_array_UNIGRID14 SWE_mask_array_UNIGRID14&lat@units = "degrees_north" SWE_mask_array_UNIGRID14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_UNIGRID14) printMinMax(SWE_mask_array_UNIGRID14,True) ;------------------------------------------------------ ;UNIGRID28smooth - SWE SWE_mask_array_UNIGRID28smooth = f_data->SWE_mask_array_UNIGRID28smooth SWE_mask_array_UNIGRID28smooth@_FillValue = 1e+36 SWE_mask_array_UNIGRID28smooth!0 = "time" SWE_mask_array_UNIGRID28smooth!1 = "lat" SWE_mask_array_UNIGRID28smooth!2 = "lon" lat1d_mask_array_UNIGRID28smooth = fspan(33,43,41) lon1d_mask_array_UNIGRID28smooth = fspan(-125,-115,41) SWE_mask_array_UNIGRID28smooth&lat = lat1d_mask_array_UNIGRID28smooth SWE_mask_array_UNIGRID28smooth&lon = lon1d_mask_array_UNIGRID28smooth SWE_mask_array_UNIGRID28smooth&lat@units = "degrees_north" SWE_mask_array_UNIGRID28smooth&lon@units = "degrees_east" printVarSummary(SWE_mask_array_UNIGRID28smooth) printMinMax(SWE_mask_array_UNIGRID28smooth,True) ;------------------------------------------------------ ;UNIGRID28rough - SWE SWE_mask_array_UNIGRID28rough = f_data->SWE_mask_array_UNIGRID28rough SWE_mask_array_UNIGRID28rough@_FillValue = 1e+36 SWE_mask_array_UNIGRID28rough!0 = "time" SWE_mask_array_UNIGRID28rough!1 = "lat" SWE_mask_array_UNIGRID28rough!2 = "lon" lat1d_mask_array_UNIGRID28rough = fspan(33,43,41) lon1d_mask_array_UNIGRID28rough = fspan(-125,-115,41) SWE_mask_array_UNIGRID28rough&lat = lat1d_mask_array_UNIGRID28rough SWE_mask_array_UNIGRID28rough&lon = lon1d_mask_array_UNIGRID28rough SWE_mask_array_UNIGRID28rough&lat@units = "degrees_north" SWE_mask_array_UNIGRID28rough&lon@units = "degrees_east" printVarSummary(SWE_mask_array_UNIGRID28rough) printMinMax(SWE_mask_array_UNIGRID28rough,True) ;------------------------------------------------------ ;VCAM - SWE SWE_mask_array_VCAM = f_data->SWE_mask_array_VCAM SWE_mask_array_VCAM@_FillValue = 1e+36 SWE_mask_array_VCAM!0 = "time" SWE_mask_array_VCAM!1 = "lat" SWE_mask_array_VCAM!2 = "lon" lat1d_mask_array_VCAM = fspan(33,43,41) lon1d_mask_array_VCAM = fspan(-125,-115,41) SWE_mask_array_VCAM&lat = lat1d_mask_array_VCAM SWE_mask_array_VCAM&lon = lon1d_mask_array_VCAM SWE_mask_array_VCAM&lat@units = "degrees_north" SWE_mask_array_VCAM&lon@units = "degrees_east" printVarSummary(SWE_mask_array_VCAM) printMinMax(SWE_mask_array_VCAM,True) ;------------------------------------------------------ ;UNIFORM28 - SWE SWE_mask_array_UNIFORM28 = f_data->SWE_mask_array_UNIFORM28 SWE_mask_array_UNIFORM28@_FillValue = 1e+36 SWE_mask_array_UNIFORM28!0 = "time" SWE_mask_array_UNIFORM28!1 = "lat" SWE_mask_array_UNIFORM28!2 = "lon" lat1d_mask_array_UNIFORM28 = fspan(33,43,41) lon1d_mask_array_UNIFORM28 = fspan(-125,-115,41) SWE_mask_array_UNIFORM28&lat = lat1d_mask_array_UNIFORM28 SWE_mask_array_UNIFORM28&lon = lon1d_mask_array_UNIFORM28 SWE_mask_array_UNIFORM28&lat@units = "degrees_north" SWE_mask_array_UNIFORM28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_UNIFORM28) printMinMax(SWE_mask_array_UNIFORM28,True) ;------------------------------------------------------ ;UNIFORM111 - SWE SWE_mask_array_UNIFORM111 = f_data->SWE_mask_array_UNIFORM111 SWE_mask_array_UNIFORM111@_FillValue = 1e+36 SWE_mask_array_UNIFORM111!0 = "time" SWE_mask_array_UNIFORM111!1 = "lat" SWE_mask_array_UNIFORM111!2 = "lon" lat1d_mask_array_UNIFORM111 = fspan(33,43,11) lon1d_mask_array_UNIFORM111 = fspan(-125,-115,11) SWE_mask_array_UNIFORM111&lat = lat1d_mask_array_UNIFORM111 SWE_mask_array_UNIFORM111&lon = lon1d_mask_array_UNIFORM111 SWE_mask_array_UNIFORM111&lat@units = "degrees_north" SWE_mask_array_UNIFORM111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_UNIFORM111) printMinMax(SWE_mask_array_UNIFORM111,True) ;------------------------------------------------------ ; WRF9 - SWE SWE_mask_array_WRF9 = f_data->SWE_mask_array_WRF9 SWE_mask_array_WRF9@_FillValue = 9.96921e+36 SWE_mask_array_WRF9!0 = "time" SWE_mask_array_WRF9!1 = "lat" SWE_mask_array_WRF9!2 = "lon" lat1d_mask_array_WRF9 = fspan(33,43,81) lon1d_mask_array_WRF9 = fspan(-125,-115,81) SWE_mask_array_WRF9&lat = lat1d_mask_array_WRF9 SWE_mask_array_WRF9&lon = lon1d_mask_array_WRF9 SWE_mask_array_WRF9&lat@units = "degrees_north" SWE_mask_array_WRF9&lon@units = "degrees_east" printVarSummary(SWE_mask_array_WRF9) printMinMax(SWE_mask_array_WRF9,True) ;------------------------------------------------------ ; WRF27 - SWE SWE_mask_array_WRF27 = f_data->SWE_mask_array_WRF27 SWE_mask_array_WRF27@_FillValue = 9.96921e+36 SWE_mask_array_WRF27!0 = "time" SWE_mask_array_WRF27!1 = "lat" SWE_mask_array_WRF27!2 = "lon" lat1d_mask_array_WRF27 = fspan(33,43,41) lon1d_mask_array_WRF27 = fspan(-125,-115,41) SWE_mask_array_WRF27&lat = lat1d_mask_array_WRF27 SWE_mask_array_WRF27&lon = lon1d_mask_array_WRF27 SWE_mask_array_WRF27&lat@units = "degrees_north" SWE_mask_array_WRF27&lon@units = "degrees_east" printVarSummary(SWE_mask_array_WRF27) printMinMax(SWE_mask_array_WRF27,True) ;------------------------------------------------------ ; DAYMET (1km,14km,28km,111km) - SWE SWE_mask_array_DAYMET1 = f_data->SWE_mask_DAYMET1 SWE_mask_array_DAYMET1@_FillValue = -9999 SWE_mask_array_DAYMET1!0 = "time" SWE_mask_array_DAYMET1!1 = "lat" SWE_mask_array_DAYMET1!2 = "lon" lat1d_mask_array_DAYMET1 = fspan(33,43,3001) lon1d_mask_array_DAYMET1 = fspan(-125,-115,1201) SWE_mask_array_DAYMET1&lat = lat1d_mask_array_DAYMET1 SWE_mask_array_DAYMET1&lon = lon1d_mask_array_DAYMET1 SWE_mask_array_DAYMET1&lat@units = "degrees_north" SWE_mask_array_DAYMET1&lon@units = "degrees_east" printVarSummary(SWE_mask_array_DAYMET1) printMinMax(SWE_mask_array_DAYMET1,True) SWE_mask_array_DAYMET14 = f_data->SWE_mask_array_DAYMET14 SWE_mask_array_DAYMET14@_FillValue = -9999 SWE_mask_array_DAYMET14!0 = "time" SWE_mask_array_DAYMET14!1 = "lat" SWE_mask_array_DAYMET14!2 = "lon" lat1d_mask_array_DAYMET14 = fspan(33,43,81) lon1d_mask_array_DAYMET14 = fspan(-125,-115,81) SWE_mask_array_DAYMET14&lat = lat1d_mask_array_DAYMET14 SWE_mask_array_DAYMET14&lon = lon1d_mask_array_DAYMET14 SWE_mask_array_DAYMET14&lat@units = "degrees_north" SWE_mask_array_DAYMET14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_DAYMET14) printMinMax(SWE_mask_array_DAYMET14,True) SWE_mask_array_DAYMET28 = f_data->SWE_mask_array_DAYMET28 SWE_mask_array_DAYMET28@_FillValue = -9999 SWE_mask_array_DAYMET28!0 = "time" SWE_mask_array_DAYMET28!1 = "lat" SWE_mask_array_DAYMET28!2 = "lon" lat1d_mask_array_DAYMET28 = fspan(33,43,41) lon1d_mask_array_DAYMET28 = fspan(-125,-115,41) SWE_mask_array_DAYMET28&lat = lat1d_mask_array_DAYMET28 SWE_mask_array_DAYMET28&lon = lon1d_mask_array_DAYMET28 SWE_mask_array_DAYMET28&lat@units = "degrees_north" SWE_mask_array_DAYMET28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_DAYMET28) printMinMax(SWE_mask_array_DAYMET28,True) SWE_mask_array_DAYMET111 = f_data->SWE_mask_array_DAYMET111 SWE_mask_array_DAYMET111@_FillValue = -9999 SWE_mask_array_DAYMET111!0 = "time" SWE_mask_array_DAYMET111!1 = "lat" SWE_mask_array_DAYMET111!2 = "lon" lat1d_mask_array_DAYMET111 = fspan(33,43,11) lon1d_mask_array_DAYMET111 = fspan(-125,-115,11) SWE_mask_array_DAYMET111&lat = lat1d_mask_array_DAYMET111 SWE_mask_array_DAYMET111&lon = lon1d_mask_array_DAYMET111 SWE_mask_array_DAYMET111&lat@units = "degrees_north" SWE_mask_array_DAYMET111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_DAYMET111) printMinMax(SWE_mask_array_DAYMET111,True) ;------------------------------------------------------ ; NLDAS VIC (14km,28km,111km) - SWE SWE_mask_array_NLDAS_VIC_14 = f_data->SWE_mask_array_NLDAS_VIC_14 SWE_mask_array_NLDAS_VIC_14@_FillValue = 1e+20 SWE_mask_array_NLDAS_VIC_14!0 = "time" SWE_mask_array_NLDAS_VIC_14!1 = "lat" SWE_mask_array_NLDAS_VIC_14!2 = "lon" lat1d_mask_array_NLDAS_VIC_14 = fspan(33,43,81) lon1d_mask_array_NLDAS_VIC_14 = fspan(-125,-115,81) SWE_mask_array_NLDAS_VIC_14&lat = lat1d_mask_array_NLDAS_VIC_14 SWE_mask_array_NLDAS_VIC_14&lon = lon1d_mask_array_NLDAS_VIC_14 SWE_mask_array_NLDAS_VIC_14&lat@units = "degrees_north" SWE_mask_array_NLDAS_VIC_14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_VIC_14) printMinMax(SWE_mask_array_NLDAS_VIC_14,True) SWE_mask_array_NLDAS_VIC_28 = f_data->SWE_mask_array_NLDAS_VIC_28 SWE_mask_array_NLDAS_VIC_28@_FillValue = 1e+20 SWE_mask_array_NLDAS_VIC_28!0 = "time" SWE_mask_array_NLDAS_VIC_28!1 = "lat" SWE_mask_array_NLDAS_VIC_28!2 = "lon" lat1d_mask_array_NLDAS_VIC_28 = fspan(33,43,41) lon1d_mask_array_NLDAS_VIC_28 = fspan(-125,-115,41) SWE_mask_array_NLDAS_VIC_28&lat = lat1d_mask_array_NLDAS_VIC_28 SWE_mask_array_NLDAS_VIC_28&lon = lon1d_mask_array_NLDAS_VIC_28 SWE_mask_array_NLDAS_VIC_28&lat@units = "degrees_north" SWE_mask_array_NLDAS_VIC_28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_VIC_28) printMinMax(SWE_mask_array_NLDAS_VIC_28,True) SWE_mask_array_NLDAS_VIC_111 = f_data->SWE_mask_array_NLDAS_VIC_111 SWE_mask_array_NLDAS_VIC_111@_FillValue = 1e+20 SWE_mask_array_NLDAS_VIC_111!0 = "time" SWE_mask_array_NLDAS_VIC_111!1 = "lat" SWE_mask_array_NLDAS_VIC_111!2 = "lon" lat1d_mask_array_NLDAS_VIC_111 = fspan(33,43,11) lon1d_mask_array_NLDAS_VIC_111 = fspan(-125,-115,11) SWE_mask_array_NLDAS_VIC_111&lat = lat1d_mask_array_NLDAS_VIC_111 SWE_mask_array_NLDAS_VIC_111&lon = lon1d_mask_array_NLDAS_VIC_111 SWE_mask_array_NLDAS_VIC_111&lat@units = "degrees_north" SWE_mask_array_NLDAS_VIC_111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_VIC_111) printMinMax(SWE_mask_array_NLDAS_VIC_111,True) ;------------------------------------------------------ ; NLDAS NOAH (14km,28km,111km) - SWE SWE_mask_array_NLDAS_NOAH_14 = f_data->SWE_mask_array_NLDAS_NOAH_14 SWE_mask_array_NLDAS_NOAH_14@_FillValue = 1e+20 SWE_mask_array_NLDAS_NOAH_14!0 = "time" SWE_mask_array_NLDAS_NOAH_14!1 = "lat" SWE_mask_array_NLDAS_NOAH_14!2 = "lon" lat1d_mask_array_NLDAS_NOAH_14 = fspan(33,43,81) lon1d_mask_array_NLDAS_NOAH_14 = fspan(-125,-115,81) SWE_mask_array_NLDAS_NOAH_14&lat = lat1d_mask_array_NLDAS_NOAH_14 SWE_mask_array_NLDAS_NOAH_14&lon = lon1d_mask_array_NLDAS_NOAH_14 SWE_mask_array_NLDAS_NOAH_14&lat@units = "degrees_north" SWE_mask_array_NLDAS_NOAH_14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_NOAH_14) printMinMax(SWE_mask_array_NLDAS_NOAH_14,True) SWE_mask_array_NLDAS_NOAH_28 = f_data->SWE_mask_array_NLDAS_NOAH_28 SWE_mask_array_NLDAS_NOAH_28@_FillValue = 1e+20 SWE_mask_array_NLDAS_NOAH_28!0 = "time" SWE_mask_array_NLDAS_NOAH_28!1 = "lat" SWE_mask_array_NLDAS_NOAH_28!2 = "lon" lat1d_mask_array_NLDAS_NOAH_28 = fspan(33,43,41) lon1d_mask_array_NLDAS_NOAH_28 = fspan(-125,-115,41) SWE_mask_array_NLDAS_NOAH_28&lat = lat1d_mask_array_NLDAS_NOAH_28 SWE_mask_array_NLDAS_NOAH_28&lon = lon1d_mask_array_NLDAS_NOAH_28 SWE_mask_array_NLDAS_NOAH_28&lat@units = "degrees_north" SWE_mask_array_NLDAS_NOAH_28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_NOAH_28) printMinMax(SWE_mask_array_NLDAS_NOAH_28,True) SWE_mask_array_NLDAS_NOAH_111 = f_data->SWE_mask_array_NLDAS_NOAH_111 SWE_mask_array_NLDAS_NOAH_111@_FillValue = 1e+20 SWE_mask_array_NLDAS_NOAH_111!0 = "time" SWE_mask_array_NLDAS_NOAH_111!1 = "lat" SWE_mask_array_NLDAS_NOAH_111!2 = "lon" lat1d_mask_array_NLDAS_NOAH_111 = fspan(33,43,11) lon1d_mask_array_NLDAS_NOAH_111 = fspan(-125,-115,11) SWE_mask_array_NLDAS_NOAH_111&lat = lat1d_mask_array_NLDAS_NOAH_111 SWE_mask_array_NLDAS_NOAH_111&lon = lon1d_mask_array_NLDAS_NOAH_111 SWE_mask_array_NLDAS_NOAH_111&lat@units = "degrees_north" SWE_mask_array_NLDAS_NOAH_111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_NOAH_111) printMinMax(SWE_mask_array_NLDAS_NOAH_111,True) ;------------------------------------------------------ ; NLDAS MOSAIC (14km,28km,111km) - SWE SWE_mask_array_NLDAS_MOSAIC_14 = f_data->SWE_mask_array_NLDAS_MOSAIC_14 SWE_mask_array_NLDAS_MOSAIC_14@_FillValue = 1e+20 SWE_mask_array_NLDAS_MOSAIC_14!0 = "time" SWE_mask_array_NLDAS_MOSAIC_14!1 = "lat" SWE_mask_array_NLDAS_MOSAIC_14!2 = "lon" lat1d_mask_array_NLDAS_MOSAIC_14 = fspan(33,43,81) lon1d_mask_array_NLDAS_MOSAIC_14 = fspan(-125,-115,81) SWE_mask_array_NLDAS_MOSAIC_14&lat = lat1d_mask_array_NLDAS_MOSAIC_14 SWE_mask_array_NLDAS_MOSAIC_14&lon = lon1d_mask_array_NLDAS_MOSAIC_14 SWE_mask_array_NLDAS_MOSAIC_14&lat@units = "degrees_north" SWE_mask_array_NLDAS_MOSAIC_14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_MOSAIC_14) printMinMax(SWE_mask_array_NLDAS_MOSAIC_14,True) SWE_mask_array_NLDAS_MOSAIC_28 = f_data->SWE_mask_array_NLDAS_MOSAIC_28 SWE_mask_array_NLDAS_MOSAIC_28@_FillValue = 1e+20 SWE_mask_array_NLDAS_MOSAIC_28!0 = "time" SWE_mask_array_NLDAS_MOSAIC_28!1 = "lat" SWE_mask_array_NLDAS_MOSAIC_28!2 = "lon" lat1d_mask_array_NLDAS_MOSAIC_28 = fspan(33,43,41) lon1d_mask_array_NLDAS_MOSAIC_28 = fspan(-125,-115,41) SWE_mask_array_NLDAS_MOSAIC_28&lat = lat1d_mask_array_NLDAS_MOSAIC_28 SWE_mask_array_NLDAS_MOSAIC_28&lon = lon1d_mask_array_NLDAS_MOSAIC_28 SWE_mask_array_NLDAS_MOSAIC_28&lat@units = "degrees_north" SWE_mask_array_NLDAS_MOSAIC_28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_MOSAIC_28) printMinMax(SWE_mask_array_NLDAS_MOSAIC_28,True) SWE_mask_array_NLDAS_MOSAIC_111 = f_data->SWE_mask_array_NLDAS_MOSAIC_111 SWE_mask_array_NLDAS_MOSAIC_111@_FillValue = 1e+20 SWE_mask_array_NLDAS_MOSAIC_111!0 = "time" SWE_mask_array_NLDAS_MOSAIC_111!1 = "lat" SWE_mask_array_NLDAS_MOSAIC_111!2 = "lon" lat1d_mask_array_NLDAS_MOSAIC_111 = fspan(33,43,11) lon1d_mask_array_NLDAS_MOSAIC_111 = fspan(-125,-115,11) SWE_mask_array_NLDAS_MOSAIC_111&lat = lat1d_mask_array_NLDAS_MOSAIC_111 SWE_mask_array_NLDAS_MOSAIC_111&lon = lon1d_mask_array_NLDAS_MOSAIC_111 SWE_mask_array_NLDAS_MOSAIC_111&lat@units = "degrees_north" SWE_mask_array_NLDAS_MOSAIC_111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_NLDAS_MOSAIC_111) printMinMax(SWE_mask_array_NLDAS_MOSAIC_111,True) ;------------------------------------------------------ ; CAL-ADAPT CCSM3 (14km,28km,111km) - SWE SWE_mask_array_CALADAPT_CCSM3_14 = f_data->SWE_mask_array_CALADAPT_CCSM3_14 SWE_mask_array_CALADAPT_CCSM3_14@_FillValue = -9999 SWE_mask_array_CALADAPT_CCSM3_14!0 = "time" SWE_mask_array_CALADAPT_CCSM3_14!1 = "lat" SWE_mask_array_CALADAPT_CCSM3_14!2 = "lon" lat1d_mask_array_CALADAPT_CCSM3_14 = fspan(33,43,81) lon1d_mask_array_CALADAPT_CCSM3_14 = fspan(-125,-115,81) SWE_mask_array_CALADAPT_CCSM3_14&lat = lat1d_mask_array_CALADAPT_CCSM3_14 SWE_mask_array_CALADAPT_CCSM3_14&lon = lon1d_mask_array_CALADAPT_CCSM3_14 SWE_mask_array_CALADAPT_CCSM3_14&lat@units = "degrees_north" SWE_mask_array_CALADAPT_CCSM3_14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_CCSM3_14) printMinMax(SWE_mask_array_CALADAPT_CCSM3_14,True) SWE_mask_array_CALADAPT_CCSM3_28 = f_data->SWE_mask_array_CALADAPT_CCSM3_28 SWE_mask_array_CALADAPT_CCSM3_28@_FillValue = -9999 SWE_mask_array_CALADAPT_CCSM3_28!0 = "time" SWE_mask_array_CALADAPT_CCSM3_28!1 = "lat" SWE_mask_array_CALADAPT_CCSM3_28!2 = "lon" lat1d_mask_array_CALADAPT_CCSM3_28 = fspan(33,43,41) lon1d_mask_array_CALADAPT_CCSM3_28 = fspan(-125,-115,41) SWE_mask_array_CALADAPT_CCSM3_28&lat = lat1d_mask_array_CALADAPT_CCSM3_28 SWE_mask_array_CALADAPT_CCSM3_28&lon = lon1d_mask_array_CALADAPT_CCSM3_28 SWE_mask_array_CALADAPT_CCSM3_28&lat@units = "degrees_north" SWE_mask_array_CALADAPT_CCSM3_28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_CCSM3_28) printMinMax(SWE_mask_array_CALADAPT_CCSM3_28,True) SWE_mask_array_CALADAPT_CCSM3_111 = f_data->SWE_mask_array_CALADAPT_CCSM3_111 SWE_mask_array_CALADAPT_CCSM3_111@_FillValue = -9999 SWE_mask_array_CALADAPT_CCSM3_111!0 = "time" SWE_mask_array_CALADAPT_CCSM3_111!1 = "lat" SWE_mask_array_CALADAPT_CCSM3_111!2 = "lon" lat1d_mask_array_CALADAPT_CCSM3_111 = fspan(33,43,11) lon1d_mask_array_CALADAPT_CCSM3_111 = fspan(-125,-115,11) SWE_mask_array_CALADAPT_CCSM3_111&lat = lat1d_mask_array_CALADAPT_CCSM3_111 SWE_mask_array_CALADAPT_CCSM3_111&lon = lon1d_mask_array_CALADAPT_CCSM3_111 SWE_mask_array_CALADAPT_CCSM3_111&lat@units = "degrees_north" SWE_mask_array_CALADAPT_CCSM3_111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_CCSM3_111) printMinMax(SWE_mask_array_CALADAPT_CCSM3_111,True) ;------------------------------------------------------ ; CAL-ADAPT CNRM (14km,28km,111km) - SWE SWE_mask_array_CALADAPT_CNRM_14 = f_data->SWE_mask_array_CALADAPT_CNRM_14 SWE_mask_array_CALADAPT_CNRM_14@_FillValue = -9999 SWE_mask_array_CALADAPT_CNRM_14!0 = "time" SWE_mask_array_CALADAPT_CNRM_14!1 = "lat" SWE_mask_array_CALADAPT_CNRM_14!2 = "lon" lat1d_mask_array_CALADAPT_CNRM_14 = fspan(33,43,81) lon1d_mask_array_CALADAPT_CNRM_14 = fspan(-125,-115,81) SWE_mask_array_CALADAPT_CNRM_14&lat = lat1d_mask_array_CALADAPT_CNRM_14 SWE_mask_array_CALADAPT_CNRM_14&lon = lon1d_mask_array_CALADAPT_CNRM_14 SWE_mask_array_CALADAPT_CNRM_14&lat@units = "degrees_north" SWE_mask_array_CALADAPT_CNRM_14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_CNRM_14) printMinMax(SWE_mask_array_CALADAPT_CNRM_14,True) SWE_mask_array_CALADAPT_CNRM_28 = f_data->SWE_mask_array_CALADAPT_CNRM_28 SWE_mask_array_CALADAPT_CNRM_28@_FillValue = -9999 SWE_mask_array_CALADAPT_CNRM_28!0 = "time" SWE_mask_array_CALADAPT_CNRM_28!1 = "lat" SWE_mask_array_CALADAPT_CNRM_28!2 = "lon" lat1d_mask_array_CALADAPT_CNRM_28 = fspan(33,43,41) lon1d_mask_array_CALADAPT_CNRM_28 = fspan(-125,-115,41) SWE_mask_array_CALADAPT_CNRM_28&lat = lat1d_mask_array_CALADAPT_CNRM_28 SWE_mask_array_CALADAPT_CNRM_28&lon = lon1d_mask_array_CALADAPT_CNRM_28 SWE_mask_array_CALADAPT_CNRM_28&lat@units = "degrees_north" SWE_mask_array_CALADAPT_CNRM_28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_CNRM_28) printMinMax(SWE_mask_array_CALADAPT_CNRM_28,True) SWE_mask_array_CALADAPT_CNRM_111 = f_data->SWE_mask_array_CALADAPT_CNRM_111 SWE_mask_array_CALADAPT_CNRM_111@_FillValue = -9999 SWE_mask_array_CALADAPT_CNRM_111!0 = "time" SWE_mask_array_CALADAPT_CNRM_111!1 = "lat" SWE_mask_array_CALADAPT_CNRM_111!2 = "lon" lat1d_mask_array_CALADAPT_CNRM_111 = fspan(33,43,11) lon1d_mask_array_CALADAPT_CNRM_111 = fspan(-125,-115,11) SWE_mask_array_CALADAPT_CNRM_111&lat = lat1d_mask_array_CALADAPT_CNRM_111 SWE_mask_array_CALADAPT_CNRM_111&lon = lon1d_mask_array_CALADAPT_CNRM_111 SWE_mask_array_CALADAPT_CNRM_111&lat@units = "degrees_north" SWE_mask_array_CALADAPT_CNRM_111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_CNRM_111) printMinMax(SWE_mask_array_CALADAPT_CNRM_111,True) ;------------------------------------------------------ ; CAL-ADAPT GFDL (14km,28km,111km) - SWE SWE_mask_array_CALADAPT_GFDL_14 = f_GFDL->SWE_mask_array_CALADAPT_GFDL_14 SWE_mask_array_CALADAPT_GFDL_14@_FillValue = -9999 SWE_mask_array_CALADAPT_GFDL_14!0 = "time" SWE_mask_array_CALADAPT_GFDL_14!1 = "lat" SWE_mask_array_CALADAPT_GFDL_14!2 = "lon" lat1d_mask_array_CALADAPT_GFDL_14 = fspan(33,43,81) lon1d_mask_array_CALADAPT_GFDL_14 = fspan(-125,-115,81) SWE_mask_array_CALADAPT_GFDL_14&lat = lat1d_mask_array_CALADAPT_GFDL_14 SWE_mask_array_CALADAPT_GFDL_14&lon = lon1d_mask_array_CALADAPT_GFDL_14 SWE_mask_array_CALADAPT_GFDL_14&lat@units = "degrees_north" SWE_mask_array_CALADAPT_GFDL_14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_GFDL_14) printMinMax(SWE_mask_array_CALADAPT_GFDL_14,True) SWE_mask_array_CALADAPT_GFDL_28 = f_GFDL->SWE_mask_array_CALADAPT_GFDL_28 SWE_mask_array_CALADAPT_GFDL_28@_FillValue = -9999 SWE_mask_array_CALADAPT_GFDL_28!0 = "time" SWE_mask_array_CALADAPT_GFDL_28!1 = "lat" SWE_mask_array_CALADAPT_GFDL_28!2 = "lon" lat1d_mask_array_CALADAPT_GFDL_28 = fspan(33,43,41) lon1d_mask_array_CALADAPT_GFDL_28 = fspan(-125,-115,41) SWE_mask_array_CALADAPT_GFDL_28&lat = lat1d_mask_array_CALADAPT_GFDL_28 SWE_mask_array_CALADAPT_GFDL_28&lon = lon1d_mask_array_CALADAPT_GFDL_28 SWE_mask_array_CALADAPT_GFDL_28&lat@units = "degrees_north" SWE_mask_array_CALADAPT_GFDL_28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_GFDL_28) printMinMax(SWE_mask_array_CALADAPT_GFDL_28,True) SWE_mask_array_CALADAPT_GFDL_111 = f_GFDL->SWE_mask_array_CALADAPT_GFDL_111 SWE_mask_array_CALADAPT_GFDL_111@_FillValue = -9999 SWE_mask_array_CALADAPT_GFDL_111!0 = "time" SWE_mask_array_CALADAPT_GFDL_111!1 = "lat" SWE_mask_array_CALADAPT_GFDL_111!2 = "lon" lat1d_mask_array_CALADAPT_GFDL_111 = fspan(33,43,11) lon1d_mask_array_CALADAPT_GFDL_111 = fspan(-125,-115,11) SWE_mask_array_CALADAPT_GFDL_111&lat = lat1d_mask_array_CALADAPT_GFDL_111 SWE_mask_array_CALADAPT_GFDL_111&lon = lon1d_mask_array_CALADAPT_GFDL_111 SWE_mask_array_CALADAPT_GFDL_111&lat@units = "degrees_north" SWE_mask_array_CALADAPT_GFDL_111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_GFDL_111) printMinMax(SWE_mask_array_CALADAPT_GFDL_111,True) ;------------------------------------------------------ ; CAL-ADAPT PCM1 (14km,28km,111km) - SWE SWE_mask_array_CALADAPT_PCM1_14 = f_data->SWE_mask_array_CALADAPT_PCM1_14 SWE_mask_array_CALADAPT_PCM1_14@_FillValue = -9999 SWE_mask_array_CALADAPT_PCM1_14!0 = "time" SWE_mask_array_CALADAPT_PCM1_14!1 = "lat" SWE_mask_array_CALADAPT_PCM1_14!2 = "lon" lat1d_mask_array_CALADAPT_PCM1_14 = fspan(33,43,81) lon1d_mask_array_CALADAPT_PCM1_14 = fspan(-125,-115,81) SWE_mask_array_CALADAPT_PCM1_14&lat = lat1d_mask_array_CALADAPT_PCM1_14 SWE_mask_array_CALADAPT_PCM1_14&lon = lon1d_mask_array_CALADAPT_PCM1_14 SWE_mask_array_CALADAPT_PCM1_14&lat@units = "degrees_north" SWE_mask_array_CALADAPT_PCM1_14&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_PCM1_14) printMinMax(SWE_mask_array_CALADAPT_PCM1_14,True) SWE_mask_array_CALADAPT_PCM1_28 = f_data->SWE_mask_array_CALADAPT_PCM1_28 SWE_mask_array_CALADAPT_PCM1_28@_FillValue = -9999 SWE_mask_array_CALADAPT_PCM1_28!0 = "time" SWE_mask_array_CALADAPT_PCM1_28!1 = "lat" SWE_mask_array_CALADAPT_PCM1_28!2 = "lon" lat1d_mask_array_CALADAPT_PCM1_28 = fspan(33,43,41) lon1d_mask_array_CALADAPT_PCM1_28 = fspan(-125,-115,41) SWE_mask_array_CALADAPT_PCM1_28&lat = lat1d_mask_array_CALADAPT_PCM1_28 SWE_mask_array_CALADAPT_PCM1_28&lon = lon1d_mask_array_CALADAPT_PCM1_28 SWE_mask_array_CALADAPT_PCM1_28&lat@units = "degrees_north" SWE_mask_array_CALADAPT_PCM1_28&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_PCM1_28) printMinMax(SWE_mask_array_CALADAPT_PCM1_28,True) SWE_mask_array_CALADAPT_PCM1_111 = f_data->SWE_mask_array_CALADAPT_PCM1_111 SWE_mask_array_CALADAPT_PCM1_111@_FillValue = -9999 SWE_mask_array_CALADAPT_PCM1_111!0 = "time" SWE_mask_array_CALADAPT_PCM1_111!1 = "lat" SWE_mask_array_CALADAPT_PCM1_111!2 = "lon" lat1d_mask_array_CALADAPT_PCM1_111 = fspan(33,43,11) lon1d_mask_array_CALADAPT_PCM1_111 = fspan(-125,-115,11) SWE_mask_array_CALADAPT_PCM1_111&lat = lat1d_mask_array_CALADAPT_PCM1_111 SWE_mask_array_CALADAPT_PCM1_111&lon = lon1d_mask_array_CALADAPT_PCM1_111 SWE_mask_array_CALADAPT_PCM1_111&lat@units = "degrees_north" SWE_mask_array_CALADAPT_PCM1_111&lon@units = "degrees_east" printVarSummary(SWE_mask_array_CALADAPT_PCM1_111) printMinMax(SWE_mask_array_CALADAPT_PCM1_111,True) ;------------------------------------------------------ ; CHANGE ARRAY TO 1D FOR STAT DISPERSION ;------------------------------------------------------ ;UNIGRID14 oned_UNIGRID14_nosum = ndtooned(SWE_mask_array_UNIGRID14) ;UNIGRID28rough oned_UNIGRID28rough_nosum = ndtooned(SWE_mask_array_UNIGRID28rough) ;UNIGRID28smooth oned_UNIGRID28smooth_nosum = ndtooned(SWE_mask_array_UNIGRID28smooth) ;VCAM oned_VCAM_nosum = ndtooned(SWE_mask_array_VCAM) ;UNIFORM28 oned_UNIFORM28_nosum = ndtooned(SWE_mask_array_UNIFORM28) ;UNIFORM111 oned_UNIFORM111_nosum = ndtooned(SWE_mask_array_UNIFORM111) ;WRF9 oned_WRF9_nosum = ndtooned(SWE_mask_array_WRF9) ;WRF27 oned_WRF27_nosum = ndtooned(SWE_mask_array_WRF27) ;DAYMET oned_DAYMET14_nosum = ndtooned(SWE_mask_array_DAYMET14) oned_DAYMET28_nosum = ndtooned(SWE_mask_array_DAYMET28) oned_DAYMET111_nosum = ndtooned(SWE_mask_array_DAYMET111) oned_DAYMET1_nosum = ndtooned(SWE_mask_array_DAYMET1) ;NLDAS VIC oned_NLDAS_VIC_14_nosum = ndtooned(SWE_mask_array_NLDAS_VIC_14) oned_NLDAS_VIC_28_nosum = ndtooned(SWE_mask_array_NLDAS_VIC_28) oned_NLDAS_VIC_111_nosum = ndtooned(SWE_mask_array_NLDAS_VIC_111) ;NLDAS NOAH oned_NLDAS_NOAH_14_nosum = ndtooned(SWE_mask_array_NLDAS_NOAH_14) oned_NLDAS_NOAH_28_nosum = ndtooned(SWE_mask_array_NLDAS_NOAH_28) oned_NLDAS_NOAH_111_nosum = ndtooned(SWE_mask_array_NLDAS_NOAH_111) ;NLDAS MOSAIC oned_NLDAS_MOSAIC_14_nosum = ndtooned(SWE_mask_array_NLDAS_MOSAIC_14) oned_NLDAS_MOSAIC_28_nosum = ndtooned(SWE_mask_array_NLDAS_MOSAIC_28) oned_NLDAS_MOSAIC_111_nosum = ndtooned(SWE_mask_array_NLDAS_MOSAIC_111) ;CALADAPT CCSM3 oned_CALADAPT_CCSM3_14_nosum = ndtooned(SWE_mask_array_CALADAPT_CCSM3_14) oned_CALADAPT_CCSM3_28_nosum = ndtooned(SWE_mask_array_CALADAPT_CCSM3_28) oned_CALADAPT_CCSM3_111_nosum = ndtooned(SWE_mask_array_CALADAPT_CCSM3_111) ;CALADAPT CNRM oned_CALADAPT_CNRM_14_nosum = ndtooned(SWE_mask_array_CALADAPT_CNRM_14) oned_CALADAPT_CNRM_28_nosum = ndtooned(SWE_mask_array_CALADAPT_CNRM_28) oned_CALADAPT_CNRM_111_nosum = ndtooned(SWE_mask_array_CALADAPT_CNRM_111) ;CALADAPT GFDL oned_CALADAPT_GFDL_14_nosum = ndtooned(SWE_mask_array_CALADAPT_GFDL_14) oned_CALADAPT_GFDL_28_nosum = ndtooned(SWE_mask_array_CALADAPT_GFDL_28) oned_CALADAPT_GFDL_111_nosum = ndtooned(SWE_mask_array_CALADAPT_GFDL_111) ;CALADAPT PCM1 oned_CALADAPT_PCM1_14_nosum = ndtooned(SWE_mask_array_CALADAPT_PCM1_14) oned_CALADAPT_PCM1_28_nosum = ndtooned(SWE_mask_array_CALADAPT_PCM1_28) oned_CALADAPT_PCM1_111_nosum = ndtooned(SWE_mask_array_CALADAPT_PCM1_111) ;------------------------------------------------------ ; STAT DISPERSION ;------------------------------------------------------ ; OPTS REQUIRED FOR STAT DISPERSION ;---------------------------------- opt = True opt@PrintStat = True ;---------------------------------- ;UNIGRID14 print("UNIGRID14 - DJF Season Sierra SWE (mm)") Sierra_Stat_UNIGRID14 = stat_dispersion(oned_UNIGRID14_nosum,opt) ;UNIGRID28rough print("UNIGRID28rough - DJF Season Sierra SWE (mm)") Sierra_Stat_UNIGRID28rough = stat_dispersion(oned_UNIGRID28rough_nosum,opt) ;UNIGRID28smooth print("UNIGRID28smooth - DJF Season Sierra SWE (mm)") Sierra_Stat_UNIGRID28smooth = stat_dispersion(oned_UNIGRID28smooth_nosum,opt) ;VCAM print("VCAM - DJF Season Sierra SWE (mm)") Sierra_Stat_VCAM = stat_dispersion(oned_VCAM_nosum,opt) ;UNIFORM28 print("UNIFORM28 - DJF Season Sierra SWE (mm)") Sierra_Stat_UNIFORM28 = stat_dispersion(oned_UNIFORM28_nosum,opt) ;UNIFORM111 print("UNIFORM111 - DJF Season Sierra SWE (mm)") Sierra_Stat_UNIFORM111 = stat_dispersion(oned_UNIFORM111_nosum,opt) ;WRF9 print("WRF9 - DJF Season Sierra SWE (mm)") Sierra_Stat_WRF9 = stat_dispersion(oned_WRF9_nosum,opt) ;WRF27 print("WRF27 - DJF Season Sierra SWE (mm)") Sierra_Stat_WRF27 = stat_dispersion(oned_WRF27_nosum,opt) ;DAYMET print("DAYMET14 - DJF Season Sierra SWE (mm)") Sierra_Stat_DAYMET14 = stat_dispersion(oned_DAYMET14_nosum,opt) print("DAYMET28 - DJF Season Sierra SWE (mm)") Sierra_Stat_DAYMET28 = stat_dispersion(oned_DAYMET28_nosum,opt) print("DAYMET111 - DJF Season Sierra SWE (mm)") Sierra_Stat_DAYMET111 = stat_dispersion(oned_DAYMET111_nosum,opt) print("DAYMET1 - DJF Season Sierra SWE (mm)") Sierra_Stat_DAYMET1 = stat_dispersion(oned_DAYMET1_nosum,opt) ;NLDAS VIC print("NLDAS_VIC_14 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_VIC_14 = stat_dispersion(oned_NLDAS_VIC_14_nosum,opt) print("NLDAS_VIC_28 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_VIC_28 = stat_dispersion(oned_NLDAS_VIC_28_nosum,opt) print("NLDAS_VIC_111 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_VIC_111 = stat_dispersion(oned_NLDAS_VIC_111_nosum,opt) ;NLDAS NOAH print("NLDAS_NOAH_14 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_NOAH_14 = stat_dispersion(oned_NLDAS_NOAH_14_nosum,opt) print("NLDAS_NOAH_28 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_NOAH_28 = stat_dispersion(oned_NLDAS_NOAH_28_nosum,opt) print("NLDAS_NOAH_111 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_NOAH_111 = stat_dispersion(oned_NLDAS_NOAH_111_nosum,opt) ;NLDAS MOSAIC print("NLDAS_MOSAIC_14 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_MOSAIC_14 = stat_dispersion(oned_NLDAS_MOSAIC_14_nosum,opt) print("NLDAS_MOSAIC_28 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_MOSAIC_28 = stat_dispersion(oned_NLDAS_MOSAIC_28_nosum,opt) print("NLDAS_MOSAIC_111 - DJF Season Sierra SWE (mm)") Sierra_Stat_NLDAS_MOSAIC_111 = stat_dispersion(oned_NLDAS_MOSAIC_111_nosum,opt) ;CALADAPT CCSM3 print("CALADAPT_CCSM3_14 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_CCSM3_14 = stat_dispersion(oned_CALADAPT_CCSM3_14_nosum,opt) print("CALADAPT_CCSM3_28 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_CCSM3_28 = stat_dispersion(oned_CALADAPT_CCSM3_28_nosum,opt) print("CALADAPT_CCSM3_111 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_CCSM3_111 = stat_dispersion(oned_CALADAPT_CCSM3_111_nosum,opt) ;CALADAPT CNRM print("CALADAPT_CNRM_14 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_CNRM_14 = stat_dispersion(oned_CALADAPT_CNRM_14_nosum,opt) print("CALADAPT_CNRM_28 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_CNRM_28 = stat_dispersion(oned_CALADAPT_CNRM_28_nosum,opt) print("CALADAPT_CNRM_111 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_CNRM_111 = stat_dispersion(oned_CALADAPT_CNRM_111_nosum,opt) ;CALADAPT GFDL print("CALADAPT_GFDL_14 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_GFDL_14 = stat_dispersion(oned_CALADAPT_GFDL_14_nosum,opt) print("CALADAPT_GFDL_28 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_GFDL_28 = stat_dispersion(oned_CALADAPT_GFDL_28_nosum,opt) print("CALADAPT_GFDL_111 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_GFDL_111 = stat_dispersion(oned_CALADAPT_GFDL_111_nosum,opt) ;CALADAPT PCM1 print("CALADAPT_PCM1_14 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_PCM1_14 = stat_dispersion(oned_CALADAPT_PCM1_14_nosum,opt) print("CALADAPT_PCM1_28 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_PCM1_28 = stat_dispersion(oned_CALADAPT_PCM1_28_nosum,opt) print("CALADAPT_PCM1_111 - DJF Season Sierra SWE (mm)") Sierra_Stat_CALADAPT_PCM1_111 = stat_dispersion(oned_CALADAPT_PCM1_111_nosum,opt) ;BOX PLOT - MIN, 25TH PERCENTILE, MEDIAN, 75TH PERCENTILE, MAX UNIGRID14_MIN = Sierra_Stat_UNIGRID14(2) UNIGRID14_25percent = Sierra_Stat_UNIGRID14(6) UNIGRID14_MEDIAN = Sierra_Stat_UNIGRID14(8) UNIGRID14_75percent = Sierra_Stat_UNIGRID14(10) UNIGRID14_MAX = Sierra_Stat_UNIGRID14(14) UNIGRID28rough_MIN = Sierra_Stat_UNIGRID28rough(2) UNIGRID28rough_25percent = Sierra_Stat_UNIGRID28rough(6) UNIGRID28rough_MEDIAN = Sierra_Stat_UNIGRID28rough(8) UNIGRID28rough_75percent = Sierra_Stat_UNIGRID28rough(10) UNIGRID28rough_MAX = Sierra_Stat_UNIGRID28rough(14) UNIGRID28smooth_MIN = Sierra_Stat_UNIGRID28smooth(2) UNIGRID28smooth_25percent = Sierra_Stat_UNIGRID28smooth(6) UNIGRID28smooth_MEDIAN = Sierra_Stat_UNIGRID28smooth(8) UNIGRID28smooth_75percent = Sierra_Stat_UNIGRID28smooth(10) UNIGRID28smooth_MAX = Sierra_Stat_UNIGRID28smooth(14) VCAM_MIN = Sierra_Stat_VCAM(2) VCAM_25percent = Sierra_Stat_VCAM(6) VCAM_MEDIAN = Sierra_Stat_VCAM(8) VCAM_75percent = Sierra_Stat_VCAM(10) VCAM_MAX = Sierra_Stat_VCAM(14) UNIFORM28_MIN = Sierra_Stat_UNIFORM28(2) UNIFORM28_25percent = Sierra_Stat_UNIFORM28(6) UNIFORM28_MEDIAN = Sierra_Stat_UNIFORM28(8) UNIFORM28_75percent = Sierra_Stat_UNIFORM28(10) UNIFORM28_MAX = Sierra_Stat_UNIFORM28(14) UNIFORM111_MIN = Sierra_Stat_UNIFORM111(2) UNIFORM111_25percent = Sierra_Stat_UNIFORM111(6) UNIFORM111_MEDIAN = Sierra_Stat_UNIFORM111(8) UNIFORM111_75percent = Sierra_Stat_UNIFORM111(10) UNIFORM111_MAX = Sierra_Stat_UNIFORM111(14) WRF9_MIN = Sierra_Stat_WRF9(2) WRF9_25percent = Sierra_Stat_WRF9(6) WRF9_MEDIAN = Sierra_Stat_WRF9(8) WRF9_75percent = Sierra_Stat_WRF9(10) WRF9_MAX = Sierra_Stat_WRF9(14) WRF27_MIN = Sierra_Stat_WRF27(2) WRF27_25percent = Sierra_Stat_WRF27(6) WRF27_MEDIAN = Sierra_Stat_WRF27(8) WRF27_75percent = Sierra_Stat_WRF27(10) WRF27_MAX = Sierra_Stat_WRF27(14) DAYMET14_MIN = Sierra_Stat_DAYMET14(2) DAYMET14_25percent = Sierra_Stat_DAYMET14(6) DAYMET14_MEDIAN = Sierra_Stat_DAYMET14(8) DAYMET14_75percent = Sierra_Stat_DAYMET14(10) DAYMET14_MAX = Sierra_Stat_DAYMET14(14) DAYMET28_MIN = Sierra_Stat_DAYMET28(2) DAYMET28_25percent = Sierra_Stat_DAYMET28(6) DAYMET28_MEDIAN = Sierra_Stat_DAYMET28(8) DAYMET28_75percent = Sierra_Stat_DAYMET28(10) DAYMET28_MAX = Sierra_Stat_DAYMET28(14) DAYMET111_MIN = Sierra_Stat_DAYMET111(2) DAYMET111_25percent = Sierra_Stat_DAYMET111(6) DAYMET111_MEDIAN = Sierra_Stat_DAYMET111(8) DAYMET111_75percent = Sierra_Stat_DAYMET111(10) DAYMET111_MAX = Sierra_Stat_DAYMET111(14) DAYMET1_MIN = Sierra_Stat_DAYMET1(2) DAYMET1_25percent = Sierra_Stat_DAYMET1(6) DAYMET1_MEDIAN = Sierra_Stat_DAYMET1(8) DAYMET1_75percent = Sierra_Stat_DAYMET1(10) DAYMET1_MAX = Sierra_Stat_DAYMET1(14) NLDAS_VIC_14_MIN = Sierra_Stat_NLDAS_VIC_14(2) NLDAS_VIC_14_25percent = Sierra_Stat_NLDAS_VIC_14(6) NLDAS_VIC_14_MEDIAN = Sierra_Stat_NLDAS_VIC_14(8) NLDAS_VIC_14_75percent = Sierra_Stat_NLDAS_VIC_14(10) NLDAS_VIC_14_MAX = Sierra_Stat_NLDAS_VIC_14(14) NLDAS_VIC_28_MIN = Sierra_Stat_NLDAS_VIC_28(2) NLDAS_VIC_28_25percent = Sierra_Stat_NLDAS_VIC_28(6) NLDAS_VIC_28_MEDIAN = Sierra_Stat_NLDAS_VIC_28(8) NLDAS_VIC_28_75percent = Sierra_Stat_NLDAS_VIC_28(10) NLDAS_VIC_28_MAX = Sierra_Stat_NLDAS_VIC_28(14) NLDAS_VIC_111_MIN = Sierra_Stat_NLDAS_VIC_111(2) NLDAS_VIC_111_25percent = Sierra_Stat_NLDAS_VIC_111(6) NLDAS_VIC_111_MEDIAN = Sierra_Stat_NLDAS_VIC_111(8) NLDAS_VIC_111_75percent = Sierra_Stat_NLDAS_VIC_111(10) NLDAS_VIC_111_MAX = Sierra_Stat_NLDAS_VIC_111(14) NLDAS_NOAH_14_MIN = Sierra_Stat_NLDAS_NOAH_14(2) NLDAS_NOAH_14_25percent = Sierra_Stat_NLDAS_NOAH_14(6) NLDAS_NOAH_14_MEDIAN = Sierra_Stat_NLDAS_NOAH_14(8) NLDAS_NOAH_14_75percent = Sierra_Stat_NLDAS_NOAH_14(10) NLDAS_NOAH_14_MAX = Sierra_Stat_NLDAS_NOAH_14(14) NLDAS_NOAH_28_MIN = Sierra_Stat_NLDAS_NOAH_28(2) NLDAS_NOAH_28_25percent = Sierra_Stat_NLDAS_NOAH_28(6) NLDAS_NOAH_28_MEDIAN = Sierra_Stat_NLDAS_NOAH_28(8) NLDAS_NOAH_28_75percent = Sierra_Stat_NLDAS_NOAH_28(10) NLDAS_NOAH_28_MAX = Sierra_Stat_NLDAS_NOAH_28(14) NLDAS_NOAH_111_MIN = Sierra_Stat_NLDAS_NOAH_111(2) NLDAS_NOAH_111_25percent = Sierra_Stat_NLDAS_NOAH_111(6) NLDAS_NOAH_111_MEDIAN = Sierra_Stat_NLDAS_NOAH_111(8) NLDAS_NOAH_111_75percent = Sierra_Stat_NLDAS_NOAH_111(10) NLDAS_NOAH_111_MAX = Sierra_Stat_NLDAS_NOAH_111(14) NLDAS_MOSAIC_14_MIN = Sierra_Stat_NLDAS_MOSAIC_14(2) NLDAS_MOSAIC_14_25percent = Sierra_Stat_NLDAS_MOSAIC_14(6) NLDAS_MOSAIC_14_MEDIAN = Sierra_Stat_NLDAS_MOSAIC_14(8) NLDAS_MOSAIC_14_75percent = Sierra_Stat_NLDAS_MOSAIC_14(10) NLDAS_MOSAIC_14_MAX = Sierra_Stat_NLDAS_MOSAIC_14(14) NLDAS_MOSAIC_28_MIN = Sierra_Stat_NLDAS_MOSAIC_28(2) NLDAS_MOSAIC_28_25percent = Sierra_Stat_NLDAS_MOSAIC_28(6) NLDAS_MOSAIC_28_MEDIAN = Sierra_Stat_NLDAS_MOSAIC_28(8) NLDAS_MOSAIC_28_75percent = Sierra_Stat_NLDAS_MOSAIC_28(10) NLDAS_MOSAIC_28_MAX = Sierra_Stat_NLDAS_MOSAIC_28(14) NLDAS_MOSAIC_111_MIN = Sierra_Stat_NLDAS_MOSAIC_111(2) NLDAS_MOSAIC_111_25percent = Sierra_Stat_NLDAS_MOSAIC_111(6) NLDAS_MOSAIC_111_MEDIAN = Sierra_Stat_NLDAS_MOSAIC_111(8) NLDAS_MOSAIC_111_75percent = Sierra_Stat_NLDAS_MOSAIC_111(10) NLDAS_MOSAIC_111_MAX = Sierra_Stat_NLDAS_MOSAIC_111(14) CALADAPT_CCSM3_14_MIN = Sierra_Stat_CALADAPT_CCSM3_14(2) CALADAPT_CCSM3_14_25percent = Sierra_Stat_CALADAPT_CCSM3_14(6) CALADAPT_CCSM3_14_MEDIAN = Sierra_Stat_CALADAPT_CCSM3_14(8) CALADAPT_CCSM3_14_75percent = Sierra_Stat_CALADAPT_CCSM3_14(10) CALADAPT_CCSM3_14_MAX = Sierra_Stat_CALADAPT_CCSM3_14(14) CALADAPT_CCSM3_28_MIN = Sierra_Stat_CALADAPT_CCSM3_28(2) CALADAPT_CCSM3_28_25percent = Sierra_Stat_CALADAPT_CCSM3_28(6) CALADAPT_CCSM3_28_MEDIAN = Sierra_Stat_CALADAPT_CCSM3_28(8) CALADAPT_CCSM3_28_75percent = Sierra_Stat_CALADAPT_CCSM3_28(10) CALADAPT_CCSM3_28_MAX = Sierra_Stat_CALADAPT_CCSM3_28(14) CALADAPT_CCSM3_111_MIN = Sierra_Stat_CALADAPT_CCSM3_111(2) CALADAPT_CCSM3_111_25percent = Sierra_Stat_CALADAPT_CCSM3_111(6) CALADAPT_CCSM3_111_MEDIAN = Sierra_Stat_CALADAPT_CCSM3_111(8) CALADAPT_CCSM3_111_75percent = Sierra_Stat_CALADAPT_CCSM3_111(10) CALADAPT_CCSM3_111_MAX = Sierra_Stat_CALADAPT_CCSM3_111(14) CALADAPT_CNRM_14_MIN = Sierra_Stat_CALADAPT_CNRM_14(2) CALADAPT_CNRM_14_25percent = Sierra_Stat_CALADAPT_CNRM_14(6) CALADAPT_CNRM_14_MEDIAN = Sierra_Stat_CALADAPT_CNRM_14(8) CALADAPT_CNRM_14_75percent = Sierra_Stat_CALADAPT_CNRM_14(10) CALADAPT_CNRM_14_MAX = Sierra_Stat_CALADAPT_CNRM_14(14) CALADAPT_CNRM_28_MIN = Sierra_Stat_CALADAPT_CNRM_28(2) CALADAPT_CNRM_28_25percent = Sierra_Stat_CALADAPT_CNRM_28(6) CALADAPT_CNRM_28_MEDIAN = Sierra_Stat_CALADAPT_CNRM_28(8) CALADAPT_CNRM_28_75percent = Sierra_Stat_CALADAPT_CNRM_28(10) CALADAPT_CNRM_28_MAX = Sierra_Stat_CALADAPT_CNRM_28(14) CALADAPT_CNRM_111_MIN = Sierra_Stat_CALADAPT_CNRM_111(2) CALADAPT_CNRM_111_25percent = Sierra_Stat_CALADAPT_CNRM_111(6) CALADAPT_CNRM_111_MEDIAN = Sierra_Stat_CALADAPT_CNRM_111(8) CALADAPT_CNRM_111_75percent = Sierra_Stat_CALADAPT_CNRM_111(10) CALADAPT_CNRM_111_MAX = Sierra_Stat_CALADAPT_CNRM_111(14) CALADAPT_GFDL_14_MIN = Sierra_Stat_CALADAPT_GFDL_14(2) CALADAPT_GFDL_14_25percent = Sierra_Stat_CALADAPT_GFDL_14(6) CALADAPT_GFDL_14_MEDIAN = Sierra_Stat_CALADAPT_GFDL_14(8) CALADAPT_GFDL_14_75percent = Sierra_Stat_CALADAPT_GFDL_14(10) CALADAPT_GFDL_14_MAX = Sierra_Stat_CALADAPT_GFDL_14(14) CALADAPT_GFDL_28_MIN = Sierra_Stat_CALADAPT_GFDL_28(2) CALADAPT_GFDL_28_25percent = Sierra_Stat_CALADAPT_GFDL_28(6) CALADAPT_GFDL_28_MEDIAN = Sierra_Stat_CALADAPT_GFDL_28(8) CALADAPT_GFDL_28_75percent = Sierra_Stat_CALADAPT_GFDL_28(10) CALADAPT_GFDL_28_MAX = Sierra_Stat_CALADAPT_GFDL_28(14) CALADAPT_GFDL_111_MIN = Sierra_Stat_CALADAPT_GFDL_111(2) CALADAPT_GFDL_111_25percent = Sierra_Stat_CALADAPT_GFDL_111(6) CALADAPT_GFDL_111_MEDIAN = Sierra_Stat_CALADAPT_GFDL_111(8) CALADAPT_GFDL_111_75percent = Sierra_Stat_CALADAPT_GFDL_111(10) CALADAPT_GFDL_111_MAX = Sierra_Stat_CALADAPT_GFDL_111(14) CALADAPT_PCM1_14_MIN = Sierra_Stat_CALADAPT_PCM1_14(2) CALADAPT_PCM1_14_25percent = Sierra_Stat_CALADAPT_PCM1_14(6) CALADAPT_PCM1_14_MEDIAN = Sierra_Stat_CALADAPT_PCM1_14(8) CALADAPT_PCM1_14_75percent = Sierra_Stat_CALADAPT_PCM1_14(10) CALADAPT_PCM1_14_MAX = Sierra_Stat_CALADAPT_PCM1_14(14) CALADAPT_PCM1_28_MIN = Sierra_Stat_CALADAPT_PCM1_28(2) CALADAPT_PCM1_28_25percent = Sierra_Stat_CALADAPT_PCM1_28(6) CALADAPT_PCM1_28_MEDIAN = Sierra_Stat_CALADAPT_PCM1_28(8) CALADAPT_PCM1_28_75percent = Sierra_Stat_CALADAPT_PCM1_28(10) CALADAPT_PCM1_28_MAX = Sierra_Stat_CALADAPT_PCM1_28(14) CALADAPT_PCM1_111_MIN = Sierra_Stat_CALADAPT_PCM1_111(2) CALADAPT_PCM1_111_25percent = Sierra_Stat_CALADAPT_PCM1_111(6) CALADAPT_PCM1_111_MEDIAN = Sierra_Stat_CALADAPT_PCM1_111(8) CALADAPT_PCM1_111_75percent = Sierra_Stat_CALADAPT_PCM1_111(10) CALADAPT_PCM1_111_MAX = Sierra_Stat_CALADAPT_PCM1_111(14) ;--------------------------------------------------- ;STORE ARRAY VALUES FOR SIERRA NEVADA MASK REGION ;BOX PLOT ARRAY MUST BE (/*,5/) ;--------------------------------------------------- array_all = new((/15,5/),float) ydata_BOX_ALL = array_all ; MAKE SURE # OF DJF YEARS ARE UNIFORM IN THE 2ND DIMENSION OF THE ARRAY ydata_BOX_ALL(6,:) = (/UNIGRID14_MIN,UNIGRID14_25percent,UNIGRID14_MEDIAN,UNIGRID14_75percent,UNIGRID14_MAX/) ydata_BOX_ALL(5,:) = (/UNIGRID28rough_MIN,UNIGRID28rough_25percent,UNIGRID28rough_MEDIAN,UNIGRID28rough_75percent,UNIGRID28rough_MAX/) ydata_BOX_ALL(4,:) = (/UNIGRID28smooth_MIN,UNIGRID28smooth_25percent,UNIGRID28smooth_MEDIAN,UNIGRID28smooth_75percent,UNIGRID28smooth_MAX/) ydata_BOX_ALL(1,:) = (/UNIFORM28_MIN,UNIFORM28_25percent,UNIFORM28_MEDIAN,UNIFORM28_75percent,UNIFORM28_MAX/) ydata_BOX_ALL(0,:) = (/UNIFORM111_MIN,UNIFORM111_25percent,UNIFORM111_MEDIAN,UNIFORM111_75percent,UNIFORM111_MAX/) ydata_BOX_ALL(3,:) = (/WRF9_MIN,WRF9_25percent,WRF9_MEDIAN,WRF9_75percent,WRF9_MAX/) ydata_BOX_ALL(2,:) = (/WRF27_MIN,WRF27_25percent,WRF27_MEDIAN,WRF27_75percent,WRF27_MAX/) ; ydata_BOX_ALL(7,:) = (/DAYMET14_MIN,DAYMET14_25percent,DAYMET14_MEDIAN,DAYMET14_75percent,DAYMET14_MAX/) ydata_BOX_ALL(7,:) = (/DAYMET1_MIN,DAYMET1_25percent,DAYMET1_MEDIAN,DAYMET1_75percent,DAYMET1_MAX/) ydata_BOX_ALL(8,:) = (/NLDAS_VIC_14_MIN,NLDAS_VIC_14_25percent,NLDAS_VIC_14_MEDIAN,NLDAS_VIC_14_75percent,NLDAS_VIC_14_MAX/) ydata_BOX_ALL(9,:) = (/NLDAS_NOAH_14_MIN,NLDAS_NOAH_14_25percent,NLDAS_NOAH_14_MEDIAN,NLDAS_NOAH_14_75percent,NLDAS_NOAH_14_MAX/) ydata_BOX_ALL(10,:) = (/NLDAS_MOSAIC_14_MIN,NLDAS_MOSAIC_14_25percent,NLDAS_MOSAIC_14_MEDIAN,NLDAS_MOSAIC_14_75percent,NLDAS_MOSAIC_14_MAX/) ydata_BOX_ALL(11,:) = (/CALADAPT_CCSM3_14_MIN,CALADAPT_CCSM3_14_25percent,CALADAPT_CCSM3_14_MEDIAN,CALADAPT_CCSM3_14_75percent,CALADAPT_CCSM3_14_MAX/) ydata_BOX_ALL(12,:) = (/CALADAPT_CNRM_14_MIN,CALADAPT_CNRM_14_25percent,CALADAPT_CNRM_14_MEDIAN,CALADAPT_CNRM_14_75percent,CALADAPT_CNRM_14_MAX/) ydata_BOX_ALL(13,:) = (/CALADAPT_GFDL_14_MIN,CALADAPT_GFDL_14_25percent,CALADAPT_GFDL_14_MEDIAN,CALADAPT_GFDL_14_75percent,CALADAPT_GFDL_14_MAX/) ydata_BOX_ALL(14,:) = (/CALADAPT_PCM1_14_MIN,CALADAPT_PCM1_14_25percent,CALADAPT_PCM1_14_MEDIAN,CALADAPT_PCM1_14_75percent,CALADAPT_PCM1_14_MAX/) print(ydata_BOX_ALL) printVarSummary(ydata_BOX_ALL) xdata_num_BOX_ALL = dimsizes(ydata_BOX_ALL(:,0)) xdata_BOX_ALL = fspan(1,xdata_num_BOX_ALL,15) ; xdata_BOX_ALL = ispan(1,4*xdata_num_BOX_ALL-1,4) print(xdata_BOX_ALL) ;--------------------------------------------------- ; Open Workstation ;--------------------------------------------------- filename = "SWE_SIERRAMASK_BOXPLOT_ALL_REGRID" wks = gsn_open_wks (type,filename) ; open workstation ;--------------------------------------------------- ; Set Plot Resources ;--------------------------------------------------- res = True ; plot mods desired ; res@gsnMaximize = True ; maximize the plot space on page - this will apply to each individual plot ; X AND Y LABELS res@tiYAxisString = "Snow Water Equivalent (mm)" ; add an Y axis title res@tiXAxisString = "" ; add an X axis title res@tiYAxisFontHeightF = 0.0100 ; Y-axis string size (X/0.025 for scaling) res@tiMainFontHeightF = 0.0100 ; Title string size (X/0.025 for scaling) res@tmXBLabelFontHeightF = 0.0100 ; Label font height (X/0.020 for scaling) res@tmYLLabelFontHeightF = 0.0100 ; Label font height (X/0.020 for scaling) ; X-AXIS LABELS ; res@tmXBLabels = (/"CESM-SE 1~S~o~N~","CESM-FV 0.25~S~o~N~","WRF-27","WRF-9","VR-CESM 0.25~S~o~N~ (s)","VR-CESM 0.25~S~o~N~ (r)","VR-CESM 0.125~S~o~N~","DAYMET 14km","NLDAS VIC 0.125~S~o~N~","NLDAS NOAH 0.125~S~o~N~","NLDAS MOSAIC 0.125~S~o~N~","CAL-ADAPT CCSM3 0.125~S~o~N~","CAL-ADAPT CNRM 0.125~S~o~N~","CAL-ADAPT GFDL 0.125~S~o~N~","CAL-ADAPT PCM1 0.125~S~o~N~"/) res@tmXBLabelAngleF = 90.0 ; Change Direction of Text in x-axis res@tmXBMajorLengthF = 0.001 ; ADD GRID TO PLOT ; res@tmXMajorGrid = True ; implement x grid ; res@tmXMajorGridThicknessF = 1.0 ; 2.0 is default ; res@tmXMajorGridLineDashPattern = 2 ; select short dash lines res@tmYMajorGrid = True ; implement y grid res@tmYMajorGridThicknessF = 1.0 ; 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 = 1250 ;BOX PLOT LINE SETTINGS (COLOR AND WIDTH) ; MUST BE SEPERATE RES SETTINGS opti = True ; opti@boxColors = (/"NavyBlue","CadetBlue","Cyan","Red","Brown","DarkGreen","Gold2","Purple","Orange"/) ; opti@boxColors = (/"DarkGreen","DarkGreen","Red","Red","NavyBlue","NavyBlue","NavyBlue","Orange","Orange","Orange","Orange","Orange","Orange","Orange","Orange"/) opti@boxWidth = 0.25 ; line thicknesses for all lines ;BOX PLOT LINE SETTINGS (THICKNESS) llres = True llres@gsLineThicknessF = 0.5 ; res@tmXBLabelDeltaF = 0.5 ; Move label away from tickmarks. ;----------------------------------------------------- ; INDIVIDUAL PLOTTING RESOURCES ;----------------------------------------------------- ; ALL resCOMP = res ; resCOMP@tiMainString = "AVG DJF SWE in Sierra Mask" resCOMP@tmXBLabels = (/"CESM-SE 1~S~o~N~","CESM-FV 0.25~S~o~N~","WRF-27","WRF-9","VR-CESM 0.25~S~o~N~ (s)","VR-CESM 0.25~S~o~N~ (r)","VR-CESM 0.125~S~o~N~","DAYMET 1km","NLDAS VIC 0.125~S~o~N~","NLDAS NOAH 0.125~S~o~N~","NLDAS MOSAIC 0.125~S~o~N~","CAL-ADAPT CCSM3 0.125~S~o~N~","CAL-ADAPT CNRM 0.125~S~o~N~","CAL-ADAPT GFDL 0.125~S~o~N~","CAL-ADAPT PCM1 0.125~S~o~N~"/) ; resCOMP@gsnPanelBottom = 0.2 ; shrink panel plot by setting bottom edge of plot ; resCOMP@gsnPanelTop = 0.9 ; shrink panel plot by setting top edge of plot ;BOX PLOT LINE SETTINGS (COLOR AND WIDTH) ; MUST BE SEPERATE RES SETTINGS opti_all = True ; opti_all@boxColors = (/"NavyBlue","CadetBlue","Cyan","Red","Brown","DarkGreen","Gold2","Purple","Orange"/) opti_all@boxColors = (/"DarkGreen","DarkGreen","Red","Red","NavyBlue","NavyBlue","NavyBlue","Orange","Purple","Purple","Purple","Cyan","Cyan","Cyan","Cyan"/) opti_all@boxWidth = 0.25 ; line thicknesses for all lines plot_box_ALL = boxplot (wks,xdata_BOX_ALL,ydata_BOX_ALL,opti_all,resCOMP,llres) draw(plot_box_ALL) frame(wks) ;------------------------------------------------------ print("Done Creating Plots") end