<div dir="ltr">Hi all, <br><div><br></div><div>I have some problems with plotting wrf output, the values are not in the correct place here is my code:</div><div><br></div><div>;; -*- mode: emacs-lisp; coding: emacs-mule; -*-;  -*- mode: ncl;-*-<br>;***************************************************************************<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br>load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" <br><br>ana=addfile("wrfout_d01_2015-08-13","r")<br><br>lat1=ana->XLAT(0,:,0)<br>lon1=ana->XLONG(0,0,:)<br>P=ana->P<br>pm10=ana->PM10<br>pp = wrf_user_getvar(ana, "pressure",-1)<br>delete(ana)<br><br>printVarSummary(pm10)<br>  printMinMax(pm10,0)<br>  print("====")<br>  printVarSummary(P)<br>  printMinMax(P,0)<br>  print("====")<br><br><br>po=(/1000,925,850,700,600,500,400,300,250,200,150,100/)<br>po!0  = "po"<br>po@long_name = "Pressure Level"<br>po@units     = "hPa"<br><br>To = int2p_n_Wrap (pp,pm10,po,1,1)<br><br><br>printVarSummary(To)<br>   printMinMax(To,0)<br>print("====")<br><br><br>lat= lat1 ;fspan(39.43722,49.86291,95)<br>lon= lon1 ;fpan(-92.1441,-75.8559,114)+360<br>lat@longname = "latitude"<br>lat@units = "degrees_north"<br>lon@longname ="longitude"<br>lon@units = "degrees_east"<br><br>To!2         = "lat"                   <br>To!3         = "lon" <br>To&lon     =  lon <br>To&lat     =  lat  <br><br>wks = gsn_open_wks("eps","pm10_")<br>;gsn_define_colormap(wks,"BlWhRe")<br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>  res                     = True               ; plot mods desired<br>  res@cnFillOn            = True              ; turn on color   <br>  res@gsnSpreadColors     = True               ; use full color map<br>  res@cnLinesOn           = False              ; no contour lines<br>  res@cnLineLabelsOn      = False              ; no line labels<br><br>  res@mpMaxLatF =  50.                         ; specify the plot domain<br>  res@mpMinLatF = 39.                         ;                         <br>  res@mpMinLonF = 267                          ;<br>  res@mpMaxLonF = 285    <br>  res@mpCenterLonF= 276 <br>  res@gsnAddCyclic        = False <br><br> res@cnLevelSelectionMode =  "AutomaticLevels"   ; AutomaticLevels<br><br>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<br>res@tiMainString        = " Test PM10"<br><br>res@gsnDraw              = False           ; Do not draw plot<br>res@gsnFrame             = False           ; Do not advance frome<br><br><br>plot= gsn_csm_contour_map(wks,To(1,0,:,:),res)<br><br>  draw (plot)<br>  frame(wks)<br>delete(wks)<br><br></div><div><br></div><div>Thanks,</div><div>Ana</div></div>