<font size="3"><span style="line-height: 24px;">Hi,all</span></font><div><font size="3"><span style="line-height: 24px;">&nbsp; I have calculated the max Qgraup of each time from the wrfout data with the following program,</span></font></div><div><font size="3"><span style="line-height: 24px;">but how can I know&nbsp;&nbsp;the (i,j)or (lat,lon) of the max&nbsp;<span style="display: none; line-height: 0px; " id="_editor_bookmark_start_1">‍</span></span></font><span style="font-size: medium; line-height: 24px; ">Qgraup of each time?<span style="display: none; line-height: 0px; " id="_editor_bookmark_start_2">‍</span></span></div><div><span style="font-size: medium; line-height: 24px; ">any information will be appreciated!</span></div><div><font size="3"><span style="line-height: 24px;"><br></span></font></div><div><font size="3"><div>&nbsp;&nbsp;ascii_file = "shishidataGRAUP1B333.txt"</div><div>&nbsp; str = new(ntimes, "string")</div><div><span style="display: none; line-height: 0px; " id="_editor_bookmark_start_4">‍</span>&nbsp; do it =0, ntimes-1,1 &nbsp; &nbsp; &nbsp; &nbsp;; TIME LOOP</div><div style="line-height: 24px; "><br></div><div style="line-height: 24px; ">&nbsp; &nbsp; print("Working on time: " + times(it) )</div><div style="line-height: 24px; ">&nbsp; &nbsp; res@TimeLabel = times(it) &nbsp; ; Set Valid time to use on plots</div><div style="line-height: 24px; "><br></div><div style="line-height: 24px; ">;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;</div><div style="line-height: 24px; ">; First get the variables we will need &nbsp; &nbsp; &nbsp; &nbsp;</div><div style="line-height: 24px; ">&nbsp; &nbsp; if(isfilevar(a,"QGRAUP"))</div><div style="line-height: 24px; ">&nbsp; &nbsp; &nbsp; qi = wrf_user_getvar(a,"QGRAUP",it)</div><div style="line-height: 24px; ">&nbsp; &nbsp; &nbsp; qi = qi*1000.</div><div style="line-height: 24px; ">&nbsp; &nbsp; &nbsp; qi@units = "g/kg"&nbsp;</div><div style="line-height: 24px; "><span class="Apple-tab-span" style="white-space:pre">        </span> &nbsp;qis=dim_max_n_Wrap(qi(13:15,:,:),0)</div><div style="line-height: 24px; ">&nbsp; &nbsp; end if</div><div style="line-height: 24px; "><span class="Apple-tab-span" style="white-space:pre">        </span>;print(qis)</div><div style="line-height: 24px; "><span class="Apple-tab-span" style="white-space:pre">        </span>printVarSummary(qis)</div><div style="line-height: 24px; ">&nbsp; &nbsp; wmax1 = dim_max_n_Wrap(qis(:,:),0)</div><div style="line-height: 24px; "><span class="Apple-tab-span" style="white-space:pre">        </span>;print(wmax1)</div><div style="line-height: 24px; "><span class="Apple-tab-span" style="white-space:pre">        </span>printVarSummary(wmax1)</div><div style="line-height: 24px; ">&nbsp; &nbsp; wmax2 = dim_max_n_Wrap(qis,0)</div><div style="line-height: 24px; "><span class="Apple-tab-span" style="white-space:pre">        </span>;print(wmax2)</div><div style="line-height: 24px; "><span class="Apple-tab-span" style="white-space:pre">        </span>printVarSummary(wmax2)</div><div style="line-height: 24px; ">&nbsp; &nbsp; wmax=max(wmax2)<span style="display: none; line-height: 0px; " id="_editor_bookmark_start_0">‍</span></div><div>&nbsp;str(it) = sprintf("%12.8f", wmax)<span style="display: none; line-height: 0px; " id="_editor_bookmark_start_3">‍</span></div></font></div>