how to know the (i,j)or (lat,lon) of the max

Hi,all I have calculated the max Qgraup of each time from the wrfout data with the following program,
but how can I know the (i,j)or (lat,lon) of the max ‍Qgraup of each time?‍
any information will be appreciated!


  ascii_file = "shishidataGRAUP1B333.txt"
  str = new(ntimes, "string")
‍ do it =0, ntimes-1,1 ; TIME LOOP


    print("Working on time: " + times(it) )
    res@TimeLabel = times(it) ; Set Valid time to use on plots


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; First get the variables we will need
    if(isfilevar(a,"QGRAUP"))
      qi = wrf_user_getvar(a,"QGRAUP",it)
      qi = qi*1000.
      qi@units = "g/kg"
          qis=dim_max_n_Wrap(qi(13:15,:,:),0)
    end if
        ;print(qis)
        printVarSummary(qis)
    wmax1 = dim_max_n_Wrap(qis(:,:),0)
        ;print(wmax1)
        printVarSummary(wmax1)
    wmax2 = dim_max_n_Wrap(qis,0)
        ;print(wmax2)
        printVarSummary(wmax2)
    wmax=max(wmax2)‍
 str(it) = sprintf("%12.8f", wmax)‍

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk



Received on Tue Jul 01 08:25:45 2014