Hi,
I am trying to interpolate WRF output to another grid (gridded observational
rain data). The interpolation works. But which lat/lon attributes do I have to
pass to the interpolated field in order to map it? The new lat/lon are one
dimensional arrays. Here is part of my script:
-----------------------------------------------------------------------------
a = addfile("wrfout_d04_2005-08-20_12_mp6.nc","r")
rain_exp = wrf_user_getvar(a,"RAINNC",it2)-wrf_user_getvar(a,"RAINNC",it1)
XLAT = wrf_user_getvar(a,"XLAT",it1)
XLONG= wrf_user_getvar(a,"XLONG",it1)
fil="2005.nc"
f=addfile(fil,"read")
RR=new((/103,241/),"float")
RR(:,:)=f->RainRate_mmh(5591,0,:,:)
do k=1,24
RR(:,:)=RR(:,:) + f->RainRate_mmh(5567+k,0,:,:)
end do
lat2dobs=new((/103,241/),"float")
lon2dobs=new((/103,241/),"float")
do k=0,102
lat2dobs(k,:)= 45.75 + (47.875-45.75)/102. * k
end do
do l=0,240
lon2dobs(:,l)=5.75 + (10.75-5.75)/240. * l
end do
RR_at_lat2d=lat2dobs
RR_at_lon2d=lon2dobs
lat1dobs = lat2dobs(:,0)
lon1dobs = lon2dobs(0,:)
rain_int = rcm2rgrid(XLONG,XLAT,rain_exp(:,:),lon1dobs,lat1dobs,1)
rain_int_at_lat2d=lat2dobs ?????
rain_int_at_lon2d=lon2dobs ????? Does not work!
-------------------------------------------------------------------------------
Thanks!
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Mar 25 2008 - 13:07:04 MDT
This archive was generated by hypermail 2.2.0 : Mon Mar 31 2008 - 09:11:46 MDT