;********************************************************************************************************************* ;This script will produce scatter plot of radiation measurements associated with Cloud, Clear, and BLSN cases for 2010 ;********************************************************************************************************************* ; These files still have to be loaded manually load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/contrib/ut_string.ncl" ;load "color_map.ncl" ;************************************************** begin ;**** For 2010, we are reading in the monthly BLSN, Cloud, and Clear LW radiation and surf. temp. per month. set_default_fillvalue("double", -999) set_default_fillvalue("float", -999) set_default_fillvalue("integer", -999) Jan_blsn_depth = (/75, 90/); Jan_blsn_pixels = (/6, 7/); Jan_blsn_domec_depth = (/-999, -999/); Jan_blsn_surflw = (/113, 113/); Jan_blsn_surftemp = (/-25.2, -26.5/); print(getFillValue(Jan_blsn_depth)) print(typeof(Jan_blsn_surftemp)) Jan_cloud_surflw = (/138, 144, 120, 185, 108/); Jan_cloud_surftemp = (/-34, -25.1, -34.9, -33.6, -24.1/); Jan_clear_surflw = (/119, 103, 106, 109/); Jan_clear_surftemp = (/-26.1, -27.5, -33.8, -40.3/); Feb_blsn_depth = (/155, 119, 180, 143, 99, 75, 93/); Feb_blsn_pixels = (/117, 26, 1, 227, 31, 8, 10/); Feb_blsn_domec_depth = (/-999, -999, -999, -999, -999, -999, -999/); Feb_blsn_surflw = (/104, 99, 110, 85, 83, 94, 91/); Feb_blsn_surftemp = (/-35, -37, -36.9, -43.6, -43.5, -53, -44.9/); Feb_cloud_surflw = (/123, 127/); Feb_cloud_surftemp = (/-38.7 ,-47.4/); Feb_clear_surflw = (/103, 93, 101, 104, 88, 82, 80, 90, 96/); Feb_clear_surftemp = (/-35.5, -33.8, -28.9, -40.4, -51.8, -53.8, -42.6, -50.4, -40.6/); Mar_blsn_depth = (/103, 109, 85, 81, 60, 105/) Mar_blsn_pixels = (/211, 394, 51, 48, 1, 2/) Mar_blsn_domec_depth = (/-999, 90, -999, -999, -999, -999/) Mar_blsn_surflw = (/81, 78, 75, 79, 86, 75/) Mar_blsn_surftemp = (/-53.1, -55.6, -53.4, -56, -55.5, -59.8/) Mar_cloud_surflw = (/119, 109, 111, 92, 73, 122, 105/) Mar_cloud_surftemp = (/-54.1, -48.9, -41.8, -48.3, -47.7, -49.5, -999/) Mar_clear_surflw = (/86, 81, 69, 63, 91/) Mar_clear_surftemp = (/-999, -55.2, -65.3, -68.5, -53.6/) Apr_blsn_depth = (/98, 75, 73/) Apr_blsn_pixels = (/4, 2, 66/) Apr_blsn_domec_depth = (/-999, -999, -999/) Apr_blsn_surflw = (/72, 66, 51/) Apr_blsn_surftemp = (/-63.2, -65.5, -76.2/) Apr_cloud_surflw = (/91, 85, 102, 64, 60, 70, -999/) Apr_cloud_surftemp = (/-56, -63.1, -52, -64.7, -72.8, -65.4, -999/) Apr_clear_surflw = (/81, 62, 58, 61, 66, 48, 68/) Apr_clear_surftemp = (/-57.7, -68.9, -999, -999, -68.5, -77.8, -68.2/) May_blsn_depth = (/77, 108, 120, 64, 70, 60, 66, 71/) May_blsn_pixels = (/7, 49, 486, 30, 6, 1, 5, 32/) May_blsn_domec_depth = (/-999, -999, 150, -999, -999, -999, -999, -999/) May_blsn_surflw = (/76, 66, 67, 60, 60, 56, 73, -999/) May_blsn_surftemp = (/-67.4, -70, -61.6, -71.9, -73.3, -76.3, -73.4, -70.4/) print(May_blsn_surflw@_FillValue) print(typeof(May_blsn_surflw)) May_cloud_surflw = (/64, 88, 60, 63, 64/) May_cloud_surftemp = (/-66.2, -50.2, -73.7, -74.9, -70/) May_clear_surflw = (/80, 67, 65, 65, 68, 64, 69/) May_clear_surftemp = (/-67.7, -999, -70.3, -70.2, -65.4, -75.1, -67.4 /) Jun_blsn_depth = (/75, 45, 60, 120/) Jun_blsn_pixels = (/2, 2, 1, 1/) Jun_blsn_domec_depth = (/-999, -999, -999, -999/) Jun_blsn_surflw = (/56, 70, 56, 60/) Jun_blsn_surftemp = (/-72, -60.2, -72.2, -999/) Jun_cloud_surflw = (/72, 85, 81, 78/) Jun_cloud_surftemp = (/-60.1, -999, -61.2, -50.4/) Jun_clear_surflw = (/64, 61, 55, 62, 59, 52, 62, 65/) Jun_clear_surftemp = (/-70.5, -71.6, -72.5, -58.2, -61.9, -76.6, -69.4, -70.8/) Jul_blsn_depth = (/47, 75, 83/) Jul_blsn_pixels = (/16, 24, 4/) Jul_blsn_domec_depth = (/-999, -999, -999/) Jul_blsn_surflw = (/50, 51, 48/) Jul_blsn_surftemp = (/-79.4, -75.7, -76.5/) Jul_cloud_surflw = (/111, 79, 75, 84, 101/) Jul_cloud_surftemp = (/-61, -66.6, -54.6, -999, -60.8/) Jul_clear_surflw = (/52, 50, 51, 68, 46, 46, 50/) Jul_clear_surftemp = (/-77.2, -82, -73, -68.8, -79, -73.3, -77.7/) Aug_cloud_surflw = (/61, 72, 59, 97, 98, 55, 159, 87, 118, 70/) Aug_cloud_surftemp = (/-61.3, -61.8, -67.2, -59.6, -58.3, -999, -48.4, -60.2, -60.2, -63.9/) Aug_clear_surflw = (/43, 45, 46, 46, 66, 66/) Aug_clear_surftemp = (/-82.1, -75.4, -70.8, -77, -66, -69.5/) Sep_blsn_depth = (/82, 80, 87, 75, 76, 76, 76, 90, 60, 63, 64/) Sep_blsn_pixels = (/43, 3, 463, 76, 94, 50, 50, 7, 9, 195, 101/) Sep_blsn_domec_depth = (/-999, -999, 210, 30, -999, -999, -999, 30, -999, 60, -999/) Sep_blsn_surflw = (/62, 60, 62, 61, 77, 65, 61, 68, 68, 58, 79 /) Sep_blsn_surftemp = (/-67.7, -75, -65.4, -62.5, -55.3, -63.8, -68.4, -65.2, -59, -61.6, -57.3/) Sep_cloud_surflw = (/76, 61, 76, 78/) Sep_cloud_surftemp = (/-999, -73.6, -56.2, -55.2/) Sep_clear_surflw = (/48/) Sep_clear_surftemp = (/-72.9/) Oct_blsn_depth = (/94, 38, 72, 68, 34/) Oct_blsn_pixels = (/14, 4, 5, 22, 31/) Oct_blsn_domec_depth = (/-999, -999, -999, 30, 30/) Oct_blsn_surflw = (/71, 67, 66, 78, 86/) Oct_blsn_surftemp = (/-57.9, -61.1, -51.5, -56.4, -52.9/) Oct_cloud_surflw = (/82, -999, -999, 119, -999/) Oct_cloud_surftemp = (/-61, -53.6, -999, -50.1, -52/) Oct_clear_surflw = (/70,66, -999, 71, 80, -999, -999/) Oct_clear_surftemp = (/-57.3, -64.9, -64.5, -53.2, -50.7, -45, -49.6/) Nov_blsn_depth = (/63, 83, 125, 70, 90, 108, 90, 120, 108/) Nov_blsn_pixels = (/19, 12, 12, 3, 1, 5, 1, 7, 5/) Nov_blsn_domec_depth = (/-999, -999, -999, -999, -999, -999, -999, -999, -999/) Nov_blsn_surflw = (/77, 67, 68, 121, 103, 95, 77, 81, 81/) Nov_blsn_surftemp = (/-51.4, -43.7, -44, -36.8, -34, -31.3, -35.6, -33.3, -33.8/) Nov_cloud_surflw = (/97, 82, 42/) Nov_cloud_surftemp = (/-45.9, -35.6, -44.9/) Nov_clear_surflw = (/79, 79, 75, 77/) Nov_clear_surftemp = (/-52.9, -45.5, -48, -38.5/) Dec_blsn_depth = (/114, 120, 101, 108, 75/) Dec_blsn_pixels = (/19, 11, 8, 20, 2/) Dec_blsn_domec_depth = (/-999, 60, -999, -999, -999/) Dec_blsn_surflw = (/78, 109, 105, 109, 118/) Dec_blsn_surftemp = (/-34.7, -24.5, -25.9, -22.8, -21.5/) Dec_cloud_surflw = (/76, 143, 121/) Dec_cloud_surftemp = (/-46.8, -31.1, -32.1/) Dec_clear_surflw = (/70, 74, 119, 119, 117, 120, 120, 110, 104/) Dec_clear_surftemp = (/-46.8, -34.7, -30.3, -21.1, -31, -24, -30.2, -31.2, -26.2/) Aug_blsn_surflw = (/-999, -999, -999/) Aug_blsn_surftemp = (/-999, -999, -999/) Jan_blsn_domec_depth@_FillValue= -999 Jan_blsn_surflw@_FillValue = -999 Jan_blsn_surftemp@FillValue = -999 Feb_blsn_surflw@_FillValue = -999 Feb_blsn_surftemp@FillValue = -999 Mar_blsn_surflw@FillValue = -999 Mar_blsn_surftemp@FillValue = -999 Apr_blsn_surflw@FillValue = -999 Apr_blsn_surftemp@FillValue = -999 May_blsn_surflw@_FillValue= -999 May_blsn_surftemp@FillValue = -999 Jun_blsn_surflw@FillValue = -999 Jun_blsn_surftemp@_FillValue = -999 Jul_blsn_surflw@FillValue = -999 Jul_blsn_surftemp@_FillValue = -999 Aug_blsn_surflw@_FillValue = -999 Aug_blsn_surftemp@_FillValue = -999 Sep_blsn_surflw@_FillValue = -999 Sep_blsn_surftemp@_FillValue = -999 Oct_blsn_surflw@_FillValue = -999 Oct_blsn_surftemp@_FillValue = -999 Nov_blsn_surflw@_FillValue = -999 Nov_blsn_surftemp@_FillValue = -999 Dec_blsn_surflw@_FillValue = -999 Dec_blsn_surftemp@_FillValue = -999 Jan_cloud_surflw@_FillValue = -999 Jan_cloud_surftemp@FillValue = -999 Feb_cloud_surflw@_FillValue = -999 Feb_cloud_surftemp@FillValue = -999 Mar_cloud_surflw@FillValue = -999 Mar_cloud_surftemp@FillValue = -999 Apr_cloud_surflw@FillValue = -999 Apr_cloud_surftemp@FillValue = -999 May_cloud_surflw@_FillValue= -999 May_cloud_surftemp@FillValue = -999 Jun_cloud_surflw@FillValue = -999 Jun_cloud_surftemp@_FillValue = -999 Jul_cloud_surflw@FillValue = -999 Jul_cloud_surftemp@_FillValue = -999 Aug_cloud_surflw@_FillValue = -999 Aug_cloud_surftemp@_FillValue = -999 Sep_cloud_surflw@_FillValue = -999 Sep_cloud_surftemp@_FillValue = -999 Oct_cloud_surflw@_FillValue = -999 Oct_cloud_surftemp@_FillValue = -999 Nov_cloud_surflw@_FillValue = -999 Nov_cloud_surftemp@_FillValue = -999 Dec_cloud_surflw@_FillValue = -999 Dec_cloud_surftemp@_FillValue = -999 Jan_clear_surflw@_FillValue = -999 Jan_clear_surftemp@FillValue = -999 Feb_clear_surflw@_FillValue = -999 Feb_clear_surftemp@FillValue = -999 Mar_clear_surflw@FillValue = -999 Mar_clear_surftemp@FillValue = -999 Apr_clear_surflw@FillValue = -999 Apr_clear_surftemp@FillValue = -999 May_clear_surflw@_FillValue= -999 May_clear_surftemp@FillValue = -999 Jun_clear_surflw@FillValue = -999 Jun_clear_surftemp@_FillValue = -999 Jul_clear_surflw@FillValue = -999 Jul_clear_surftemp@_FillValue = -999 Aug_clear_surflw@_FillValue = -999 Aug_clear_surftemp@_FillValue = -999 Sep_clear_surflw@_FillValue = -999 Sep_clear_surftemp@_FillValue = -999 Oct_clear_surflw@_FillValue = -999 Oct_clear_surftemp@_FillValue = -999 Nov_clear_surflw@_FillValue = -999 Nov_clear_surftemp@_FillValue = -999 Dec_clear_surflw@_FillValue = -999 Dec_clear_surftemp@_FillValue = -999 blsn_surflw = (/dim_avg(Jan_blsn_surflw), dim_avg(Feb_blsn_surflw), dim_avg(Mar_blsn_surflw), dim_avg(Apr_blsn_surflw), dim_avg(May_blsn_surflw), dim_avg(Jun_blsn_surflw), dim_avg(Jul_blsn_surflw), dim_avg(Aug_blsn_surflw), dim_avg(Sep_blsn_surflw), dim_avg(Oct_blsn_surflw), dim_avg(Nov_blsn_surflw), dim_avg(Dec_blsn_surflw)/) blsn_surftemp = (/dim_avg(Jan_blsn_surftemp),dim_avg(Feb_blsn_surftemp),dim_avg(Mar_blsn_surftemp), dim_avg(Apr_blsn_surftemp), dim_avg(May_blsn_surftemp), dim_avg(Jun_blsn_surftemp), dim_avg(Jul_blsn_surftemp), dim_avg(Aug_blsn_surftemp), dim_avg(Sep_blsn_surftemp), dim_avg(Oct_blsn_surftemp), dim_avg(Nov_blsn_surftemp), dim_avg(Dec_blsn_surftemp)/) print(blsn_surflw) print(blsn_surftemp) cloud_surflw = (/dim_avg(Jan_cloud_surflw), dim_avg(Feb_cloud_surflw), dim_avg(Mar_cloud_surflw), dim_avg(Apr_cloud_surflw), dim_avg(May_cloud_surflw), dim_avg(Jun_cloud_surflw), dim_avg(Jul_cloud_surflw), dim_avg(Aug_cloud_surflw), dim_avg(Sep_cloud_surflw), dim_avg(Oct_cloud_surflw), dim_avg(Nov_cloud_surflw), dim_avg(Dec_cloud_surflw)/) cloud_surftemp = (/dim_avg(Jan_cloud_surftemp),dim_avg(Feb_cloud_surftemp),dim_avg(Mar_cloud_surftemp), dim_avg(Apr_cloud_surftemp), dim_avg(May_cloud_surftemp), dim_avg(Jun_cloud_surftemp), dim_avg(Jul_cloud_surftemp), dim_avg(Aug_cloud_surftemp), dim_avg(Sep_cloud_surftemp), dim_avg(Oct_cloud_surftemp), dim_avg(Nov_cloud_surftemp), dim_avg(Dec_cloud_surftemp)/) print(cloud_surflw) print(cloud_surftemp) clear_surflw = (/dim_avg(Jan_clear_surflw), dim_avg(Feb_clear_surflw), dim_avg(Mar_clear_surflw), dim_avg(Apr_clear_surflw), dim_avg(May_clear_surflw), dim_avg(Jun_clear_surflw), dim_avg(Jul_clear_surflw), dim_avg(Aug_clear_surflw), dim_avg(Sep_clear_surflw), dim_avg(Oct_clear_surflw), dim_avg(Nov_clear_surflw), dim_avg(Dec_clear_surflw)/) clear_surftemp = (/dim_avg(Jan_clear_surftemp),dim_avg(Feb_clear_surftemp),dim_avg(Mar_clear_surftemp), dim_avg(Apr_clear_surftemp), dim_avg(May_clear_surftemp), dim_avg(Jun_clear_surftemp), dim_avg(Jul_clear_surftemp), dim_avg(Aug_clear_surftemp), dim_avg(Sep_clear_surftemp), dim_avg(Oct_clear_surftemp), dim_avg(Nov_clear_surftemp), dim_avg(Dec_clear_surftemp)/) print(clear_surflw) print(clear_surftemp) month = (/1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12/) ;************************************************ ; plotting parameters ;************************************************ wks = gsn_open_wks("x11","Plot_Monthly") ; send graphics to PNG file res = True ; plot mods desired res@tiMainString = "Monthy plot" ; add title res@xyMarkLineMode = "Markers" ; choose to use markers res@xyMarkers = 16 ; choose type of marker res@xyMarkerColor = "gold2" ; Marker color res@xyMarkerSizeF = 0.015 ; Marker size (default 0.01) res@trYMinF = -80 ; min value on y-axis res@trYMaxF = 120 ; max value on y-axis ; res@trXMinF = 40.0 ; min value on x-axis ; res@trXMaxF = 150.0 res@gsnDraw = False res@gsnFrame = False ; plot = gsn_csm_xy(wks,Nov_blsn_surflw,Nov_blsn_surftemp,res) ; create plot plot = gsn_csm_xy(wks,month,blsn_surflw,res) ; create plot ; res@xyMarkerColor = "green" ; plot2 = gsn_csm_xy(wks,Nov_cloud_surflw,Nov_cloud_surftemp,res) ; create plot res@xyMarkers = 3 res@xyMarkerThicknessF = 4 plot2 = gsn_csm_xy(wks,month,blsn_surftemp,res) ; create plot ;res@xyMarkerColor = "blue" ;plot3 = gsn_csm_xy(wks,Nov_clear_surflw,Nov_clear_surftemp,res) ; create plot res@xyMarkers = 16 res@xyMarkerThicknessF = 1 res@xyMarkerColor = "green" plot3 = gsn_csm_xy(wks,month,cloud_surflw,res) ; create plot res@xyMarkers = 3 res@xyMarkerThicknessF = 4 plot4 = gsn_csm_xy(wks,month,cloud_surftemp,res) ; create plot res@xyMarkers = 16 res@xyMarkerThicknessF = 1 res@xyMarkerColor = "blue" plot5 = gsn_csm_xy(wks,month,clear_surflw,res) ; create plot res@xyMarkers = 3 res@xyMarkerThicknessF = 4 plot6 = gsn_csm_xy(wks,month,clear_surftemp,res) ; create plot overlay(plot,plot2) ; overlay the U-wind plot on the temperature plot overlay(plot,plot3) overlay(plot,plot4) overlay(plot,plot5) overlay(plot,plot6) draw(plot) ; draw the temperature plot (with the U-wind plot overlaid) frame(wks) end