<div dir="ltr">I am working on zooming in on Texas/ New Mexico for a case study
that I'm performing with the WRF. I have the actual map zoomed in to the
region that I want, but the plotting of the qnifa variable doesn't
scale and stays as if the map projected is that of the CONUS. Can anyone
tell me what I'm doing wrong? I have attached two plots to illustrate the issue. The plot of the CONUS is correct in its placement and areas of qnifa, the zoomed plot is not.<br><br>Thanks in advance for the help!<br><br>-Andrew<br><br>load "$CSM/gsn_code.ncl"<br>load "$CSM/gsn_csm.ncl"<br>load "$CSM/contributed.ncl"<br><br>begin<br> res = True<br> res@cnLinesOn = False<br> res@cnFillOn = True ; color plot desired<br> res@cnLineLabelsOn = False ; turn off contour lines<br> res@cnInfoLabelOn = False ; turn off cn info label<br><br> res@mpGeophysicalLineColor= "black" ; color of continental outlines<br> res@mpPerimOn = True ; draw box around map<br> res@mpPerimDrawOrder = "PostDraw"<br> res@mpGridLineDashPattern = 2 ; lat/lon lines as dashed<br> res@mpOutlineBoundarySets = "GeophysicalAndUSStates"<br> res@mpUSStateLineColor = "black"<br><br> cmap = read_colormap_file("CBR_<wbr>drywet")<br><br> ff = addfile("<a href="http://wrfout.nc" target="_blank">wrfout.nc</a>","r")<br> if (ff@MAP_PROJ .eq. 1) then<br> res@mpProjection = "LambertConformal"<br> end if<br> res@mpLambertMeridianF = ff@CEN_LON<br> res@mpLambertParallel1F = ff@TRUELAT1<br> res@mpLambertParallel2F = ff@TRUELAT2<br> res@gsnAddCyclic = False<br> res@tfDoNDCOverlay = True<br> res@mpLimitMode = "Corners"<br> xlats = ff->XLAT<br> xlons = ff->XLONG<br> ijdim = dimsizes(xlats)<br> imax = ijdim(1)<br> jmax = ijdim(2)<br> res@mpLeftCornerLatF = 37<br> res@mpLeftCornerLonF = -110<br> res@mpRightCornerLatF = 28<br> res@mpRightCornerLonF = -96<br> res@mpFillOn <wbr> = True<br><br> res@cnLevelSelectionMode = "ExplicitLevels"<br> res@cnLevels = (/4E6,8E6,16E6,32E6,64E6/)<br> res@cnFillColors = (/7,6,5,4,3,2/)<br> res@cnMissingValFillColor = 13<br><br> hgt = ff->QNIFA(0,0,:,:)<br><br>
res@tiMainString = "Ice-Friendly Aerosol Concentration at Lowest Model
Level ~C~ W/ IC_BC and -30% SMOIS 19:00Z February 24, 2007"<br> res@tiMainFontHeightF = 0.015<br> res@gsnRightString = " "<br> res@gsnLeftString = " "<br> res@gsnMaximize = True<br> res@gsnFrame = False ; Don't advance the frame<br><br> wks = gsn_open_wks("png","qnifa_<wbr>plot")<br> gsn_define_colormap(wks, "CBR_drywet")<br><br> map = gsn_csm_contour_map(wks,hgt,<wbr>res)<br><br> frame(wks)<br><br>end</div>