<div dir="ltr"><div class="gmail_default" style="font-family:georgia,serif">Hi,</div><div class="gmail_default" style="font-family:georgia,serif"><br></div><div class="gmail_default" style="font-family:georgia,serif">I'm trying to create a panel contour plot of cloud liquid water, and I particularly want to have as many contour levels as possible, i.e. I need to be able to show the smallest amount of cloud water. In order to do so, I've set:</div><div class="gmail_default" style="font-family:georgia,serif">res@ContourParameters = (/0.0, max(qc_plane), 0.01 /)</div><div class="gmail_default" style="font-family:georgia,serif"><br></div><div class="gmail_default" style="font-family:georgia,serif">However, the plots that I get don't have any designated color for 0 value.</div><div class="gmail_default" style="font-family:georgia,serif">I use a similar code to produce one such plot (instead of a panel plot), and it works fine, in the sense that I can divide the panel bar to the number of contour levels.<br></div><div class="gmail_default" style="font-family:georgia,serif"><br></div><div class="gmail_default" style="font-family:georgia,serif">Here is my NCL code and the plot that I get: (NCL.6.4.0)<br></div><div class="gmail_default" style="font-family:georgia,serif"><br clear="all"></div> type = "png"<br><br> wks = gsn_open_wks(type,"Cloud_panel")<br><br> times = wrf_user_list_times(a)<br> time = 166<br><br> p = wrf_user_getvar(a, "pressure", time) ; Total Pressure<br> z = wrf_user_getvar(a,"z",time)<br> qc = wrf_user_getvar(a,"QCLOUD", time)<br><br> qc = qc*1000.<br> qc@units = "g/kg"<br><br> res = True<br> res@NoHeaderFooter = True<br><br> pltres = True<br> pltres@PanelPlot = True<br> mpres = True<br><br> ; Generate Contours<br><br> height_levels = (/750., 1000., 1200, 1400./) ; heigth levels to plot<br> nlevels = dimsizes(height_levels) ; number of height levels<br> nplots = 4<br> plots = new(nplots,"graphic")<br> iplot = 0<br><br> do level = 0,nlevels-1<br> res@cnFillOn = True<br> res@cnFillPalette = "WhViBlGrYeOrRe"<br><div style="font-family:georgia,serif" class="gmail_default"> res@lbLabelBarOn = False<br> res@cnLinesOn = False<br> res@cnLevelSpacingF = 0.1<br><br> height = height_levels(level)<br> qc_plane = wrf_user_intrp3d(qc,z,"h",height,0.,False)<br><br> res@ContourParameters = (/0.0, max(qc_plane), 0.01 /)<br><br> contour = wrf_contour(a,wks,qc_plane,res)<br> pltres@CommonTitle = True<br> pltres@PlotTitle = height+"m"<br> plots(iplot) = wrf_map_overlays(a, wks, (/contour/), pltres, mpres)<br> delete(contour)<br> iplot = iplot + 1<br> end do<br><br> ; Panel the WRF Plots<br> pnlres = True<br> pnlres@txString = qc@description + " (" + qc@units + ")"+ "at" + " at "+ times(time)<br> pnlres@gsnPanelYWhiteSpacePercent = 10 ; Add white space b/w plots.<br> pnlres@gsnPanelLabelBar = True ; Turn on common labelbar<br><br> pnlres@lbLabelAlignment = "BoxCenters"<br> pnlres@lbLabelsOn = True<br> pnlres@lbPerimOn = False<br> pnlres@lbTitleOn = False</div><div style="font-family:georgia,serif" class="gmail_default"><br></div><div style="font-family:georgia,serif" class="gmail_default"> pnlres@pmLabelBarWidthF = 0.6<br><br> gsn_panel(wks,(/plots/),(/2,2/),pnlres)<br><br>end<br></div><br>-- <br><div class="gmail_signature"><div dir="ltr"><font face="georgia, serif"><i>Melissa Kazemi Rad</i></font><div><font face="georgia, serif"><i>Atmospheric Sciences Department</i></font></div><div><font face="georgia, serif"><i>Rutgers University</i></font></div></div></div>
</div>