Re: interpolate wrf data to lat/lon points

After posting I saw and error. Please change the print

    print("n="+(n+1)+" XLAT="+ sprintf("%9.3f", XLAT(ij(n,0),ij(n,1)) \
                               +" XLON="+ sprintf("%9.3f",
XLONG(ij(n,0),ij(n,1)) \
                               +" rain="+ sprintf("%9.3f",
rain_exp(ij(n,0),ij(n,1)) \
                                + RAIN="+ sprintf("%9.3f", RAIN(n) )

Dennis Shea wrote:
> [1] The documentation for linint2 and linint2_points states that the
> input x and y arrays be one-dimensional *monotonically
> {in/de}creasing* arrays. You made the input one dimensional
> but the arrays are not monotonic.
>
> [2] Probably best for you to use a WRF provided function:
>
> http://www.ncl.ucar.edu/Document/Functions/wrf.shtml
>
> Specifically,
>
>
> http://www.ncl.ucar.edu/Document/Functions/Built-in/wrf_latlon_to_ij.shtml
>
>
> [3] Please note: wrfhelp_at_ucar.edu if the official WRF site for
> help. They may be better able to help.
>
> Completely untested! If the variables have a time dimension, you will
> have to
> change accordingly.
>
>
> a = addfile(dir+fil,"r")
> rain_exp = wrf_user_getvar(a,"RAINNC",it2)
> XLAT = wrf_user_getvar(a,"XLAT",it1)
> XLONG= wrf_user_getvar(a,"XLONG",it1)
>
> RAIN = asciiread("RAIN_OUT",(/283,3/),"float")
> dimsR = dimsizes(RAIN)
> nobs = dimR(0)
>
> lat=RAIN(:,0)
> lon=RAIN(:,1)
> ij =wrf_latlon_to_ij(XLAT,XLONG,lat,lon)
>
> do n=0,nobs-1
> print("n="+(n+1)+" XLAT="+ sprintf("%9.3f",
> XLAT(ij(n,0),ij(n,1)) \
> +" XLON="+ sprintf("%9.3f",
> XLONG(ij(n,0),ij(n,1)) \
> +" rain="+ sprintf("%9.3f",
> rain_exp(ij(n,0),ij(n,1)) \
> + RAIN="+ sprintf("%9.3f",
> RAIN(ij(n,0),ij(n,1)) )
> end do
> Wolfgang Langhans wrote:
>> Hi,
>> I try to interploate bilinearly from the wrf output to certain lat/lon
>> observation points. Therefor I convert the 2-dimensional XLAT and
>> XLONG to
>> 1-dimensional coordinates XLAT_1d and XLONG_1d. Then interpolating to
>> the
>> observation points does not work, because now the dimensions of
>> rain_exp (WRF
>> raindata) are different as the dimensions of XLAT_1d, XLONG_1d.
>> Could you please help me?
>> Thanks!
>>
>> This is my script:
>>
>> a = addfile(dir+fil,"r")
>> rain_exp = wrf_user_getvar(a,"RAINNC",it2)
>> XLAT = wrf_user_getvar(a,"XLAT",it1)
>> XLONG= wrf_user_getvar(a,"XLONG",it1)
>> XLAT_1d=ndtooned(XLAT)
>> XLONG_1d=ndtooned(XLONG)
>>
>> RAIN = asciiread("RAIN_OUT",(/283,3/),"float")
>>
>> lat=RAIN(:,0)
>> lon=RAIN(:,1)
>> fo = linint2_points(XLONG_1d,XLAT_1d,rain_exp,False,lon,lat,0)
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk_at_ucar.edu
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
>

-- 
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk