<div dir="ltr">Hello Everyone<div><br></div><div>I use ncl 6.4.0 under Cygwin on Windows, when I deal with wrfout_* files generated by WRF, I have this problem:</div><div>I just want to plot a cross section through a specific longitude and pressure levels, so I decide to use qrain2 = wrf_user_vert_interp(a,qrain,"pressure",level,opts)</div><div>to interpolate the model data to pressure level and then use qrain3 =  rcm2rgrid_Wrap(lat2d,lon2d,qrain2,lat,lon,1) to interpolate the data to longitude-latitude grids.</div><div><br></div><div>But I got this fatal message:</div><div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">fatal:wrf_vintrp: The pres and field arrays must have the same dimensionality<br>fatal:["Execute.c":8640]:Execute: Error occurred at or near line 2348 in file /usr/local/ncl6.4.0/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl<br>fatal:["Execute.c":8640]:Execute: Error occurred at or near line 83 in file wrf_cross_q_pl2.ncl</blockquote></div><div><br></div><div>Here is part of my code:</div><div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">lon = fspan(110,125,15)<br>lat = fspan(25,55,30)<br>level = ispan(200,1000,100)<br>lon@units = "degrees_east"<br>lat@units = "degrees_noeth"<br>level@units = "hPa"<br>a = addfile("wrfout_d01","r")<br>times = wrf_user_getvar(a,"times",-1)<br>lat2d = wrf_user_getvar(a,"XLAT",0)<br>lon2d = wrf_user_getvar(a,"XLONG",0)<br>opts = True<br>do it = 0,ntimes-1,12<br>qrain = wrf_user_getvar(a,"QRAIN",it)<br>qrain2 = wrf_user_vert_interp(a,qrain,"theta",level,opts)<br>qrain3 = rcm2rgrid_Wrap(lat2d,lon2d,qrain2,lat,lon,1)<br>qrain3!0 = "pressure"<br>qrain3!1 = "lat"<br>qrain3!2 = "lon"<br>qrain3&pressure = level<br>qrain3&lat = lat<br>qrain3&lon = lon<br>plot(ii) = gsn_csm_pres_hgt(wks,qrain3(:,:,{115}),res)<br>end do</blockquote></div><div><div><br></div></div><div>What's wrong with my code? And is this a correct method to plot a cross section through a specific longitude and pressure levels?</div><div><br></div><div>Fangda</div></div>