<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Hello,</div><div><br></div><div>A few questions/comments/suggestions:</div><div><br></div><div>[1] <br></div><div>I am not sure where you get 10m wind 1st guess. Adding "<b>.grb2</b>" on the command line and looking at the <b>lv_HGT2 </b>variable:<br></div><div><br></div><div><b>%></b> <b>ncl_filedump</b> gfs.t18z.pgrb2.0p25.anl<b>.<span style="color:rgb(0,0,255)">grb2</span></b><span style="color:rgb(0,0,255)"> <b>-v lv_HTGL2</b></span> | less<br></div><div><b><br></b></div><div><b>GRD_P0_L103_GLL0 </b>(<b> <span style="color:rgb(0,0,255)">lv_HTGL2</span></b>, lat_0, lon_0 ) ; <= 3D<br></div><div><br></div><div><b>l<span style="color:rgb(0,0,255)">v_HTGL2</span></b>: 80.00 100.0 <== meters<br></div><div><br></div><div>[2] Apparently, you have created a 10 m wind: EG<br></div><div> g = addfile (DirAnData+"<b><a href="http://wind10m.nc">wind10m.nc</a>"</b>,"r")<br></div><div> I have no idea how this was created.<br></div><div><br></div><div>[3] Inverting and the latitude and longitude coordinates does *NOT* result in the actual grid being reordered<br></div><div><br></div><div>;invert the latitude order<br>lati=lat_g(::-1)<br>;convert to west latitude<br>loni=lon_g-360.0</div><div><br></div><div>[4] Please <b>carefully </b>examine the <b>coordinates</b> after each operation<br></div><div><br></div><div> dirg = "./"<br> filg = "gfs.t18z.pgrb2.0p25.anl"<br> pthg = dirg+filg+".grb2" ; add .grb extension<br> fg = <b>addfile</b>(pthg, "r")<br> ;print(fg) ; ncl_filedump gfs.t18z.pgrb2.0p25.anl.grb2| less<br><br> uh = fg->UGRD_P0_L103_GLL0 ; ( lv_HTGL2, lat_0, lon_0)<br> <b>printVarSummary</b>(uh) <br> print("---")<br><br> uh = uh(:,<b>::-1</b>,:) ; Invert the grid (S->N) via NCL syntax<br> printVarSummary(uh) <br> print("---")<br><br> uh =<b> lonFlip</b>(uh) ; -180 to 180<br> printVarSummary(uh) <br> print("---")<br><br> lat = uh<b>&</b>lat_0 ; extract the latitudes <b>associated with the variable</b> <br> lon = uh<b>&</b>lon_0 ; longitudes</div><div> print(lat) ; look at values<br></div><div> print(lon)</div><div><br></div><div>If, after you makes changes, you still have problens<br></div><div>I will take a look. However, I would need the files:</div><div><br></div><div>VarMetar1H.txt</div><div>VieMetar1H.txt</div></div></div></div></div></div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Mar 5, 2019 at 7:55 AM Gerardo Montoya <<a href="mailto:gemonga@gmail.com">gemonga@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><p class="MsoNormal" style="text-align:justify;line-height:12.65pt;background-image:initial;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial;margin:0in 0in 10pt;font-size:11pt;font-family:Calibri,sans-serif"><span style="font-size:12pt;font-family:"Times New Roman",serif"> </span><span style="font-family:"Times New Roman",serif;font-size:12pt">Dear NCL talk
colleagues</span></p>
<pre style="margin:0in 0in 0.0001pt;font-size:10pt;font-family:"Courier New""><span style="font-size:12pt;font-family:"Times New Roman",serif">I’m interpolating wind, extracted from METAR, using obj_anal_ic function. </span></pre><pre style="margin:0in 0in 0.0001pt;font-size:10pt;font-family:"Courier New""><span style="font-size:12pt;font-family:"Times New Roman",serif">I take as first guess the wind at 10m (UGRD_P0_L103_GLL0 and VGRD_P0_L103_GLL0 variables) from GFS analysis at 2019-02-28:18Z (</span><span style="color:black"><a href="https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2019022818/gfs.t18z.pgrb2.0p25.anl" target="_blank"><span style="color:blue">gfs.t18z.pgrb2.0p25.anl</span></a>)</span><span style="font-size:12pt;font-family:"Times New Roman",serif">and use a rscan=(/0.1,0.01/). </span></pre><pre style="margin:0in 0in 0.0001pt;font-size:10pt;font-family:"Courier New""><span style="font-size:12pt;font-family:"Times New Roman",serif">The METAR information is also taken for this data and time. </span></pre><pre style="margin:0in 0in 0.0001pt;font-size:10pt;font-family:"Courier New""><span style="font-size:12pt;font-family:"Times New Roman",serif"> </span></pre><pre style="margin:0in 0in 0.0001pt;font-size:10pt;font-family:"Courier New""><span style="font-size:12pt;font-family:"Times New Roman",serif">Fig. </span><span style="font-family:"Times New Roman",serif">Vgfs_A20190228_1800png</span><span style="font-size:12pt;font-family:"Times New Roman",serif"> shows the GFS analysis wind and Fig. </span><span style="font-family:"Times New Roman",serif">CressmanInterpolation19-02-28_18Z</span><span style="font-size:12pt;font-family:"Times New Roman",serif">.png, the interpolated wind. </span></pre><pre style="margin:0in 0in 0.0001pt;font-size:10pt;font-family:"Courier New""><span style="font-size:12pt;font-family:"Times New Roman",serif">These Figs, also contain wind barbs (plotted in black) at airports. </span></pre><pre style="margin:0in 0in 0.0001pt;font-size:10pt;font-family:"Courier New""><span style="color:black"> </span></pre>
<p class="MsoNormal" style="text-align:justify;line-height:12.65pt;background-image:initial;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial;margin:0in 0in 10pt;font-size:11pt;font-family:Calibri,sans-serif"><span style="font-size:12pt;font-family:"Times New Roman",serif">Comparing these two
Figs., it can be seen that, in some regions and in the vicinity of airports,
The interpolation function works as expected: the interpolated wind deviates
from the large scale flow (trade winds) and tends to curve in the direction
indicated by wind barbs (see the areas highlighted with green arrow). Although,
in some areas, specially over complex terrain (central Colombia) and Cuba, the interpolated wind is extremely
strong (see the areas highlighted with red arrow). It there is the possibility
I made a mistake in my code. However, I can’t find it. I also noted that
Cressman interpolation is very sensible to the rscan values. I’m attaching the
code. The NCL 6.5 version is used.</span></p>
<p class="MsoNormal" style="text-align:justify;line-height:12.65pt;background-image:initial;background-position:initial;background-size:initial;background-repeat:initial;background-origin:initial;background-clip:initial;margin:0in 0in 10pt;font-size:11pt;font-family:Calibri,sans-serif"><span style="font-size:12pt;font-family:"Times New Roman",serif">Any suggestion will be
greatly appreciated.</span></p>
<span style="font-size:12pt;line-height:115%;font-family:"Times New Roman",serif">Gerardo Montoya, full
professor (retired) at Universidad Nacional de Colombia </span> <br></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>