<DIV>&nbsp;Hi:</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; All,I am trying to&nbsp;find out&nbsp;the lat and lon location of the max graup of wrfout data,I have used the&nbsp;funtion " ind_resolve","maxind" and "wrf_user_ij_to_ll" according to the web examples.</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; But after I compared the output data with the contour plot that I ploted ,them are not <SPAN closure_uid_407902488="7">corresponding.</SPAN></DIV>
<DIV><SPAN closure_uid_407902488="7">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I thought maybe I confused the lat and lon ,but after I switch it ,the output data are not <SPAN closure_uid_407902488="7">correspond to the plots too.</SPAN></SPAN></DIV>
<DIV><SPAN closure_uid_407902488="7"><SPAN closure_uid_407902488="7">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; I can not figure out where the problem is? I think maybe it is because of the zoom ?</SPAN></SPAN></DIV>
<DIV><SPAN closure_uid_407902488="7"><SPAN closure_uid_407902488="7">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;But I do not know how to fix it,can you guys fix this?</SPAN></SPAN></DIV>
<DIV>&nbsp;</DIV>
<DIV><BR></DIV>
<DIV>This is my script:</DIV>
<DIV>
<DIV>;&nbsp;&nbsp; Example script to produce plots for a WRF real-data run,<BR>;&nbsp;&nbsp; with the ARW coordinate dynamics option.</DIV>
<DIV>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/wrf/WRFUserARW.ncl"</DIV>
<DIV>begin<BR>;<BR>; The WRF ARW input file.&nbsp; <BR>; This needs to have a ".nc" appended, so just do it.<BR>&nbsp; a = addfile("/backup/hl/z_20min/wrfout_d03_2010-08-11_00:00:00"+".nc","r")&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </DIV>
<DIV><BR>; We generate plots, but what kind do we prefer?<BR>;&nbsp; type = "eps"<BR>&nbsp;type = "pdf"<BR>; type = "ps"<BR>; type = "ncgm"<BR>&nbsp; wks = gsn_open_wks(type,"plt_Cloudgraup_20min_20100811")<BR>&nbsp;&nbsp; ;colors = (/"white","black","white","royal blue","light sky blue",\<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;&nbsp; "powder blue","light sea green","pale green","wheat","brown",\<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; "pink"/)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;gsn_define_colormap(wks,"BlAqGrYeOrRevi200")&nbsp;&nbsp; ; overwrite the .hluresfile color map<BR><BR>; Set some basic resources<BR>&nbsp; res = True<BR>&nbsp; <A href="mailto:res@MainTitle">res@MainTitle</A> = "REAL-TIME WRF"</DIV>
<DIV>&nbsp; mpres&nbsp; = True&nbsp; ; Map resources<BR>&nbsp; <A href="mailto:mpres@mpOutlineOn">mpres@mpOutlineOn</A> = False&nbsp; ; Turn off map outlines<BR>&nbsp; <A href="mailto:mpres@mpFillOn">mpres@mpFillOn</A>&nbsp;&nbsp;&nbsp; = False&nbsp; ; Turn off map fill<BR>&nbsp; <A href="mailto:mpres@mpGridAndLimbOn">mpres@mpGridAndLimbOn</A> = True<BR>&nbsp;;res@mpProjection&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = "Lambert"<BR>&nbsp; pltres = True ; Plot resources<BR>&nbsp; <A href="mailto:pltres@PanelPlot">pltres@PanelPlot</A>&nbsp; = True&nbsp;&nbsp; ; Tells wrf_map_overlays not to remove overlays</DIV>
<DIV><BR>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<BR>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</DIV>
<DIV>; What times and how many time steps are in the data set?<BR>&nbsp; times = wrf_user_getvar(a,"times",-1)&nbsp; ; get all times in the file<BR>&nbsp; ntimes = dimsizes(times)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; number of times in the file<BR>&nbsp;lat = new(ntimes,float)&nbsp; ;creat new variable to store the lat of the max qi<BR>&nbsp;lon = new(ntimes,float)<BR>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</DIV>
<DIV>&nbsp; do it =0, ntimes-1,1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; TIME LOOP</DIV>
<DIV>&nbsp;&nbsp;&nbsp; print("Working on time: " + times(it) )<BR>&nbsp;&nbsp;&nbsp; <A href="mailto:res@TimeLabel">res@TimeLabel</A> = times(it)&nbsp;&nbsp; ; Set Valid time to use on plots<BR>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;<BR>; First get the variables we will need&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp; if(isfilevar(a,"QGRAUP"))<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qi = wrf_user_getvar(a,"QGRAUP",it )<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qi = qi*1000.<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:qi@units">qi@units</A> = "g/kg"&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; qiav = dim_avg_n_Wrap(qi, 0)<BR>&nbsp;&nbsp;&nbsp; end if<BR>&nbsp;&nbsp; tes=True<BR>&nbsp;&nbsp; <A href="mailto:tes@returnInt">tes@returnInt</A> = False<BR>&nbsp;&nbsp; loc1=wrf_user_ll_to_ij(a,105.0,32.0,tes)<BR>&nbsp;&nbsp; ;print("X/Y location is: "+ loc1)<BR>&nbsp;&nbsp; <BR>&nbsp;&nbsp; loc2=wrf_user_ll_to_ij(a,112.0,38.0,tes)<BR>&nbsp; ; print("X/Y location is: "+ loc2)<BR>&nbsp;&nbsp; x_start = 15<BR>&nbsp;&nbsp; x_end&nbsp;&nbsp; = 227<BR>&nbsp;&nbsp; y_start = 71<BR>&nbsp;&nbsp; y_end&nbsp;&nbsp; = 274</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp; level = 13</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; display_level = level + 2<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; opts = res<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:opts@cnFillOn">opts@cnFillOn</A>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = True<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:opts@gsnSpreadColors">opts@gsnSpreadColors</A>&nbsp; = False<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:opts@ContourParameters">opts@ContourParameters</A>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = (/ 1, 5, 0.5 /)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:opts@PlotLevelID">opts@PlotLevelID</A>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = " Level -10 - -20 ~S~o~N~C " <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:opts@gsnDraw">opts@gsnDraw</A>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; False&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:opts@gsnFrame">opts@gsnFrame</A>&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; False</DIV>
<DIV><BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (isvar("qi"))<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qis&nbsp; = qi(level + (/0,1,2,3,4,5,6/),:,:)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qisum = dim_sum_n_Wrap(qis, 0)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; mpres1 = True<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:mpres1@ZoomIn">mpres1@ZoomIn</A> = True<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:mpres1@Xstart">mpres1@Xstart</A> = x_start<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:mpres1@Ystart">mpres1@Ystart</A> = y_start<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:mpres1@Xend">mpres1@Xend</A>&nbsp;&nbsp; = x_end<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <A href="mailto:mpres1@Yend">mpres1@Yend</A>&nbsp;&nbsp; = y_end<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;&nbsp; printVarSummary(qisum)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qi_zoom = qisum(y_start:y_end,x_start:x_end)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; printVarSummary(qi_zoom)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;;;;;;;convert to 1D;;;;;;;;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; qi_zoom1D&nbsp; = ndtooned(qi_zoom)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dsizes_a = dimsizes(qi_zoom)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;;;;;;resolve the 1D indices back to their original 2D arry.;;;<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; indices = ind_resolve(maxind(qi_zoom1D),dsizes_a)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; print(indices)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ilatlon = wrf_user_ij_to_ll(a, indices(0,1),indices(0,0),True)&nbsp; ; select the latitude index where the X array is at its' maximum <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lon(it)= ilatlon(0)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lat(it)= ilatlon(1)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; contour = wrf_contour(a,wks,qi_zoom,opts)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres1)</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; delete(contour)<BR>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; end if<BR>;&gt;============================================================&lt;<BR>;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; add&nbsp; map<BR>;&gt;------------------------------------------------------------&lt;</DIV>
<DIV>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<BR>&nbsp; shp_name2&nbsp;&nbsp;&nbsp; = "/public/home/huanglei/map/shaanxi_city_l.shp"</DIV>
<DIV>&nbsp; prres=True<BR>&nbsp; <A href="mailto:prres@gsLineThicknessF">prres@gsLineThicknessF</A> = 1.0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp; <A href="mailto:prres@gsLineColor">prres@gsLineColor</A> = "black"<BR>&nbsp; plotcn3 = gsn_add_shapefile_polylines(wks,plot,shp_name2,prres)</DIV>
<DIV>&nbsp;&nbsp; txres2&nbsp; = True<BR>&nbsp;&nbsp; <A href="mailto:txres2@txFont">txres2@txFont</A>&nbsp; = 10<BR>&nbsp;&nbsp; <A href="mailto:txres2@txFontHeightF">txres2@txFontHeightF</A> =0.01<BR>&nbsp;&nbsp; <A href="mailto:txres2@txFontColor">txres2@txFontColor</A> = "Blue"<BR>&nbsp;&nbsp; txdum1 =gsn_add_text(wks, plot, "Xian", 108.93,34.27, txres2)</DIV>
<DIV>&nbsp; draw(plot)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; This will draw the map and the shapefile outlines.<BR>&nbsp; frame(wks)<BR>&nbsp;&nbsp; delete(opts)</DIV>
<DIV>&nbsp;</DIV>
<DIV>;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</DIV>
<DIV>&nbsp; end do&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ; END OF TIME LOOP&nbsp;&nbsp;&nbsp;&nbsp; <BR>&nbsp;<BR>&nbsp;output_file = "/public/home/huanglei/z_data/20100811graup.txt"<BR>&nbsp;write_table(output_file,"w",[/times,lon,lat/],"%s%9.4f%9.4f")<BR>end<BR></DIV></DIV>
<DIV><BR></DIV>
<DIV><SPAN style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'">&nbsp;How can I slove the problem?</SPAN></DIV>
<DIV></DIV>