interpolate wrf data to lat/lon points

From: Wolfgang Langhans <Wolfgang.Langhans_at_nyahnyahspammersnyahnyah>
Date: Mon, 7 Apr 2008 18:07:03 +0200

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
Received on Mon Apr 07 2008 - 10:07:03 MDT

This archive was generated by hypermail 2.2.0 : Fri Apr 11 2008 - 11:01:25 MDT