; Example script to produce plots for a WRF real-data run, ; with the ARW coordinate dynamics option. ; Plot data on a cross section ; This script will plot data at a set angle through a specified point ; This script adds lon/lat info along X-axis 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" begin ; ----------------------- ; Type of plot to generate ; wks_type = "x11" ; wks_type = "ncgm" ; wks_type = "ps" ; wks_type@wkOrientation = "portrait" ; wks_type@wkPaperWidthF = 8.5 ; in inches ; wks_type@wkPaperHeightF = 11.0 ; in inches wks_type = "pdf" wks_type@wkPaperSize = "A4" ; wks_type = "png" ; wks_type@wkWidth = 1000 ; wks_type@wkHeight = 800 ; ----------------------- wks = gsn_open_wks(wks_type,"wrf.hrz.cross") ; ----------------------- ; gsn_define_colormap (wks,"gui_default") ; gsn_define_colormap (wks,"wgne15") ; gsn_define_colormap (wks,"ViReOrYeGrAqBl") ; gsn_define_colormap (wks,"WhViBlGrYeOrReWh") gsn_define_colormap (wks,"BlAqGrYeOrReVi200") ; gsn_reverse_colormap(wks) ; Reverse the color map. ; ----------------------- ; Set some basic resources res = True ; res@MainTitle = "REAL-TIME WRF" res = True res@NoTitles = False res@Footer = False res@Header = False res@NoHeaderFooter = True mpres = True mpres@NoTitles = False mpres@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries mpres@mpDataSetName = "Earth..4" mpres@mpFillBoundarySets = "National" mpres@mpGeophysicalLineColor = "Black" mpres@mpNationalLineColor = "Black" mpres@mpUSStateLineColor = "Black" ; mpres@mpGeophysicalLineColor = "Maroon" ; mpres@mpNationalLineColor = "Maroon" ; mpres@mpUSStateLineColor = "Maroon" ; mpres@mpGeophysicalLineThicknessF = 2 ; mpres@mpNationalLineThicknessF = 2 ; mpres@mpUSStateLineThicknessF = 2 ; mpres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; mpres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; mpres@tmXBMinorOn = True ; no lon minor tickmarks mpres@Footer = False mpres@Header = False mpres@NoHeaderFooter = True res@gsnMaximize = True ; maximize plot in frame pltres = True pltres@NoTitles = False ; ----------------------- FirstTime = True FirstTimeMap = True ; ------------------------------------------ ; First set of input files. wrffiles1 = systemfunc("ls input.files/wrfout_d03*") numFiles1 = dimsizes(wrffiles1) do i = 0, numFiles1 -1 wrffiles1(i) = wrffiles1(i)+".nc" end do inpFiles1 = addfiles(wrffiles1,"r") ; ------------------------------------------ ; Second set of input files. ; wrffiles2 = systemfunc("ls run2/wrfout_d03*") ; numFiles2 = dimsizes(wrffiles2) ; do i = 0, numFiles2 -1 ; wrffiles2(i) = wrffiles2(i)+".nc" ; end do ; inpFiles2 = addfiles(wrffiles2,"r") ; ------------------------------------------ ; LOOP OVER FILES do ifile = 0, numFiles1-1,3 a = inpFiles1[ifile] ; b = inpFiles2[ifile] times1 = wrf_user_list_times(a) ; get times in the file ntimes1 = dimsizes(times1) ; number of times in the file ; if(ifile.eq.0) then ; times1 = wrf_user_list_times(a) ; get times in the file ; ntimes1 = dimsizes(times1) ; number of times in the file ; istrt = 3 ; Model run started at 12PM. Start processing at 2 PM. ; iend = 133 ; iend = 5 ; end if ; if(ifile.eq.1) then ; times2 = wrf_user_list_times(a) ; get times in the file ; ntimes2 = dimsizes(times2) ; number of times in the file ; istrt = 13 ; iend = 133 ; iend = 17 ; end if irpt = 3 ; every third datafile ; time = a->Times(:,0) ; times = chartostring(a->Times) ; ntime = dimsizes(time) if(ifile.eq.0) then u10_1 = wrf_user_getvar(a,"U10",-1) ; u at 10 m, mass point v10_1 = wrf_user_getvar(a,"V10",-1) ; v at 10 m, mass point u10_1@units = "" u10_1@description = "" v10_1@units = "" v10_1@description = "" t2_1 = wrf_user_getvar(a,"T2",-1) ; theta at 2m only (K) t2_1@units = "" t2_1@description = "" t2_1@description = "" tc1 = wrf_user_getvar(a,"tc",-1) ; 3D temperature (C) tc1@units = "" tc1@description = "" qv1 = wrf_user_getvar(a,"QVAPOR",-1)*1000. ; mixing ratio (g/kg) qv1@units = "" qv1@description = "" rh1 = wrf_user_getvar(a,"rh",-1) ; relative humidity (%) rh1@units = "" rh1@description = "" xlat = wrf_user_getvar(a,"XLAT",0) xlon = wrf_user_getvar(a,"XLONG",0) xlatmin = min(xlat) xlatmax = max(xlat) xlonmin = min(xlon) xlonmax = max(xlon) ; print("xlatmin = "+xlatmin) ; print("xlatmax = "+xlatmax) ; print("xlonmin = "+xlonmin) ; print("xlonmax = "+xlonmax) ; print("") ; printMinMax(xlat,0) ; printMinMax(xlon,0) else u10_2 = wrf_user_getvar(a,"U10",-1) ; u at 10 m, mass point v10_2 = wrf_user_getvar(a,"V10",-1) ; v at 10 m, mass point u10_2@units = "" u10_2@description = "" v10_2@units = "" v10_2@description = "" t2_2 = wrf_user_getvar(a,"T2",-1) ; theta at 2m only (K) t2_2@units = "" t2_2@description = "" tc2 = wrf_user_getvar(a,"tc",-1) ; 3D temperature (C) tc2@units = "" tc2@description = "" qv2 = wrf_user_getvar(a,"QVAPOR",-1)*1000. ; mixing ratio (g/kg) qv2@units = "" qv2@description = "" rh2 = wrf_user_getvar(a,"rh",-1) ; relative humidity (%) rh2@units = "" rh2@description = "" xlat = wrf_user_getvar(a,"XLAT",0) xlon = wrf_user_getvar(a,"XLONG",0) xlatmin = min(xlat) xlatmax = max(xlat) xlonmin = min(xlon) xlonmax = max(xlon) ; print("xlatmin = "+xlatmin) ; print("xlatmax = "+xlatmax) ; print("xlonmin = "+xlonmin) ; print("xlonmax = "+xlonmax) ; print("") ; printMinMax(xlat,0) ; printMinMax(xlon,0) end if ; printVarSummary(u10_1) ; umin = min(u10_1) ; umax = max(u10_1) ; print(umin) ; print(umax) ; printVarSummary(t2) ; tmin = min(t2) ; tmax = max(t2) ; print(tmin) ; print(tmax) ; u10 = u10*1.94386 ; Convert m/s into knots ; v10 = v10*1.94386 ; u10@units = "kts" ; v10@units = "kts" ; tc1 = wrf_user_getvar(a,"tc",-1) ; temperature ; tc2 = wrf_user_getvar(b,"tc",-1) ; temperature ; rh1 = wrf_user_getvar(a,"rh",-1) ; relative humidity ; rh2 = wrf_user_getvar(b,"rh",-1) ; relative humidity ; do it = 0,ntimes1-1,3 ; TIME LOOP ; do it = 0,ntimes1-1,irpt ; TIME LOOP do it = 1,3,2 ; TIME LOOP ; do it = 0,0 ; TIME LOOP ; if(ifile.eq.0) then print("Working on time = "+times1(it)) ; else ; print("Working on time = "+times2(it)) ; end if ;******************************************** ; fancy curly vector plot ;******************************************** if(ifile.eq.0) then u10 = u10_1(it,:,:) v10 = v10_1(it,:,:) else u10 = u10_2(it,:,:) v10 = v10_2(it,:,:) end if spd = sqrt(u10^2+v10^2) gsn_define_colormap (wks,"BlAqGrYeOrReVi200") ; gsn_reverse_colormap(wks) ; Reverse the color map. vecres1 = True ; plot mods desired vecres1@tfDoNDCOverlay = True vecres1@gsnMaximize = True ; maximize plot in frame vecres1@mpShapeMode = "FreeAspect" vecres1@vpWidthF = 0.80 vecres1@vpHeightF = 0.75 vecres1@mpLeftCornerLatF = xlatmin vecres1@mpRightCornerLatF = xlatmax vecres1@mpLeftCornerLonF = xlonmin vecres1@mpRightCornerLonF = xlonmax vecres1@mpMinLatF = xlatmin vecres1@mpMaxLatF = xlatmax vecres1@mpMinLonF = xlonmin vecres1@mpMaxLonF = xlonmax vecres1@mpOutlineDrawOrder = "PostDraw" ; draw continental outline last ; vecres1@mpOutLineOn = True vecres1@pmTickMarkDisplayMode = "Always" vecres1@mpFillOn = False ; turn off map fill vecres1@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries vecres1@mpDataSetName = "Earth..4" vecres1@mpFillBoundarySets = "National" vecres1@mpGridLineColor = "Black" vecres1@mpPerimLineColor = "Black" vecres1@mpGeophysicalLineColor = "Black" vecres1@mpNationalLineColor = "Black" vecres1@mpUSStateLineColor = "Black" vecres1@mpLimbLineColor = "Black" vecres1@mpPerimLineColor = "Black" vecres1@gsnAddCyclic = False ; regional data ; ---- Just arrows for wind vectors ------------------------------------------------------- ; vecres1@vcGlyphStyle = "LineArrow" ; vecres1@vcLineArrowThicknessF = 1 ; vecres1@vcMinDistanceF = 0.015 ; vecres1@vcRefLengthF = 0.035 ; vecres1@vcRefAnnoOn = True ; vecres1@vcRefMagnitudeF = 10 ; vecres1@vcRefAnnoString1 = "10" ; vecres1@vcRefAnnoSide = "Top" ; vecres1@vcRefAnnoString2On = False ; vecres1@vcRefAnnoPerimOn = True ; vecres1@vcRefAnnoOrthogonalPosF = -0.12 ; move ref vector up ; vecres1@vcRefAnnoParallelPosF = 0.999 ; vecres1@vcRefAnnoBackgroundColor = "Yellow" ; vecres1@vcVectorDrawOrder = "PostDraw" ; ---- Curly vectors ---------------------------------------------------------------------- vecres1@vcGlyphStyle = "CurlyVector" ; turn on curly vectors vecres1@vcRefMagnitudeF = 1. ; define vector ref mag vecres1@vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector down vecres1@vcRefAnnoArrowLineColor = "red" ; change ref vector color vecres1@vcRefLengthF = 0.045 ; define length of vec ref vecres1@vcMinDistanceF = 0.017 ; thin vectors vecres1@vcMinDistanceF = 0.02 ; larger means sparser vecres1@vcLineArrowHeadMaxSizeF = 0.0075 ; vecres1@vcLineArrowHeadMaxSizeF = 0.01 ; ; ---------------------------------------------------------------------------------------- vecres1@tiMainOn = True vecres1@tiMainString = "WRF - Wind Vectors (m/s) - "+times1(it) vecres1@tiMainStringOrthogonalPosF = 0.02 vecres1@tiMainFont = "helvetica" vecres1@tiMainFontHeightF = 0.016 vecres1@tiMainFontColor = "black" vecres1@tmXBLabelFontHeightF = 0.016 vecres1@tmYLLabelFontHeightF = 0.016 vecres1@mpGeophysicalLineThicknessF = 2 vecres1@mpNationalLineThicknessF = 2 vecres1@mpUSStateLineThicknessF = 2 vecres1@vcLevelSelectionMode = "ManualLevels" vecres1@vcMinLevelValF = 0. vecres1@vcMaxLevelValF = 20. vecres1@vcLevelSpacingF = 1. ; vecres1@gsnSpreadColorEnd = -2 vecres1@lbLabelAutoStride = True vecres1@lbAutoManage = False ; we control label bar vecres1@pmLabelBarDisplayMode = "Always" ; turns on label bar vecres1@lbOrientation = "Horizontal" ; ncl default is vertical vecres1@pmLabelBarSide = "Bottom" ; default is right vecres1@lbLabelStride = 2 ; skip every other label vecres1@pmLabelBarWidthF = 0.80 ; default is shorter vecres1@pmLabelBarHeightF = 0.10 ; default is taller vecres1@lbLabelFontHeightF = 0.016 ; default is HUGE vecres1@lbTitleFontHeightF = 0.016 vecres1@lbPerimOn = False ; default has box plot = gsn_csm_vector_scalar_map(wks,u10,v10,spd,vecres1) delete(vecres1) ;******************************************** ; Plotting options for Potential Temperature ;******************************************** if(ifile.eq.0) then tk = t2_1(it,:,:) else tk = t2_2(it,:,:) end if gsn_define_colormap (wks,"BlAqGrYeOrReVi200") ; gsn_reverse_colormap(wks) ; Reverse the color map. opts_tk = res opts_tk@gsnMaximize = True ; maximize plot in frame opts_tk@Footer = False opts_tk@Header = False opts_tk@NoHeaderFooter = True ; opts_tk@vpXF = 0.1 ; opts_tk@vpYF = 0.95 opts_tk@mpShapeMode = "FreeAspect" opts_tk@vpWidthF = 0.80 opts_tk@vpHeightF = 0.75 ; mpres@vpXF = 0.1 ; mpres@vpYF = 0.97 ; mpres@mpShapeMode = "FreeAspect" ; mpres@vpWidthF = 0.85 ; mpres@vpHeightF = 0.68 opts_tk@mpLeftCornerLatF = xlatmin opts_tk@mpRightCornerLatF = xlatmax opts_tk@mpLeftCornerLonF = xlonmin opts_tk@mpRightCornerLonF = xlonmax ; opts_tk@mpMinLatF = xlatmin ; opts_tk@mpMaxLatF = xlatmax ; opts_tk@mpMinLonF = xlonmin ; opts_tk@mpMaxLonF = xlonmax ; opts_tk@mpLeftCornerLatF = 31. ; opts_tk@mpRightCornerLatF = 43. ; opts_tk@mpLeftCornerLonF = -128.5 ; opts_tk@mpRightCornerLonF = -113.5 ; opts_tk@mpMinLatF = 31. ; opts_tk@mpMaxLatF = 43. ; opts_tk@mpMinLonF = -128.5 ; opts_tk@mpMaxLonF = -113.5 ; opts_tk@tiMainOn = True ; opts_tk@tiMainString = "Potential Temperature (K) - "+times1(it) ; opts_tk@tiMainFont = "helvetica" ; opts_tk@tiMainFontHeightF = 0.012 ; opts_tk@tiMainFontColor = "black" ; opts_tk@FieldTitle = "Potential Temperature (K) - "+times1(it) ; opts_tk@FieldTitleFontFont = "helvetica" ; opts_tk@FieldTitleFontHeightF = 0.012 ; opts_tk@FieldTitleFontColor = "black" opts_tk@tmXBLabelFontHeightF = 0.013 opts_tk@tmYLLabelFontHeightF = 0.013 mpres@tmXBLabelFontHeightF = 0.013 mpres@tmYLLabelFontHeightF = 0.013 mpres@tiMainOn = True mpres@tiMainString = "WRF - Potential Temperature (K) - "+times1(it) mpres@tiMainStringOrthogonalPosF = 0.02 mpres@tiMainFont = "helvetica" mpres@tiMainFontHeightF = 0.013 mpres@tiMainFontColor = "black" ; mpres@gsnLeftString = "WRF - Potential Temperature (K)" ; mpres@gsnRightString = times1(it) ; mpres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; mpres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; mpres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; pltres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; pltres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@tiMainOn = True ; pltres@tiMainFont = "helvetica" ; pltres@tiMainFontHeightF = 0.012 ; pltres@tiMainFontColor = "black" ; pltres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; pltres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; pltres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@tiMainString = "Potential Temperature (K) - "+times1(it); overwrite Field Title ; opts_tk@gsnLeftString = "Potential Temperature (K)" ; opts_tk@gsnLeftStringFontHeightF = 0.012 ; opts_tk@gsnCenterString = "Potential Temperature (K)"+times1(it) ; opts_tk@gsnCenterStringFontHeightF = 0.013 ; opts_tk@gsnCenterStringOrthogonalPosF = 0.013 ; opts_tk@gsnRightString = times1(it) ; opts_tk@gsnRightStringFontHeightF = 0.012 ; if(ifile.eq.0) then ; opts_tk@tiMainString = "Potential Temperature (K) - "+times1(it); overwrite Field Title ; opts_tk@gsnLeftString = "Potential Temperature (K) - "+times1(it) ; opts_tk@FieldTitle = "Potential Temperature (K) "+times1(it) ; else ; opts_tk@tiMainString = "Potential Temperature (K) - "+times2(it); overwrite Field Title ; opts_tk@gsnLeftString = "Potential Temperature (K) - "+times2(it) ; opts_tk@FieldTitle = "Potential Temperature (K) "+times2(it) ; end if ; cmap = read_colormap_file("BlAqGrYeOrReVi200") ; opts_tk@vcLevelPalette = cmap(6:193,:) ; opts_tk@pmTickMarkDisplayMode = "Always" ; opts_tk@mpFillOn = False ; turn off map fill opts_tk@mpOutlineDrawOrder = "PostDraw" ; draw continental outline last ; opts_tk@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_tk@mpDataSetName = "Earth..4" ; opts_tk@mpGeophysicalLineColor = "Red" ; opts_tk@mpGeophysicalLineThicknessF = 2 ; opts_tk@mpNationalLineThicknessF = 2 ; opts_tk@NationalLineColor = "Red" ; opts_tk@mpOutLineOn = True ; opts_tk@mpGridSpacingF = 2.5 ; opts_tk@mpGridLineThicknessF = 2.0 ; opts_tk@mpGridLineColor = "Red" ; opts_tk@mpPerimLineColor = "Red" ; opts_tk@mpOutLineOn = True ; opts_tk@pmTickMarkDisplayMode = "Always" ; opts_tk@mpFillOn = False ; turn off map fill ; opts_tk@mpFillOn = True ; turn on map fill ; opts_tk@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_tk@mpDataSetName = "Earth..4" ; opts_tk@mpFillBoundarySets = "National" ; opts_tk@mpGridLineColor = "Blue" ; opts_tk@mpPerimLineColor = "Blue" ; opts_tk@mpLimbLineColor = "Blue" ; opts_tk@mpPerimLineColor = "Blue" ; opts_tk@mpGeophysicalLineColor = "Blue" ; opts_tk@mpNationalLineColor = "Blue" ; opts_tk@mpUSStateLineColor = "Blue" ; opts_tk@mpGridLineThicknessF = 2 opts_tk@mpGeophysicalLineThicknessF = 2 opts_tk@mpNationalLineThicknessF = 2 opts_tk@mpUSStateLineThicknessF = 2 ; opts_tk@mpPerimOn = False ; opts_tk@mpGridAndLimbOn = False ; opts_tk@gsnTickMarksOn = False ; opts_tk@mpDataSetName = "Earth..3" ; opts_tk@mpDataBaseVersion = "MediumRes" ; opts_tk@mpOutlineOn = True ; opts_tk@mpOutlineSpecifiers = "Conterminous US: Counties" ; opts_tk@mpOFillAreaSpecifiers = "Conterminous US: Counties" ; opts_tk@mpOutlineBoundarySets = "USStates" ; opts_tk@mpFillOn = True ; opts_tk@mpSpecifiedFillColors = "red" ; opts_tk@mpFillAreaSpecifiers = "California : 05" ; opts_tk@mpOceanFillColor = "transparent" ; opts_tk@mpLandFillColor = "transparent" ; opts_tk@mpInlandWaterFillColor = "transparent" ; opts_tk@mpUSStateLineColor = "red" ; opts_tk@mpNationalLineColor = "red" ; opts_tk@mpGeophysicalLineColor = "red" opts_tk@cnFillOn = True ; opts_tk@cnLevelSelectionMode = "AutomaticLevels" opts_tk@cnLevelSelectionMode = "ManualLevels" opts_tk@ContourParameters = (/ 270., 310., 2./) ; opts_tk@gsnSpreadColorEnd = -2 opts_tk@lbLabelAutoStride = True opts_tk@lbAutoManage = False ; we control label bar opts_tk@pmLabelBarDisplayMode = "Always" ; turns on label bar opts_tk@lbOrientation = "Horizontal" ; ncl default is vertical opts_tk@pmLabelBarSide = "Bottom" ; default is right opts_tk@lbLabelStride = 2 ; skip every other label opts_tk@pmLabelBarWidthF = 0.80 ; default is shorter opts_tk@pmLabelBarHeightF = 0.30 ; default is taller opts_tk@lbLabelFontHeightF = 0.015 ; default is HUGE opts_tk@lbTitleFontHeightF = 0.015 ; opts_tk@lbTitleOn = True ; turn on title ; opts_tk@lbTitleString = "Potential Temperature (K)" ; opts_tk@pmLabelBarOrthogonalPosF = 0.0001 ; move whole thing down opts_tk@lbPerimOn = False ; default has box contour_tk = wrf_contour(a,wks,tk,opts_tk) plot = wrf_map_overlays(a,wks,(/contour_tk/),pltres,mpres) delete(opts_tk) ; map_cont_tk = wrf_map_overlays(a,wks,tk,True,mpres) ;******************************************** ; Plotting options for TC ;******************************************** if(ifile.eq.0) then tc = tc1(it,0,:,:) else tc = tc2(it,0,:,:) end if gsn_define_colormap (wks,"BlAqGrYeOrReVi200") ; gsn_reverse_colormap(wks) ; Reverse the color map. opts_tc = res opts_tc@gsnMaximize = True ; maximize plot in frame opts_tc@Footer = False opts_tc@Header = False opts_tc@NoHeaderFooter = True ; opts_tc@tiMainOn = True ; opts_tc@tiMainString = "Temperature (C) - "+times1(it) ; opts_tc@tiMainFontHeightF = 0.012 ; opts_tc@tiMainFontColor = "black" ; opts_tc@tiMainFont = "helvetica" ; opts_tc@FieldTitle = "Potential Temperature (K) - "+times1(it) ; opts_tc@FieldTitleFontFont = "helvetica" ; opts_tc@FieldTitleFontHeightF = 0.012 ; opts_tc@FieldTitleFontColor = "black" ; opts_tc@vpXF = 0.1 ; opts_tc@vpYF = 0.95 opts_tc@mpShapeMode = "FreeAspect" opts_tc@vpWidthF = 0.8 opts_tc@vpHeightF = 0.75 ; mpres@vpXF = 0.1 ; mpres@vpYF = 0.97 ; mpres@mpShapeMode = "FreeAspect" ; mpres@vpWidthF = 0.85 ; mpres@vpHeightF = 0.68 opts_tc@mpLeftCornerLatF = xlatmin opts_tc@mpRightCornerLatF = xlatmax opts_tc@mpLeftCornerLonF = xlonmin opts_tc@mpRightCornerLonF = xlonmax ; opts_tc@mpMinLatF = xlatmin ; opts_tc@mpMaxLatF = xlatmax ; opts_tc@mpMinLonF = xlonmin ; opts_tc@mpMaxLonF = xlonmax ; opts_tc@mpLeftCornerLatF = 31. ; opts_tc@mpRightCornerLatF = 43. ; opts_tc@mpLeftCornerLonF = -128.5 ; opts_tc@mpRightCornerLonF = -113.5 opts_tc@tmXBLabelFontHeightF = 0.013 opts_tc@tmYLLabelFontHeightF = 0.013 mpres@tmXBLabelFontHeightF = 0.013 mpres@tmYLLabelFontHeightF = 0.013 ; opts_tc@mpMinLatF = 31. ; opts_tc@mpMaxLatF = 43. ; opts_tc@mpMinLonF = -128.5 ; opts_tc@mpMaxLonF = -113.5 mpres@tiMainOn = True mpres@tiMainString = "WRF - Temperature (C) - "+times1(it) mpres@tiMainStringOrthogonalPosF = 0.02 mpres@tiMainFontHeightF = 0.013 mpres@tiMainFontColor = "black" mpres@tiMainFont = "helvetica" ; mpres@gsnLeftString = "WRF - Temperature (C)" ; mpres@gsnRightString = times1(it) ; pltres@tiMainOn = True ; pltres@tiMainFontHeightF = 0.012 ; pltres@tiMainFontColor = "black" ; pltres@tiMainFont = "helvetica" ; mpres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; mpres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; mpres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; pltres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; pltres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@tiMainString = "Temperature (C) - "+times1(it); overwrite Field Title ; opts_tc@gsnLeftString = "Temperature (C)" ; opts_tc@gsnLeftStringFontHeightF = 0.012 ; opts_tc@gsnCenterString = "Temperature (C)"+times1(it) ; opts_tc@gsnCenterStringFontHeightF = 0.013 ; opts_tc@gsnCenterStringOrthogonalPosF = 0.013 ; opts_tc@gsnRightString = times1(it) ; opts_tc@gsnRightStringFontHeightF = 0.012 ; if(ifile.eq.0) then ; opts_tc@tiMainString = "Temperature (C) - "+times1(it); overwrite Field Title ; opts_tc@gsnLeftString = "Temperature (C)~C~Wind Vectors - "+times1(it) ; opts_tc@gsnLeftString = "Temperature (C) - "+times1(it) ; opts_tc@FieldTitle = "Temperature (C) "+times1(it) ; else ; opts_tc@tiMainString = "Temperature (C) - "+times2(it); overwrite Field Title ; opts_tc@gsnLeftString = "Temperature (C)~C~Wind Vectors - "+times2(it) ; opts_tc@gsnLeftString = "Temperature (C) - "+times2(it) ; opts_tc@FieldTitle = "Temperature (C) "+times2(it) ; end if ; cmap = read_colormap_file("BlAqGrYeOrReVi200") ; opts_tc@vcLevelPalette = cmap(6:193,:) ; opts_tc@pmTickMarkDisplayMode = "Always" ; opts_tc@mpFillOn = False ; turn off map fill opts_tc@mpOutlineDrawOrder = "PostDraw" ; draw continental outline last ; opts_tc@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_tc@mpDataSetName = "Earth..4" ; opts_tc@mpGeophysicalLineColor = "Red" ; opts_tc@mpGeophysicalLineThicknessF = 2 ; opts_tc@mpNationalLineThicknessF = 2 ; opts_tc@NationalLineColor = "Red" ; opts_tc@mpOutLineOn = True ; opts_tc@mpGridSpacingF = 2.5 ; opts_tc@mpGridLineThicknessF = 2.0 ; opts_tc@mpGridLineColor = "Red" ; opts_tc@mpPerimLineColor = "Red" ; opts_tc@mpOutLineOn = True ; opts_tc@pmTickMarkDisplayMode = "Always" ; opts_tc@mpFillOn = False ; turn off map fill ; opts_tc@mpFillOn = True ; turn on map fill ; opts_tc@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_tc@mpDataSetName = "Earth..4" ; opts_tc@mpFillBoundarySets = "National" ; opts_tc@mpGridLineColor = "Blue" ; opts_tc@mpPerimLineColor = "Blue" ; opts_tc@mpLimbLineColor = "Blue" ; opts_tc@mpPerimLineColor = "Blue" ; opts_tc@mpGeophysicalLineColor = "Blue" ; opts_tc@mpNationalLineColor = "Blue" ; opts_tc@mpUSStateLineColor = "Blue" ; opts_tc@mpGridLineThicknessF = 2 opts_tc@mpGeophysicalLineThicknessF = 2 opts_tc@mpNationalLineThicknessF = 2 opts_tc@mpUSStateLineThicknessF = 2 ; opts_tc@mpOutLineOn = True ; opts_tc@pmTickMarkDisplayMode = "Always" ; opts_tc@mpFillOn = False ; turn off map fill ; opts_tc@mpFillOn = True ; turn on map fill ; opts_tc@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_tc@mpDataSetName = "Earth..4" ; opts_tc@mpFillBoundarySets = "National" ; opts_tc@mpPerimOn = False ; opts_tc@mpGridAndLimbOn = False ; opts_tc@gsnTickMarksOn = False ; opts_tc@mpDataSetName = "Earth..3" ; opts_tc@mpDataBaseVersion = "MediumRes" ; opts_tc@mpOutlineOn = True ; opts_tc@mpOutlineSpecifiers = "Conterminous US: Counties" ; opts_tc@mpOFillAreaSpecifiers = "Conterminous US: Counties" ; opts_tc@mpOutlineBoundarySets = "USStates" ; opts_tc@mpFillOn = True ; opts_tc@mpSpecifiedFillColors = "red" ; opts_tc@mpFillAreaSpecifiers = "California : 05" ; opts_tc@mpOceanFillColor = "transparent" ; opts_tc@mpLandFillColor = "transparent" ; opts_tc@mpInlandWaterFillColor = "transparent" ; opts_tc@mpUSStateLineColor = "red" ; opts_tc@mpNationalLineColor = "red" ; opts_tc@mpGeophysicalLineColor = "red" ; opts_tc@cnFillOn = True ; turn on color ; opts_tc@cnFillPalette = "gui_default" ; opts_tc@cnFillColors = "Red" ; opts_tc@cnFillMode = "RasterFill" ; opts_tc@cnRasterSmoothingOn = True opts_tc@cnFillOn = True ; opts_tc@cnLevelSelectionMode = "AutomaticLevels" opts_tc@cnLevelSelectionMode = "ManualLevels" opts_tc@ContourParameters = (/ 0., 30., 1./) ; opts_tc@gsnSpreadColorEnd = -2 opts_tc@lbLabelAutoStride = True opts_tc@lbAutoManage = False ; we control label bar opts_tc@pmLabelBarDisplayMode = "Always" ; turns on label bar opts_tc@lbOrientation = "Horizontal" ; ncl default is vertical opts_tc@pmLabelBarSide = "Bottom" ; default is right opts_tc@lbLabelStride = 2 ; skip every other label opts_tc@pmLabelBarWidthF = 0.80 ; default is shorter opts_tc@pmLabelBarHeightF = 0.30 ; default is taller opts_tc@lbLabelFontHeightF = 0.015 ; default is HUGE opts_tc@lbTitleFontHeightF = 0.015 ; opts_tc@lbTitleOn = True ; turn on title ; opts_tc@lbTitleString = "Temperature (K)" ; opts_tc@pmLabelBarOrthogonalPosF = 0.0001 ; move whole thing down opts_tc@lbPerimOn = False ; default has box contour_tc = wrf_contour(a,wks,tc,opts_tc) plot = wrf_map_overlays(a,wks,(/contour_tc/),pltres,mpres) delete(opts_tc) ; map_cont_tc = wrf_map_overlays(a,wks,tc,True,mpres) ;******************************************** ; Plotting options for mixing ratio ;******************************************** if(ifile.eq.0) then qv = qv1(it,0,:,:) else qv = qv2(it,0,:,:) end if gsn_define_colormap (wks,"BlAqGrYeOrReVi200") ; gsn_reverse_colormap(wks) ; Reverse the color map. opts_qv = res opts_qv@gsnMaximize = True ; maximize plot in frame opts_qv@Footer = False opts_qv@Header = False opts_qv@NoHeaderFooter = True ; opts_qv@vpXF = 0.1 ; opts_qv@vpYF = 0.95 opts_qv@mpShapeMode = "FreeAspect" opts_qv@vpWidthF = 0.8 opts_qv@vpHeightF = 0.75 ; mpres@vpXF = 0.1 ; mpres@vpYF = 0.97 ; mpres@mpShapeMode = "FreeAspect" ; mpres@vpWidthF = 0.85 ; mpres@vpHeightF = 0.68 opts_qv@mpLeftCornerLatF = xlatmin opts_qv@mpRightCornerLatF = xlatmax opts_qv@mpLeftCornerLonF = xlonmin opts_qv@mpRightCornerLonF = xlonmax ; opts_qv@mpMinLatF = xlatmin ; opts_qv@mpMaxLatF = xlatmax ; opts_qv@mpMinLonF = xlonmin ; opts_qv@mpMaxLonF = xlonmax ; opts_qv@mpLeftCornerLatF = 31. ; opts_qv@mpRightCornerLatF = 43. ; opts_qv@mpLeftCornerLonF = -128.5 ; opts_qv@mpRightCornerLonF = -113.5 ; opts_qv@mpMinLatF = 31. ; opts_qv@mpMaxLatF = 43. ; opts_qv@mpMinLonF = -128.5 ; opts_qv@mpMaxLonF = -113.5 ; opts_qv@tiMainOn = True ; opts_qv@tiMainString = "Mixing Ratio (g/kg) - "+times1(it) ; opts_qv@tiMainFont = "helvetica" ; opts_qv@tiMainFontHeightF = 0.012 ; opts_qv@tiMainFontColor = "black" ; opts_qv@FieldTitle = "Potential Temperature (K) - "+times1(it) ; opts_qv@FieldTitleFontFont = "helvetica" ; opts_qv@FieldTitleFontHeightF = 0.012 ; opts_qv@FieldTitleFontColor = "black" opts_qv@tmXBLabelFontHeightF = 0.013 opts_qv@tmYLLabelFontHeightF = 0.013 mpres@tmXBLabelFontHeightF = 0.013 mpres@tmYLLabelFontHeightF = 0.013 mpres@tiMainOn = True mpres@tiMainString = "WRF - Mixing Ratio (g/kg) - "+times1(it) mpres@tiMainStringOrthogonalPosF = 0.02 mpres@tiMainFont = "helvetica" mpres@tiMainFontHeightF = 0.013 mpres@tiMainFontColor = "black" ; mpres@gsnLeftString = "WRF - Mixing Ratio (g/kg)" ; mpres@gsnRightString = times1(it) ; mpres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; mpres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; mpres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; pltres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; pltres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@tiMainOn = True ; pltres@tiMainFont = "helvetica" ; pltres@tiMainFontHeightF = 0.012 ; pltres@tiMainFontColor = "black" ; pltres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; pltres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; pltres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@tiMainString = "Mixing Ratio (g/kg) - "+times1(it); overwrite Field Title ; opts_qv@gsnLeftString = "Mixing Ratio (g/kg)" ; opts_qv@gsnLeftStringFontHeightF = 0.012 ; opts_qv@gsnCenterString = "Mixing Ratio (g/kg)"+times1(it) ; opts_qv@gsnCenterStringFontHeightF = 0.013 ; opts_qv@gsnCenterStringOrthogonalPosF = 0.013 ; opts_qv@gsnRightString = times1(it) ; opts_qv@gsnRightStringFontHeightF = 0.012 ; if(ifile.eq.0) then ; opts_qv@tiMainString = "Mixing Ratio (g/kg) - "+times1(it); overwrite Field Title ; opts_qv@gsnLeftString = "Mixing Ratio (g/kg) - "+times1(it) ; opts_qv@FieldTitle = "Mixing Ratio (g/kg) "+times1(it) ; else ; opts_qv@tiMainString = "Mixing Ratio (g/kg) - "+times2(it); overwrite Field Title ; opts_qv@gsnLeftString = "Mixing Ratio (g/kg) - "+times2(it) ; opts_qv@FieldTitle = "Mixing Ratio (g/kg) "+times2(it) ; end if ; cmap = read_colormap_file("BlAqGrYeOrReVi200") ; opts_qv@vcLevelPalette = cmap(6:193,:) ; opts_qv@pmTickMarkDisplayMode = "Always" ; opts_qv@mpFillOn = False ; turn off map fill opts_qv@mpOutlineDrawOrder = "PostDraw" ; draw continental outline last ; opts_qv@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_qv@mpDataSetName = "Earth..4" ; opts_qv@mpGeophysicalLineColor = "Red" ; opts_qv@mpGeophysicalLineThicknessF = 2 ; opts_qv@mpNationalLineThicknessF = 2 ; opts_qv@NationalLineColor = "Red" ; opts_qv@mpOutLineOn = True ; opts_qv@mpGridSpacingF = 2.5 ; opts_qv@mpGridLineThicknessF = 2.0 ; opts_qv@mpGridLineColor = "Red" ; opts_qv@mpPerimLineColor = "Red" ; opts_qv@mpOutLineOn = True ; opts_qv@pmTickMarkDisplayMode = "Always" ; opts_qv@mpFillOn = False ; turn off map fill ; opts_qv@mpFillOn = True ; turn on map fill ; opts_qv@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_qv@mpDataSetName = "Earth..4" ; opts_qv@mpFillBoundarySets = "National" ; opts_qv@mpGridLineColor = "Blue" ; opts_qv@mpPerimLineColor = "Blue" ; opts_qv@mpLimbLineColor = "Blue" ; opts_qv@mpPerimLineColor = "Blue" ; opts_qv@mpGeophysicalLineColor = "Blue" ; opts_qv@mpNationalLineColor = "Blue" ; opts_qv@mpUSStateLineColor = "Blue" ; opts_qv@mpGridLineThicknessF = 2 opts_qv@mpGeophysicalLineThicknessF = 2 opts_qv@mpNationalLineThicknessF = 2 opts_qv@mpUSStateLineThicknessF = 2 ; opts_qv@mpPerimOn = False ; opts_qv@mpGridAndLimbOn = False ; opts_qv@gsnTickMarksOn = False ; opts_qv@mpDataSetName = "Earth..3" ; opts_qv@mpDataBaseVersion = "MediumRes" ; opts_qv@mpOutlineOn = True ; opts_qv@mpOutlineSpecifiers = "Conterminous US: Counties" ; opts_qv@mpOFillAreaSpecifiers = "Conterminous US: Counties" ; opts_qv@mpOutlineBoundarySets = "USStates" ; opts_qv@mpFillOn = True ; opts_qv@mpSpecifiedFillColors = "red" ; opts_qv@mpFillAreaSpecifiers = "California : 05" ; opts_qv@mpOceanFillColor = "transparent" ; opts_qv@mpLandFillColor = "transparent" ; opts_qv@mpInlandWaterFillColor = "transparent" ; opts_qv@mpUSStateLineColor = "red" ; opts_qv@mpNationalLineColor = "red" ; opts_qv@mpGeophysicalLineColor = "red" opts_qv@cnFillOn = True ; opts_qv@cnLevelSelectionMode = "AutomaticLevels" opts_qv@cnLevelSelectionMode = "ManualLevels" opts_qv@ContourParameters = (/ 0., 12., 0.5/) ; opts_qv@gsnSpreadColorEnd = -2 opts_qv@lbLabelAutoStride = True opts_qv@lbAutoManage = False ; we control label bar opts_qv@pmLabelBarDisplayMode = "Always" ; turns on label bar opts_qv@lbOrientation = "Horizontal" ; ncl default is vertical opts_qv@pmLabelBarSide = "Bottom" ; default is right opts_qv@lbLabelStride = 2 ; skip every other label opts_qv@pmLabelBarWidthF = 0.80 ; default is shorter opts_qv@pmLabelBarHeightF = 0.30 ; default is taller opts_qv@lbLabelFontHeightF = 0.015 ; default is HUGE opts_qv@lbTitleFontHeightF = 0.015 ; opts_qv@lbTitleOn = True ; turn on title ; opts_qv@lbTitleString = "Mixing Ratio (g/kg)" ; opts_qv@pmLabelBarOrthogonalPosF = 0.0001 ; move whole thing down opts_qv@lbPerimOn = False ; default has box contour_qv = wrf_contour(a,wks,qv,opts_qv) plot = wrf_map_overlays(a,wks,(/contour_qv/),pltres,mpres) delete(opts_qv) ; map_cont_qv = wrf_map_overlays(a,wks,qv,True,mpres) ;******************************************** ; Plotting options for relative humidity ;******************************************** if(ifile.eq.0) then rh = rh1(it,0,:,:) else rh = rh2(it,0,:,:) end if gsn_define_colormap (wks,"BlAqGrYeOrReVi200") gsn_reverse_colormap(wks) ; Reverse the color map. opts_rh = res opts_rh@Footer = False opts_rh@Header = False opts_rh@NoHeaderFooter = True opts_rh@gsnMaximize = True ; maximize plot in frame ; opts_rh@vpXF = 0.1 ; opts_rh@vpYF = 0.95 opts_rh@mpShapeMode = "FreeAspect" opts_rh@vpWidthF = 0.8 opts_rh@vpHeightF = 0.75 ; mpres@vpXF = 0.1 ; mpres@vpYF = 0.97 ; mpres@mpShapeMode = "FreeAspect" ; mpres@vpWidthF = 0.85 ; mpres@vpHeightF = 0.68 opts_rh@mpLeftCornerLatF = xlatmin opts_rh@mpRightCornerLatF = xlatmax opts_rh@mpLeftCornerLonF = xlonmin opts_rh@mpRightCornerLonF = xlonmax ; opts_rh@mpMinLatF = xlatmin ; opts_rh@mpMaxLatF = xlatmax ; opts_rh@mpMinLonF = xlonmin ; opts_rh@mpMaxLonF = xlonmax ; opts_rh@mpLeftCornerLatF = 31. ; opts_rh@mpRightCornerLatF = 43. ; opts_rh@mpLeftCornerLonF = -128.5 ; opts_rh@mpRightCornerLonF = -113.5 ; opts_rh@mpMinLatF = 31. ; opts_rh@mpMaxLatF = 43. ; opts_rh@mpMinLonF = -128.5 ; opts_rh@mpMaxLonF = -113.5 ; opts_rh@tiMainOn = True ; opts_rh@tiMainString = "Relative Humidity (%) - "+times1(it) ; opts_rh@tiMainFontHeightF = 0.012 ; opts_rh@tiMainFontColor = "black" ; opts_rh@tiMainFont = "helvetica" ; opts_rh@FieldTitle = "Potential Temperature (K) - "+times1(it) ; opts_rh@FieldTitleFontFont = "helvetica" ; opts_rh@FieldTitleFontHeightF = 0.012 ; opts_rh@FieldTitleFontColor = "black" opts_rh@tmXBLabelFontHeightF = 0.013 opts_rh@tmYLLabelFontHeightF = 0.013 mpres@tmXBLabelFontHeightF = 0.013 mpres@tmYLLabelFontHeightF = 0.013 mpres@tiMainOn = True mpres@tiMainString = "WRF - Relative Humidity (%) - "+times1(it) mpres@tiMainStringOrthogonalPosF = 0.02 mpres@tiMainFontHeightF = 0.013 mpres@tiMainFontColor = "black" mpres@tiMainFont = "helvetica" ; mpres@gsnLeftString = "WRF - Relative Humidity (%)" ; mpres@gsnRightString = times1(it) ; pltres@tiMainOn = True ; pltres@tiMainFontHeightF = 0.012 ; pltres@tiMainFontColor = "black" ; pltres@tiMainFont = "helvetica" ; mpres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; mpres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; mpres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@gsnMajorLatSpacing = 5 ; change maj lat tm spacing ; pltres@gsnMajorLonSpacing = 5 ; change maj lon tm spacing ; pltres@tmXBMinorOn = True ; no lon minor tickmarks ; pltres@tiMainString = "Relative Humidity (%) - "+times1(it); overwrite Field Title ; opts_rh@gsnLeftString = "Relative Humidity (%)" ; opts_rh@gsnLeftStringFontHeightF = 0.012 ; opts_rh@gsnCenterString = "Relative Humidity (%)"+times1(it) ; opts_rh@gsnCenterStringFontHeightF = 0.013 ; opts_rh@gsnCenterStringOrthogonalPosF = 0.013 ; opts_rh@gsnRightString = times1(it) ; opts_rh@gsnRightStringFontHeightF = 0.012 ; if(ifile.eq.0) then ; opts_rh@tiMainString = "Relative Humidity (%) - "+times1(it); overwrite Field Title ; opts_rh@gsnLeftString = "Relative Humidity (%)~C~Wind Vectors - "+times1(it) ; opts_rh@gsnLeftString = "Relative Humidity (%) - "+times1(it) ; opts_rh@FieldTitle = "Relative Humidity (%) "+times1(it) ; else ; opts_rh@tiMainString = "Relative Humidity (%) - "+times2(it); overwrite Field Title ; opts_rh@gsnLeftString = "Relative Humidity (%)~C~Wind Vectors - "+times2(it) ; opts_rh@gsnLeftString = "Relative Humidity (%) - "+times2(it) ; opts_rh@FieldTitle = "Relative Humidity (%) "+times2(it) ; end if ; cmap = read_colormap_file("BlAqGrYeOrReVi200") ; opts_rh@vcLevelPalette = cmap(6:193,:) ; opts_rh@pmTickMarkDisplayMode = "Always" ; opts_rh@mpFillOn = False ; turn off map fill opts_rh@mpOutlineDrawOrder = "PostDraw" ; draw continental outline last ; opts_rh@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_rh@mpDataSetName = "Earth..4" ; opts_rh@mpGeophysicalLineColor = "Red" ; opts_rh@mpGeophysicalLineThicknessF = 2 ; opts_rh@mpNationalLineThicknessF = 2 ; opts_rh@NationalLineColor = "Red" ; opts_rh@mpOutLineOn = True ; opts_rh@mpGridSpacingF = 2.5 ; opts_rh@mpGridLineThicknessF = 2.0 ; opts_rh@mpGridLineColor = "Red" ; opts_rh@mpPerimLineColor = "Red" ; opts_rh@mpOutLineOn = True ; opts_rh@pmTickMarkDisplayMode = "Always" ; opts_rh@mpFillOn = False ; turn off map fill ; opts_rh@mpFillOn = True ; turn on map fill ; opts_rh@mpOutlineBoundarySets = "GeophysicalAndUSStates" ; state boundaries ; opts_rh@mpDataSetName = "Earth..4" ; opts_rh@mpFillBoundarySets = "National" ; opts_rh@mpGridLineColor = "Blue" ; opts_rh@mpPerimLineColor = "Blue" ; opts_rh@mpLimbLineColor = "Blue" ; opts_rh@mpPerimLineColor = "Blue" ; opts_rh@mpGeophysicalLineColor = "Blue" ; opts_rh@mpNationalLineColor = "Blue" ; opts_rh@mpUSStateLineColor = "Blue" ; opts_rh@mpGridLineThicknessF = 2 opts_rh@mpGeophysicalLineThicknessF = 2 opts_rh@mpNationalLineThicknessF = 2 opts_rh@mpUSStateLineThicknessF = 2 ; opts_rh@mpPerimOn = False ; opts_rh@mpGridAndLimbOn = False ; opts_rh@gsnTickMarksOn = False ; opts_rh@mpDataSetName = "Earth..3" ; opts_rh@mpDataBaseVersion = "MediumRes" ; opts_rh@mpOutlineOn = True ; opts_rh@mpOutlineSpecifiers = "Conterminous US: Counties" ; opts_rh@mpOFillAreaSpecifiers = "Conterminous US: Counties" ; opts_rh@mpOutlineBoundarySets = "USStates" ; opts_rh@mpFillOn = True ; opts_rh@mpSpecifiedFillColors = "red" ; opts_rh@mpFillAreaSpecifiers = "California : 05" ; opts_rh@mpOceanFillColor = "transparent" ; opts_rh@mpLandFillColor = "transparent" ; opts_rh@mpInlandWaterFillColor = "transparent" ; opts_rh@mpUSStateLineColor = "red" ; opts_rh@mpNationalLineColor = "red" ; opts_rh@mpGeophysicalLineColor = "red" ; opts_rh@cnFillOn = True ; turn on color ; opts_rh@cnFillPalette = "gui_default" ; opts_rh@cnFillColors = "Red" ; opts_rh@cnFillMode = "RasterFill" ; opts_rh@cnRasterSmoothingOn = True opts_rh@cnFillOn = True ; opts_rh@cnLevelSelectionMode = "AutomaticLevels" opts_rh@cnLevelSelectionMode = "ManualLevels" opts_rh@ContourParameters = (/ 0., 100., 5./) ; opts_rh@gsnSpreadColorEnd = -2 opts_rh@lbLabelAutoStride = True opts_rh@lbAutoManage = False ; we control label bar opts_rh@pmLabelBarDisplayMode = "Always" ; turns on label bar opts_rh@lbOrientation = "Horizontal" ; ncl default is vertical opts_rh@pmLabelBarSide = "Bottom" ; default is right opts_rh@lbLabelStride = 2 ; skip every other label opts_rh@pmLabelBarWidthF = 0.80 ; default is shorter opts_rh@pmLabelBarHeightF = 0.30 ; default is taller opts_rh@lbLabelFontHeightF = 0.015 ; default is HUGE opts_rh@lbTitleFontHeightF = 0.015 ; opts_rh@lbTitleOn = True ; turn on title ; opts_rh@lbTitleString = "Relative Humidity (%)" ; opts_rh@pmLabelBarOrthogonalPosF = 0.0001 ; move whole thing down opts_rh@lbPerimOn = False ; default has box contour_rh = wrf_contour(a,wks,rh,opts_rh) plot = wrf_map_overlays(a,wks,(/contour_rh/),pltres,mpres) delete(opts_rh) ; map_cont_rh = wrf_map_overlays(a,wks,rh,True,mpres) end do ; end of time loop end do ; end of file loop end