<div dir="ltr"><div dir="ltr">
<h1 class="gmail-western" style="font-family:"Liberation Sans",sans-serif;font-size:14pt;font-weight:normal;margin-bottom:0.08in;direction:ltr;color:rgb(0,0,10);line-height:115%;text-align:left"><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">Thanks Den</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">n</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">is
for your com</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">m</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">ents:</span></font></font></font></h1>
<h1 class="gmail-western" style="font-family:"Liberation Sans",sans-serif;font-size:14pt;font-weight:normal;margin-bottom:0.08in;direction:ltr;color:rgb(0,0,10);line-height:115%;text-align:left"><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">I
checked that, the variable lv_HTGL2 in the GFS analysis,</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US"> effectively</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">
has two values corresponding to 80m and 100 respectively. </span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">Here
was my error</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">.
I also found that this level variable, in the GFS forecast it has 3
values corresponding to the levels 2, 80 and 100m respectively, as it
can be seen for the following ncl_filedump:</span></font></font></font></h1>
<h1 class="gmail-western" style="font-family:"Liberation Sans",sans-serif;font-size:14pt;font-weight:normal;margin-bottom:0.08in;direction:ltr;color:rgb(0,0,10);line-height:115%;text-align:left"><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">ncl_filedump
gfs.t18z.pgrb2.0p25.f001.grb2 -v lv_HTGL2 | less</span></font></font></font></h1>
<h1 class="gmail-western" style="font-family:"Liberation Sans",sans-serif;font-size:14pt;font-weight:normal;margin-bottom:0.08in;direction:ltr;color:rgb(0,0,10);line-height:115%;text-align:left"><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">Variable
No. 0</span></font></font></font></h1>
<h1 class="gmail-western" style="font-family:"Liberation Sans",sans-serif;font-size:14pt;font-weight:normal;margin-bottom:0.08in;direction:ltr;color:rgb(0,0,10);line-height:115%;text-align:left"><font color="#00000a"> <font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">lv_HTGL2:</span></font></font></font></h1>
<h1 class="gmail-western" style="font-family:"Liberation Sans",sans-serif;font-size:14pt;font-weight:normal;margin-bottom:0.08in;direction:ltr;color:rgb(0,0,10);line-height:115%;text-align:left"><font color="#00000a"> <font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">2.000
80.00 100.0 </span></font></font></font>
</h1>
<h1 class="gmail-western" style="font-family:"Liberation Sans",sans-serif;font-size:14pt;font-weight:normal;margin-bottom:0.08in;direction:ltr;color:rgb(0,0,10);line-height:115%;text-align:left"><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">However,
</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">even
</span></font></font></font><font color="#00000a"><font face="Liberation Sans, serif"><font style="font-size:12pt" size="3"><span lang="en-US">I
don't know the name of the level variable (instead of lv_HTGL2) I should use
in order to download the 10m wind. Thanks for your help.</span></font></font></font></h1><p>Gerardo<br></p>
</div></div><br><div class="gmail_quote"><div class="gmail_attr" dir="ltr">El mar., 5 mar. 2019 a las 22:22, Dennis Shea (<<a href="mailto:shea@ucar.edu">shea@ucar.edu</a>>) escribió:<br></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><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 target="_blank" 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 class="gmail_attr" dir="ltr">On Tue, Mar 5, 2019 at 7:55 AM Gerardo Montoya <<a target="_blank" href="mailto:gemonga@gmail.com">gemonga@gmail.com</a>> wrote:<br></div><blockquote style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex" class="gmail_quote"><div dir="ltr"><p 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" class="MsoNormal"><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 target="_blank" href="https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2019022818/gfs.t18z.pgrb2.0p25.anl"><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 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" class="MsoNormal"><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 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" class="MsoNormal"><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 target="_blank" href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a target="_blank" rel="noreferrer" href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div>