; load libraries, twice for good measure 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/wrf/WRFUserARW.ncl" ;load "/home/torn/ncl/wrf_map_info.ncl" ;load "/home/torn/ncl/array_stat.ncl" ;load "/home/torn/ncl/image_conv.ncl" 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/wrf/WRFUserARW.ncl" ;load "/home/torn/ncl/wrf_map_info.ncl" ;load "/home/torn/ncl/array_stat.ncl" ;load "/home/torn/ncl/image_conv.ncl" ;load "wrf_dart_util.ncl" begin wks_type = "png" wks_type@wkWidth = 2048 wks_type@wkHeight = 2048 wks = gsn_open_wks(wks_type,"header") gsn_define_colormap(wks,"WhiteBlueGreenYellowRed") gres = True tres = True tres@txFont = "helvetica-bold" mres = True mres@gsMarkerIndex = 16 ; marker style (filled circle) ; add temperature stuff gres@gsFillColor = 190 plat = (/.90, .96, .96, .9, .9/) plon = (/.05, .05, .15, .15, .05/) gsn_polygon_ndc(wks,plon,plat,gres) mres@gsMarkerSizeF = 61. ; marker size mres@gsMarkerColor = 190 ; marker color gsn_polymarker_ndc(wks,.15,.93,mres) gsn_polymarker_ndc(wks,.05,.93,mres) tres@txFontHeightF = 0.0105 gsn_text_ndc(wks,"Temperature",-79.55,45.19,tres) gsn_text_ndc(wks,"Anomaly",-79.55,45.06,tres) gsn_text_ndc(wks,"________",-79.55,45.05,tres) mres@gsMarkerSizeF = 50.0 ; marker size mres@gsMarkerColor = 200 ; marker color gsn_polymarker_ndc(wks,-79.55,44.75,mres) gsn_text_ndc(wks,"1",-79.55,44.85,tres) gsn_text_ndc(wks,"Day",-79.55,44.71,tres) mres@gsMarkerColor = 210 ; marker color gsn_polymarker_ndc(wks,-79.55,44.3,mres) gsn_text_ndc(wks,"1",-79.55,44.40,tres) gsn_text_ndc(wks,"Week",-79.55,44.26,tres) mres@gsMarkerColor = 220 ; marker color gsn_polymarker_ndc(wks,-79.55,43.85,mres) gsn_text_ndc(wks,"1",-79.55,43.95,tres) gsn_text_ndc(wks,"Month",-79.55,43.81,tres) ;;;; Precipitation Section gres@gsFillColor = 130 plat = (/45, 45.3, 45.3, 45, 45/) plon = (/-78.9, -78.9, -78.1, -78.1, -78.9/) gsn_polygon_ndc(wks,plon,plat,gres) mres@gsMarkerSizeF = 36.5 ; marker size mres@gsMarkerColor = 130 ; marker color gsn_polymarker_ndc(wks,-78.9,45.15,mres) gsn_polymarker_ndc(wks,-78.1,45.15,mres) gsn_text_ndc(wks,"Precipitation",-78.5,45.19,tres) gsn_text_ndc(wks," Total Anomaly",-78.49,45.06,tres) gsn_text_ndc(wks,"____",-78.76,45.05,tres) gsn_text_ndc(wks,"________",-78.3,45.05,tres) mres@gsMarkerSizeF = 50.0 ; marker size mres@gsMarkerColor = 120 ; marker color gsn_polymarker_ndc(wks,-78.75,44.75,mres) gsn_text_ndc(wks,"1",-78.75,44.85,tres) gsn_text_ndc(wks,"Day",-78.75,44.71,tres) mres@gsMarkerColor = 110 ; marker color gsn_polymarker_ndc(wks,-78.75,44.3,mres) gsn_text_ndc(wks,"1",-78.75,44.40,tres) gsn_text_ndc(wks,"Week",-78.75,44.26,tres) mres@gsMarkerColor = 100 ; marker color gsn_polymarker_ndc(wks,-78.75,43.85,mres) gsn_text_ndc(wks,"1",-78.75,43.95,tres) gsn_text_ndc(wks,"Month",-78.75,43.81,tres) mres@gsMarkerColor = 50 ; marker color gsn_polymarker_ndc(wks,-78.25,44.75,mres) gsn_text_ndc(wks,"N/A",-78.25,44.75,tres) mres@gsMarkerColor = 60 ; marker color gsn_polymarker_ndc(wks,-78.25,44.3,mres) gsn_text_ndc(wks,"1",-78.25,44.40,tres) gsn_text_ndc(wks,"Week",-78.25,44.26,tres) mres@gsMarkerColor = 70 ; marker color gsn_polymarker_ndc(wks,-78.25,43.85,mres) gsn_text_ndc(wks,"1",-78.25,43.95,tres) gsn_text_ndc(wks,"Month",-78.25,43.81,tres) ; add solar section gres@gsFillColor = 150 plat = (/45, 45.3, 45.3, 45, 45/) plon = (/-77.7, -77.7, -77.2, -77.2, -77.7/) gsn_polygon_ndc(wks,plon,plat,gres) mres@gsMarkerSizeF = 36.5 ; marker size mres@gsMarkerColor = 150 ; marker color gsn_polymarker_ndc(wks,-77.7,45.15,mres) gsn_polymarker_ndc(wks,-77.2,45.15,mres) gsn_text_ndc(wks,"Insolation",-77.45,45.19,tres) gsn_text_ndc(wks,"Total",-77.45,45.068,tres) gsn_text_ndc(wks,"_____",-77.45,45.05,tres) mres@gsMarkerSizeF = 50.0 ; marker size mres@gsMarkerColor = 160 ; marker color gsn_polymarker_ndc(wks,-77.45,44.75,mres) gsn_text_ndc(wks,"1",-77.45,44.85,tres) gsn_text_ndc(wks,"Day",-77.45,44.71,tres) mres@gsMarkerColor = 170 ; marker color gsn_polymarker_ndc(wks,-77.45,44.3,mres) gsn_text_ndc(wks,"1",-77.45,44.40,tres) gsn_text_ndc(wks,"Week",-77.45,44.26,tres) mres@gsMarkerColor = 180 ; marker color gsn_polymarker_ndc(wks,-77.45,43.85,mres) gsn_text_ndc(wks,"1",-77.45,43.95,tres) gsn_text_ndc(wks,"Month",-77.45,43.81,tres) frame(wks) end