[ncl-talk] Unrealistic (strange) wind values using obj_anal_ic interpolation function

Gerardo Montoya gemonga at gmail.com
Thu Mar 7 22:25:41 MST 2019


 Thanks Dennis for your comments: I checked that, the variable lv_HTGL2 in
the GFS analysis, effectively has two values corresponding to 80m and 100
respectively. Here was my error. 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: ncl_filedump
gfs.t18z.pgrb2.0p25.f001.grb2 -v lv_HTGL2 | less Variable No. 0 lv_HTGL2: 2.000
80.00 100.0 However, even 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.

Gerardo

El mar., 5 mar. 2019 a las 22:22, Dennis Shea (<shea at ucar.edu>) escribió:

> Hello,
>
> A few questions/comments/suggestions:
>
> [1]
> I am not sure where you get 10m wind 1st guess.  Adding "*.grb2*" on the
> command line and looking at the *lv_HGT2 *variable:
>
> *%>*  *ncl_filedump*   gfs.t18z.pgrb2.0p25.anl*.grb2*    *-v lv_HTGL2*
> | less
>
> *GRD_P0_L103_GLL0 *(* lv_HTGL2*, lat_0, lon_0 )       ; <= 3D
>
> *lv_HTGL2*:    80.00       100.0      <== meters
>
> [2]  Apparently, you have created a 10 m wind: EG
>           g = addfile (DirAnData+"*wind10m.nc <http://wind10m.nc>"*,"r")
>       I have no idea how this was created.
>
> [3]  Inverting and the latitude and longitude coordinates does *NOT*
> result in the actual grid being reordered
>
> ;invert the latitude order
> lati=lat_g(::-1)
> ;convert to west latitude
> loni=lon_g-360.0
>
> [4] Please *carefully *examine the *coordinates* after each operation
>
>   dirg = "./"
>   filg = "gfs.t18z.pgrb2.0p25.anl"
>   pthg = dirg+filg+".grb2"            ; add .grb extension
>   fg   = *addfile*(pthg, "r")
>  ;print(fg)                           ; ncl_filedump
> gfs.t18z.pgrb2.0p25.anl.grb2| less
>
>   uh   = fg->UGRD_P0_L103_GLL0        ; ( lv_HTGL2, lat_0, lon_0)
>   *printVarSummary*(uh)
>   print("---")
>
>   uh   = uh(:,*::-1*,:)                 ; Invert the grid (S->N)  via NCL
> syntax
>   printVarSummary(uh)
>   print("---")
>
>   uh   =* lonFlip*(uh)                  ; -180 to 180
>   printVarSummary(uh)
>   print("---")
>
>   lat  = uh*&*lat_0                     ; extract the latitudes *associated
> with the variable*
>   lon  = uh*&*lon_0                    ;             longitudes
>   print(lat)                              ; look at values
>   print(lon)
>
> If, after you makes changes, you still have problens
> I will take a look. However, I would need the files:
>
> VarMetar1H.txt
> VieMetar1H.txt
>
> On Tue, Mar 5, 2019 at 7:55 AM Gerardo Montoya <gemonga at gmail.com> wrote:
>
>>  Dear NCL talk colleagues
>>
>> I’m interpolating wind, extracted from METAR, using  obj_anal_ic function.
>>
>> 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 (gfs.t18z.pgrb2.0p25.anl <https://nomads.ncep.noaa.gov/pub/data/nccf/com/gfs/prod/gfs.2019022818/gfs.t18z.pgrb2.0p25.anl>)and use a rscan=(/0.1,0.01/).
>>
>> The METAR information is also taken for this data and time.
>>
>>
>>
>> Fig. Vgfs_A20190228_1800png  shows the GFS analysis wind and Fig. CressmanInterpolation19-02-28_18Z.png, the interpolated wind.
>>
>> These Figs, also contain wind barbs (plotted in black) at airports.
>>
>>
>>
>> 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.
>>
>> Any suggestion will be greatly appreciated.
>> Gerardo Montoya, full professor (retired) at Universidad Nacional de
>> Colombia
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190308/9b0258d6/attachment.html>


More information about the ncl-talk mailing list