[ncl-talk] Plot ind_resolve output

[Wang, Jiali]

Hi Barry,

I would say I need more help with converting the ‘i’ to map coordinate. Right now the output of ind_resolve is organized by ‘time’ dimension, but different grid cell would have different ‘time’ (first day) for having the precip greater than a threshold. Once I have a normal variable with X and Y dimension, I have no problem with making maps.

Thank you!
Jiali

From: Barry Lynn <barry.h.lynn at gmail.com>
Date: Sunday, December 13, 2020 at 1:01 AM
To: Jiali Wang <jialiwang at anl.gov>
Cc: "ncl-talk at mailman.ucar.edu" <ncl-talk at mailman.ucar.edu>, "ncl-talk at ucar.edu" <ncl-talk at ucar.edu>
Subject: Re: [ncl-talk] Plot ind_resolve output

Hi Jiali:

It would be helpful if you could be more specific.

Are you in need of help making a map, or converting the information you have in "i" to map coordinate information?

Barry

On Sun, Dec 13, 2020 at 7:44 AM Wang, Jiali via ncl-talk <ncl-talk at mailman.ucar.edu<mailto:ncl-talk at mailman.ucar.edu>> wrote:
Dear NCL users,

I am trying to find out the first and last day of precip greater than 0.1inch for all the grid cells over a study area. My data is at daily scale. I use ind_resolve and was able to print out the (time, lat, lon) information for all the grid cells that meet the threshold.

However, I am not sure how to move forward.. basically I’d like to make a map, showing the first (or last) day of 0.1’’ over all the grid cells.

Here is my short script (data is too large to share here). Can anyone please help?

Thanks!
Jiali
========================================
a = addfile("rain1995.nc.CA.wy","r")
rain1 = wrf_user_getvar(a,"RAINC",-1)
rain2 = wrf_user_getvar(a,"RAINNC",-1)
rain  = rain1 + rain2
printVarSummary(rain)

a1D      = ndtooned(rain)
dsizes_a = dimsizes(rain)
i  = ind_resolve(ind(a1D.ge.2.54),dsizes_a)
print(i(:,0) +"  "+ i(:,1) +"  "+i(:,2))

_______________________________________________
ncl-talk mailing list
ncl-talk at mailman.ucar.edu<mailto:ncl-talk at mailman.ucar.edu>
List instructions, subscriber options, unsubscribe:
https://mailman.ucar.edu/mailman/listinfo/ncl-talk


--
Guest Editor of Special Issue “Extreme Hydroclimatic Events and Prediction”
LINK:
https://www.mdpi.com/journal/water/special_issues/Extreme_Hydroclimatic_Events

Barry H. Lynn, Ph.D
Senior Associate Scientist, Lecturer,
The Institute of the Earth Science,
The Hebrew University of Jerusalem,
Givat Ram, Jerusalem 91904, Israel
Tel: 972 547 231 170
Fax: (972)-25662581

C.E.O, Weather It Is, LTD
Weather and Climate Focus
http://weather-it-is.com
Jerusalem, Israel
Local: 02 930 9525
Cell: 054 7 231 170
Int-IS: x972 2 930 9525
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20201213/73073c55/attachment-0001.html>