<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">Hi Bill,</div>
<div class=""><br class="">
</div>
<div class="">Thanks a lot for your answer. I tried the first one and it gave an error because of the lack of xlat variable, then tried the second routine with giving the parameters, and it is working now.</div>
<div class=""><br class="">
</div>
<div class="">I am now wondering if i will be able to plot my wrf outputs with wrf-python, because I noticed that getvar routine doesn’t work due to missing lat-lon variables, do you have any advice for this?</div>
<div class=""><br class="">
</div>
<div class="">All the best,</div>
<div class=""><br class="">
</div>
<div class="">Huseyin</div>
<div class=""><br class="">
</div>
<br class="">
<div>
<blockquote type="cite" class="">
<div class="">Bill Ladwig <<a href="mailto:ladwig@ucar.edu" class="">ladwig@ucar.edu</a>> şunları yazdı (29 Eki 2018 17:40):</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div dir="ltr" class="">
<div dir="ltr" class="">Hi Huseyin,
<div class=""><br class="">
</div>
<div class="">Without that information in the files, you're going to have to calculate the latitude and longitude values. Since you're using Python, you can use wrf-python which contains a xy_to_ll routine that converts grid points to latitude and longitude
points, as long as you're using one of the standard WRF projections, which it looks like you are (<a href="https://wrf-python.readthedocs.io/en/latest/user_api/generated/wrf.xy_to_ll.html#wrf.xy_to_ll" class="">https://wrf-python.readthedocs.io/en/latest/user_api/generated/wrf.xy_to_ll.html#wrf.xy_to_ll</a>)</div>
<div class=""><br class="">
</div>
<div class="">Note that this will only work if you have the map projection global parameters still in your file. If not, there is an xy_to_ll_proj version where you can specify the parameters yourself. (<a href="https://wrf-python.readthedocs.io/en/latest/user_api/generated/wrf.ll_to_xy_proj.html#wrf.ll_to_xy_proj" class="">https://wrf-python.readthedocs.io/en/latest/user_api/generated/wrf.ll_to_xy_proj.html#wrf.ll_to_xy_proj</a>). </div>
<div class=""><br class="">
</div>
<div class="">Hope this helps,</div>
<div class=""><br class="">
</div>
<div class="">Bill</div>
</div>
</div>
</div>
<br class="">
<div class="gmail_quote">
<div dir="ltr" class="">On Mon, Oct 29, 2018 at 7:17 AM Huseyin Ozdemir via pyngl-talk <<a href="mailto:pyngl-talk@ucar.edu" class="">pyngl-talk@ucar.edu</a>> wrote:<br class="">
</div>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi everyone,<br class="">
<br class="">
I have a problem with reading wrf output file (Projection is Polar Stereographic); i want to compare wrf model output with observations at specific latitude-longitude coordinates, but output netcdf file has no latitude and longitude variables.
<br class="">
<br class="">
I am using my python environment for reading this netcdf file, where PyNGL and PyNIO are already installed. As an example, temp variable in this output has attributes below:<br class="">
<br class="">
short T2 [ Time|1, south_north|150, west_east|150 ] <br class="">
FieldType : 104<br class="">
MemoryOrder : XY <br class="">
description : TEMP at 2 M<br class="">
units : K<br class="">
stagger : <br class="">
coordinates : XLONG XLAT XTIME<br class="">
add_offset : 273.8086<br class="">
scale_factor : -0.0005068767<br class="">
<br class="">
Here, T2 variable has “Time", “south_north" and “west_east” dimensions and “south_north" and "west _east" dimensions have both 150 grids. There is no XLAT and XLONG variables in the netcdf output, so how can i get the temperatures at any lat-lon points with
just grid numbers? <br class="">
<br class="">
Thanks in advance,<br class="">
<br class="">
Huseyin<br class="">
<br class="">
_______________________________________________<br class="">
pyngl-talk mailing list<br class="">
List instructions, subscriber options, unsubscribe:<br class="">
<a href="http://mailman.ucar.edu/mailman/listinfo/pyngl-talk" rel="noreferrer" target="_blank" class="">http://mailman.ucar.edu/mailman/listinfo/pyngl-talk</a><br class="">
</blockquote>
</div>
</div>
</blockquote>
</div>
<br class="">
</body>
</html>